ANISOTROPIC POLARIZATION TENSORS FOR ELLIPSES AND ELLIPSOIDS

In this paper we present a systematic way of computing the polarization tensors, anisotropic as well as isotropic, based on the boundary integral method. We then use this method to compute the anisotropic polarization tensor for ellipses and ellipsoids. The computation reveals the pair of anisotropy and ellipses which produce the same polarization tensors.     

一种基于小波分析的各向异性图象去噪方法

本文利用非线性各向异性扩散方程结合小波变换提出一种图象去噪的方法。首先对图像进行离散小波变换,然后对其各个分量分别用各向异性的方法实现去噪。实验结果表明,该方法能够较好的去除噪声的同时,很好的保留边缘信息。     

A combined finite volume–finite element scheme for the discretization of strongly nonlinear convection–diffusion–reaction problems on nonmatching grids

We propose and analyze in this paper a numerical scheme for nonlinear degenerate parabolic convection–diffusion–reaction equations in two or three space dimensions. We discretize the time evolution, convection, reaction, and source terms on a given grid, which can be nonmatching and can contain nonconvex elements, by means of the cell‐centered finite volume method. To discretize the diffusion term, we construct a conforming simplicial mesh with the vertices given by the original grid and use the conforming piecewise linear finite element method. In this way, the scheme is fully consistent and the discrete solution is naturally continuous across the interfaces between the subdomains with nonmatching grids, without introducing any supplementary equations and unknowns or using any interpolation at the interfaces. We allow for general inhomogeneous and anisotropic diffusion–dispersion tensors, propose two variants corresponding respectively to arithmetic and harmonic averaging, and use the local Péclet upstream weighting in order to only add the minimal numerical diffusion necessary to avoid spurious oscillations in the convection‐dominated case. The scheme is robust, efficient since it leads to positive definite matrices and one unknown per element, locally conservative, and satisfies the discrete maximum principle under the conditions on the simplicial mesh and the diffusion tensor usual in the finite element method. We prove its convergence using a priori estimates and the Kolmogorov relative compactness theorem and illustrate its behavior on a numerical experiment. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Flow-induced anisotropic viscosity in short fiber reinforced polymers

The commonly used flow models for fiber reinforced polymers often neglect the flow induced mechanical anisotropy of the suspension. With an increasing fiber volume fraction, this plays, however, an important role. There are some models which count on this effect, they are, however, phenomenological and require a fitted model parameter. In this paper, a micromechanically based constitutive law is proposed which considers the flow induced anisotropic viscosity of the fiber suspension. The introduced viscosity tensor can handle arbitrary anisotropy of the fluid-fiber mixture depending on the actual fiber orientation distribution. A homogenization method for unidirectional structures in contribution with orientation averaging is used to determine the effective viscosity tensor. The motion of rigid ellipsoidal fibers induced by the flow of the matrix material is described by Jeffery's equation. A numerical implementation of the introduced model is applied to representative flow modes. The calculated stress values are analyzed in transient and stationary flow cases. They show a less pronounced anisotropic viscous behaviour in every investigated case compared to the results obtained by the use of the Dinh-Armstrong constitutive law. The reason for the qualitative difference is that the presented model depends on the complete orientation information, while the other one is linear in the fourth-order fiber orientation tensor. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

一个求解约束非线性优化问题的微分方程方法

本文构造的求解非线性优化问题的微分方程方法包括两个微分方程系统,第一个系统基于问题函数的一阶信息,第二个系统基于二阶信息．这两个系统具有性质:非线性优化问题的局部最优解是它们的渐近稳定的平衡点,并且初始点是可行点时,解轨迹都落于可行域中．我们证明了两个微分方程系统的离散迭代格式的收敛性定理和基于第二个系统的离散迭代格式的局部二次收敛性质．还给出了基于两个系统的离散迭代方法的数值算例,数值结果表明基于二阶信息的微分方程方法速度更快．     

A phenomenological method of calculating the residual stresses and plastic deformations in a hollow surface-hardened cylindrical sample

