An alternating direction implicit （ADI） numerical model for two－dimensional hydrodynamic equations

A two-dimensional computational model is develope for the calulation of tides, storm surges and otherlong-period waves in coastal and shelf waters. The Partial differental equations are approximated by two sets of difference equations on a space-staggered grid system. Both sets are explicit with one set for water level and x-component velocity, and another for water level and y-component velocity. These two sets are used successively for stepby-step solution in time. An analytical investigation on the linearized sets of the difference equations indicates that thecomputational scheme is unconditionally stable. The model is of second order accuracy both in space and in time andconserves mass and momentum. Simulations of surface elevation caused by periodic forcing in one-opening rectangularbasin with flat topography and by steady wind stress in the basin with flat or slope topography show that the computed results are in excellent agreement with the corresponding analytic solutions. The steady-tate wind-induced setupin a ofed basin with discontinuous topography computed with the present model are also in excellent agreement withthe results from Leendertse's model. Finally, the model is applied to hindcast a storm surge in the South China Seaand reproduces the surge elevation satisfactorily.

Advancing Hierarchical Equations of Motion for Efficient Evaluation of Coherent Two-dimensional Spectroscopy

To advance hierarchical equations of motion as a standard theory for quantum dissipative dynamics, we put forward a mixed Heisenberg-SchrSdinger scheme with block-matrix implementation on efficient evaluation of nonlinear optical response function. The new approach is also integrated with optimized hierarchical theory and numerical filtering algorithm. Different configurations of coherent two-dimensional spectroscopy of model excitonic dimer systems are investigated, with focusing on the effects of intermolecular transfer coupling and bi-exciton interaction.

Normalization of Hydrodynamic Coefficients in Morison Equation

The hydrodynamic coefficients C-d and C-m are not only dependent on the size of slender cylinder, its location in water, KC number and Re number, but also vary with environmental conditions, i.e., in regular waves or in irregular waves, in pure waves or in wave-current coexisting field. In this paper, the normalization of hydrodynamic coefficients for various environmental conditions is discussed. When a proper definition of KC number and proper characteristic values of irregular waves are used, a unified relationship between C-d, C-m and KC number for regular waves, irregular waves, pure waves and wave-current coexisting field can be obtained.

TWO-DIMENSIONAL RIEMANN PROBLEMS:FROM SCALAR CONSERVATION LAWS TO COMPRESSIBLE EULER EQUATIONS

In this paper we survey the authors＇ and related work on two-dimensional Riemann problems for hyperbolic conservation laws, mainly those related to the compressible Euler equations in gas dynamics. It contains four sections： 1. Historical review. 2. Scalar conservation laws. 3. Euler equations. 4. Simplified models.

Rational and Periodic Wave Solutions of Two-Dimensional Boussinesq Equation

Two new exact, rational and periodic wave solutions are derived for the two-dimensional Boussinesq equation. For the first solution it is obtained by performing an appropriate limiting procedure on the soliton solutions obtained by Hirota bilinear method. The second one in terms of Riemann theta function is explicitly presented by virtue of Hirota bilinear method and its asymptotic property is also analyzed in detail. Moreover, it is of interest to note that classical soliton solutions can be reduced from the periodic wave solutions.

Multi-Symplectic Splitting Method for Two-Dimensional Nonlinear Schrodinger Equation

Using the idea of splitting numerical methods and the multi-symplectic methods, we propose a multisymplectic splitting （MSS） method to solve the two-dimensional nonlinear Schrodinger equation （2D-NLSE） in this paper. It is further shown that the method constructed in this way preserve the global symplectieity exactly. Numerical experiments for the plane wave solution and singular solution of the 2D-NLSE show the accuracy and effectiveness of the proposed method.

Symmetry Reduction of Two-Dimensional Damped Kuramoto-Sivashinsky Equation

