Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functions

This paper addresses the convergence properties of implicit numerical solution algorithms for nonlinear hyperbolic transport problems. It is shown that the Newton–Raphson (NR) method converges for any time step size, if the flux function is convex, concave, or linear, which is, in general, the case for CFD problems. In some problems, e.g., multiphase flow in porous media, the nonlinear flux function is S-shaped (not uniformly convex or concave); as a result, a standard NR iteration can diverge for large time steps, even if an implicit discretization scheme is used to solve the nonlinear system of equations. In practice, when such convergence difficulties are encountered, the current time step is cut, previous iterations are discarded, a smaller time step size is tried, and the NR process is repeated. The criteria for time step cutting and selection are usually based on heuristics that limit the allowable change in the solution over a time step and/or NR iteration. Here, we propose a simple modification to the NR iteration scheme for conservation laws with S-shaped flux functions that converges for any time step size. The new scheme allows one to choose the time step size based on accuracy consideration only without worrying about the convergence behavior of the nonlinear solver. The proposed method can be implemented in an existing simulator, e.g., for CO₂ sequestration or reservoir flow modeling, quite easily. The numerical analysis is confirmed with simulation studies using various test cases of nonlinear multiphase transport in porous media. The analysis and numerical experiments demonstrate that the modified scheme allows for the use of arbitrarily large time steps for this class of problems.     

Reply to comment on “On the controversy around Daganzo’s requiem for and Aw-Rascle’s resurrection of second-order traffic flow models” by H.M. Zhang

This contribution compares several different approaches allowing one to derive macroscopic traffic equation directly from microscopic car-following models. While it is shown that some conventional approaches lead to theoretical problems, it is proposed to use an approach reminding of smoothed particle hydrodynamics to avoid gradient expansions. The derivation circumvents approximations and, therefore, demonstrates the large range of validity of macroscopic traffic equations, without the need of averaging over many vehicles. It also gives an expression for the “traffic pressure”, which generalizes previously used formulas. Furthermore, the method avoids theoretical inconsistencies of macroscopic traffic models, which have been criticized in the past by Daganzo and others.     

Synchronization in evolving snowdrift game model

The interaction between the evolution of the game and the underlying network structure with evolving snowdrift game model is investigated. The constructed network follows a power-law degree distribution typically showing scale-free feature. The topological features of average path length, clustering coefficient, degree-degree correlations and the dynamical feature of synchronizability are studied. The synchronizability of the constructed networks changes by the interaction. It will converge to a certain value when sufficient new nodes are added. It is found that initial payoffs of nodes greatly affect the synchronizability. When initial payoffs for players are equal, low common initial payoffs may lead to more heterogeneity of the network and good synchronizability. When initial payoffs follow certain distributions, better synchronizability is obtained compared to equal initial payoff. The result is also true for phase synchronization of nonidentical oscillators.     

Elasticity transition and loop formation in vibrated bead chains: A simulation of polymer chains

By measuring the distribution function of the end-to-end distance, we find that strongly shaken bead chains exhibit many properties, such as the rigid-rod-to-Gaussian chain transition, scaling, fast drop of loop formation probability in the short-chain regime, and enhancement of loop formation probability for kinked chains, of long-chain polymers. Though there is difference in local details between our chains and the worm-like chains, our results are consistent with recent calculations based on the worm-like chain model in many respects.     

Multi-velocity and multi-temperature model of the mixture of polyatomic gases issuing from kinetic theory

In this paper, we consider Euler-like balance laws for mixture components that involve macroscopic velocities and temperatures for each different species. These laws are not conservation laws due to mutual interaction between species. In particular, source terms that describe the rate of change of momentum and energy of the constituents appear. These source terms are computed with the help of kinetic theory for mixtures of polyatomic gases. Moreover, if we restrict the attention to processes which occur in the neighborhood of the average velocity and temperature of the mixture, the phenomenological coefficients of extended thermodynamics can be determined from the computed source terms.     

The gauge theory of dislocations: conservation and balance laws

We derive conservation and balance laws for the translational gauge theory of dislocations by applying Noether's theorem. We present an improved translational gauge theory of dislocations including the dislocation density tensor and the dislocation current tensor. The invariance of the variational principle under the continuous group of transformations is studied. Through Lie's infinitesimal invariance criterion we obtain conserved translational and rotational currents for the total Lagrangian made up of an elastic and dislocation part. We calculate the broken scaling current. Looking only on one part of the whole system, the conservation laws are changed into balance laws. Because of the lack of translational, rotational and dilatation invariance for each part, a configurational force, moment and power appears. The corresponding J , L and M integrals are obtained. Only isotropic and homogeneous materials are considered and we restrict ourselves to a linear theory. We choose constitutive laws for the most general linear form of material isotropy. Also we give the conservation and balance laws corresponding to the gauge symmetry and the addition of solutions. From the addition of solutions we derive a reciprocity theorem for the gauge theory of dislocations. Also, we derive the conservation laws for stress-free states of dislocations.     

CONSERVATION LAWS FOR HEATED LAMINAR RADIAL LIQUID AND FREE JETS

The conserved quantities for the heated radial liquid jet and the heated radial free jet are established by using conservation laws. The flow in a heated radial jet is described by Prandtl's momentum boundary layer equation, the continuity equation and the energy equation. Viscous dissipation is neglected. The multiplier approach is used to derive the conservation laws for the system of three equations for the velocity components and the temperature and three conserved vectors are obtained. The conservation laws for the system of two partial differential equations for the stream function formulation are also computed by the multiplier approach and three conserved vectors are obtained. One of these is a non-local conserved vector for the system. The conserved quantities for the heated radial liquid jet and the heated radial free jet, emitted into a stationary fluid of uniform temperature θ_∞, are derived by integrating the conservation laws across the jet.     

On symmetries,reductions, conservation laws and conserved quantities of optical solitons with inter-modal dispersion

This paper studies the compressional dispersive Alfvén (CDA) waves where Noether symmetries will be calculated from which the corresponding conservation laws will be obtained via Noether's theorem. Furthermore, one case of double reduction is performed via the association of a conserved vector with a Noether symmetry (with zero gauge). The conserved quantities of optical solitons in the presence of intermodal dispersion that is governed by the perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity. The invariance-multiplier method is adopted to carry out the analysis, from which the conserved densities are then retrieved. Finally, the conserved quantities are obtained using the 1-soliton solution of the governing equation.