Structured condition numbers and statistical condition estimation for the LDU factorization

In this article, we consider the structured condition numbers for LDU, factorization by using the modified matrix-vector approach and the differential calculus, which can be represented by sets of parameters. By setting the specific norms and weight parameters, we present the expressions of the structured normwise, mixed, componentwise condition numbers and the corresponding results for unstructured ones. In addition, we investigate the statistical estimation of condition numbers of LDU factorization using the probabilistic spectral norm estimator and the small-sample statistical condition estimation method, and devise three algorithms. Finally, we compare the structured condition numbers with the corresponding unstructured ones in numerical experiments.

On the Preconditioning Properties of RHSS Preconditioner for Saddle-Point Linear Systems

In this paper,for the regularized Hermitian and skew-Hermitian splitting(RHSS)preconditioner introduced by Bai and Benzi(BIT Numer Math 57:287–311,2017)for the solution of saddle-point linear systems,we analyze the spectral properties of the preconditioned matrix when the regularization matrix is a special Hermitian positive semidefinite matrix which depends on certain parameters.We accurately describe the numbers of eigenvalues clustered at(0,0)and(2,0),if the iteration parameter is close to 0.An estimate about the condition number of the corresponding eigenvector matrix,which partly determines the convergence rate of the RHSS-preconditioned Krylov subspace method,is also studied in this work.

On Strong Law of Large Numbers for Random Sequence

这笔记被奉献介绍有条件地统治的随机的变量的一个新概念。在下面合适限制条件，为任意的连续随机的变量的大数字的一条一般强壮的法律被获得。

2D Minimum Compliance Topology Optimization Based on a Region Partitioning Strategy

This paper presents an extended sequential element rejection and admission(SERA)topology optimizationmethod with a region partitioning strategy.Based on the partitioning of a design domain into solid regions and weak regions,the proposed optimizationmethod sequentially implements finite element analysis(FEA)in these regions.After standard FEA in the solid regions,the boundary displacement of the weak regions is constrained using the numerical solution of the solid regions as Dirichlet boundary conditions.This treatment can alleviate the negative effect of the material interpolation model of the topology optimization method in the weak regions,such as the condition number of the structural global stiffness matrix.For optimization,in which the forward problem requires nonlinear structural analysis,a linear solver can be applied in weak regions to avoid numerical singularities caused by the over-deformedmesh.To enhance the robustness of the proposedmethod,the nonmanifold point and island are identified and handled separately.The performance of the proposed method is verified by three 2D minimum compliance examples.

On Iterative Algorithm and Perturbation Analysis for the Nonlinear Matrix Equation

In this study,an iterative algorithm is proposed to solve the nonlinear matrix equation X+A∗eXA=In.Explicit expressions for mixed and componentwise condition numbers with their upper bounds are derived to measure the sensitivity of the considered nonlinear matrix equation.Comparative analysis for the derived condition numbers and the proposed algorithm are presented.The proposed iterative algorithm reduces the number of iterations significantly when incorporated with exact line searches.Componentwise condition number seems more reliable to detect the sensitivity of the considered equation than mixed condition number as validated by numerical examples.

A complete solution of an improved universal 3D coordinate similarity transformation model

When linearizing three-dimensional(3 D)coordinate similarity transformation model with large rotations,we usually encounter the ill-posed normal matrix which may aggravate the instability of solutions.To alleviate the problem,a series of conversions are contributed to the 3 D coordinate similarity transformation model in this paper.We deduced a complete solution for the 3 D coordinate similarity transformation at any rotation with the nonlinear adjustment methodology,which involves the errors of the common and the non-common points.Furthermore,as the large condition number of the normal matrix resulted in an intractable form,we introduced the bary-centralization technique and a surrogate process for deterministic element of the normal matrix,and proved its benefit for alleviating the condition number.The experimental results show that our approach can obtain the smaller condition number to stabilize the convergence of the interested parameters.Especially,our approach can be implemented for considering the errors of the common and the non-common points,thus the accuracy of the transformed coordinates improves.

Domain Decomposition for Wavelet Single Layer on Geometries with Patches

