Search

In this paper, an efficient multiscale finite element method via local defect-correction technique is developed. This method is used to solve the Schrödinger eigenvalue problem with three-dimensional domain. First, this paper considers a three-dimensional bounded spherical region, which is the truncation of a three-dimensional unbounded region. Using polar coordinate transformation, we successfully transform the three-dimensional problem into a series of one-dimensional eigenvalue problems. These one-dimensional eigenvalue problems also bring singularity. Second, using local refinement technique, we establish a new multiscale finite element discretization method. The scheme can correct the defects repeatedly on the local refinement grid, which can solve the singularity problem efficiently. Finally, the error estimates of eigenvalues and eigenfunctions are also proved. Numerical examples show that our numerical method can significantly improve the accuracy of eigenvalues.

Diffusion equations with Riemann–Liouville fractional derivatives are Volterra integro-partial differential equations with weakly singular kernels and present fundamental challenges for numerical computation. In this paper, we make a convergence analysis of the Schwarz waveform relaxation (SWR) algorithms with Robin transmission conditions (TCs) for these problems. We focus on deriving good choice of the parameter involved in the Robin TCs, at the continuous and fully discretized level. Particularly, at the space-time continuous level, we show that the derived Robin parameter is much better than the one predicted by the well-understood equioscillation principle. At the fully discretized level, the problem of determining a good Robin parameter is studied in the convolution quadrature framework, which permits us to precisely capture the effects of different temporal discretization methods on the convergence rate of the SWR algorithms. The results obtained in this paper will be preliminary preparations for our further study of the SWR algorithms for integro-partial differential equations.

In order to solve the complicated process and low efficiency and low accuracy of solving a class of matrix equations, this paper introduces the linear saturated system model neural network architecture to solve the bisymmetric solution of a class of matrix equations. Firstly, a class of matrix equations is constructed to determine the key problems of solving the equations. Secondly, the linear saturated system model neural network structure is constructed to determine the characteristic parameters in the process of bisymmetric solution. Then, the matrix equations is solved by using backpropagation neural network topology. Finally, the class normalization is realized by using the objective function of bisymmetric solution, and the bisymmetric solution of a class of matrix equations is realized. In order to verify the solving effect of the method in this paper, three indexes (accuracy, correction accuracy, and solving time) are designed in the experiment. The experimental results show that the proposed method can effectively reduce the solving time, can improve the accuracy and correction effect of the bisymmetric solution, and has high practicability.

This paper presents some new formulas and properties of the generalized fractional derivative with non-singular kernel that covers various types of fractional derivatives such as the Caputo–Fabrizio fractional derivative, the Atangana–Baleanu fractional derivative, and the weighted Atangana–Baleanu fractional derivative. These new properties extend many recent results existing in the literature. Furthermore, the paper proposes some interesting inequalities that estimate the generalized fractional derivatives of some specific functions. These inequalities can be used to construct Lyapunov functions with the aim to study the global asymptotic stability of several fractional-order systems arising from diverse fields of science and engineering.

This paper proposes a new class of difference methods with intrinsic parallelism for solving the Burgers–Fisher equation. A new class of parallel difference schemes of pure alternating segment explicit-implicit (PASE-I) and pure alternating segment implicit-explicit (PASI-E) are constructed by taking simple classical explicit and implicit schemes, combined with the alternating segment technique. The existence, uniqueness, linear absolute stability, and convergence for the solutions of PASE-I and PASI-E schemes are well illustrated. Both theoretical analysis and numerical experiments show that PASE-I and PASI-E schemes are linearly absolute stable, with 2-order time accuracy and 2-order spatial accuracy. Compared with the implicit scheme and the Crank–Nicolson (C-N) scheme, the computational efficiency of the PASE-I (PASI-E) scheme is greatly improved. The PASE-I and PASI-E schemes have obvious parallel computing properties, which show that the difference methods with intrinsic parallelism in this paper are feasible to solve the Burgers–Fisher equation.

