A multi-moment finite volume formulation for shallow water equations on unstructured mesh

A novel and accurate finite volume method has been presented to solve the shallow water equations on unstructured grid in plane geometry. In addition to the volume integrated average (VIA moment) for each mesh cell, the point values (PV moment) defined on cell boundary are also treated as the model variables. The volume integrated average is updated via a finite volume formulation, and thus is numerically conserved, while the point value is computed by a point-wise Riemann solver. The cell-wise local interpolation reconstruction is built based on both the VIA and the PV moments, which results in a scheme of almost third order accuracy. Efforts have also been made to formulate the source term of the bottom topography in a way to balance the numerical flux function to satisfy the so-called C-property. The proposed numerical model is validated by numerical tests in comparison with other methods reported in the literature. (C) 2010 Elsevier Inc. All rights reserved.

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.

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.

Hot Nuclear Matter Equation of State and Finite Temperature Kaon Condensation

We perform a systematic calculation of the equation of state of asymmetric nuclear matter at finite temperature within the framework of the Brueckner-Hartree-Fock approach with a microscopic three-body force. When applying it to the study of hotka on condensed matter, we find that the thermal effect is more profound in comparison with normal matter, in particular around the threshold density. Also, the increase of temperature makes the equation of state slightly stiffer through suppression of kaon condensation.

3D quench simulation using finite element method for CR superconducting magnet's strand

The CR superconducting magnet is a dipole of the FAIR project of GSI in Germany. The quench of the strand is simulated using FEM software ANSYS. From the simulation, the quench propagation can be visualized. Programming with APDL, the value of propagation velocity of normal zone is calculated. Also the voltage increasing over time of the strand is computed and pictured. Furthermore, the Minimum Propagation Zone (MPZ) is studied. At last, the relation between the current and the propagation velocity of normal zone, and the influence of initial temperature on quench propagation are studied.

Photon production in a chemically equilibrating quark-gluon plasma at finite Baryon density: Complete leading order results

We investigate hard photon production of the near-collinear bremsstrahlung and a new process called the inelastic pair annihilation, fully including the LPM effect, in a chemically equilibrating quark-gluon plasma at finite baryon density, and find that the effect of the system evolution on the photon production and large contribution of the bremsstrahlung make the total photon yield of the two processes as a strongly increasing function of the initial quark chemical potential.

Photon emission from a chemically non-equilibrated parton gas at finite baryon density

We study hard photon production from a chemically non-equilibrated quark-gluon plasma with finite baryon density on the basis of Juttner distribution of partons of the system. We find that the photon production is ruled by early times, main contributions are given by rapidities y <= 6, and photon yield is a strongly increasing function of the initial quark chemical potential. In addition, we note that contribution from bremsstrahlung and Compton process qg -> q gamma dominates.

Critical behavior of multifragmentation in finite nuclei and extraction of critical exponents

Using a phenomenological asymmetric nuclear equation of state, we obtained pressure-density isotherms of the finite nucleus Sn-112 simulated in r-space and in p-space and constructed the nuclear fragments by using the coalescence model. After correlatively analysing the fragments, the signal of critical behavior has been found and critical exponents were also extracted.