Sparse-grid finite-volume multigrid for 3D-problems

We introduce a multigrid algorithm for the solution of a second order elliptic equation in three dimensions. For the approximation of the solution we use a partially ordered hierarchy of finite-volume discretisations. We show that there is a relation with semicoarsening and approximation by more-dimensional Haar wavelets. By taking a proper subset of all possible meshes in the hierarchy, a sparse grid finite-volume discretisation can be constructed.The multigrid algorithm consists of a simple damped point-Jacobi relaxation as the smoothing procedure, while the coarse grid correction is made by interpolation from several coarser grid levels.The combination of sparse grids and multigrid with semi-coarsening leads to a relatively small number of degrees of freedom,N, to obtain an accurate approximation, together with anO(N) method for the solution. The algorithm is symmetric with respect to the three coordinate directions and it is fit for combination with adaptive techniques.To analyse the convergence of the multigrid algorithm we develop the necessary Fourier analysis tools. All techniques, designed for 3D-problems, can also be applied for the 2D case, and — for simplicity — we apply the tools to study the convergence behaviour for the anisotropic Poisson equation for this 2D case.     

长柱形CsI(Tl)闪烁晶体探测器的能量重建

提出了用于TEXONO实验中双端读出的长柱形CsI(Tl)闪烁晶体探测器的两种粒子能量重建方法———“算术平均值能量重建法”和“几何平均值能量重建法”,讨论了两种算法的理论依据和计算方法,并对两种算法进行了比较     

Discrete duality finite volume schemes for Leray−Lions−type elliptic problems on general 2D meshes

Discrete duality finite volume schemes on general meshes, introduced by Hermeline and Domelevo and Omnès for the Laplace equation, are proposed for nonlinear diffusion problems in 2D with nonhomogeneous Dirichlet boundary condition. This approach allows the discretization of non linear fluxes in such a way that the discrete operator inherits the key properties of the continuous one. Furthermore, it is well adapted to very general meshes including the case of nonconformal locally refined meshes. We show that the approximate solution exists and is unique, which is not obvious since the scheme is nonlinear. We prove that, for general W^?1,p′(Ω) source term and W^1‐(1/p),p(?Ω) boundary data, the approximate solution and its discrete gradient converge strongly towards the exact solution and its gradient, respectively, in appropriate Lebesgue spaces. Finally, error estimates are given in the case where the solution is assumed to be in W^2,p(Ω). Numerical examples are given, including those on locally refined meshes. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

An invariance of geometric mean with respect to Lagrangian means

The invariance of the geometric mean G with respect to the Lagrangian mean-type mapping (Lf,Lg), i.e. the equation G○(Lf,Lg)=G, is considered. We show that the functions f and g must be of high class regularity. This fact allows to reduce the problem to a differential equation and determine the second derivatives of the generators f and g.     

17N+197Au反应中碎片与中子的符合测量

在3?3?.4MeV/u ¹⁷N束轰击¹⁹⁷Au靶产生的反应中,利用放置于不同角度组合的17个中子探测器(4°—83°)和14个半导体望远镜(2.3°—9.0°)对反应产物碎片与中子进行了符合测量.经对所得实验角分布积分得到Z=3—6元素的同位素产额分布.在参加者-旁观者模型框架下,采用¹⁷N原子核内部的不同密度分布计算了同位素产额分布并与实验数据做比较.     

Horizontal Connection and Horizontal Mean Curvature in Carnot Groups

In this paper we give a geometric interpretation of the notion of the horizontal mean curvature which is introduced by Danielli Garofalo-Nhieu and Pauls who recently introduced sub- Riemannian minimal surfaces in Carnot groups. This will be done by introducing a natural nonholonomic connection which is the restriction （projection） of the natural Riemannian connection on the horizontal bundle. For this nonholonomic connection and （intrinsic） regular hypersurfaces we introduce the notions of the horizontal second fundamental form and the horizontal shape operator. It turns out that the horizontal mean curvature is the trace of the horizontal shape operator.     

On the uniqueness of the foliation of spheres of constant mean curvature in asymptotically flat 3-manifolds

In this note we study constant mean curvature surfaces in asymptotically flat 3-manifolds. We prove that, outside a given compact subset in an asymptotically flat 3-manifold with positive mass, stable spheres of given constant mean curvature are unique. Therefore we are able to conclude that the foliation of stable spheres of constant mean curvature in an asymptotically flat 3-manifold with positive mass outside a given compact subset is unique.

     

Error estimate for the approximation of nonlinear conservation laws on bounded domains by the finite volume method

In this paper we derive a priori and a posteriori error estimates for cell centered finite volume approximations of nonlinear conservation laws on polygonal bounded domains. Numerical experiments show the applicability of the a posteriori result for the derivation of local adaptive solution strategies.

     

On the Isoperimetric Type Problems with Current Hyperplanes

Analytical approach is presented to the isoperimetric type problems of finding convex figures separated by current hyperplanes.