Impacting and friction are always encountered in the mechanical system and make the system experience nonlinearities in both force-position and force-velocity characteristics. In this paper, we propose one effective n...
详细信息
Impacting and friction are always encountered in the mechanical system and make the system experience nonlinearities in both force-position and force-velocity characteristics. In this paper, we propose one effective numerical algorithm for a two-dimensional frictional system involving the impacting by employing the closed form solution to yield the slip responses and extend the closed form solutions for curve length to predict the impacting location and time with high accuracy. By performing the accuracy analysis and robustness analysis through several representative scenarios, we reveal that, (1) the impacting dynamic responses can be predicted through our new algorithm with high accuracy, (2) the accumulative errors of dynamic responses are insensitive to the complex boundary, (3) the new algorithm can keep better robustness when involving repeated impacting.
The object of this article is to present the computational solution of the time-fractional Schrodinger equation subject to given constraint condition based on the generalized Taylor series formula in the Caputo sense....
详细信息
The object of this article is to present the computational solution of the time-fractional Schrodinger equation subject to given constraint condition based on the generalized Taylor series formula in the Caputo sense. The algorithm methodology is based on construct a multiple fractional power series solution in the form of a rabidly convergent series with minimum size of calculations using symbolic computation software. The proposed technique is fully compatible with the complexity of this problem and obtained results are highly encouraging. Efficacious computational experiments are provided to guarantee the procedure and to illustrate the theoretical statements of the present algorithm in order to show its potentiality, generality, and superiority for solving such fractional equation. Graphical results and numerical comparisons are presented and discussed quantitatively to illustrate the solution.
For several applications, it is important to know the location of the singularities of a complex function: just for example, the rightmost singularity of a Laplace Transform is related to the exponential order of its ...
详细信息
For several applications, it is important to know the location of the singularities of a complex function: just for example, the rightmost singularity of a Laplace Transform is related to the exponential order of its inverse function. We discuss a numerical method to approximate, within an input accuracy tolerance, a finite sequence of Laurent coefficients of a function by means of the Discrete Fourier Transform (DFT) of its samples along an input circle. The circle may also enclose some singularities, since the method works with the Laurent expansion. The DFT is computed by the FFT algorithm so that, from a computational point of view, the efficiency is guaranteed. The function samples may be obtained by solving a numerical problem such as, for example, a differential problem. We derive, as consequences of the method, some new outcomes able to detect those singularities which are close to the circle and to discover if the singularities are all external or internal to the circle so that the Laurent expansion reduces to its regular or singular part, respectively. Other singularities may be located by means of a repeated application of the method, as well as an analytic continuation. Some examples and results, obtained by a first implementation, are reported.
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 performances. To solve this proble...
详细信息
ISBN:
(纸本)9781509046577
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 performances. To 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 reported. These benchmark problems and the solutions are constituted as a framework, and the numerical algorithms for solving the FODEs could be compared in the same framework. The comparing result could assess which algorithm is better to a concrete FODE.
Demands for high-power converters are increasing rapidly. Nevertheless, due to their complex structures and multi-time-scale features, accurate and efficient simulation of such systems has long been a huge obstacle. A...
详细信息
ISBN:
(纸本)9781728103952
Demands for high-power converters are increasing rapidly. Nevertheless, due to their complex structures and multi-time-scale features, accurate and efficient simulation of such systems has long been a huge obstacle. Aiming at improving the unacceptably slow simulation speed of high-power converters, this paper improves Discrete State Event-Driven (DSED) approach by implementing a decoupling method. Taking a 10kV 1MW electric energy router as an example, this paper evaluates the effectiveness of DSED with the proposed decoupling method in high-power applications. In the studied case, compared with commercially available tool for power electronics simulation, PLECS (R), DSED can achieve an up to 1000-fold improvement in simulation speed under the same accuracy: PLECS (R) takes 6 hours while DSED only takes 18 seconds. DSED provides a powerful simulation tool for analysis, design and control of high-power converters.
The periodic QZ algorithm (pQZ) is required in many applications, including periodic linear systems, cyclic matrices and matrix pencils, and some structured eigenproblems. The implicit double pQZ step is the essential...
详细信息
ISBN:
(纸本)9781728106991
The periodic QZ algorithm (pQZ) is required in many applications, including periodic linear systems, cyclic matrices and matrix pencils, and some structured eigenproblems. The implicit double pQZ step is the essential part of this algorithm, since it determines its convergence rate. The pQZ step acts on a formal matrix product already reduced to the Hessenberg-triangular form. Each pQZ step starts with an initial transformation aimed to make the first column of the identity matrix and of the Wilkinson double-shift polynomial parallel. This transformation is found in an implicit manner, without evaluating the formal matrix product and the polynomial. This paper presents in detail the computation of the initial transformation in the periodic QZ step, and summarizes the benefits obtained using a new such approach for solving skew-Hamiltonian/Hamiltonian eigenproblems. The previous convergence failures are avoided and the number of iterations necessary to converge can be significantly reduced.
Vibration problems of pipelines made of composite materials conveying pulsating flow of gas and fluid are investigated in the paper. A dynamic model of motion of pipelines conveying pulsating fluid flow supported by a...
详细信息
Vibration problems of pipelines made of composite materials conveying pulsating flow of gas and fluid are investigated in the paper. A dynamic model of motion of pipelines conveying pulsating fluid flow supported by a Hetenyi's base is developed taking into account the viscosity properties of the structure material, axial forces, internal pressure and Winkler's viscoelastic base. To describe the processes of viscoelastic material strain, the Boltzmann-Volterra integral model with weakly singular hereditary kernels is used. Using the Bubnov-Galerkin method, the problem is reduced to the study of a system of ordinary integro-differential equations (IDE). A computational algorithm is developed based on the elimination of the features of IDE with weakly singular kernels, followed by the use of quadrature formulas. The effect of rheological parameters of the pipeline material, flow rate and base parameters on the vibration of a viscoelastic pipeline conveying pulsating fluid is analyzed. The convergence analysis of the approximate solution of the Bubnov-Galerkin method is carried out. It was revealed that the viscosity parameters of the material and the pipeline base lead to a significant change in the critical flow rate. It was stated that an increase in excitation coefficient of pulsating flow and the parameter of internal pressure leads to a decrease in the critical flow rate. It is shown that an increase in the singularity parameter, the Winkler base parameter, the rigidity parameter of the continuous base layer and the Reynolds number increases the critical flow rate.
Simulating the steam Injection process of pressurized decompression tanks is a challenging engineering problem due to the large amount of computational resource requirements and its complex physical model. This study ...
详细信息
ISBN:
(纸本)9781728109459
Simulating the steam Injection process of pressurized decompression tanks is a challenging engineering problem due to the large amount of computational resource requirements and its complex physical model. This study is a necessary algorithmic study to perform a 3D numerical simulation from the PRT (Pressurizer Relief Tank) to the transient dumping process in a chamber containing a tank. Parallel simulation was performed by the open source CFD tool OpenFOAM based on the numerical algorithm we discussed. The results show that the simulation results are in accordance with the basic physical laws, and the results reflect the pressure and temperature of the entire pressure relief tank in detail, which can provide guidance for the safe design of the reactor coolant system.
The next generation of the public domain Storm Water Management Model (SWMM) code needs to be able to effectively use parallelization of 1000 + threads for high-resolution simulation of large urban drainage systems. W...
详细信息
ISBN:
(纸本)9783319998671;9783319998664
The next generation of the public domain Storm Water Management Model (SWMM) code needs to be able to effectively use parallelization of 1000 + threads for high-resolution simulation of large urban drainage systems. We propose a new explicit solution technique using a "no-neighbour" finite-volume algorithm. A no-neighbour algorithm is strictly limited in the design of its time/space discretization stencil to limit the communication bandwidth issues that typically cause codes to lose efficiency when applied over large numbers of parallel threads. Returning SWMM to an explicit solver creates a number of challenges for implementing the solution technique, but also raises some new possibilities for how the model can be used. In this work, we briefly describe the proposed algorithm and examine its performance in simulating a steep urban stream.
This paper is devoted to the numerical scheme for the impulsive differential equations. The main idea of this method is, for the first time, to establish a broken reproducing kernel space that can be used in pulse mod...
详细信息
This paper is devoted to the numerical scheme for the impulsive differential equations. The main idea of this method is, for the first time, to establish a broken reproducing kernel space that can be used in pulse models. Then the uniform convergence of the numerical solution is proved, and the time consuming Schmidt orthogonalization process is avoided. The proposed method is proved to be stable and have the second-order convergence. The algorithm is proved to be feasible and effective through some numerical examples.
暂无评论