Aggregation Processes with Catalysis-Driven Decomposition

We propose a three-species aggregation model with catalysis-driven decomposition. Based on the mean-field rate equations, we investigate the evoIution behavior of the system with the size-dependent catalysis-driven decomposition rate J（i; j; k） = Jijk^v and the constant aggregation rates. The results show that the cluster size distribution of the species without decomposition can always obey the conventional scaling law in the case of 0 ≤v ≤ 1, while the kinetic evolution of the decomposed species depends crucially on the index v. Moreover, the total size of the species without decomposition can keep a nonzero value at large times, while the total size of the decomposed species decreases exponentially with time and vanishes finally.     

Influence of the halo upon angular distributions for elastic scattering and breakup

The angular distributions for elastic scattering and breakup of halo nuclei are analysed using a near-side/far-side decomposition within the framework of the dynamical eikonal approximation. This analysis is performed for ¹¹Be impinging on Pb at 69 MeV/nucleon. These distributions exhibit very similar features. In particular they are both near-side dominated, as expected from Coulomb-dominated reactions. The general shape of these distributions is sensitive mostly to the projectile-target interactions, but is also affected by the extension of the halo. This suggests the elastic scattering not to be affected by a loss of flux towards the breakup channel.     

Upwind-biased FORCE schemes with applications to free-surface shallow flows

In this paper we develop numerical fluxes of the centred type for one-step schemes in conservative form for solving general systems of conservation laws in multiple-space dimensions on structured meshes. The proposed method is an extension of the multidimensional FORCE flux developed by Toro et al. (2009) [14]. Here we introduce upwind bias by modifying the shape of the staggered mesh of the original FORCE method. The upwind bias is evaluated using an estimate of the largest eigenvalue, which in any case is needed for selecting a time step. The resulting basic flux is first-order accurate and monotone. For the linear advection equation, the proposed UFORCE method reproduces exactly the upwind Godunov method. Extension to non-linear systems has been done empirically via the two-dimensional inviscid shallow water equations. Second order of accuracy in space and time on structured meshes is obtained in the framework of finite volume methods. The proposed method improves the accuracy of the solution for small Courant numbers and intermediate waves associated with linearly degenerate fields (contact discontinuities, shear waves and material interfaces). It achieves comparable accuracy to that of upwind methods with approximate Riemann solvers, though retaining the simplicity and efficiency of centred methods. The performance of the schemes is assessed on a suite of test problems for the two-dimensional shallow water equations.     

A cell functional minimization scheme for parabolic problem

We are interested in a robust and accurate finite volume scheme for 2-D parabolic problems derived from the cell functional minimization approach. The scheme has a local stencil, is locally conservative, treats discontinuity rigorously and leads to a symmetric positive definite linear system. Since the scheme has both cell centered unknowns and cell edge unknowns, the computational cost is an issue and a parallel algorithm is then suggested based on nonoverlapping domain decomposition approach. The interface condition is of the Dirichlet–Robin type and has a parameter λ. By choosing this parameter properly, the convergence of the iteration process could be sped up. Numerical results for linear and nonlinear problems demonstrate the good performance of the cell functional minimization scheme and its parallel version on distorted meshes.     

Early stages of decomposition of alloys (spinodal or nucleation)

Progress in theory and experiment on the early stage of decomposition of alloys is reviewed. In the first half of this paper, theoretical treatments on decomposition of alloys by spinodal and nucleation mechanism and also continuum theories and the discrete lattice theories are comparatively reviewed. In the latter half of this paper, experimental results on Al-Zn, Ni-Al, Cu-Co and Al-Ag alloys are discussed.     

A covariant representation of the Ball–Chiu vertex

