There are many numerical algorithms for solving the fractional-order differential equations(FODEs).The numerical algorithms are very different,and it is difficult to compare their *** solve this problem,some different...
详细信息
There are many numerical algorithms for solving the fractional-order differential equations(FODEs).The numerical algorithms are very different,and it is difficult to compare their *** solve this problem,some different FODEs with the known analytical solution are designed and proposed,they could be used as the benchmark problems for testing the numerical algorithms.A Simulink block diagram scheme is proposed for solving these benchmark problems,and the computing errors and the running times are *** benchmark problems and the solutions are constituted as a framework,and the numerical algorithms for solving the FODEs could be compared in the same *** comparing result could assess which algorithm is better to a concrete FODE.
The paper deals with the numerical analysis of large elastic-plastic deformation behavior of anisotropically damaged ductile solids based on a generalized macroscopic theory within the framework of nonlinear continuum...
详细信息
The paper deals with the numerical analysis of large elastic-plastic deformation behavior of anisotropically damaged ductile solids based on a generalized macroscopic theory within the framework of nonlinear continuum damage mechanics. Estimates of the stress and strain histories are obtained from a straightforward numerical integration algorithm based on operator split methodology which employs an inelastic (damage plastic) predictor followed by an elastic corrector step. The finite element method is used to approximate the linearized variational problem. Furthermore, identification of material parameters is discussed. numerical simulation of the elastic-plastic deformation behavior of damaged tension specimens demonstrate the efficiency of the formulation. (C) 2003 Elsevier B.V. All rights reserved.
In this paper the authors present the idea for reconstructing one-photon states. Reconstructing a quantum state means measuring the probability distribution P that allows one to write the density operator for the anal...
详细信息
In this paper the authors present the idea for reconstructing one-photon states. Reconstructing a quantum state means measuring the probability distribution P that allows one to write the density operator for the analysed state. The most commonly known approach for the quantum reconstruction is the quantum tomography. Our alternative method assumes that the analysed field is coupled with the reference field which is described by the parameters settled during a measurement. In the proposed gedankenexperiment the degree of linear polarisation (DOLP) of this combined beam is measured using a rotating linear polariser. We state that it is possible to obtain the P-function by changing the parameters of reference beams and by preparing the series of DOLP measurements. This series of data leads to the system of equations. The values of the P-function at chosen points are the unknowns of this system. This article focuses on the numerical algorithm for solving these equations.
For automated fiber placement (AFP), it is crucial to research the contact characteristics between tow and roller during the tow feeding process. However, it has not been widely explored. In this paper, a theoretical ...
详细信息
For automated fiber placement (AFP), it is crucial to research the contact characteristics between tow and roller during the tow feeding process. However, it has not been widely explored. In this paper, a theoretical model of contact characteristics for the tow feeding process is established by analyzing the distribution of normal stress, tangential stress, adhesion zone, and slip zone, respectively. The model proposed a numerical algorithm for calculating the contact characteristic. Then the contact characteristic of the tow feeding system can be acquired and validated by the experiment and finite element method. The model that provides theoretical support for AFP optimization. Furthermore, proper pressure (P = 0.6Mpa) and tension (T = 3 N) for towpreg feeding under experimental conditions were obtained.
Introduces a statistical quantity, termed the extended Hadamard variance, for the spectral analysis of the frequency fluctuations nu (t) of an oscillator. The samples of the time series nu (t) are the gated averages n...
详细信息
Introduces a statistical quantity, termed the extended Hadamard variance, for the spectral analysis of the frequency fluctuations nu (t) of an oscillator. The samples of the time series nu (t) are the gated averages nu tau or counts of a digital frequency counter spaced by the dead time tau d. The analysis technique and the numerical algorithm combine the flexibility of the extended two-sample variance with a much higher resolution spectral analysis capability, as shown by experimental results. The limitations are also indicated.
This paper presents a mass flow model that includes gravity force, material stresses, drag force and topography effects solving a set of hyperbolic partial differential equations by using a so-called depth-averaged te...
详细信息
This paper presents a mass flow model that includes gravity force, material stresses, drag force and topography effects solving a set of hyperbolic partial differential equations by using a so-called depth-averaged technique. The model is non-linear and general enough to tackle various problems of interest for geophysics and environmental engineering, such as the dynamic evolution of flow-like avalanches, the dam break problem (involving only water flow) and the generation of tsunami waves by landslides. The model is based on a Eulerian fluid solver, using a second-order central scheme with a minmod-like limiter;is tested against a number of typical benchmark cases, including analytical solutions and experimental laboratory data;and also compared with other numerical codes. Through this model, we study a historical tsunamigenic event occurred in 1783 in Scilla, Italy, that resulted to be catastrophic with a toll exceeding 1500 fatalities. The landslide is reconstructed by a mixture debris flow, and results are compared with the observational data and other numerical simulations.
We analyze a mathematical model of a tumor that includes the interactions among tumor cells (live or dead), macrophages, endothelial cells, and the cytokines based on a system of partial differential equations which d...
详细信息
We analyze a mathematical model of a tumor that includes the interactions among tumor cells (live or dead), macrophages, endothelial cells, and the cytokines based on a system of partial differential equations which describes inhibition of tumor growth by granulocyte-macrophage colony-stimulating factor (GM-CSF) treatment. The aim of the model is to predict the growth of the tumor volume under partial blocking of hypoxia-inducible factor (HIF)-1 or stabilization of HIF-2, with injection of GM-CSF It was assumed that tumor is a spherical and the treatment by GM-CSF is not optimized. In this paper, we do not assume that tumor is a spherical model and we develop a dual dynamic programming approach to formulate a verification condition of an approximate treatment for an injection of GM-CSF. Next, we analyze numerically an epsilon-optimal treatment.
The dynamic process of frictionless contact between a viscoelastic body and a reactive foundation is modelled, analyzed, and simulated. The contact is adhesive and it is described by introducing an internal variable, ...
详细信息
The dynamic process of frictionless contact between a viscoelastic body and a reactive foundation is modelled, analyzed, and simulated. The contact is adhesive and it is described by introducing an internal variable, the bonding field beta, which measures the fractional density of active bonds. The evolution of beta is described by an ordinary differential equation that depends on the process history, taking into account possible adhesive degradation during cycles of debonding and rebonding. The existence of the unique weak solution of the model is proved by using arguments of nonlinear evolutionary equations with monotone operators and a fixed-point theorem. A fully discrete numerical scheme is proposed for the model and implemented in a computer code. numerical simulations of one- and two-dimensional examples are presented. (C) 2003 Elsevier Science Ltd. All rights reserved.
It is known that the generating function of a sequence of Toeplitz matrices may not describe the asymptotic distribution of the eigenvalues of the considered matrix sequence in the non-Hermitian setting. In a recent w...
详细信息
It is known that the generating function of a sequence of Toeplitz matrices may not describe the asymptotic distribution of the eigenvalues of the considered matrix sequence in the non-Hermitian setting. In a recent work, under the assumption that the eigenvalues are real, admitting an asymptotic expansion whose first term is the distribution function, fast algorithms computing all the spectra were proposed in different settings. In the current work, we extend this idea to non-Hermitian Toeplitz matrices with complex eigenvalues, in the case where the range of the generating function does not disconnect the complex field or the limiting set of the spectra, as the matrix-size tends to infinity, has one nonclosed analytic arc. For a generating function having a power singularity, we prove the existence of an asymptotic expansion, that can be used as a theoretical base for the respective numerical algorithm. Different generating functions are explored, highlighting different numerical and theoretical aspects;for example, non-Hermitian and complex symmetric matrix sequences, the reconstruction of the generating function, a consistent eigenvalue ordering, the requirements of high-precision data types. Several numerical experiments are reported and critically discussed, and avenues of possible future research are presented.
numerical methods are proposed for constructing Nash and Stackelberg solutions in a two-player linear non-zero-sum positional differential game with terminal cost functionals and geometric constraints on the players...
详细信息
numerical methods are proposed for constructing Nash and Stackelberg solutions in a two-player linear non-zero-sum positional differential game with terminal cost functionals and geometric constraints on the players' controls. The formalization of the players' strategies and of the motions generated by them is based on the formalization and results from the theory of positional zero-sum differential games developed by N.N. Krasovskii and his school. It is assumed that the game is reduced to a planar game and the constraints on the players' controls are given in the form of convex polygons. The problem of finding solutions of the game may be reduced to solving nonstandard optimal control problems. Several computational geometry algorithms are used to construct approximate trajectories in these problems, in particular, algorithms for constructing the convex hull as well as the union, intersection, and algebraic sum of polygons.
暂无评论