In this paper, by using the partial order method, the existence and uniqueness of a solution for systems of a class of abstract operator equations in Banach spaces are discussed. The result obtained in this paper improves and unifies many recent results. Two applications to the system of nonlinear differential equations and the systems of nonlinear differential equations in Banach spaces are given, and the unique solution and interactive sequences which converge the unique solution and the error estimation are obtained.

Linearized model of the system is often used in control design. It is generally believed that we can obtain the linearized model as long as the Taylor expansion method is used for the nonlinear model. This paper points out that the Taylor expansion method is only applicable to the linearization of the original nonlinear function. If the Taylor expansion is used for the derived nonlinear equation, wrong results are often obtained. Taking the linearization model of the maglev system as an example, it is shown that the linearization should be carried out with the process of equation derivation. The model is verified by nonlinear system simulation in Simulink. The method in this paper is helpful to write the linearized equation of the control system correctly.

An adaptive high-order difference solution about a 2D nonlinear degenerate singular reaction-diffusion equation with a convection term is initially proposed in the paper. After the first and the second central difference operator approximating the first-order and the second-order spatial derivative, respectively, the higher-order spatial derivatives are discretized by applying the Taylor series rule and the temporal derivative is discretized by using the Crank–Nicolson (CN) difference scheme. An alternating direction implicit (ADI) scheme with a nonuniform grid is built in this way. Meanwhile, accuracy analysis declares the second order in time and the fourth order in space under certain conditions. Sequentially, the high-order scheme is performed on an adaptive mesh to demonstrate quenching behaviors of the singular parabolic equation and analyse the influence of combustion chamber size on quenching. The paper displays rationally that the proposed scheme is practicable for solving the 2D quenching-type problem.

This paper considers some multiscale radial basis function collocation methods for solving the two-dimensional Monge–Ampère equation with Dirichlet boundary. We discuss and study the performance of the three kinds of multiscale methods. The first method is the cascadic meshfree method, which was proposed by Liu and He (2013). The second method is the stationary multilevel method, which was proposed by Floater and Iske (1996), and is used to solve the fully nonlinear partial differential equation in the paper for the first time. The third is the hierarchical radial basis function method, which is constructed by employing successive refinement scattered data sets and scaled compactly supported radial basis functions with varying support radii. Compared with the first two methods, the hierarchical radial basis function method can not only solve the present problem on a single level with higher accuracy and lower computational cost but also produce highly sparse nonlinear discrete system. These observations are obtained by taking the direct approach of numerical experimentation.

This research paper is about the existence and uniqueness of the coupled system of nonlinear fractional Langevin equations with multipoint and nonlocal integral boundary conditions. The Caputo fractional derivative is used to formulate the fractional differential equations, and the fractional integrals mentioned in the boundary conditions are due to Atangana–Baleanu and Katugampola. The existence of solution has been proven by two main fixed-point theorems: O’Regan’s fixed-point theorem and Krasnoselskii’s fixed-point theorem. By applying Banach’s fixed-point theorem, we proved the uniqueness result for the concerned problem. This research paper highlights the examples related with theorems that have already been proven.

The paper investigates the dynamical behaviors of a two-species discrete predator-prey system with Crowley–Martin functional response incorporating prey refuge proportional to prey density. The existence of equilibrium points, stability of three fixed points, period-doubling bifurcation, Neimark–Sacker bifurcation, Marottos chaos, and Control Chaos are analyzed for the discrete-time domain. The time graphs, phase portraits, and bifurcation diagrams are obtained for different parameters of the model. Numerical simulations and graphics show that the discrete model exhibits rich dynamics, which also present that the system is a chaotic and complex one. This paper attempts to present a feedback control method which can stabilize chaotic orbits at an unstable equilibrium point.

The well-known Stokes’ problems are reexamined by applying the Adomian decomposition method (ADM) associated with other mathematical techniques in this paper. Both the finite-depth (bounded) and infinite-depth (unbounded) cases are analyzed. The present paper raises and deals with two major concerns. The first one is that, for Stokes’ problems, it lacks one boundary condition at the expansion point to fully determine all coefficients of the ADM solution in which an unknown function appears. This unknown function which is dependent on the transformed variable will be determined by the boundary condition at the far end. The second concern is that the derived solution begins to deviate from the exact solution as the spatial variable grows for the unbounded problems. This can be greatly improved by introducing the Padé approximant to satisfy the boundary condition at the far end. For the second problems, the derived ADM solution can be easily separated into the steady-state and the transient parts for a deeper comprehension of the flow. The present result shows an excellent agreement with the exact solution. The ADM is therefore verified to be a reliable mathematical method to analyze Stokes’ problems of finite and infinite depths.

