A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

A posteriori H1 error estimation for a singularly perturbed reaction diffusion problem on anisotropic meshes

The paper deals with a singularly perturbed reaction diffusionmodel problem. The focus is on reliable a posteriori error estimatorsfor the H¹ seminorm that can be applied to anisotropic finiteelement meshes. A residual error estimator and a local problemerror estimator are proposed and rigorously analysed. They arelocally equivalent, and both bound the error reliably. Threemodifications of these estimators are introduced and discussed. Much attention is given to the performance of the error estimatorin numerical experiments. This helps to identify those estimatorsthat are suitable for practical applications.     

Axisymmetric frictionless contact of functionally graded materials

The main interest of this study is a new method to solve the axisymmetric frictionless contact problem of functionally graded materials (FGMs). Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into a series of sub-layers with shear modulus varying linearly in each sub-layer and continuous at the sub-interfaces. With this model, the axisymmetric frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Hankel integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically. An erratum to this article can be found at     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

On a problem of Zygmund

It is proved that there exists an integrable function on [0, 1]² whose integral is nondifferentiable in each direction belonging to a set everywhere dense in [0, 2π] but is strongly differentiable. Translated fromMatematicheskie Zametki, Vol. 64, No. 5, pp. 749–762, November, 1998.     

Minimum energy configurations on a toric lattice as a quadratic assignment problem

We consider three known bounds for the quadratic assignment problem (QAP): an eigenvalue, a convex quadratic programming (CQP), and a semidefinite programming (SDP) bound. Since the last two bounds were not compared directly before, we prove that the SDP bound is stronger than the CQP bound. We then apply these to improve known bounds on a discrete energy minimization problem, reformulated as a QAP, which aims to minimize the potential energy between repulsive particles on a toric grid. Thus we are able to prove optimality for several configurations of particles and grid sizes, complementing earlier results by Bouman et al. (2013). The semidefinite programs in question are too large to solve without pre-processing, and we use a symmetry reduction method by Permenter and Parrilo (2020) to make computation of the SDP bounds possible.     

Necessary and sufficient conditions for the existence of positive solutions to algebraic Riccati equations with indefinite quadratic term

Necessary and sufficient conditions are established in this paper for the existence of positive- and/or negative-definite solutions to the algebraic Riccati equation with indefinite coefficient. An iterative procedure is also given for computing such a solution.Project supported by the National Science Foundation of China and by the special program of the State Education Commission of China under grant 9033507.     

Picture Fuzzy Threshold Graphs with Application in Medicine Replenishment

In this study, a novel concept of picture fuzzy threshold graph (PFTG) is introduced. It has been shown that PFTGs are free from alternating 4-cycle and it can be constructed by repeatedly adding a dominating or an isolated node. Several properties about PFTGs are discussed and obtained the results that every picture fuzzy graph (PFG) is equivalent to a PFTG under certain conditions. Also, the underlying crisp graph (UCG) of PFTG is a split graph (SG), and conversely, a given SG can be applied to constitute a PFTG. A PFTG can be decomposed in a unique way and it generates three distinct fuzzy threshold graphs (FTGs). Furthermore, two important parameters i.e., picture fuzzy (PF) threshold dimension (TD) and PF partition number (PN) of PFGs are defined. Several properties on TD and PN have also been discussed. Lastly, an application of these developed results are presented in controlling medicine resources.     

Exact solution of plane isolated crack normal to a bimaterial interface of infinite extent

This paper presents an exact solution for the transverse interface crack in the plane strain case. The crack is perpendicular to the interface and in one material. The exact complex stress functions are first obtained with some unknown constants. The satisfactions of all boundary conditions are then checked, the condition at infinity is considered and the unknown constants are determined. Further study may focus on the case with different shear moduli and the influence of the large deformation.The English text was polished by Keren Wang.     

A Bregman-Split-Based Compressive Sensing Method for Dynamic Harmonic Estimation

In order to overcome the spectral interference of the conventional Fourier transform in the International Electrotechnical Commission framework, this paper introduces a Bregman-split-based compressive sensing (BSCS) method to estimate the Taylor–Fourier coefficients in a multi-frequency dynamic phasor model. Considering the DDC component estimation, this paper transforms the phasor problem into a compressive sensing model based on the regularity and sparsity of the dynamic harmonic signal distribution. It then derives an optimized hybrid regularization algorithm with the Bregman split method to reconstruct the dynamic phasor estimation. The accuracy of the model was verified by using the cross entropy to measure the distribution differences of values. Composite tests derived from the dynamic phasor test conditions were then used to verify the potentialities of the BSCS method. Simulation results show that the algorithm can alleviate the impact of dynamic signals on phasor estimation and significantly improve the estimation accuracy, which provides a theoretical basis for P-class phasor measurement units (PMUs).