The current scientific work examines a pseudo-plastic Williamson fluid over a homogenous porous media. Momentum and energy equations refer to the executive branch partial differential equations (PDE). A new transforma...
详细信息
Focusing on hybrid diffusions in which continuous dynamics and discrete events coexist, this work is concerned with approximation of solutions for hybrid stochastic differential equations with a state-dependent switch...
详细信息
Focusing on hybrid diffusions in which continuous dynamics and discrete events coexist, this work is concerned with approximation of solutions for hybrid stochastic differential equations with a state-dependent switching process. Iterative algorithms are developed. The continuous-state-dependent switching process presents added difficulties in analyzing the numerical procedures. Weak convergence of the algorithms is established by a martingale problem formulation first. This weak convergence result is then used as a bridge to obtain strong convergence. In this process, the existence and uniqueness of the solution of the switching diffusions with continuous-state-dependent switching are obtained. In contrast to existing results of solutions of stochastic differential equations in which the Picard iterations are utilized, Euler's numerical schemes are considered here. Moreover, decreasing-stepsize algorithms together with their weak convergence are given. numerical experiments are also provided for demonstration.
For k + 1 power series a(0)(z),...,a(k)(z), we present a new iterative, look-ahead algorithm for numerically computing Pade-Hermite systems and simultaneous Pade systems along a diagonal of the associated Pade tables....
详细信息
For k + 1 power series a(0)(z),...,a(k)(z), we present a new iterative, look-ahead algorithm for numerically computing Pade-Hermite systems and simultaneous Pade systems along a diagonal of the associated Pade tables. The algorithm computes the systems at all those points along the diagonal at which the associated striped Sylvester and mosaic Sylvester matrices are well conditioned. The operation and the stability of the algorithm is controlled by a single parameter tau which serves as a threshold in deciding if the Sylvester matrices at a point are sufficiently well conditioned. We show that the algorithm is weakly stable and provide bounds for the error in the computed solutions as a function of tau. Experimental results are given which show that the bounds reflect the actual behavior of the error. The algorithm requires O(\\n\\(2) + s(3)\\n\\) Operations to compute Pade-Hermite and simultaneous Fade systems of type n = [n(0),...,n(k)], where \\n\\ = n(0)+...+n(k) and s is the largest step-size taken along the diagonal. An additional application of the algorithm is the stable inversion of striped and mosaic Sylvester matrices.
The complex technical object usually can be described as a set of technological components (blocks). Its functional diagram can be represented by an oriented multi-graph. The vertices of the graph correspond to separa...
详细信息
ISBN:
(纸本)9783031734168;9783031734175
The complex technical object usually can be described as a set of technological components (blocks). Its functional diagram can be represented by an oriented multi-graph. The vertices of the graph correspond to separate components (blocks) of a technical object. The blocks are characterized by design parameters. The graph may contain feedback loops. Blocks' output parameters can be determined from the values of their input parameters. The software developed in C++ language serves for calculating functional diagrams of steam boilers. The software determines also the values of derived output parameters from the design parameters values. It can be used to improve the technical characteristics of an object and determine control parameters for various modes of its operation. The calculation of the thermal scheme is based on the use of a model heat exchanger. In this case, the calculation of the thermal scheme is reduced to an iterative solution of a system of linear equations. A brief description of the software is provided.
A new method is developed to calculate characteristics of contaminant transport (including non-classical regimes) in statistically homogeneous sharply contrasting media. A transport integro-differential equation in th...
详细信息
A new method is developed to calculate characteristics of contaminant transport (including non-classical regimes) in statistically homogeneous sharply contrasting media. A transport integro-differential equation in the space-time representation is formulated on the basis of the model earlier proposed by one of the authors (L. M.). Analytical expressions for transport characteristics in limiting time intervals in the one-dimensional case are derived. An interpolation form is proposed for the integral kernel of the transport equation. On a basis of this expression, an algorithm is developed for numerical modelling the contaminant transport in statistically homogeneous sharply contrasting media. Trial numerical 1D calculations are performed based on this algorithm. Good agreement was found between the numerical simulation results and the asymptotic analytical expressions.
Lithography simulators have been playing an indispensable role in process optimization and design for manufacturability (DFM). The ever smaller feature sizes demand higher numerical accuracy and faster runtime on thes...
详细信息
Lithography simulators have been playing an indispensable role in process optimization and design for manufacturability (DFM). The ever smaller feature sizes demand higher numerical accuracy and faster runtime on these lithography simulators. Aerial image simulation is the first key step in lithography simulation, and the method using transmission cross coefficient (TCC), which is a two-dimensional integral, is the most commonly used technique for full-chip aerial image simulation. In this paper, we present a very accurate, yet efficient and extensible aerial image simulator, ELIAS. We find that the majority of the numerical error during the TCC computation is due to the discontinuous boundaries of the support of the TCC integrand. We reduce the error dramatically by using a recursive integration algorithm. Because TCC is usually computed on uniform grids, we further speed up the algorithm without increasing the errors. Given the same accuracy, our new algorithm can speed up the runtime by 100 x -1000 X. Our algorithm also provides smooth tradeoff between accuracy and runtime. It can be used to benchmark other lithography aerial simulators. In addition, ELIAS provides an open-source, flexible software framework to incorporate different lithography settings.
A linear two-point boundary value problem for a system of loaded differential equations with impulse effect is investigated. Values in the previous impulse points are taken into consideration in the conditions of impu...
详细信息
A linear two-point boundary value problem for a system of loaded differential equations with impulse effect is investigated. Values in the previous impulse points are taken into consideration in the conditions of impulse effect. The considered problem is reduced to an equivalent multi-point boundary value problem for the system of ordinary differential equations with parameters. A numerical implementation of parametrization method is offered using the Runge-Kutta method of 4th-order accuracy for solving the Cauchy problems for ordinary differential equations. The constructed numerical algorithms are illustrated by examples.
Fragmentation is a well-known method of the parallelization of numerical algorithms and programs. algorithm fragmentation allows creating fragmented parallel programs that can be executed on parallel computers of diff...
详细信息
Fragmentation is a well-known method of the parallelization of numerical algorithms and programs. algorithm fragmentation allows creating fragmented parallel programs that can be executed on parallel computers of different types (multiprocessors and/or multicomputers) and can be dynamically tuned to all the available resources. Fragmentation of the often used numerical algorithms, their representation for inclusion into the library of parallel numerical subroutines and properties of the runtime system are considered.
Plasma dispersion function is an important parameter in the ionospheric physics, fast and accurate computation of this function is extremely valuable in practical use. Many numerical algorithms have been proposed, and...
详细信息
Plasma dispersion function is an important parameter in the ionospheric physics, fast and accurate computation of this function is extremely valuable in practical use. Many numerical algorithms have been proposed, and it is very necessary to evaluate their performances. In present paper, we first introduce an alternative method to derive the plasma dispersion function and its first derivative using the differential and integral calculus method, and then compare the accuracy and efficiency of three well-known numerical algorithms (algorithm of Steven G. Johnson, algorithm 916, and algorithm of Abrarov and Quine) in Matlab environment under the uniform and non-uniform distribution of the grid-points, and finally analyze the variations of the real and imaginary part of the plasma dispersion function and its first derivative preliminarily. The results show that Abrarov and Quine's algorithm performs better than other two numerical algorithms from the comprehensive viewpoint of accuracy and efficiency. (C) 2019 Elsevier Ltd. All rights reserved.
An antenna array composed of one transmitter and multi-receivers and dedicated to measuring terminal trajectory of the target of interest, which is supposed to be in uniform rectilinear motion, is set up. On the basis...
详细信息
ISBN:
(纸本)9783540741701
An antenna array composed of one transmitter and multi-receivers and dedicated to measuring terminal trajectory of the target of interest, which is supposed to be in uniform rectilinear motion, is set up. On the basis of the model, the Vector parameter that can uniquely determine the terminal trajectory of target is introduced, and the measurement equations which describe the respective relationships between the Vector parameter and the instantaneous Doppler frequency and the phase differences extracted from echoes are established. Taking advantage of the measurement equations, we propose an algorithm of estimating the Vector parameter without resolving the phase difference ambiguity;furthermore, the detailed steps in estimating the Vector parameter using numerical optimization techniques are put forward. The Monte Carlo simulation results demonstrate the effectiveness and reliability of our numerical algorithm compared with the traditional method.
暂无评论