Finding exact solutions of nonlinear equations plays an important role in nonlinear science, especially in engineering and mathematical physics. In this paper, we employed the complex method to get eight exact solutions of the modified BBM equation for the first time, including two elliptic function solutions, two simply periodic solutions, and four rational function solutions. We used the exp − ϕ z -expansion methods to get fourteen forms of solutions of the modified BBM equation. We also used the sine-cosine method to obtain eight styles’ exact solutions of the modified BBM equation. Only the complex method can obtain elliptic function solutions. We believe that the complex method presented in this paper can be more effectively applied to seek solutions of other nonlinear evolution equations.

In this paper, the dependent competing risks’ model is considered in the constant-stress accelerated life test under the adaptive type-II progressive hybrid-censored scheme. The dependency between failure causes is modeled by Marshall–Olkin bivariate Gompertz distribution. The scale and shape parameters in the model both change with the stress levels, and the failure causes of some test units are unknown. Then, the maximum likelihood estimations and approximation confidence intervals of the unknown parameters are considered. And, the necessary and sufficient condition is established for the existence and uniqueness of the maximum likelihood estimations for unknown parameters. The Bayes approach is also employed to estimate the unknown parameters under suitable prior distributions. The Bayes estimations and highest posterior credible intervals of the unknown parameters are obtained. Finally, a simulation experiment has been performed to illustrate the methods proposed in this paper.

This paper focuses on a singular boundary value (SBV) problem of nonlinear fractional differential (NFD) equation defined as follows: D 0 + β υ τ + f τ , υ τ = 0 , τ ∈ 0,1 , υ 0 = υ ′ 0 = υ ″ 0 = υ ″ 1 = 0 , where 3 < β ≤ 4 , D 0 + β is the standard Riemann–Liouville fractional (RLF) derivative. The nonlinear function f τ , υ τ might be singular on the spatial and temporal variables. This paper proves that a positive solution to the SBV problem exists and is unique, taking advantage of Green’s function through a fixed-point (FP) theory on cones and mixed monotone operators.

Solving nonlinear equation systems for engineering applications is one of the broadest and most essential numerical studies. Several methods and combinations were developed to solve such problems by either finding their roots mathematically or formalizing such problems as an optimization task to obtain the optimal solution using a predetermined objective function. This paper proposes a new algorithm for solving square and nonsquare nonlinear systems combining the genetic algorithm (GA) and the homotopy analysis method (HAM). First, the GA is applied to find out the solution. If it is realized, the algorithm is terminated at this stage as the target solution is determined. Otherwise, the HAM is initiated based on the GA stage’s computed initial guess and linear operator. Moreover, the GA is utilized to calculate the optimum value of the convergence control parameter (h) algebraically without plotting the h-curves or identifying the valid region. Four test functions are examined in this paper to verify the proposed algorithm’s accuracy and efficiency. The results are compared to the Newton HAM (NHAM) and Newton homotopy differential equation (NHDE). The results corroborated the superiority of the proposed algorithm in solving nonlinear equation systems efficiently.

This paper is concerned with the quantized identification of rational systems, where the systems’ output is quantized by a logarithmic quantizer. Under the assumption that the systems’ input is periodic, the identification procedure is categorized into two steps. The first step is to identify the noise-free output of systems based on the quantized data. The second is to identify the unknown parameter based on the input and the estimation of the noise-free output. The identification algorithm is also summarized. Asymptotic convergence of the estimators is analyzed in detail, which shows that the estimators are convergent almost everywhere. A numerical example is given to illustrate the results obtained in this paper.

