Search results

The purpose of this paper is to both determine the effects of the nonlinearity on the wave dynamics and assess the temporal and spatial accuracy of five finite difference methods for the solution of the inviscid generalized regularized long-wave (GRLW) equation subject to initial Gaussian conditions.

Two implicit second- and fourth-order accurate finite difference methods and three Runge-Kutta procedures are introduced. The methods employ a new dependent variable which contains the wave amplitude and its second-order spatial derivative. Numerical experiments are reported for several temporal and spatial step sizes in order to assess their accuracy and the preservation of the first two invariants of the inviscid GRLW equation as functions of the spatial and temporal orders of accuracy, and thus determine the conditions under which grid-independent results are obtained.

It has been found that the steepening of the wave increase as the nonlinearity exponent is increased and that the accuracy of the fourth-order Runge-Kutta method is comparable to that of a second-order implicit procedure for time steps smaller than 100th, and that only the fourth-order compact method is almost grid-independent if the time step is on the order of 1,000th and more than 5,000 grid points are used, because of the initial steepening of the initial profile, wave breakup and solitary wave propagation.

This is the first study where an accuracy assessment of wave breakup of the inviscid GRLW equation subject to initial Gaussian conditions is reported.

The purpose of this paper is to determine both analytically and numerically the existence of smooth, cusped and sharp shock wave solutions to a one-dimensional model of microfluidic droplet ensembles, water flow in unsaturated flows, infiltration, etc., as functions of the powers of the convection and diffusion fluxes and upstream boundary condition; to study numerically the evolution of the wave for two different initial conditions; and to assess the accuracy of several finite difference methods for the solution of the degenerate, nonlinear, advection--diffusion equation that governs the model.

The theory of ordinary differential equations and several explicit, finite difference methods that use first- and second-order, accurate upwind, central and compact discretizations for the convection terms are used to determine the analytical solution for steadily propagating waves and the evolution of the wave fronts from hyperbolic tangent and piecewise linear initial conditions to steadily propagating waves, respectively. The amplitude and phase errors of the semi-discrete schemes are determined analytically and the accuracy of the discrete methods is assessed.

For non-zero upstream boundary conditions, it has been found both analytically and numerically that the shock wave is smooth and its steepness increases as the power of the diffusion term is increased and as the upstream boundary value is decreased. For zero upstream boundary conditions, smooth, cusped and sharp shock waves may be encountered depending on the powers of the convection and diffusion terms. For a linear diffusion flux, the shock wave is smooth, whereas, for a quadratic diffusion flux, the wave exhibits a cusped front whose left spatial derivative decreases as the power of the convection term is increased. For higher nonlinear diffusion fluxes, a sharp shock wave is observed. The wave speed decreases as the powers of both the convection and the diffusion terms are increased. The evolution of the solution from hyperbolic tangent and piecewise linear initial conditions shows that the wave back adapts rapidly to its final steady value, whereas the wave front takes much longer, especially for piecewise linear initial conditions, but the steady wave profile and speed are independent of the initial conditions. It is also shown that discretization of the nonlinear diffusion flux plays a more important role in the accuracy of first- and second-order upwind discretizations of the convection term than either a conservative or a non-conservative discretization of the latter. Second-order upwind and compact discretizations of the convection terms are shown to exhibit oscillations at the foot of the wave’s front where the solution is nil but its left spatial derivative is largest. The results obtained with a conservative, centered second--order accurate finite difference method are found to be in good agreement with those of the second-order accurate, central-upwind Kurganov--Tadmor method which is a non-oscillatory high-resolution shock-capturing procedure, but differ greatly from those obtained with a non-conservative, centered, second-order accurate scheme, where the gradients are largest.

A new, one-dimensional model for microfluidic droplet transport, water flow in unsaturated flows, infiltration, etc., that includes high-order convection fluxes and degenerate diffusion, is proposed and studied both analytically and numerically. Its smooth, cusped and sharp shock wave solutions have been determined analytically as functions of the powers of the nonlinear convection and diffusion fluxes and the boundary conditions. These solutions are used to assess the accuracy of several finite difference methods that use different orders of accuracy in space, and different discretizations of the convection and diffusion fluxes, and can be used to assess the accuracy of other numerical procedures for one-dimensional, degenerate, convection--diffusion equations.

The purpose of this paper is to evaluate the performance of a numerical method for the solution to shallow-water equations on a staggered grid, in simulations for shear instabilities at two convective Froude numbers.

The simulations start from a small perturbation to a base flow with a hyperbolic-tangent velocity profile. The subsequent development of the shear instabilities is studied from the simulations using a number of flux-limiting schemes, including the second-order MINMOD, the third-order ULTRA-QUICK and the fifth-order WENO schemes for the spatial interpolation of the nonlinear fluxes. The fourth-order Runge-Kutta method advances the simulation in time.

The simulations determine two parameters: the fractional growth rate of the linear instabilities; and the vorticity thickness of the first nonlinear peak. Grid refinement using 32, 64, 128, 256 and 512 nodes over one wave length determines the exact values by extrapolation and the computational error for the parameters. It also determines the overall order of convergence for each of the flux-limiting schemes used in the numerical simulations.

The four-digit accuracy of the numerical simulations presented in this paper are comparable to analytical solutions. The development of this reliable numerical simulation method has paved the way for further study of the instabilities in shear flows that radiate waves.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

The purpose of this paper is to study the new (3 + 1)-dimensional integrable fourth-order nonlinear equation which is used to model the shallow water waves.

