• Coupled systems mechanics
• /
• v.1 no.4
• /
• pp.325-343
• /
• 2012

We present different numerical methods for solving the shallow shelf equations with basal drag (SSAB). An alternative approach of splitting the SSAB equation into a Laplacian and diagonal shift operator is discussed with respect to the underlying eigenvalue problem. First, we solve the equations using standard methods. Then, the coupled equations are decomposed into operators for membranes stresses, basal shear stress and driving stress. Applying reasonable parameter values, we demonstrate that the operator of the membrane stresses is much stiffer than the operator of the basal shear stress. Here, we could apply a new splitting method, which alternates between the iteration on the membrane-stress operator and the basal-shear operator, with a more frequent iteration on the operator of the membrane stresses. We show that this splitting accelerates and stabilize the computational performance of the numerical method, although an appropriate choice of the standard method used to solve for all operators in one step speeds up the scheme as well.

• Journal of Environmental Science International
• /
• v.11 no.9
• /
• pp.907-914
• /
• 2002

Various Eulerian-Lagrangian models for the one-dimensional longitudinal dispersion equation in nonuniform flow were studied comparatively. In the models studied, the transport equation was decoupled into two component parts by the operator-splitting approach; one part is governing advection and the other is governing dispersion. The advection equation has been solved by using the method of characteristics following fluid particles along the characteristic line and the results were interpolated onto an Eulerian grid on which the dispersion equation was solved by Crank-Nicholson type finite difference method. In the solution of the advection equation, Lagrange fifth, cubic spline, Hermite third and fifth interpolating polynomials were tested by numerical experiment and theoretical error analysis. Among these, Hermite interpolating polynomials are generally superior to Lagrange and cubic spline interpolating polynomials in reducing both dissipation and dispersion errors.

• Korean Journal of Hydrosciences
• /
• v.6
• /
• pp.51-66
• /
• 1995

Various Eulerian-Lagerangian numerical models for the one-dimensional longtudinal dispersion equation are studied comparatively. In the models studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing advection and the other dispersion. The advection equation has been solved using the method of characteristics following flud particles along the characteristic line and the result are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpo;ation po;ynomials are superor to Lagrange interpolation polynomials in reducing both dissipation and dispersion errors.

• Water for future
• /
• v.27 no.2
• /
• pp.155-166
• /
• 1994

Various Eulerian-Lagrangian numerical models for the one-dimensional longitudinal dispersion equation are studied comparatively. In the model studied, the transport equation is decoupled into two component parts by the operator-splitting approach ; one part governing adveciton and the other dispersion. The advection equation has been solved using the method of characteristics following fluid particles along the characteristic line and the results are interpolated onto an Eulerian grid on which the dispersion equation is solved by Crank-Nicholson type finite difference method. In solving the advection equation, various interpolation schemes are tested. Among those, Hermite interpolation polynomials are superior to Lagrange interpolation polynomials in reducing dissipation and dispersion errors in the simulation.

• Water for future
• /
• v.28 no.5
• /
• pp.151-162
• /
• 1995

A two-dimensional stream tube dispersion model is developed to simulate accurately transverse dispersion processes of pollutants in natural channels. Two distinct features of the stream tube dispersion model derived herein are that it employs the transverse cumulative discharge as an independent variable replacing the transverse distance and that it is developed in a natural coordinate system which follows the general direction of the channel flow. In the model studied, Eulerian-Lagrangian method is used to solve the stream tube dispersion equation. The stream tube dispersion equation is decoupled into two components by the operator-splitting approach; one is governing advection and the other is governing dispersion. The advection equation has been solved using the method of characteristics and the results are interpolated onto Eulerian grid on which the dispersion equation is solved by centered difference method. In solving the advection equation, cubic spline interpolating polynomials is used. In the present study, the results of the application of this model to a natural channel are compared with a steady-state flow measurements. Simulation results are in good accordance with measured data.