Recently, the development of neural network method for solving differential equations has made a remarkable progress for solving fractional differential equations. In this paper, a neural network method is employed to solve time-fractional telegraph equation. The loss function containing initial/boundary conditions with adjustable parameters (weights and biases) is constructed. Also, in this paper, a time-fractional telegraph equation was formulated as an optimization problem. Numerical examples with known analytic solutions including numerical results, their graphs, weights, and biases were also discussed to confirm the accuracy of the method used. Also, the graphical and tabular results were analyzed thoroughly. The mean square errors for different choices of neurons and epochs have been presented in tables along with graphical presentations.

This paper uses a merit function derived from the Fishcher–Burmeister function and formulates box-constrained stochastic variational inequality problems as an optimization problem that minimizes this merit function. A sufficient condition for the existence of a solution to the optimization problem is suggested. Finally, this paper proposes a Monte Carlo sampling method for solving the problem. Under some moderate conditions, comprehensive convergence analysis is included as well.

This paper studies and analyzes a model describing the flow of contaminated brines through the porous media under severe thermal conditions caused by the radioactive contaminants. The problem is approximated based on combining the mixed finite element method with the modified method of characteristics. In order to solve the resulting algebraic nonlinear equations efficiently, a two-grid method is presented and discussed in this paper. This approach includes a small nonlinear system on a coarse grid with size H and a linear system on a fine grid with size h . It follows from error estimates that asymptotically optimal accuracy can be obtained as long as the mesh sizes satisfy H = O h 1 / 3 .

In this paper, we present a smoothing Newton method for solving the monotone weighted complementarity problem (WCP). In each iteration of our method, the iterative direction is achieved by solving a system of linear equations and the iterative step length is achieved by adopting a line search. A feature of the line search criteria used in this paper is that monotone and nonmonotone line search are mixed used. The proposed method is new even when the WCP reduces to the standard complementarity problem. Particularly, the proposed method is proved to possess the global convergence under a weak assumption. The preliminary experimental results show the effectiveness and robustness of the proposed method for solving the concerned WCP.

In this paper, we deal with the global behavior of the positive solutions of the system of k -difference equations u n + 1 1 = α 1 u n − 1 1 / β 1 + α 1 u n − 2 2 r 1 , u n + 1 2 = α 2 u n − 1 2 / β 2 + α 2 u n − 2 3 r 2 , … , u n + 1 k = α k u n − 1 k / β k + α k u n − 2 1 r k , n ∈ ℕ 0 , where the initial conditions u − l i l = 0,1,2 are nonnegative real numbers and the parameters α i , β i , γ i , and r i are positive real numbers for i = 1,2 , … , k , by extending some results in the literature. By the end of the paper, we give three numerical examples to support our theoretical results related to the system with some restrictions on the parameters.

In this paper, we investigate pattern dynamics with multivariable by using the method of matrix analysis and obtain a condition under which the system loses stability and Turing bifurcation occurs. In addition, we also derive the amplitude equation with multivariable. This is an effective tool to investigate multivariate pattern dynamics. The example and simulation used in this paper validate our theoretical results. The method presented is a novel approach to the investigation of specific real systems based on the model developed in this paper.

In this paper, the delayed doubly stochastic linear quadratic optimal control problem is discussed. It deduces the expression of the optimal control for the general delayed doubly stochastic control system which contained time delay both in the state variable and in the control variable at the same time and proves its uniqueness by using the classical parallelogram rule. The paper is concerned with the generalized matrix value Riccati equation for a special delayed doubly stochastic linear quadratic control system and aims to give the expression of optimal control and value function by the solution of the Riccati equation.

In this paper, a novel chaotic new three-dimensional system has been studied by Zhang et al. in 2012. In the system, there are three control parameters and three different nonlinear terms which governed equations. Zhang et al. studied elementary (preliminary) dynamic properties of the chaotic new three-dimensional system by means of bifurcation diagram, maximum Lyapunov exponent, phase portraits, dynamics behaviors by changing some parameters etc., using all possible theoretical analysis and numerical simulation. In this paper, we have demonstrated its complete synchronization. The proposed results are verified by numerical simulations.

