Analytical solution for the time-fractional heat conduction equation in spherical coordinate system by the method of variable separation

In this paper,using the fractional Fourier law,we obtain the fractional heat conduction equation with a time-fractional derivative in the spherical coordinate system.The method of variable separation is used to solve the timefractional heat conduction equation.The Caputo fractional derivative of the order 0 〈 α≤ 1 is used.The solution is presented in terms of the Mittag-Leffler functions.Numerical results are illustrated graphically for various values of fractional derivative.

Mathematical Algorithm for Calculating the Velocity Vectors of Fluid by CTA in Spherical Coordinates

This paper presents a mathematical algorithm that determines the fluid flow velocity vector （direction, intensity and orientation）, based on measured voltages on multi-channel hot-wire anemometer. As the voltage on Constant Temperature hot-wire Anemometer （CTA） is non-linear function of velocity and angle of the fluid, inverse function is also non-linear and has several mathematically correct solutions. In the Laboratory of Non-linear Mechanics at the Faculty of Mechanical Engineering in Ljubljana, the authors have decided to try developing multi-charmel hot-wire anemometer with constant temperature at which it is possible to select physically correct solutions from several mathematically correct solutions. The mathematical algorithm works correctly if the range of instrument operation is limited for the value of spherical angles ｜φ｜≤ 60°and ｜ψ｜＇1 ≤ 58°.

SPHERICAL-SYMMETRIC STEADY-STATE RESPONSE OF PIEZOELECTRIC SPHERICAL SHELL UNDER EXTERNAL EXCITATION

Spherical-symmetric steady-state response problem of piezoelectric spherical shell in the absence of body force and free charges is discussed. The steady-state response solutions of mechanical displacement, stresses, strains, potential and electric displacement were derived Sram constitutive relations, geometric and motion equations for the piezoelectric medium under external excitation (i.e. applied surface traction and potential) in spherical coordinate system. As an application of the? general solutions, the problem of an elastic spherical shell with piezoelectric actuator and sensor layers was solved. The results could provide good theoretical basis for the spherical-symmetric dynamic control problem of piezoelectric intelligent structure. Furthermore, the solutions can serve as reference for the research of general dynamic control problem.

Coupling dynamic analysis of a liquid-filled spherical container subject to arbitrary excitation

Using spherical coordinates, the coupling nonlinear dynamic system of a liquid-filled spherical tank, which can be excited discretionarily, is deduced by the H-O varia- tional principle, and the viscous damping is introduced via the liquid dissipation function. The kinetic equations of the coupling system are deduced by the relationship between the velocity of liquid particles and the disturbed liquid surface equation. Normal differential equations are obtained through the Galerkin method. An equivalent mechanical model is developed for liquid sloshing in a spherical tank subject to arbitrary excitation. The fixed and slosh masses, as well as the spring and damping constants, are determined in such a way as to satisfy the principle of equivalence. Numerical simulations illustrate the theoretical results in this paper as well.

A spherical higher-order finite-difference time-domain algorithm with perfectly matched layer

A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.

Using three-dimensional discrete spherical Fourier descriptors based on surface curvature voxels for pollen particle recognition

This paper presents a new method for extract three-dimensional （3D） discrete spherical Fourier descriptors based on surface curvature voxels for pollen particle recognition. In order to reduce the high amount of pollen information and noise disturbance, the geometric normalized curvature voxels with the principal curvedness are first extracted to represent the intrinsic pollen volumetric data. Then the curvature voxels are decomposed into radial and angular components with spherical harmonic transform in spherical coordinates. Finally the 3D discrete Fourier transform is applied to the decomposed curvature voxels to obtain the 3D spherical Fourier descriptors for pollen recognition. Experimental results show that the presented descriptors are invariant to different pollen particle geometric transformations, such as pose change and spatial rotation, and can obtain high recognition accuracy and speed simultaneously.

A Semi-Analytical Approach to Green＇s Functions for Problems in Multiply-Connected Regions on a Spherical Surface

This work aims at potential fields generated by point sources in conductive perforated fragments of spherical shells. Such fields are interpreted as profiles of Green＇s functions of relevant boundary-value problems stated in multiply-connected regions for Laplace equation written in geographical coordinates. Those are efficiently computed by a modification of the method of functional equations, with closed analytical forms preliminary obtained for Green＇s functions for the corresponding simply-connected regions.

Bald Eagle Search Optimization Algorithm Combined with Spherical Random Shrinkage Mechanism and Its Application