A phenomenological method is proposed for calculating the residual stress and plastic deformation fields in a hollow surface-hardened cylindrical sample. Versions of the hardening are considered that lead to isotropy and anisotropy in the plastic deformations in the surface layer. A hardening anisotropy parameter is introduced that relates the axial and circumferential components of the residual plastic deformation tensor. The experimentally determined axial and/or circumferential components of the residual plastic stress tensor are used as the initial information. The tensor fields of the residual stresses and deformations are constructed assuming the hypothesis of surface hardening anisotropy and the absence of secondary plastic compression deformations and that the tangential components of the residual stress tensor and the plastic incompressibility of the material are small. A technique is developed for identifying the parameters of the proposed method. The adequacy is checked using experimental data for test pieces of type 45 and 12X18H10T steels hardened by hydro-shot blasting treatment and of type 45 steel hardened by treatment with a roller. Good agreement is observed between the calculated and experimental results. It is noted that the anisotropic hardening procedure leads to a substantial difference between the circumferential and axial components of the residual stresses in the hardened layer, unlike the case of isotropic hardening where they are practically identical.     

Riemannian exponential and quantization

This article continues and completes the previous one [18]. First of all, we present two methods of quantization associated with a linear connection given on a differentiable manifold, one of them being the one presented in [18]. The two methods allow quantization of functions that come from covariant tensor fields. The equivalence of both is demonstrated as a consequence of a remarkable property of the Riemannian exponential (Theorem 5.1) that, as far as we know, is new to the literature. In addition, we provide a characterization of the Schrödinger operators as the only ones that by quantization correspond to classical mechanical systems. Finally, it is shown that the extension of the above quantization to functions of a very broad type can be carried out by generalizing the method of [18] in terms of fields of distributions.     

基于方向多尺度变换和统计建模的纹理分类方法

纹理是图像分析和识别中经常使用的关键特征, 而小波变换则是图像纹理表示和分类中的常用工具. 然而, 基于小波变换的纹理分类方法常常忽略了小波低频子带信息, 并且无法提取图像纹理的块状奇异信息. 本文提出小波子带系数的局部能量直方图建模方法、轮廓波特征的Poisson 混合模型建模方法和基于轮廓波子带系数聚类的特征提取方法, 并将其应用于图像纹理分类上. 基于局部能量直方图的纹理分类方法解决了小波低频子带的建模难题, 基于Poisson 混合模型的纹理分类方法则首次将Poisson 混合模型用于轮廓子带特征的建模, 而基于轮廓波域聚类的纹理分类方法是一种快速的分类方法. 实验结果显示, 本文所提出的三类方法都超过了当前典型的纹理分类方法.     

A combined finite volume–nonconforming/mixed-hybrid finite element scheme for degenerate parabolic problems

We propose and analyze a numerical scheme for nonlinear degenerate parabolic convection–diffusion–reaction equations in two or three space dimensions. We discretize the diffusion term, which generally involves an inhomogeneous and anisotropic diffusion tensor, over an unstructured simplicial mesh of the space domain by means of the piecewise linear nonconforming (Crouzeix–Raviart) finite element method, or using the stiffness matrix of the hybridization of the lowest-order Raviart–Thomas mixed finite element method. The other terms are discretized by means of a cell-centered finite volume scheme on a dual mesh, where the dual volumes are constructed around the sides of the original mesh. Checking the local Péclet number, we set up the exact necessary amount of upstream weighting to avoid spurious oscillations in the convection-dominated case. This technique also ensures the validity of the discrete maximum principle under some conditions on the mesh and the diffusion tensor. We prove the convergence of the scheme, only supposing the shape regularity condition for the original mesh. We use a priori estimates and the Kolmogorov relative compactness theorem for this purpose. The proposed scheme is robust, only 5-point (7-point in space dimension three), locally conservative, efficient, and stable, which is confirmed by numerical experiments.This work was supported by the GdR MoMaS, CNRS-2439, ANDRA, BRGM, CEA, EdF, France.     

Robust iterative methods for elliptic problems with highly varying coefficients in thin substructures

Summary. In this paper we introduce a class of robust multilevel interface solvers for two-dimensional finite element discrete elliptic problems with highly varying coefficients corresponding to geometric decompositions by a tensor product of strongly non-uniform meshes. The global iterations convergence rate is shown to be of the order with respect to the number of degrees of freedom on the single subdomain boundaries, uniformly upon the coarse and fine mesh sizes, jumps in the coefficients and aspect ratios of substructures. As the first approach, we adapt the frequency filtering techniques [28] to construct robust smoothers on the highly non-uniform coarse grid. As an alternative, a multilevel averaging procedure for successive coarse grid correction is proposed and analyzed. The resultant multilevel coarse grid preconditioner is shown to have (in a two level case) the condition number independent of the coarse mesh grading and jumps in the coefficients related to the coarsest refinement level. The proposed technique exhibited high serial and parallel performance in the skin diffusion processes modelling [20] where the high dimensional coarse mesh problem inherits a strong geometrical and coefficients anisotropy. The approach may be also applied to magnetostatics problems as well as in some composite materials simulation. Received December 27, 1994     