This paper is intended to study the limit theorem of Markov chain function in the environment of single infinite Markovian systems. Moreover, the problem of the strong law of large numbers in the infinite environment is presented by means of constructing martingale differential sequence for the measurement under some different sufficient conditions. If the sequence of even functions gnx,n≥0 satisfies different conditions when the value ranges of x are different, we have obtained SLLN for function of Markov chain in the environment of single infinite Markovian systems. In addition, the paper studies the strong convergence of the weighted sums of function for finite state Markov Chains in single infinitely Markovian environments. Although the similar conclusions have been carried out, the difference results performed by previous scholars are that we give weaker different sufficient conditions of the strong convergence of weighted sums compared with the previous conclusions.

One of the most popular statistical models is a low-order polynomial response surface model, i.e., a polynomial of first order or second order. These polynomials can be used for global metamodels in weakly nonlinear simulation to approximate their global tendency and local metamodels in response surface methodology (RSM), which has been studied in various applications in engineering design and analysis. The order of the selected polynomial determines the number of sampling points (input combinations) and the resulting accuracy (validity, adequacy). This paper derives a novel method to obtain an accurate high-order polynomial while requiring fewer sampling points. This method uses a two-stage procedure such that the second stage modifies the low-order polynomial estimated in the first stage; this second stage does not require new points. This paper evaluates the performance of the method numerically by using several test functions. These numerical results show that the proposed method can provide more accurate predictions than the traditional method.

Current density distribution in electric arcs inside low voltage circuit breakers is a crucial parameter for us to understand the complex physical behavior during the arcing process. In this paper, we investigate the inverse problem of reconstructing the current density distribution in arcs by inverting the magnetic fields. A simplified 2D arc chamber is considered. The aim of this paper is the computational side of the regularization method, regularization parameter selection strategies, and the estimation of systematic error. To address the ill-posedness of the inverse problem, Tikhonov regularization is analyzed, with the regularization parameter chosen by Morozov’s discrepancy principle, the L-curve, the generalized cross-validation, and the quasi-optimality criteria. The provided range of regularization parameter selection strategies is much wider than in the previous works. Effects of several features on the performance of these criteria have been investigated, including the signal-to-noise ratio, dimension of measurement space, and the measurement distance. The numerical simulations show that the generalized cross-validation and quasi-optimality criteria provide a more satisfactory performance on the robustness and accuracy. Moreover, an optimal measurement distance can be expected when using a planner sensor array to perform magnetic measurements.

By using the semidiscrete method of differential equations, a new version of discrete analogue of stochastic fuzzy BAM neural networks was formulated, which gives a more accurate characterization for continuous-time stochastic neural networks than that by the Euler scheme. Firstly, the existence of the 2p-th mean almost periodic sequence solution of the discrete-time stochastic fuzzy BAM neural networks is investigated with the help of Minkowski inequality, Hölder inequality, and Krasnoselskii’s fixed point theorem. Secondly, the 2p-th moment global exponential stability of the discrete-time stochastic fuzzy BAM neural networks is also studied by using some analytical skills in stochastic theory. Finally, two examples with computer simulations are given to demonstrate that our results are feasible. The main results obtained in this paper are completely new, and the methods used in this paper provide a possible technique to study 2p-th mean almost periodic sequence solution and 2p-th moment global exponential stability of semidiscrete stochastic fuzzy models.

This paper first provides a concept of almost sure stability for uncertain delay differential equations and analyzes this new sort of stability. In addition, this paper derives three sufficient conditions for uncertain delay differential equations being stable almost surely. Finally, the relationship between almost sure stability and stability in measure for uncertain delay differential equations is discussed.

This paper presents a modified Precise Integration Method (PIM) for long-time duration dynamic analysis. The fundamental solution and loading operator matrices in the first time substep are numerically computed employing a known unconditionally stable method (referred to as original method in this paper). By using efficient recursive algorithms to evaluate these matrices in the first time-step, the same numerical results as those using the original method with small time-step are obtained. The proposed method avoids the need of matrix inversion and numerical quadrature formulae, while the particular solution obtained has the same accuracy as that of the homogeneous solution. Through setting a proper value of the time substep, satisfactory accuracy and numerical dissipation can be achieved.

