The discontinuous Galerkin (DG) or local discontinuous Galerkin (LDG) method is a spatial discretization procedure for convection-diffusion equations, which employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes and limiters. The LaxWendroff time discretization procedure is an alternative method for time d...

The objective of this problem is to study ant foraging lines in one and two dimensions in order to understand their dynamics and examine the fitness of ant colonies. A simple model to describe ant traffic is formulated using conservation laws and a behavioral model. In this work, the behavioral model assumes that the concentration of pheromones is simply proportional to the density of ants. Thi...

This paper is devoted to the construction of numerical fluxes for hyperbolic systems. We first present a GFORCE numerical flux, which is a weighted average of the Lax-Friedrichs and Lax-Wendroff fluxes. For the linear advection equation with constant coefficient, the new flux reduces identically to that of the Godunov first order upwind method. Then we incorporate GFORCE in the framework of the...

This paper introduces the use of moving least-squares (MLS) approximations for the development of high-order finite volume discretizations on unstructured grids. The field variables and their successive derivatives can be accurately reconstructed using this mesh-free technique in a general nodal arrangement. The methodology proposed is used in the construction of two numerical schemes for the s...

We present a novel implicit scheme for the numerical solution of time-dependent conservation laws. The core idea presented method is to exploit and approximate mixed spatial-temporal derivative that occurs naturally when deriving some second order accurate schemes in time. Such an approach introduced context Lax-Wendroff (or Cauchy-Kowalevski) procedure time not completely replaced by space der...

The Kreiss matrix theorem asserts that a family of N X N matrices is L,-stable if and only if either a resolvent condition (R) or a Hennitian norm condition (H) is satisfied. We give a direct, considerahly shorter proof of the power-houndedness of an N X N matrix satisfying (R), sharpening former results by showing that powerhoundedness depends, at most, linearly on the dimension M. We also sho...

Harmonic wave excitation in a semi-infinite incompressible hyperelastic 1D rod with the Mooney–Rivlin equation of state reveals formation and propagation shock fronts arising between faster slower moving parts initially harmonic wave. The observed result collapse being absorbed by parts; hence, to attenuation kinetic elastic strain energy corresponding heat generation. Both geometrically physic...

Many of the recently developed high-resolution schemes for hyperbolic conservation laws are based on upwind di erencing. The building block of these schemes is the averaging of an approximate Godunov solver; its time consuming part involves the eld-byeld decomposition which is required in order to identify the \direction of the wind." Instead, we propose to use as a building block the more robu...

We study the stability of mesh refinement in space and time for several different interface equations and finite-difference approximations. First, we derive a root condition which implies stability for the initial-boundary value problem for this type of interface. From the root condition, we prove the stability of several interface equations using the maximum principle. In some cases, the final...