Robust multilevel methods for quadratic finite element anisotropic elliptic problems

This paper discusses a class of multilevel preconditioners based on approximate block factorization for conforming finite element methods employing quadratic trial and test functions. The main focus is on diffusion problems governed by a scalar elliptic partial differential equation with a strongly anisotropic coefficient tensor. The proposed method provides a high robustness with respect to non‐grid‐aligned anisotropy, which is achieved by the interaction of the following components: (i) an additive Schur complement approximation to construct the coarse‐grid operator; (ii) a global block (Jacobi or Gauss–Seidel) smoother complementing the coarse‐grid correction based on (i); and (iii) utilization of an augmented coarse grid, which enhances the efficiency of the interplay between (i) and (ii). The performed analysis indicates the high robustness of the resulting two‐level method. Moreover, numerical tests with a nonlinear algebraic multilevel iteration method demonstrate that the presented two‐level method can be applied successfully in the recursive construction of uniform multilevel preconditioners of optimal or nearly optimal order of computational complexity. Copyright © 2013 John Wiley & Sons, Ltd.     

MIXED METHODS FOR COMPRESSIBLE MISCIBLE DISPLACEMENT WITH THE EFFECT OF MOLECULAR DISPERSION

ＭＩＸＥＤＭＥＴＨＯＤＳＦＯＲＣＯＭＰＲＥＳＳＩＢＬＥＭＩＳＣＩＢＬＥＤＩＳＰＬＡＣＥＭＥＮＴＷＩＴＨＴＨＥＥＦＦＥＣＴＯＦＭＯＬＥＣＵＬＡＲＤＩＳＰＥＲＳＩＯＮＬＩＱＩＡＮ（李潜）（ＤｅｐａｒｔｍｅｎｔｏｆＭａｔｈｅｍａｔｉｃｓ，ＳｈａｎｄｏｎｇＮ...     

Tensor structures on a differentiable manifold

Summary A tensor structure is a class of equivalent G-structures and it is defined by a special tensor field. Such fields are characterized by the existence of a linear connection relative to which they have covariant derivative zero. Two tensor structures may admit a common subordinate structure. Exaples of such subordinate stuctures are given and some cases, when one stucture is a Riemannian metric, are considered. To Enrico Bompiani on his scientific Jubilee     

Propagation of spiral waves in anisotropic media: from waves to stripes

The effects of anisotropic diffusion on the propagation of spiral waves in excitable media are studied numerically by means of time-linearized methods in very refined meshes and with very small time steps. It is shown that the anisotropy of the inhibitor's diffusivity tensor does not play as important role on wave propagation as that of the activator. It is also shown that the off-diagonal components of the activator's diffusivity tensor cause stretching of the activator's concentration along the principal directions of this tensor, and that a large difference between the diffusion coefficients in the x- and y-directions may result in the spiral wave annihilation and stripe formation.     

Multidimensional upwind convection schemes for flow in porous media on structured and unstructured quadrilateral grids

Standard reservoir simulation schemes employ first order upwind schemes for approximation of the convective fluxes when multiple phases or components are present. These convective flux schemes rely upon upwind information that is determined according to grid geometry. As a consequence directional diffusion is introduced into the solution that is grid dependent. The effect can be particularly important for cases where the flow is across grid coordinate lines and is known as cross-wind diffusion.Truly higher dimensional upwind schemes that minimize cross-wind diffusion are presented for convective flow approximation on quadrilateral unstructured grids. The schemes are locally conservative and yield improved results that are essentially free of spurious oscillations. The higher dimensional schemes are coupled with full tensor Darcy flux approximations.The benefits of the resulting schemes are demonstrated for classical test problems in reservoir simulation including cases with full tensor permeability fields. The test cases involve a range of structured and unstructured grids with variations in orientation and permeability that lead to flow fields that are poorly resolved by standard simulation methods. The higher dimensional formulations are shown to effectively reduce the numerical cross-wind diffusion effect, leading to improved resolution of concentration and saturation fronts.     

A positive cell vertex Godunov scheme for a Beeler–Reuter based model of cardiac electrical activity