The purpose of this paper is to explore modeling mechanism of a nonhomogeneous multivariable grey prediction NMGM(1, m, kα) model and its application. Although multi-variable grey prediction MGM(1, m) model has been employed in many fields, its prediction results are not always satisfactory. Traditional MGM(1, m) model is constructed on the hypothesis that original data sequences are in accord with homogeneous index trend; however, the nonhomogeneous index data sequences are the most common data existing in all systems, and how to handle multivariable nonhomogeneous index data sequences is an urgent problem. This paper proposes a novel nonhomogeneous multivariable grey prediction model termed NMGM(1, m, kα) to deal with those data sequences that are not in accord with homogeneous index trend. Based on grey prediction theory, by least square method and solutions of differential equations, the modeling mechanism and time response function of the proposed model are expounded. A case study demonstrates that the novel model provides preferable prediction performance compared with traditional MGM(1, m) model. This work is an extension of the multivariable grey prediction model and enriches the study of grey prediction theory.

This paper investigates the consistency definition and the weight-deriving method for additive interval fuzzy preference relations (IFPRs) using a particular characterization based on logarithms. In a recently published paper, a new approach with a parameter is developed to obtain priority weights from fuzzy preference relations (FPRs), then a new consistency definition for the additive IFPRs is defined, and finally linear programming models for deriving interval weights from consistent and inconsistent IFPRs are proposed. However, the discussion of the parameter value is not adequate and the weights obtained by the linear models for inconsistent IFPRs are dependent on alternative labels and not robust to permutations of the decision makers’ judgments. In this paper, we first investigate the value of the parameter more thoroughly and give the closed form solution for the parameter. Then, we design a numerical example to illustrate the drawback of the linear models. Finally, we construct a linear model to derive interval weights from IFPRs based on the additive transitivity based consistency definition. To demonstrate the effectiveness of our proposed method, we compare our method to the existing method on three numerical examples. The results show that our method performs better on both consistent and inconsistent IFPRs.

The concept of diagonally invariant exponential stability (DIES) was originally introduced for single-model linear systems and subsequently expanded in the study of linear systems with interval-type uncertainties and linear systems with arbitrary switching. The results presented in this article refer to new approaches to DIES characterization for arbitrary switching systems, which exploit mathematical tools completely different from earlier work. The previous papers are based on the properties of matrix norms and measures applied to the constituent matrices defining the switching system, while the present paper uses the eigenvalues and eigenvectors of the column and row representatives built for a set of matrices derived from the constituent matrices of the switching system. The applicability of previous and new results, respectively, is illustrated by case studies (in both continuous- and discrete-time) that lead to relevant comparisons between the two classes of analysis methods.

The complexity of cancer has motivated the development of different approaches to understand the dynamics of this large group of diseases. One that may allow us to better comprehend the behavior of cancer cells, in both short- and long-term, is mathematical modelling through ordinary differential equations. Several ODE mathematical models concerning tumor evolution and immune response have been formulated through the years, but only a few may exhibit chaotic attractors and oscillations such as stable limit cycles and periodic orbits; these dynamics are not that common among cancer systems. In this paper, we apply the Localization of Compact Invariant Sets (LCIS) method and Lyapunov stability theory to investigate the global dynamics and the main factors involved in tumor growth and immune response for a chaotic-cancer system presented by Itik and Banks in 2010. The LCIS method allows us to compute what we define as the localizing domain, which is formulated by the intersection of all lower and upper bounds of each cells population in the nonnegative octant, R+,03. Bounds of this domain are given by inequalities in terms of the system parameters. Then, we apply Lyapunov stability theory and LaSalle’s invariance principle to establish existence conditions of a global attractor. The latter implies that given any nonnegative initial condition, all trajectories will go to the largest compact invariant set (a stable equilibrium point, limit cycles, periodic orbits, or a chaotic attractor) located either inside or at the boundaries of the localizing domain. In order to complement our analysis, numerical simulations are performed throughout the paper to illustrate all mathematical results and to better explain their biological implications.