Over the last two decades,stochastic optimization algorithms have proved to be a very promising approach to solving a variety of complex optimization problems.Bald eagle search optimization(BES)as a new stochastic optimization algorithm with fast convergence speed has the ability of prominent optimization and the defect of collapsing in the local best.To avoid BES collapse at local optima,inspired by the fact that the volume of the sphere is the largest when the surface area is certain,an improved bald eagle search optimization algorithm(INMBES)integrating the random shrinkage mechanism of the sphere is proposed.Firstly,the INMBES embeds spherical coordinates to design a more accurate parameter update method to modify the coverage and dispersion of the population.Secondly,the population splits into elite and non-elite groups and the Bernoulli chaos is applied to elite group to tap around potential solutions of the INMBES.The non-elite group is redistributed again and the Nelder-Mead simplex strategy is applied to each group to accelerate the evolution of the worst individual and the convergence process of the INMBES.The results of Friedman and Wilcoxon rank sum tests of CEC2017 in 10,30,50,and 100 dimensions numerical optimization confirm that the INMBES has superior performance in convergence accuracy and avoiding falling into local optimization compared with other potential improved algorithms but inferior to the champion algorithm and ranking third.The three engineering constraint optimization problems and 26 real world problems and the problem of extracting the best feature subset by encapsulated feature selection method verify that the INMBES’s performance ranks first and has achieved satisfactory accuracy in solving practical problems.

EXPERIMENTS OF A REDUCED GRID IN LASG/IAP WORLD OCEAN GENERAL CIRCULATION MODELS (OGCMs)

Due to the decrease in grid size associated with the convergence of meridians toward the poles inspherical coordinates, the time steps in many global climate models with finite-difference method are restrictedto be unpleasantly small. To overcome the problem, a reduced grid is introduced to LASG/IAP world oceangeneral circulation models. The reduced grid is implemented successfully in the coarser resolutions versionmodel L30T63 at first. Then, it is carried out in the improved version model LICOM with finer resolutions. Inthe experiment with model L30T63, under time step unchanged though, execution time per single model run isshortened significantly owing to the decrease of grid number and filtering execution in high latitudes. Resultsfrom additional experiments with L30T63 show that the time step of integration can be quadrupled at most inreduced grid with refinement ratio 3. In the experiment with model LICOM and with the model’s original timestep unchanged, the model covered area is extended to the whole globe from its original case with the grid pointof North Pole considered as an isolated island and the results of experiment are shown to be acceptable.

A Global Spectral Element Model for Poisson Equations and Advective Flow over a Sphere

A global spherical Fourier-Legendre spectral element method is proposed to solve Poisson equations and advective flow over a sphere. In the meridional direction, Legendre polynomials are used and the region is divided into several elements. In order to avoid coordinate singularities at the north and south poles in the meridional direction, Legendre-Gauss-Radau points are chosen at the elements involving the two poles. Fourier polynomials are applied in the zonal direction for its periodicity, with only one element. Then, the partial differential equations are solved on the longitude-latitude meshes without coordinate transformation between spherical and Cartesian coordinates. For verification of the proposed method, a few Poisson equations and advective flows are tested. Firstly, the method is found to be valid for test cases with smooth solution. The results of the Poisson equations demonstrate that the present method exhibits high accuracy and exponential convergence. High- precision solutions are also obtained with near negligible numerical diffusion during the time evolution for advective flow with smooth shape. Secondly, the results of advective flow with non-smooth shape and deformational flow are also shown to be reasonable and effective. As a result, the present method is proved to be capable of solving flow through different types of elements, and thereby a desirable method with reliability and high accuracy for solving partial differential equations over a sphere.

APPLICATIONS OF FRACTIONAL EXTERIOR DIFFERENTIAL IN THREE-DIMENSIONAL SPACE

A brief survey of fractional calculus and fractional differential forms was firstly given.The fractional exterior transition to curvilinear coordinate at the origin were discussed and the two coordinate transformations for the fractional differentials for three-dimensional Cartesian coordinates to spherical and cylindrical coordinates are obtained, respectively. In particular, for v=m=1 ,the usual exterior transformations, between the spherical coordinate and Cartesian coordinate, as well as the cylindrical coordinate and Cartesian coordinate, are found respectively, from fractional exterior transformation.

Strong/weak limits and strong/weak differentials of multi-variable indeterminate forms. III: Transformations

Various transforms of the indeterminate forms are presented in this part, which include simplification in spherical coordinates, origin translation, axis alteration, transformation of limit conservation and application of Xh?K0. Fundamental factors for numerical simplification are provided respectively for bi-variable indeterminate forms, tri-variable indeterminate forms and the universal extending multiplier.

On the Axisymmetric Steady Incompressible Beltrami Flows

In this paper, Beltrami vector fields in several orthogonal coordinate systems are obtained analytically and numerically. Specifically, axisymmetric incompressible inviscid steady state Beltrami (Trkalian) fluid flows are obtained with the motivation to model flows that have been hypothesized to occur in tornadic flows. The studied coordinate systems include those that appear amenable to modeling such flows: the cylindrical, spherical, paraboloidal, and prolate and oblate spheroidal systems. The usual Euler equations are reformulated using the Bragg-Hawthorne equation for the stream function of the flow, which is solved analytically or numerically in each coordinate system under the assumption of separability of variables. Many of the obtained flows are visualized via contour plots of their stream functions in the rz-plane. Finally, the results are combined to provide a qualitative quasi-static model for a progression of tornado-like flows that develop as swirl increases. The results in this paper are equally applicable in electromagnetics, where the equivalent concept is that of a force-free magnetic field.

