Search results

The purpose of this paper is to introduce a new homotopy perturbation method (NHPM) to solve Cauchy problem of unidimensional non‐linear diffusion equation.

In this paper a modified version of HPM, which the authors call NHPM, has been presented; this technique performs much better than the HPM. HPM and NHPM start by considering a homotopy, and the solution of the problem under study is assumed to be as the summation of a power series in p, the difference between two methods starts from the form of initial approximation of the solution.

In this article, the authors have applied the NHPM for solving nonlinear Cauchy diffusion equation. In comparison with the homotopy perturbation method (HPM), in the present method, the authors achieve exact solutions while HPM does not lead to exact solutions. The authors believe that the new method is a promising technique in finding the exact solutions for a wide variety of mathematical problems.

The basic idea described in this paper is expected to be further employed to solve other functional equations.

To reduce the heat load of a gas turbine blade, its surface is covered with an outer layer of ceramics with high thermal resistance. The purpose of this paper is the selection of ceramics with such a low heat conduction coefficient and thickness, so that the permissible metal temperature is not exceeded on the metal-ceramics interface due to the loss ofmechanical properties.

Therefore, for given temperature changes over time on the metal-ceramics interface, temperature changes over time on the inner side of the blade and the assumed initial temperature, the temperature change over time on the outer surface of the ceramics should be determined. The problem presented in this way is a Cauchy type problem. When analyzing the problem, it is taken into account that thermophysical properties of metal and ceramics may depend on temperature. Due to the thin layer of ceramics in relation to the wall thickness, the problem is considered in the area in the flat layer. Thus, a one-dimensional non-stationary heat flow is considered.

The range of stability of the Cauchy problem as a function of time step, thickness of ceramics and thermophysical properties of metal and ceramics are examined. The numerical computations also involved the influence of disturbances in the temperature on metal-ceramics interface on the solution to the inverse problem.

The computational model can be used to analyze the heat flow in gas turbine blades with thermal barrier.

A number of inverse problems of the type considered in the paper are presented in the literature. Inverse problems, especially those Cauchy-type, are ill-conditioned numerically, which means that a small change in the inputs may result in significant errors of the solution. In such a case, regularization of the inverse problem is needed. However, the Cauchy problem presented in the paper does not require regularization.

We present an abstract mathematical and numerical analysis for Drift‐Diffusion equation of heterojunction semiconductor devices with Fermi‐Dirac statistic. For the approximation, a mixed finite element method is considered. This can be profitably used in the investigation of the current through the device structure. A peculiar feature of this mixed formulation is that the electric displacement D and the current densities jn and jp for electrons and holes, are taken as unknowns, together with the potential φ and quas‐Fermi levels φn and φp. This enably D, jn and jp to be determined directly and accurately. For decoupled system, existence, uniqueness, regularity and stability results of the approximate solution are given. A priori and a posteriori error estimates are also presented. A nonlinear implicit scheme with local time steps is used. This algorithm appears to be efficient and gives satisfactory results. Numerical results for an heterojunction bipolar transistor, In two dimension, are presented.

To provide a new proof of convergence of the Adomian decomposition series for solving nonlinear ordinary and partial differential equations based upon a thorough examination of the historical milieu preceding the Adomian decomposition method.

Develops a theoretical background of the Adomian decomposition method under the auspices of the Cauchy‐Kovalevskaya theorem of existence and uniqueness for solution of differential equations. Beginning from the concepts of a parametrized Taylor expansion series as previously introduced in the Murray‐Miller theorem based on analytic parameters, and the Banach‐space analog of the Taylor expansion series about a function instead of a constant as briefly discussed by Cherruault et al., the Adomian decompositions series and the series of Adomian polynomials are found to be a uniformly convergent series of analytic functions for the solution u and the nonlinear composite function f(u). To derive the unifying formula for the family of classes of Adomian polynomials, the author develops the novel notion of a sequence of parametrized partial sums as defined by truncation operators, acting upon infinite series, which induce these parametrized sums for simple discard rules and appropriate decomposition parameters. Thus, the defining algorithm of the Adomian polynomials is the difference of these consecutive parametrized partial sums.

The four classes of Adomian polynomials are shown to belong to a common family of decomposition series, which admit solution by recursion, and are derived from one unifying formula. The series of Adomian polynomials and hence the solution as computed as an Adomian decomposition series are shown to be uniformly convergent. Furthermore, the limiting value of the mth Adomian polynomial approaches zero as the index m approaches infinity for the prerequisites of the Cauchy‐Kovalevskaya theorem. The novel truncation operators as governed by discard rules are analogous to an ideal low‐pass filter, where the decomposition parameters represent the cut‐off frequency for rearranging a uniformly convergent series so as to induce the parametrized partial sums.

