A global shallow water model using high order multi-moment constrained finite volume method and icosahedral grid

A novel accurate numerical model for shallow water equations on sphere have been developed by implementing the high order multi-moment constrained finite volume (MCV) method on the icosahedral geodesic grid. High order reconstructions are conducted cell-wisely by making use of the point values as the unknowns distributed within each triangular cell element. The time evolution equations to update the unknowns are derived from a set of constrained conditions for two types of moments, i.e. the point values on the cell boundary edges and the cell-integrated average. The numerical conservation is rigorously guaranteed. in the present model, all unknowns or computational variables are point values and no numerical quadrature is involved, which particularly benefits the computational accuracy and efficiency in handling the spherical geometry, such as coordinate transformation and curved surface. Numerical formulations of third and fourth order accuracy are presented in detail. The proposed numerical model has been validated by widely used benchmark tests and competitive results are obtained. The present numerical framework provides a promising and practical base for further development of atmospheric and oceanic general circulation models. (C) 2009 Elsevier Inc. All rights reserved.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Comparative on causes and accumulation of selenium in the tree-rings ambient high-selenium coal combustion area from Yutangba；Hubei；China

IEECAS SKLLQG

Magnets for HIRFL-CSR rings

The magnet design, fabrication, and measurement of HIRFL-CSR (Heavy Ion Research Facility in Lanzhou Cooling Storage Ring) are presented. All magnets will be laminated And welded with an armor-coated surface between two big endplates made of sticking glue 0.5 mm-thick sheets. The dipole of CSRm was chosen an H type with an air circle on the pole to improve the field uniformity. The dipole of CSRe was chosen the C type with an air circle and two air slots on the pole to improve the field homogeneity. Its reproducibility of magnet to magnet was adjusted with inserting small laminating pieces before demountable pole ends to reach less than +/- 2 x 10(-4) at optimized field level. CSRm quadrupoles diameter is 170 mm and has two different lengths, and its endplates were made with punching pieces after coating with epoxy glue, there is chamfered directly on the pole ends to reduce 12th-order contribution of field and without the demountable pole ends. CSRe main quadrupoles diameter is 240 mm and has two different lengths, and its endplates were also made with punching pieces coated with epoxy glue, there is also chamfered directly on the pole ends to reduce 12th-order contribution of field like CSRm.

Symmetry energy and isovector giant dipole resonance in finite nuclei

We study the relationship between the properties of the isovector giant dipole resonance of finite nuclei and the symmetry energy in the framework of the relativistic mean field theory with six different parameter sets of nonlinear effective Lagrangian. A strong linear correlation of excited energies of the dipole resonance in finite nuclei and symmetry energy at and below the saturation density is found. This linear correlation leads to the symmetry energy at the saturation density at the interval 33.0MeV <= S(po) <= 37.0 MeV. The comparison to the present experimental data in the soft dipole mode of (132) Sn constrains approximately the symmetry energy at p = 0.1 fm(-3) at the interval 21.2MeV similar to 22.5 MeV. It is proposed that a precise measurement of the soft dipole mode in neutron rich nuclei could set up an important constraint on the equation of state for asymmetric nuclear matter.