Spectral Matrix Conditioning in Cylindrical and Spherical Elliptic Equations

In the spectral solution of 3-D Poisson equations in cylindrical and spherical coordinates including the axis or the center,it is convenient to employ radial basis functions that depend on the Fourier wavenumber or on the latitudinal mode.This idea has been adopted by Matsushima and Marcus and by Verkley for planar problems and pursued by the present authors for spherical ones.For the Dirichlet boundary value problem in both geometries,original bases have been introduced built upon Jacobi polynomials which lead to a purely diagonal representation of the radial second-order differential operator of all spectral modes.This note details the origin of such a diagonalization which extends to cylindrical and spherical regions the properties of the Legendre basis introduced by Jie Shen for Cartesian domains.Closed form expressions are derived for the diagonal elements of the stiffness matrices as well as for the elements of the tridiagonal mass matrices occurring in evolutionary problems.Furthermore,the bound on the condition number of the spectral matrices associated with the Helmholtz equation are determined,proving in a rigorous way one of the main advantages of the proposed radial bases.

Error Inhomogeneity in the Computation of Spherical Mean Displacement

The traditional method for computing the mean displacement in latitude-longitude coordinates is a spherical meridional-zonal resultant displacement method （MRDM）, which regards the displacement as the resultant vector of the meridional and zonal displacement components. However, there are inhomogeneity and singularity in the computation error of the MRDM, especially at high latitudes. Using the NCEP/NCAR long-term monthly mean wind and idealized wind fields, the inhomogeneity in the MRDM was accessed by using a great circle displacement computing method （GCDM） for non-iterative cases. The MRDM and GCDM were also compared for iteration cases by taking the trajectories from a three-time level reference method as the real trajectories. In the horizontal direction, the GCDM assumes that an air particle moves along its locating great circle and that the magnitude of the displacement equals the arc length of the great circle. The inhomogeneity of the MRDM is evaluated in terms of the horizontal dis- tance error from the products of wind speed, lapse time, and angle difference from the GCDM displacement orient. The non-iterative results show that the mean horizontal displacement computed through the MRDM has both compu- tational and analytical errors. The displacement error of the MRDM depends on the wind speed, wind direction, and the departure latitude of the air particle. It increases with the wind speed and the departure latitude. The displacement magnitude error has a four-wave pattern and the displacement direction error has a two-wave feature in the definition range of the wind direction. The iterative result shows that the displacement magnitude error and angle error of the MRDM and GCDM with respect to the reference method increase with the lapse time and have similar distribution patterns. The mean magnitude error and the angle error of the MRDM are nearly twice as large as those of the GCDM.

An FFT Based Fast Poisson Solver on Spherical Shells

We present a fast Poisson solver on spherical shells.With a special change of variable,the radial part of the Laplacian transforms to a constant coefficient differential operator.As a result,the Fast Fourier Transform can be applied to solve the Poisson equation with O(N^(3) logN)operations.Numerical examples have confirmed the accuracy and robustness of the new scheme.

Higher-Order Compact Scheme for the Incompressible Navier-Stokes Equations in Spherical Geometry

A higher-order compact scheme on the nine point 2-D stencil is developed for the steady stream-function vorticity form of the incompressible Navier-Stokes(NS)equations in spherical polar coordinates,which was used earlier only for the cartesian and cylindrical geometries.The steady,incompressible,viscous and axially symmetric flow past a sphere is used as a model problem.The non-linearity in the N-S equations is handled in a comprehensive manner avoiding complications in calculations.The scheme is combined with the multigrid method to enhance the convergence rate.The solutions are obtained over a non-uniform grid generated using the transformation r=ex while maintaining a uniform grid in the computational plane.The superiority of the higher order compact scheme is clearly illustrated in comparison with upwind scheme and defect correction technique at high Reynolds numbers by taking a large domain.This is a pioneering effort,because for the first time,the fourth order accurate solutions for the problem of viscous flow past a sphere are presented here.The drag coefficient and surface pressures are calculated and compared with available experimental and theoretical results.It is observed that these values simulated over coarser grids using the present scheme aremore accuratewhen compared to other conventional schemes.It has also been observed that the flow separation initially occurred at Re=21.

The ONAD method for solving the SH-wave equation and simulation of the SH-wave propagation in the Earth's mantle

The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD method can effectively suppress the numerical dispersion when coarse grids are used.In this paper,the ONAD method is extended to solve the 2-dimensional SH-wave equation in the spherical coordinates.To investigate the accuracy and the efficiency of the ONAD method,we compare the numerical results calculated by the ONAD method and other methods for both the homogeneous model and the inhomogeneous IASP91 model.The comparisons indicate that the ONAD method for solving the SH-wave equation in the spherical coordinates has the advantages of less numerical dispersion,small memory requirement for computer codes,and fast calculation.As an application,we use the ONAD method to simulate the SH-wave propagating between the Earth's surface and the core-mantle boundary(CMB).Meanwhile,we investigate the SH-wave propagating in the mantle through analyzing the wave-field snapshots in different times and synthetic seismograms.