In this paper, the problem of determining the largest possible set of symmetries for an important nonlinear dynamical system： the two-dimensional damped Kuramoto-Sivashinsky （（21）） DKS ） equation is studied. By applying the basic Lie symmetry method for the （217）） DKS equation, the classical Lie point symmetry operators are obtained. Also, the optimal system of one-dimensional subalgebras of the equation is constructed. The Lie invariants as well as similarity reduced equations corresponding to infinitesimal symmetries are obtained. The nonclassicaJ symmetries of the （2D） DKS equation are also investigated.

Two-dimensional equations for thin-films of ionic conductors

A theoretical model is developed for predicting both conduction and diffusion in thin-film ionic conductors or cables. With the linearized Poisson-Nernst-Planck(PNP)theory, the two-dimensional(2D) equations for thin ionic conductor films are obtained from the three-dimensional(3D) equations by power series expansions in the film thickness coordinate, retaining the lower-order equations. The thin-film equations for ionic conductors are combined with similar equations for one thin dielectric film to derive the 2D equations of thin sandwich films composed of a dielectric layer and two ionic conductor layers. A sandwich film in the literature, as an ionic cable, is analyzed as an example of the equations obtained in this paper. The numerical results show the effect of diffusion in addition to the conduction treated in the literature. The obtained theoretical model including both conduction and diffusion phenomena can be used to investigate the performance of ionic-conductor devices with any frequency.

Adaptive Moving Mesh Central-Upwind Schemes for Hyperbolic System of PDEs:Applications to Compressible Euler Equations and Granular Hydrodynamics

We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.

High-Order Local Discontinuous Galerkin Algorithm with Time Second-Order Schemes for the Two-Dimensional Nonlinear Fractional Diffusion Equation

In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.

LOGARITHMICALLY IMPROVED REGULARITY CRITERION FOR THE 3D GENERALIZED MAGNETO-HYDRODYNAMIC EQUATIONS

This article proves the logarithmically improved Serrin＇s criterion for solutions of the 3D generalized magneto-hydrodynamic equations in terms of the gradient of the velocity field, which can be regarded as improvement of results in [10] （Luo Y W. On the regularity of generalized MHD equations. J Math Anal Appl, 2010, 365： 806-808） and [18] （Zhang Z J. Remarks on the regularity criteria for generalized MHD equations. J Math Anal Appl, 2011, 375：799 802）.

Numerical Analysis of Upwind Difference Schemes for Two-Dimensional First-Order Hyperbolic Equations with Variable Coefficients

In this paper, we consider the initial-boundary value problem of two-dimensional first-order linear hyperbolic equation with variable coefficients. By using the upwind difference method to discretize the spatial derivative term and the forward and backward Euler method to discretize the time derivative term, the explicit and implicit upwind difference schemes are obtained respectively. It is proved that the explicit upwind scheme is conditionally stable and the implicit upwind scheme is unconditionally stable. Then the convergence of the schemes is derived. Numerical examples verify the results of theoretical analysis.

Monolithic Coupling of the Pressure and Rigid Body Motion Equations in Computational Marine Hydrodynamics

