A mass conservative well‐balanced reconstruction at wet/dry interfaces for the Godunov‐type shallow water model

This paper presents a novel mass conservative, positivity preserving wetting and drying treatment for Godunov‐type shallow water models with second‐order bed elevation discretization. The novel method allows to compute water depths equal to machine accuracy without any restrictions on the time step or any threshold that defines whether the finite volume cell is considered to be wet or dry. The resulting scheme is second‐order accurate in space and keeps the C‐property condition at fully flooded area and also at the wet/dry interface. For the time integration, a second‐order accurate Runge–Kutta method is used. The method is tested in two well‐known computational benchmarks for which an analytical solution can be derived, a C‐property benchmark and in an additional example where the experimental results are reproduced. Overall, the presented scheme shows very good agreement with the reference solutions. The method can also be used in the discontinuous Galerkin method. Copyright © 2016 John Wiley & Sons, Ltd.     

A three‐dimensional Cartesian cut cell method for incompressible viscous flow with irregular domains

A three‐dimensional Cartesion cut cell method is presented for the simulations of incompressible viscous flows with irregular domains. A new model (referred to as ‘6+N’ model) is proposed to describe arbitrarily shaped cut cells and treat all the cells as polyhedrons with 6+N faces. The finite volume discretization of the Navier–Stokes equation is then implemented by using the ‘6+N’ model to separate the surface flux integrals into two parts, that is, the fluxes through the basic face of the hexahedron and those through the cutting surfaces. The previously proposed Kitta Cube algorithm and volume computer‐aided design platform (J. Comput. Aided. Des. 2005; 37(4): 1509–1520. Doi:10.1016/j.cad.2005.03.006) are adopted to generate cut cells and provide shape data and physical attributes for the numerical analysis. A modified SIMPLE‐based smoothing pressure correction scheme is applied to suppress checkerboard pressure oscillations caused by the collocated arrangement of velocities and pressure. The calculation accuracy of the numerical method expressed by L₁ and L _∞ norm errors is first demonstrated by the simulation of a pipe flow. Then its feasibility, efficiency, and potential in engineering applications are verified by applying it to solve natural convections between concentric spheres and between eccentric spheres. The heat transfer patterns in eccentric spheres are also obtained by using the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

Implementation of ADER scheme for a bore on an unsaturated permeable slope

The paper details the implementation of the Godunov‐type finite volume Arbitrary high order schemes using Derivatives (ADER) scheme for the case of a large source term in the continuity equation of the nonlinear shallow water equations. The particular application is the movement of a bore on a highly permeable slope. The large source term is caused by the infiltration into the initially unsaturated slope material. Infiltration is modelled as vertical downwards piston‐like flow with Forchheimer quadratic parameterisation of the resistance law. The corresponding ODE is solved using the fourth‐order Runge–Kutta method. The surface and subsurface flow models have been tested by comparison with analytical solutions. Example predictions of surface bore propagation and wetting front propagation are presented for a range of slope permeabilities. The effects of permeability on bore run‐up, water depths and velocities are illustrated. The ADER scheme is capable of handling the source term, including the extreme case when this term dominates the volume balance. Copyright © 2011 John Wiley & Sons, Ltd.     

Approximant crystals of icosahedral Mg–(Ga,Al)–Zn alloys

The zero field cooled (ZFC) and field cooled (FC) low-field magnetic moment m of a dense frozen ferrofluid containing Fe₅₅Co₄₅ particles of size 4.6nm in hexane exhibits irreversibility at temperatures T?T _b≈ 30?K. FC in μ ₀ H ≤ 1?T gives rise to shifted minor hysteresis loops below T _b. At T _c≈ 10?K, sharp peaks of m ^ZFC and of the ac susceptibility χ ′, a kink of the thermoremanent magnetic moment m ^TRM, a sizeable reduction of the coercive field H _c, and the appearance of a spontaneous moment m ^SFM indicate a phase transition with near mean-field critical behaviour of both m ^SFM and χ ′ . These features are explained within a core-shell model of nanoparticles, whose strongly disordered shells gradually become blocked below T _b, while their soft ferromagnetic cores couple dipolarly and become superferromagnetic (SFM) below T _c.     

A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

Accurate double many‐body expansion potential energy surface by extrapolation to the complete basis set limit and dynamics calculations for ground state of NH2

An accurate single‐sheeted double many‐body expansion potential energy surface is reported for the title system. A switching function formalism has been used to warrant the correct behavior at the and dissociation channels involving nitrogen in the ground and first excited states. The topographical features of the novel global potential energy surface are examined in detail, and found to be in good agreement with those calculated directly from the raw ab initio energies, as well as previous calculations available in the literature. The novel surface can be using to treat well the Renner–Teller degeneracy of the and states of . Such a work can both be recommended for dynamics studies of the reaction and as building blocks for constructing the double many‐body expansion potential energy surface of larger nitrogen/hydrogen‐containing systems. In turn, a test theoretical study of the reaction has been carried out with the method of quantum wave packet on the new potential energy surface. Reaction probabilities, integral cross sections, and differential cross sections have been calculated. Threshold exists because of the energy barrier (68.5 meV) along the minimum energy path. On the curve of reaction probability for total angular momentum J = 0, there are two sharp peaks just above threshold. The value of integral cross section increases quickly from zero to maximum with the increase of collision energy, and then stays stable with small oscillations. The differential cross section result shows that the reaction is a typical forward and backward scatter in agreement with experimental measurement result. © 2013 Wiley Periodicals, Inc.