By means of the Cole–Hopf transform, the bilinear form of the studied equation is extracted. Then the ansatz function method combined with the symbolic computation is implemented to construct the breather, multiwave and the interaction wave solutions. In addition, the subequation method tis also used to search for the diverse travelling wave solutions.

The breather, multiwave and the interaction wave solutions and other wave solutions like the singular periodic wave structure and dark wave structure are obtained. To the author’s knowledge, the solutions obtained are all new and have never been reported before.

The solutions obtained in this work have never appeared in other literature and can be regarded as an extension of the solutions for the new (3 + 1)-dimensional integrable fourth-order nonlinear equation.

This paper aims to develop a new 3-level implicit numerical method of order 2 in time and 4 in space based on half-step cubic polynomial approximations for the solution of 1D quasi-linear hyperbolic partial differential equations. The method is derived directly from the consistency condition of spline function which is fourth-order accurate. The method is directly applied to hyperbolic equations, irrespective of coordinate system, and fourth-order nonlinear hyperbolic equation, which is main advantage of the work.

In this method, three grid points for the unknown function w(x,t) and two half-step points for the known variable x in spatial direction are used. The methodology followed in this paper is construction of a cubic spline polynomial and using its continuity properties to obtain fourth-order consistency condition. The proposed method, when applied to a linear equation is shown to be unconditionally stable. The technique is extended to solve system of quasi-linear hyperbolic equations. To assess the validity and accuracy, the method is applied to solve several benchmark problems, and numerical results are provided to demonstrate the usefulness of the method.

The paper provides a fourth-order numerical scheme obtained directly from fourth-order consistency condition. In earlier methods, consistency conditions were only second-order accurate. This brings an edge over other past methods. In addition, the method is directly applicable to physical problems involving singular coefficients. Therefore, no modification in the method is required at singular points. This saves CPU time, as well as computational costs.

There are no limitations. Obtaining a fourth-order method directly from consistency condition is a new work. In addition, being an implicit method, this method is unconditionally stable for a linear test equation.

Physical problems with singular and nonsingular coefficients are directly solved by this method.

The paper develops a new fourth-order implicit method which is original and has substantial value because many benchmark problems of physical significance are solved in this method.

The purpose of this paper is to determine both analytically and numerically the solution to a new one-dimensional equation for the propagation of small-amplitude waves in shallow waters that accounts for linear and nonlinear drift, diffusive attenuation, viscosity and dispersion, its dependence on the initial conditions, and its linear stability.

An implicit, finite difference method valid for both parabolic and second-order hyperbolic equations has been used to solve the equation in a truncated domain for five different initial conditions, a nil initial first-order time derivative and relaxation times linearly proportional to the viscosity coefficient.

A fast transition that depends on the coefficient of the linear drift, the diffusive attenuation and the power of the nonlinear drift are found for initial conditions corresponding to the exact solution of the generalized regularized long-wave equation. For initial Gaussian, rectangular and triangular conditions, the wave’s amplitude and speed increase as both the amplitude and the width of these conditions increase and decrease, respectively; wide initial conditions evolve into a narrow leading traveling wave of the pulse type and a train of slower oscillatory secondary ones. For the same initial mass and amplitude, rectangular initial conditions result in larger amplitude and velocity waves of the pulse type than Gaussian and triangular ones. The wave’s kinetic, potential and stretching energies undergo large changes in an initial layer whose thickness is on the order of the diffusive attenuation coefficient.

A new, one-dimensional equation for the propagation of small-amplitude waves in shallow waters is proposed and studied analytically and numerically. The equation may also be used to study the displacement of porous media subject to seismic effects, the dispersion of sound in tunnels, the attenuation of sound because of viscosity and/or heat and mass diffusion, the dynamics of second-order, viscoelastic fluids, etc., by appropriate choices of the parameters that appear in it.

This paper aims to introduce a new (3 + 1)-dimensional fourth-order integrable equation characterized by second-order derivative in time t. The new equation models both right- and left-going waves in a like manner to the Boussinesq equation.

This formally uses the simplified Hirota’s method and lump schemes for determining multiple soliton solutions and lump solutions, which are rationally localized in all directions in space.

This paper confirms the complete integrability of the newly developed (3 + 1)-dimensional model in the Painevé sense.

This paper addresses the integrability features of this model via using the Painlevé analysis.

This paper presents a variety of lump solutions via using a variety of numerical values of the included parameters.

This work formally furnishes useful algorithms for extending integrable equations and for the determination of lump solutions.

To the best of the author’s knowledge, this paper introduces an original work with newly developed integrable equation and shows useful findings of solitons and lump solutions.

The purpose of this paper is to apply the exp‐function method to construct exact solutions of nonlinear wave equations. The proposed technique is tested on the (2+1) and (3+1) dimensional extended shallow water wave equations. These equations play a very important role in mathematical physics and engineering sciences.

In this paper, the authors apply the exp‐function method to construct exact solutions of nonlinear wave equations.

In total, four forms of the extended shallow water wave equation have been studied, from the point of view of its exact solutions using computational method. Exp‐function method was employed to achieve the goal set for this work. The applied method will be used in further works to establish more entirely new solutions for other kinds of nonlinear wave equations. Finally, it is worthwhile to mention that the proposed method is straightforward, concise, and it is a promising and powerful new method for other nonlinear wave equations in mathematical physics.

The algorithm suggested in the paper is quite efficient and is practically well suited for use in these problems. The method is straightforward and concise, and its applications are promising.