This paper presents a numerical solution to the inverse problem consisting of recovering time-dependent thermal conductivity and  heat source coefficients  in the one-dimensional  parabolic heat equation.   This  mathematical  formulation  ensures that the inverse problem  has a unique  solution.   However, the problem  is still  ill-posed since small errors  in the input data lead to a drastic  amount  of errors in the output coefficients.  The  finite  difference method  with  the Crank-Nicolson  scheme is adopted  as a direct  solver of the problem in a fixed domain.   The inverse problem is solved sub

A non-polynomial spline (NPS) is an approximation method that relies on the triangular and polynomial parts, so the method has infinite derivatives of the triangular part of the NPS to compensate for the loss of smoothness inherited by the polynomial. In this paper, we propose polynomial-free linear and quadratic spline types to solve fuzzy Volterra integral equations (FVIE) of the 2nd kind with the weakly singular kernel (FVIEWSK) and Abel's type kernel. The linear type algorithm gives four parameters to form a linear spline. In comparison, the quadratic type algorithm gives five parameters to create a quadratic spline, which is more of a credit for the exact solution. These algorithms process kernel singularities with a simple techniqu

We study the physics of flow due to the interaction between a viscous dipole and boundaries that permit slip. This includes partial and free slip, and interactions near corners. The problem is investigated by using a two relaxation time lattice Boltzmann equation with moment-based boundary conditions. Navier-slip conditions, which involve gradients of the velocity, are formulated and applied locally. The implementation of free-slip conditions with the moment-based approach is discussed. Collision angles of 0°, 30°, and 45° are investigated. Stable simulations are shown for Reynolds numbers between 625 and 10 000 and various slip lengths. Vorticity generation on the wall is shown to be affected by slip length, angle of incidence,

The simulation of passively Q-switching is four non – linear first order differential equations. The optimization of passively Q-switching simulation was carried out using the constrained Rosenbrock technique. The maximization option in this technique was utilized to the fourth equation as an objective function; the parameters, γa, γc and β as were dealt with as decision variables. A FORTRAN program was written to determine the optimum values of the decision variables through the simulation of the four coupled equations, for ruby laser Q–switched by Dy +2: CaF2.For different Dy +2:CaF2 molecules number, the values of decision variables was predicted using our written program. The relaxation time of Dy +2: CaF2, used with ruby was

     This paper investigates the simultaneous recovery for two time-dependent coefficients for heat equation under Neumann boundary condition. This problem is considered under extra conditions of nonlocal type. The main issue with this problem is the solution unstable to small contamination of noise in the input data. The Crank-Nicolson finite difference method is utilized to solve the direct problem whilst the inverse problem is viewed as nonlinear optimization problem. The later problem is solved numerically using optimization toolbox from MATLAB. We found that the numerical results are accurate and stable.

In this work, we employ a new normalization Bernstein basis for solving linear Freadholm of fractional integro-differential equations  nonhomogeneous  of the second type (LFFIDE_s). We adopt Petrov-Galerkian method (PGM) to approximate solution of the (LFFIDE_s) via normalization Bernstein basis that yields linear system. Some examples are given and their results are shown in tables and figures, the Petrov-Galerkian method (PGM) is very effective and convenient and overcome the difficulty of traditional methods. We solve this problem (LFFIDE_s) by the assistance of Matlab10.   

A new procedure of depth estimation to the apex of dyke-like sources from
magnetic data has been achieved through the application of a derived equation. The
procedure consists of applying a simple filtering technique to the total magnetic
intensity data profiles resulting from dyke-like bodies, having various depths, widths
and inclination angles. A background trending line is drawn for the filtered profile
and the output profile is considered for further calculations.
Two straight lines are drawn along the maximum slopes of the filtered profile
flanks. Then, the horizontal distances between the two lines at various amplitude
levels are measured and plotted against the amplitudes and the resulted relation is a