Fast algorithms for placing large entries along the diagonal of a sparse matrix

Solving a sparse system of linear equations Ax=b is one of the most fundamental operations inside any circuit simulator. The equations/rows in the matrix A are often rearranged/permuted before factorization and applying direct or iterative methods to obtain the solution. Permuting the rows of the matrix A so that the entries with large absolute values lie on the diagonal has several advantages like better numerical stability for direct methods (e.g., Gaussian elimination) and faster convergence for indirect methods (such as the Jacobi method). Duff (2009) [3] has formulated this as a weighted bipartite matching problem (the MC64 algorithm). In this paper we improve the performance of the MC64 algorithm with a new labeling technique which improves the asymptotic complexity of updating dual variables from O(|V|+|E|) to O(|V|), where |V| is the order of the matrix A and |E| is the number of non-zeros. Experimental results from using the new algorithm, when benchmarked with both industry benchmarks and UFL sparse matrix collection, are very promising. Our algorithm is more than 60 times faster (than Duff’s algorithm) for sparse matrices with at least a million non-zeros.     

Peristaltic flow of viscoelastic fluid with fractional Maxwell model through a channel

The paper presents the transportation of viscoelastic fluid with fractional Maxwell model by peristalsis through a channel under long wavelength and low Reynolds number approximations. The propagation of wall of channel is taken as sinusoidal wave propagation (contraction and relaxation). Homotopy perturbation method (HPM) and Adomian decomposition method (ADM) are used to obtain the analytical approximate solutions of the problem. The expressions of axial velocity, volume flow rate and pressure gradient are obtained. The effects of fractional parameters (α), relaxation time (λ₁) and amplitude (?) on the pressure difference and friction force across one wavelength are calculated numerically for different particular cases and depicted through graphs.     

稀土焦绿石化合物R2Fe4/3W2/3O7的高温高压稳定性研究和X射线相分析

 本文利用X射线粉末衍射和位敏探测技术，研究了R₂Fe_4/3W_2/3O₇（R=Er、Yb、Dy）化合物经高温高压处理后的变化情况。在3.7 GPa，1 200 ℃条件下，六方相R₂Fe_4/3W_2/3O₇化合物按两种方式分解，而直接由R₂O₃，Fe₂O₃和WO₃原料出发，经上述同样的高温高压条件合成所得的产物与六方相高温高压分解产物相同，均为R₂WO₆、RFeO₃、WO₃和Fe₂O₃的多相聚合物。同时给出了R₂Fe_4/3W_2/3O₇六方相高温高压下的稳定区范围。     

具有时变滞后的随机系统的时滞依赖鲁棒稳定性与H_∞分析

研究具有时变滞后与Markov跳跃的随机系统的鲁棒随机稳定性与H_∞分析.通过引进随机Liapunov-Krasovskii泛函,使用自由权矩阵技术(即何技术),得到时滞依赖鲁棒随机稳定性与H_∞扰动衰减的线性矩阵不等式判据.3个数值例子说明提供的方法是有效的,并且与相关文献中存在的一些结果相比是低保守的.     

单函数(A,Γ)小波的存在性问题

小波分析是除了Fourier分析和Gabor分析以外一种新的时频分析工具．它被应用在信号处理、图像处理以及许多其他领域．小波分析的—个基本问题是什么样的(A,Γ)对,使得存在单函数(A,Γ)小波.本文填补了Ionascu,Yang Wang关于单函数小波存在性问题在二维情形论证中的漏洞.     

Toeplitz-Bezout矩阵的若干性质

从定义出发,利用矩阵生成函数的方法来研究Toeplitz-Bezout若干基本性质,同时利用极限的思想将其对角约化.     

Adaptive digital twins of combustion systems using sparse sensing strategies

This work proposes to implement a sparse sensing framework to build a hybrid numerical-experimental Digital Twin of a practical combustion system. The goal is to find the optimal sensor placement that minimizes the prediction error, and to predict the distribution of reacting scalars using few measurements. Three-dimensional CFD simulations with detailed chemistry were used to build the design space by varying the fuel composition (from pure methane to pure hydrogen), the equivalence ratio (from 0.7 to 1) and the air velocity. The Proper Orthogonal Decomposition (POD) was applied to the numerical data to find a tailored basis for dimensionality reduction. Then, the QR decomposition with column pivoting was applied to the tailored basis to find the optimal sensor placement. Finally, the model was employed to predict the three-dimensional temperature distribution in the unexplored part of the design space, using the experimental samples as input. The optimal placement of the sensors provides valuable information on the key locations and features, which can then be used in the design of reactor network models, for example. Also, the results show that the hybrid Digital Twin could predict an adjusted temperature distribution which reduces the error with the experimental measurements, when compared to the original CFD temperature distribution.     

Perturbation models

A Markov chain (with a discrete state space and a continuous parameter) is perturbed by forcing a chain to return to “permissible” states whenever it happens to enter “forbidden” states, with returns governed by a replacement distribution.The compensation method is employed to obtain the distribution for the modified chain, in terms of the original chain and the perturbation mechanism.Emphasis is placed on ergodic chains, and interpretation of results in terms of perturbation theory of semi-groups and the ergodic potential theory (based on the fundamental matrix of a chain) is mentioned.     

The generalized Kolmogorov-Petrovskii-Piskunov equation

Nonlinear nonlocal equations of mathematical physics such as the K.P.P. equation, the generalized nonlinear Schrödinger equation, the Witham equation for water waves et al. are solved by decomposition.     

Conductance and Statistical Properties of Chaotic and Integrable Electron Waveguides

We show that the S-matrix for electrons propagating in a waveguide has different statistical properties depending on whether the waveguide cavity shape gives rise to chaotic or integrable behavior classically. We obtain distributions of energy level spacings for integrable and chaotic billiards shaped like the waveguide cavity. We also obtain distributions for Wigner delay times and resonance widths for the waveguide, for integrable and chaotic cavity geometries. Our results, obtained by direct numerical calculation of the electron wave function, are consistent with the predictions of random matrix theory.