We focus on the single layer formulation which provides an integral equation of the first kind that is very badly conditioned. The condition number of the unpreconditioned system increases exponentially with the multiscale levels. A remedy utilizing overlapping domain decompositions applied to the Boundary Element Method by means of wavelets is examined. The width of the overlapping of the subdomains plays an important role in the estimation of the eigenvalues as well as the condition number of the additive domain decomposition operator. We examine the convergence analysis of the domain decomposition method which depends on the wavelet levels and on the size of the subdomain overlaps. Our theoretical results related to the additive Schwarz method are corroborated by numerical outputs.

ON THE SENSITIVITY OF THE LDU FACTORIZATION

In this paper,we give the sensitivity analyses by two approaches for L,D,U in factorization A=LDU of general perturbations in A which sufficiently small in norm. By the matrix vector equation approach,we derive the sharp condition number for L,D and U factors .By the matrix equation approach we derive corresponding condition estimates. When A is a symmetric matrix,the corresponding results can be obtained for LDL T factorization.

新型旋叶式汽车空调压缩机受力分析

本文针对为汽车空调系统而开发的新型旋叶式压缩机的结构特点，建立了受力分析的数学模型。并且在此基础上开发出了新型旋叶式汽车空调压缩机分析软件。应用这一软件，详细分析了新型旋叶式汽车空调压缩机中各种作用力的变化规律及其对压缩机性能的影响。同时详细分析了旋叶式压缩机中叶片数这一重要参数对压缩机受力的影响。

弱不连续问题高阶有限元离散系统的GAMG法

弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。

Re ned rigorous perturbation bounds for the SR decomposition

In this article,some new rigorous perturbation bounds for the SR decomposition un-der normwise or componentwise perturbations for a given matrix are derived.Also,the explicit expressions for the mixed and componentwise condition numbers are presented by utilizing the block matrix-vector equation approach.Hypothetical and trial results demonstrate that these new bounds are constantly more tightly than the comparing ones in the literature.

Perturbation Analysis for the Matrix-Scaled Total Least Squares Problem

In this paper, we extend matrix scaled total least squares (MSTLS) problem with a single right-hand side to the case of multiple right-hand sides. Firstly, under some mild conditions, this paper gives an explicit expression of the minimum norm solution of MSTLS problem with multiple right-hand sides. Then, we present the Kronecker-product-based formulae for the normwise, mixed and componentwise condition numbers of the MSTLS problem. For easy estimation, we also exhibit Kronecker-product-free upper bounds for these condition numbers. All these results can reduce to those of the total least squares (TLS) problem which were given by Zheng et al. Finally, two numerical experiments are performed to illustrate our results.

窗式空调器平行流冷凝器流程和扁管的优化分析

选择百叶窗肋片及小通道内流体传热流动的经验关联式,采用效能-传热单元数法建立了平行流冷凝器的稳态仿真模型。对不同流程数和扁管数方案的冷凝热和制冷剂压降进行了模拟分析,得出了扁管总数为39时的优化方案。

一类数列收敛的充分条件

分析了用于定义Euler-Mascheroni常数的数列，并在此基础上给出了其收敛的两个充分条件而非必要条件，这些条件较已有的条件更为简洁。

5mm管房间空调器流路的优化设计

将目前房间空调系统中室外机和室内机常用的直径7mm及以上的内螺纹铜管改为5mm内螺纹管,可以达到节省材料与提高传热效率的效果。本文提出了空调系统换热管由7mm及以上铜管改为5mm铜管的设计方法和改进措施,如调整换热器纵向和横向换热管数目,制冷剂流路的分路数和压缩机更换等。通过一个分体式房间空调系统设计算例,具体说明将空调系统由7mm及以上管径改为5mm管径的具体的改进措施。设计结果表明在换热器具有相同换热性能的条件下,室外机由9mm管改为5mm管时,需同时增加换热器纵向和横向换热管数目,并调整分路数;室内机由7mm管改为5mm管室内机时,仅需要改变分路数。同时,采用高效压缩机,可进一步提高5mm管空调系统的性能。

A Well-Conditioned, Nonconforming Nitsche’s Extended Finite Element Method for Elliptic Interface Problems

