On the numerical solution of ill‐conditioned linear systems by regularization and iteration

We propose to reduce the (spectral) condition number of a given linear system by adding a suitable diagonal matrix to the system matrix, in particular by shifting its spectrum. Iterative procedures are then adopted to recover the solution of the original system. The case of real symmetric positive definite matrices is considered in particular, and several numerical examples are given. This approach has some close relations with Riley's method and with Tikhonov regularization. Moreover, we identify approximately the aforementioned procedure with a true action of preconditioning.     

Development of an improved divergence‐free‐condition compensated coupled framework to solve flow problems with time‐varying geometries

In this article a coupled version of the improved divergence‐free‐condition compensated method will be proposed to simulate time‐varying geometries by direct forcing immersed boundary method. The proposed method can be seen as a quasi‐multi‐moment framework due to the fact that the momentum equations are discretized by both cell‐centered and cell‐face velocity. For simulating time‐varying geometries, a semi‐implicit iterative method is proposed for calculating the direct forcing terms. Treatments for suppressing spurious force oscillations, calculating drag/lift forces, and evaluating velocity and pressure for freshly cells will also be addressed. In order to show the applicability and accuracy, analytical as well as benchmark problems will be investigated by the present framework and compared with other numerical and experimental results.     

模糊选择函数半序合理的充分条件

在G?del t-模下,研究了模糊选择函数的半序合理性.首先给出了模糊选择函数的合理性条件FA1.然后研究了该条件与模糊选择函数半序合理性之间的关系,得到了半序合理的一个充分条件.     

Imposing accurate wall boundary conditions in corrective‐matrix‐based moving particle semi‐implicit method for free surface flow

Corrective matrix that is derived to restore consistency of discretization schemes can significantly enhance accuracy for the inside particles in the Moving Particle Semi‐implicit method. In this situation, the error due to free surface and wall boundaries becomes dominant. Based on the recent study on Neumann boundary condition (Matsunaga et al, CMAME, 2020), the corrective matrix schemes in MPS are generalized to straightforwardly and accurately impose Neumann boundary condition. However, the new schemes can still easily trigger instability at free surface because of the biased error caused by the incomplete/biased neighbor support. Therefore, the existing stable schemes based on virtual particles and conservative gradient models are applied to free surface and nearby particles to produce a stable transitional layer at free surface. The new corrective matrix schemes are only applied to the particles under the stable transitional layer for improving the wall boundary conditions. Three numerical examples of free surface flows demonstrate that the proposed method can help to reduce the pressure/velocity fluctuations and hence enhance accuracy further.     

Optimization of Transducer Location for Novel Non-Intrusive Methodologies of Diagnosis in Diesel Engines

The health monitoring has been studied to ensure integrity of design of engine structure by detection, quantification, and prediction of damages. Early detection of faults may allow the downtime of maintenance to be rescheduled, thus preventing sudden shutdown of machines. In cylinder pressure developed, vibrations and noise emissions data provide a rich source of information about condition of engines. Monitoring of vibrations and noise emissions are novel non-intrusive methodologies for which positioning of various transducers are important issue. The presented work shows applicability of these diagnosis methodologies adopted in case of diesel engines. The effects of changing various fuel injection parameters was analyzed. Scope of using non-intrusive technique has been analyzed by changing locations of microphone. Novelty of this worklies in exploring signal processing methods for various locations around the engine test set up. Various frequency ranges of contributing noise and vibration sources were identified. Time-Frequency analysis showed the onset of various cyclic. Based on the identification of various frequency bands, it is possible to device suitable filters in order to extract more information.     

Convergence rates of the blocked Gibbs sampler with random scan in the Wasserstein metric

ABSTRACT

This paper establishes explicit estimates of convergence rates for the blocked Gibbs sampler with random scan under the Dobrushin conditions. The estimates of convergence in the Wasserstein metric are obtained by taking purely analytic approaches.     

核电结构土-结相互作用分析分区混合计算方法

土-结构相互作用分析是核电结构进行抗震设计和安全评估的重要环节.在核电结构的土-结相互作用分析中,阻尼和非线性是影响结构反应的重要因素. 若采用频域分析,可以方便考虑阻尼,但需通过等效线性化来考虑非线性,不适合于强震作用下的土体非线性.若采用时域分析的逐步积分方法,适合于考虑非线性,但材料阻尼一般采用瑞利阻尼模型,除了紧靠指定阻尼比的少数几个振型外,其他振型的反应将受到瑞利阻尼模型所确定的大阻尼所抑制,造成地震反应与真实情形有较大差异.若采用时域分析的模态叠加法,可合理计入阻尼效应,但模态叠加法不能考虑非线性.因此,如何合理考虑阻尼和非线性是核电结构土-结相互作用分析需要关注的问题.基于此,本文提出一种模态叠加和时步积分结合的土-结相互作用分区算法.其中,出于安全性考虑,地震作用下核电主体结构一般不允许进入非线性,因此结构可采用模态叠加方法,以便合理考虑结构阻尼;土体和基础采用显式时步积分法,可考虑土体非线性;通过人工边界条件考虑无限域的影响 (辐射阻尼).通过简单算例对该方法进行了验证,并用于CAP1400核电结构的土-结相互作用分析中,对比分析了采用模态阻尼和瑞利阻尼时核电结构和场地反应的差异,结果表明结构阻尼模型对场地的反应影响不大,但对结构反应影响明显,在实际工程中应合理选取阻尼模型.      

基于Hopf-Cole变换的Burgers方程初边值问题的Sinc-Galerkin法

利用Sinc-Galerkin法数值求解Burgers方程的初边值问题。首先,用Hopf-Cole变换将二阶非线性的Burgers方程变换为二阶线性方程,同时把第一类边界条件变为第二类边界条件。时间上的导数采用θ加权格式离散,空间导数采用Sinc-Galerkin法离散,端点处分别引入权函数处理变换后的第二类边界条件。最后,通过数值算例验证了Sinc-Galerkin法的指数收敛性,与精确解相比,本文构造的数值格式精度高,能够有效捕捉激波等物理现象。