In Fluid Structure Interaction(FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned manner by solving the rigid body motion equations once per nonlinear correction loop, updating the position of the body and solving the fluid flow equations in the new configuration. The partitioned approach requires a large number of nonlinear iteration loops per time–step. In order to enhance the coupling, a monolithic approach is proposed in Finite Volume(FV) framework,where the pressure equation and the rigid body motion equations are solved in a single linear system. The coupling is resolved by solving the rigid body motion equations once per linear solver iteration of the pressure equation, where updated pressure field is used to calculate new forces acting on the body, and by introducing the updated rigid body boundary velocity in to the pressure equation. In this paper the monolithic coupling is validated on a simple 2D heave decay case. Additionally, the method is compared to the traditional partitioned approach(i.e. "strongly coupled" approach) in terms of computational efficiency and accuracy. The comparison is performed on a seakeeping case in regular head waves, and it shows that the monolithic approach achieves similar accuracy with fewer nonlinear correctors per time–step. Hence, significant savings in computational time can be achieved while retaining the same level of accuracy.

ON THE ENERGY EQUATION FOR A SPATIALLY DEVELOPING INSTABILITY WAVE IN THE THEORY OF HYDRODYNAMIC STABILITY

In this article,it is shown that the energy equation for a spatially developing disturbance used in all the literatures dealing with the problem of hydrodynamic stability suffers from a small,but crucial error.

General solution of the modified Korteweg-de-Vries equation in the lattice hydrodynamic model

Traffic congestion is related to various density waves, which might be described by the nonlinear wave equations, such as the Burgers, Korteweg-de-Vries （KdV） and modified Korteweg-de-Vries （mKdV） equations. In this paper, the mKdV equations of four different versions of lattice hydrodynamic models, which describe the kink-antikink soliton waves are derived by nonlinear analysis. Furthermore, the general solution is given, which is applied to solving a new model -- the lattice hydrodynamic model with bidirectional pedestrian flow. The result shows that this general solution is consistent with that given by previous work.

GLOBAL WEAK SOLUTIONS TO ONE-DIMENSIONAL COMPRESSIBLE VISCOUS HYDRODYNAMIC EQUATIONS

In this article, we are concerned with the global weak solutions to the 1D com- pressible viscous hydrodynamic equations with dispersion correction δ2ρ（（φ（ρ））xxφ′（ρ））x with φ（ρ） = ρα. The model consists of viscous stabilizations because of quantum Fokker-Planck operator in the Wigner equation and is supplemented with periodic boundary and initial con- ditions. The diffusion term εuxx in the momentum equation may be interpreted as a classical conservative friction term because of particle interactions. We extend the existence result in [1] （α=1/2） to 0 〈 α ≤ 1. In addition, we perform the limit ε→0 with respect to 0 〈 α ≤1/2.

Time-dependent Ginzburg-Landau equation for lattice hydrodynamic model describing pedestrian flow

A thermodynamic theory is formulated to describe the phase transition and critical phenomena in pedestrian flow. Based on the extended lattice hydrodynamic pedestrian model taking the interaction of the next-nearest-neighbor persons into account, the time-dependent Ginzburg-Landau （TDGL） equation is derived to describe the pedestrian flow near the critical point through the nonlinear analysis method. The corresponding two solutions, the uniform and the kink solutions, are given. The coexisting curve, spinodal line, and critical point are obtained by the first and second derivatives of the thermodynamic potential.

Numerical Method for Solving Relativistic Hydrodynamic Equation Including Quark-Fragment Effect

In this paper, a numerical method is given to solve relativistic hydrodynamic equations with source terms by conservative finite difference scheme. In calculation, QGP (quark gluon plasma) phase transition is also considered. The numerical experiments have verified the effectiveness of the numerical method and some computational results are illustrated.

THE CHEBYSHEV PSEUDOSPECTRAL DOMAIN DECOMPO SITION METHOD FOR SOLVING TWO-DIMENSIONAL ELLIPTIC EQUATION

This paper ix devoted to establishment of the Chebyshev pseudospectral domain de-composition scheme for solving two-dimensional elliptic equation. By the generalized equivalent variatiunal form, we can get the stability and convergence of this new scheme.

Numerical solutions of two-dimensional nonlinear integral equations via Laguerre Wavelet method with convergence analysis

In this paper,the approximate solutions for two different type of two-dimensional nonlinear integral equations:two-dimensional nonlinear Volterra-Fredholm integral equations and the nonlinear mixed Volterra-Fredholm integral equations are obtained using the Laguerre wavelet method.To do this,these two-dimensional nonlinear integral equations are transformed into a system of nonlinear algebraic equations in matrix form.By solving these systems,unknown coefficients are obtained.Also,some theorems are proved for convergence analysis.Some numerical examples are presented and results are compared with the analytical solution to demonstrate the validity and applicability of the proposed method.