This paper unifies the notion of the family of Adomian polynomials for solving nonlinear differential equations. Further it presents the new notion of parametrized partial sums as a tool for rearranging a uniformly convergent series. It offers a deeper understanding of the elegant and powerful Adomian decomposition method for solving nonlinear ordinary and partial differential equations, which are of paramount importance in modeling natural phenomena and man‐made device performance parameters.

The purpose of the current flow configurations is to bring to attention the thermophysical aspects of magnetohydrodynamics (MHD) Williamson nanofluid flow under the effects of Joule heating, nonlinear thermal radiation, variable thermal coefficient and activation energy past a rotating stretchable surface.

A mathematical model is examined to study the heat and mass transport analysis of steady MHD Williamson fluid flow past a rotating stretchable surface. Impact of activation energy with newly introduced variable diffusion coefficient at the mass equation is considered. The transport phenomenon is modeled by using highly nonlinear PDEs which are then reduced into dimensionless form by using similarity transformation. The resulting equations are then solved with the aid of fifth-order Fehlberg method.

The rotating fluid, heat and mass transport effects are analyzed for different values of parameters on velocity, energy and diffusion distributions. Parameters like the rotation parameter, Hartmann number and Weissenberg number control the flow field. In addition, the solar radiation, Joule heating, Prandtl number, thermal conductivity, concentration diffusion coefficient and activation energy control the temperature and concentration profiles inside the stretching surface. It can be analyzed that for higher values of thermal conductivity, Eckret number and solar radiation parameter the temperature profile increases, whereas opposite behavior is noticed for Prandtl number. Moreover, for increasing values of temperature difference parameter and thermal diffusion coefficient, the concentration profile shows reducing behavior.

This paper is useful for researchers working in mathematical and theoretical physics. Moreover, numerical results are very useful in industry and daily-use processes.

This study aims to apply a numerical meshless method, namely, the boundary knot method (BKM) combined with the meshless analog equation method (MAEM) in space and use a semi-implicit scheme in time for finding a new numerical solution of the advection–reaction–diffusion and reaction–diffusion systems in two-dimensional spaces, which arise in biology.

First, the BKM is applied to approximate the spatial variables of the studied mathematical models. Then, this study derives fully discrete scheme of the studied models using a semi-implicit scheme based on Crank–Nicolson idea, which gives a linear system of algebraic equations with a non-square matrix per time step that is solved by the singular value decomposition. The proposed approach approximates the solution of a given partial differential equation using particular and homogeneous solutions and without considering the fundamental solutions of the proposed equations.

This study reports some numerical simulations for showing the ability of the presented technique in solving the studied mathematical models arising in biology. The obtained results by the developed numerical scheme are in good agreement with the results reported in the literature. Besides, a simulation of the proposed model is done on buttery shape domain in two-dimensional space.

This study develops the BKM combined with MAEM for solving the coupled systems of (advection) reaction–diffusion equations in two-dimensional spaces. Besides, it does not need the fundamental solution of the mathematical models studied here, which omits any difficulties.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

The purpose of this paper is to develop a mesh-free algorithm based on the least square support vector machines method for numerical simulation of the modified Helmholtz equations.

The proposed method deals with a Cauchy problem for the modified Helmholtz equations. The algorithm converts the problem into a quadratic programming. It can be divided into three steps. First, some training points are allocated. Then, an approximate function is constructed. Finally, the shape parameters are estimated.

The proposed method's stability is discussed. Numerical experiments are conducted to check the efficiency of the algorithm. The proposed method is found to feasible for the ill-posed problems of the modified Helmholtz equations.

The originality lies in that the proposed method is applied to solve the modified Helmholtz equations for the first time, and the expected results are obtained.

Presents a finite element formulation of hot isostatic pressing (HIP) based on a continuum approach using thermal‐elastoviscoplastic constitutive equations with compressibility. The formulation takes into consideration dependence of the viscoplastic part on the porosity. Also takes into account the thermomechanical response, including nonlinear effects in both the thermal and mechanical analyses. Implements the material model in an implicit finite element code. Presents experimental procedures for evaluating the inelastic behaviour of metal powders during densification and experimental data. Chooses the simulation of the dilatometer measurement of a cylindrical component during HIP and manufacturing simulation of a turbine component to near net shape (NNS) as a demonstrator example. Both components are made of a hot isostatically pressed hot‐working martensitic steel. Compares the result of the simulation in the form of the final geometry of the container with the geometry of a real component produced by HIP. Makes a comparison between the calculated and measured deformations during the HIP process for the cylindrical component. Measures the final geometry of the turbine component by means of a computer controlled measuring machine (CMM). Performs the complete process from design and simulation to geometry verification within a computer‐aided concurrent engineering (CACE) system.