In this paper,we introduce a nonconforming Nitsche’s extended finite element method(NXFEM)for elliptic interface problems on unfitted triangulation elements.The solution on each side of the interface is separately expanded in the standard nonconforming piecewise linear polynomials with the edge averages as degrees of freedom.The jump conditions on the interface and the discontinuities on the cut edges(the segment of edges cut by the interface)are weakly enforced by the Nitsche’s approach.In the method,the harmonic weighted fluxes are used and the extra stabilization terms on the interface edges and cut edges are added to guarantee the stability and the well conditioning.We prove that the convergence order of the errors in energy and L 2 norms are optimal.Moreover,the errors are independent of the position of the interface relative to the mesh and the ratio of the discontinuous coefficients.Furthermore,we prove that the condition number of the system matrix is independent of the interface position.Numerical examples are given to confirm the theoretical results.

WELL-CONDITIONED FRAMES FOR HIGH ORDER FINITE ELEMENT METHODS

The purpose of this paper is to discuss representations of high order C^(0)finite element spaces on simplicial meshes in any dimension.When computing with high order piecewise polynomials the conditioning of the basis is likely to be important.The main result of this paper is a construction of representations by frames such that the associated L^(2)condition number is bounded independently of the polynomial degree.To our knowledge,such a representation has not been presented earlier.The main tools we will use for the construction is the bubble transform,introduced previously in[1],and properties of Jacobi polynomials on simplexes in higher dimensions.We also include a brief discussion of preconditioned iterative methods for the finite element systems in the setting of representations by frames.

Spectral Matrix Conditioning in Cylindrical and Spherical Elliptic Equations

In the spectral solution of 3-D Poisson equations in cylindrical and spherical coordinates including the axis or the center,it is convenient to employ radial basis functions that depend on the Fourier wavenumber or on the latitudinal mode.This idea has been adopted by Matsushima and Marcus and by Verkley for planar problems and pursued by the present authors for spherical ones.For the Dirichlet boundary value problem in both geometries,original bases have been introduced built upon Jacobi polynomials which lead to a purely diagonal representation of the radial second-order differential operator of all spectral modes.This note details the origin of such a diagonalization which extends to cylindrical and spherical regions the properties of the Legendre basis introduced by Jie Shen for Cartesian domains.Closed form expressions are derived for the diagonal elements of the stiffness matrices as well as for the elements of the tridiagonal mass matrices occurring in evolutionary problems.Furthermore,the bound on the condition number of the spectral matrices associated with the Helmholtz equation are determined,proving in a rigorous way one of the main advantages of the proposed radial bases.

Restricted Additive Schwarz Preconditioner for Elliptic Equations with Jump Coefficients

This paper provides a proof of robustness of the restricted additive Schwarz preconditioner with harmonic overlap(RASHO)for the second order elliptic problems with jump coefficients.By analyzing the eigenvalue distribution of the RASHO preconditioner,we prove that the convergence rate of preconditioned conjugate gradient method with RASHO preconditioner is uniform with respect to the large jump and meshsize.

On hp refinements of independent cover numerical manifold method——some strategies and observations

In this paper,strategies are provided for a powerful numerical manifold method(NMM)with h and p refinement in analyses of elasticity and elasto-plasticity.The new NMM is based on the concept of the independent cover,which gets rid of NMM's important defect of rank deficiency when using higher-order local approximation functions.Several techniques are presented.In terms of mesh generation,a relationship between the quadtree structure and the mathematical mesh is established to allow a robust h-refinement.As to the condition number,a scaling based on the physical patch is much better than the classical scaling based on the mathematical patch;an overlapping width of 1%–10%can ensure a good condition number for 2nd,3rd,and 4th order local approximation functions;the small element issue can be overcome after the local approximation on small patch is replaced by that on a regular patch.On numerical accuracy,local approximation using complete polynomials is necessary for the optimal convergence rate.Two issues that may damage the convergence rate should be prevented.The first is to approximate the curved boundary of a higher-order element by overly few straight lines,and the second is excessive overlapping width.Finally,several refinement strategies are verified by numerical examples.