The monodomain model is a widely used model in electrocardiology to simulate the propagation of electrical potential in the myocardium. In this paper, we investigate a positive nonlinear control volume finite element scheme, based on Godunov's flux approximation of the diffusion term, for the monodomain model coupled to a physiological ionic model (the Beeler–Reuter model) and using an anisotropic diffusion tensor. In this scheme, degrees of freedom are assigned to vertices of a primal triangular mesh, as in conforming finite element methods. The diffusion term which involves an anisotropic tensor is discretized on a dual mesh using the diffusion fluxes provided by the conforming finite element reconstruction on the primal mesh and the other terms are discretized by means of an upwind finite volume method on the dual mesh. The scheme ensures the validity of the discrete maximum principle without any restriction on the transmissibility coefficients. By using a compactness argument, we obtain the convergence of the discrete solution and as a consequence, we get the existence of a weak solution of the original model. Finally, we illustrate by numerical simulations that the proposed scheme successfully removes nonphysical oscillations in the propagation of the wavefront and maintains conduction velocity close to physiological values.     

Directional Frames for Image Recovery: Multi-scale Discrete Gabor Frames

Sparsity-driven image recovery methods assume that images of interest can be sparsely approximated under some suitable system. As discontinuities of 2D images often show geometrical regularities along image edges with different orientations, an effective sparsifying system should have high orientation selectivity. There have been enduring efforts on constructing discrete frames and tight frames for improving the orientation selectivity of tensor product real-valued wavelet bases/frames. In this paper, we studied the general theory of discrete Gabor frames for finite signals, and constructed a class of discrete 2D Gabor frames with optimal orientation selectivity for sparse image approximation. Besides high orientation selectivity, the proposed multi-scale discrete 2D Gabor frames also allow us to simultaneously exploit sparsity prior of cartoon image regions in spatial domain and the sparsity prior of textural image regions in local frequency domain. Using a composite sparse image model, we showed the advantages of the proposed discrete Gabor frames over the existing wavelet frames in several image recovery experiments.     

Anomalous scaling in statistical models of passively advected vector fields

We use the methods of the renormalization group and the operator product expansion to consider the problem of the stochastic advection of a passive vector field with the most general form of the nonlinear term allowed by the Galilean symmetry. The external velocity field satisfies the Navier-Stokes equation. We show that the correlation functions have anomalous scaling in the inertial range. The corresponding anomalous exponents are determined by the critical dimensions of tensor composite fields (operators) built from only the fields themselves. We calculate the anomalous dimensions in the leading order of the expansion in the exponent in the correlator of the external force in the Navier-Stokes equation (the oneloop approximation of the renormalization group). The anomalous exponents exhibit a hierarchy related to the anisotropy degree: the lower the rank of the tensor operator is, the lower its dimension. The leading asymptotic terms are determined by the scalar operators in both the isotropic and the anisotropic cases, which completely agrees with Kolmogorov’s hypothesis of local isotropy restoration.     

一个基于张量火车分解的张量填充方法及在图像恢复中的应用

低秩张量填充在数据恢复中有广泛应用, 基于张量火车(TT) 分解的张量填充模型在彩色图像和视频以及互联网数据恢复中应用效果良好。本文提出一个基于三阶张量TT分解的填充模型。在模型中, 引入稀疏正则项与时空正则项, 分别刻画核张量的稀疏性和数据固有的块相似性。根据问题的结构特点, 引入辅助变量将原模型等价转化成可分离形式, 并采用临近交替极小化(PAM) 与交替方向乘子法(ADMM) 相结合的方法求解模型。数值实验表明, 两正则项的引入有利于提高数据恢复的稳定性和实际效果, 所提出方法优于其他方法。在采样率较低或图像出现结构性缺失时, 其方法效果较为显著。     

一个基于张量火车分解的张量填充方法及在图像恢复中的应用

低秩张量填充在数据恢复中有广泛应用, 基于张量火车(TT) 分解的张量填充模型在彩色图像和视频以及互联网数据恢复中应用效果良好。本文提出一个基于三阶张量TT分解的填充模型。在模型中, 引入稀疏正则项与时空正则项, 分别刻画核张量的稀疏性和数据固有的块相似性。根据问题的结构特点, 引入辅助变量将原模型等价转化成可分离形式, 并采用临近交替极小化(PAM) 与交替方向乘子法(ADMM) 相结合的方法求解模型。数值实验表明, 两正则项的引入有利于提高数据恢复的稳定性和实际效果, 所提出方法优于其他方法。在采样率较低或图像出现结构性缺失时, 其方法效果较为显著。