In nonabelian gauge theory the three-gluon vertex function contains important structural information, in particular on infrared divergences, and is also an essential ingredient in the Schwinger–Dyson equations. Much effort has gone into analyzing its general structure, and at the one-loop level also a number of explicit computations have been done, using various approaches. Here we use the string-inspired formalism to unify the calculations of the scalar, spinor and gluon loop contributions to the one-loop vertex, leading to an extremely compact representation in all cases. The vertex is computed fully off-shell and in dimensionally continued form, so that it can be used as a building block for higher-loop calculations. We find that the Bern–Kosower loop replacement rules, originally derived for the on-shell case, hold off-shell as well. We explain the relation of the structure of this representation to the low-energy effective action, and establish the precise connection with the standard Ball–Chiu decomposition of the vertex. This allows us also to predict that the vanishing of the completely antisymmetric coefficient function S   of this decomposition is not a one-loop accident, but persists at higher-loop orders. The sum rule found by Binger and Brodsky, which leads to the vanishing of the one-loop vertex in N=4$N=4$ SYM theory, in the present approach relates to worldline supersymmetry.     

求分块鳞状因子循环矩阵逆矩阵的一种快速算法

给出了一种计算分块鳞状因子循环矩阵逆矩阵的快速算法,该算法主要利用了离散傅立叶变换和对角块矩阵求逆的递归算法,与标准的利用LU分解法求逆的算法相比,在计算复杂性上有很大的优势.     

Spatiotemporal and statistical symmetries

The notion of symmetries, either statistical or deterministic, can be useful for the characterization of complex systems and their bifurcations. In this paper, we investigate the connection between the (microscopic) spatiotemporal symmetries of a space-time functionu(x, t), on the one hand, and the (macroscopic) symmetries of statistical quantities such as the spatial (resp. temporal) two-point correlations and the spatial (resp. temporal) average, on the other hand. We show, how, under certain conditions, these symmetries are related to the symmetries of the orbits described byu(x, t) in the characteristic (phase) spaces. We also determine the largest group of spatiotemporal symmetries (in the sense introduced in our earlier work) satisfied by a given space-time functionu(x, t) and indicate how to extract the subgroups of point symmetries, namely those directly implemented on the space and time variables. Conversely, we determine all the functions invariant by a given space-time symmetry group. Finally, we illustrate all the previous points with specific examples.     

Optimal time advancing dispersion relation preserving schemes

In this paper we examine the constrained optimization of explicit Runge–Kutta (RK) schemes coupled with central spatial discretization schemes to solve the one-dimensional convection equation. The constraints are defined with respect to the correct error propagation equation which goes beyond the traditional von Neumann analysis   developed in Sengupta et al. [T.K. Sengupta, A. Dipankar, P. Sagaut, Error dynamics: beyond von Neumann analysis, J. Comput. Phys. 226 (2007) 1211–1218]. The efficiency of these optimal schemes is demonstrated for the one-dimensional convection problem and also by solving the Navier–Stokes equations for a two-dimensional lid-driven cavity (LDC) problem. For the LDC problem, results for Re=1000$\mathit{Re}=1000$ are compared with results using spectral methods in Botella and Peyret [O. Botella, R. Peyret, Benchmark spectral results on the lid-driven cavity flow, Comput. Fluids 27 (1998) 421–433] to calibrate the method in solving the steady state problem. We also report the results of the same flow at Re=10,000$\mathit{Re}=10\text{,}000$ and compare them with some recent results to establish the correctness and accuracy of the scheme for solving unsteady flow problems. Finally, we also compare our results for a wave-packet propagation problem with another method developed for computational aeroacoustics.     

The influence of cell geometry on the Godunov scheme applied to the linear wave equation

By studying the structure of the discrete kernel of the linear acoustic operator discretized with a Godunov scheme, we clearly explain why the behaviour of the Godunov scheme applied to the linear wave equation deeply depends on the space dimension and, especially, on the type of mesh. This approach allows us to explain why, in the periodic case, the Godunov scheme applied to the resolution of the compressible Euler or Navier–Stokes system is accurate at low Mach number when the mesh is triangular or tetrahedral and is not accurate when the mesh is a 2D (or 3D) cartesian mesh. This approach confirms also the fact that a Godunov scheme remains accurate when it is modified by simply centering the discretization of the pressure gradient.