EFFICIENT NUMERICAL METHODS FOR THE KDV EQUATION

In this work, the dynamic response of a cable-stayed bridge that consists of a simply supported four-cable-stayed deck beam and two rigid towers, subjected to a moving vehicle load, is studied. The vehicle is modeled as a mass–spring–damper system moving at a constant velocity, which is assumed to keep contact with the deck beam at all times. Convective velocity and acceleration terms associated with the moving vehicle as it traverses along the vibrating deck beam are taken into consideration, as well as geometric nonlinearities of stay cables. The nonlinear response of the cable-stayed bridge is obtained by solving nonlinear and linear partial differential equations which govern transverse and longitudinal vibrations of stay cables and transverse vibrations of segments of the deck beam, respectively, along with their boundary and matching conditions. Orthogonality relations of exact mode shapes of the linearized undamped cable-stayed bridge model are employed to convert the coupled nonlinear partial differential equations of the original nonlinear cable-stayed bridge model to a set of ordinary differential equations by using the Galerkin method. The dynamic response of the cable-stayed bridge is calculated using the Runge–Kutta–Fehlberg method in MATLAB. Convergence of the dynamic response from the Galerkin method is investigated for two cases in which the velocities or masses of the moving vehicle are different. Results show that an accurate calculation of the dynamic response of the cable-stayed bridge needs use of a large number of modes of the linearized undamped cable-stayed bridge model, and needs many more modes for the deck beam than stay cables. Moreover, effects of the velocity and mass of the moving vehicle and the convective terms on the dynamic response of the cable-stayed bridge are studied with convergent Galerkin truncation.

The standard approach to applying IRK methods in the solution of two-point boundary value problems involves the solution of a non-linear system ofn×s equations in order to calculate the stages of the method, wheren is the number of differential equations ands is the number of stages of the implicit Runge-Kutta method. For two-point boundary value problems, we can select a subset of the implicit Runge-Kutta methods that do not require us to solve a non-linear system; the calculation of the stages can be done explicitly, as is the case for explicit Runge-Kutta methods. However, these methods have better stability properties than the explicit Runge-Kutta methods. We have called these new formulas two-point explicit Runge-Kutta (TPERK) methods. Their most important property is that, because their stages can be computed explicity, the solution of a two-point boundary value problem can be computed more efficiently than is possible using an implicit Runge-Kutta method. We have also developed a symmetric subclass of the TPERK methods, called ATPERK methods, which exhibit a number of useful properties.

This paper investigates a class of explicit pseudo three-step Runge-Kutta methods for arbitrarily high order nonstiff initial value problems for systems of first-order differential equations. By using collocation techniques and by suitably choosing collocation points we can obtain a stable \(s\)-stage explicit pseudo three-step Runge-Kutta method (EPThRK method) of order \(p=2s\) requiring only one effective sequential f-evaluation per step on \(s\)-processor parallel computers. By a few widely-used test problems, we show the superiority of the new EPThRK methods proposed in this paper over red well-known parallel PIRK codes and efficient sequential ODEX, DOPRI5 and DOP853 codes available in the literature.

In this paper, an application of symplectic integrators for solving an important problem of non-linear dynamics is considered. Numerical algorithms for solving six-body problem are described, built on explicit second order Runge-Kutta method, implicit Euler method, symplectic Verlet method and D-method. The examined six-body problem time-domain solution is given, as well as the system's full energy dynamics for the implemented numerical algorithm. Conclusions of symplectic integrators efficiency are made based on simulation results.

We consider dispersive shock waves of the focusing nonlinear Schrödinger equation generated by discontinuous initial conditions which are periodic or quasiperiodic on the left semiaxis and zero on the right semiaxis. As an initial function, we use a finite-gap potential of the Dirac operator given in an explicit form through hyperelliptic theta-functions. The aim of this paper is to study the long-time asymptotics of the solution of this problem in a vicinity of the leading edge, where a train of asymptotic solitons are generated. Such a problem was studied in the work of Kotlyarov and Khruslov [Teor. Mat. Fiz. 68(2), 751–761 (1986)] and Kotlyarov {Mat. Zametki 49(2), 84–94 (1991) [Math. Notes 49(1-2), 172–180 (1991)]} using Marchenko’s inverse scattering techniques. We investigate this problem exceptionally using the Riemann-Hilbert (RH) problem techniques that allow us to obtain explicit formulas for asymptotic solitons themselves in contrast with the cited papers where asymptotic formulas are obtained only for the square of the absolute value of solution. Using transformations of the main RH problems, we arrive at a model problem corresponding to the parametrix at the end points of the continuous spectrum of the Zakharov-Shabat spectral problem. The parametrix problem is effectively solved in terms of the generalized Laguerre polynomials, which naturally appeared after appropriate scaling of the Riemann-Hilbert problem in small neighborhoods of the end points of the continuous spectrum. Further asymptotic analysis gives an explicit formula for solitons at the edge of dispersive waves. Thus, we give the complete description of the train of asymptotic solitons: not only bearing the envelope of each asymptotic soliton, but its oscillating structure is found explicitly. Besides, the second term of asymptotics describing an interaction between these solitons and oscillating background is also found. This gives the fine structure of the edge of dispersive shock waves.

It is possible to construct fully implicit Runge-Kutta methods like Gaus-Legendre, Radau-IA, Radau-IIA, Lobatto-IIIA, -IIIB, and -IIIC methods of arbitrary high order of convergence. The aim of this paper is to find a new adaptive time stepping for these classes which is based on the embedding technique. Adaptive time step control with embedding is well-known for Runge-Kutta methods, and therefore new embedded methods of order s-1 for the above classes of fully implicit Runge-Kutta methods are constructed. Since these fully implicit methods need the solution of a huge non-linear system of equations different approaches for non-linear equations are discussed and compared. It can be observed that non-linear solvers like the usually used simplified Newton method have a step size restriction if they are applied on higher order methods. We apply our new methods on some lower dimensional ODEs to show that our approach leads to an efficient method.

Many optical and other nonlinear media are governed by dispersive, or diffractive, wave equations, for which initial jump discontinuities are resolved into a dispersive shock wave. The dispersive shock wave smooths the initial discontinuity and is a modulated wavetrain consisting of solitary waves at its leading edge and linear waves at its trailing edge. For integrable equations the dispersive shock wave solution can be found using Whitham modulation theory. For nonlinear wave equations which are hyperbolic outside the dispersive shock region, the amplitudes of the solitary waves at the leading edge and the linear waves at the trailing edge of the dispersive shock can be determined. In this paper an approximate method is presented for calculating the amplitude of the lead solitary waves of a dispersive shock for general nonlinear wave equations, even if these equations are not hyperbolic in the dispersionless limit. The approximate method is validated using known dispersive shock solutions and then applied to calculate approximate dispersive shock solutions for equations governing nonlinear optical media, such as nematic liquid crystals, thermal glasses and colloids. These approximate solutions are compared with numerical results and excellent comparisons are obtained.

In the previous chapter, we introduced explicit Runge-Kutta (ERK) methods and demonstrated how they can be implemented as a hierarchy of Python classes. For most ODE systems, replacing the simple forward Euler method with a higher-order ERK method will significantly reduce the number of time steps needed to reach a specified accuracy. Furthermore, it often leads to reduced computation time, since the additional cost per time step is outweighed by the reduced number of steps. However, there exists a class of ODEs known as stiff systems, where all the ERK methods require very small time steps, and any attempt to increase the time step leads to spurious oscillations and possible divergence of the solution. Stiff ODE systems pose a challenge for explicit methods, and they are better addressed by implicit solvers such as implicit Runge-Kutta (IRK) methods. IRK methods are well-suited for stiff problems and can offer substantial reductions in computation time when tackling challenging problems.

We present an algorithm for determining the stepsize in an explicit Runge-Kutta method that is suitable when solving moderately stiff differential equations. The algorithm has a geometric character, and is based on a pair of semicircles that enclose the boundary of the stability region in the left half of the complex plane. The algorithm includes an error control device. We describe a vectorized form of the algorithm, and present a corresponding MATLAB code. Numerical examples for Runge-Kutta methods of third and fourth order demonstrate the properties and capabilities of the algorithm.

Two algorithms used to correct the semi-implicit numerical scheme for time integration of the conservation equations incorporated into the KORSAR/GP computer code two-liquid model are presented. The KORSAR/GP code has been developed jointly by specialists of the Federal State Unitary Enterprise Aleksandrov Research Institute of Technology and Gidropress Experimental Design Office; in 2009, the code was certified at the Federal Service for Environmental, Technological, and Nuclear Supervision (Rostekhnadzor) as applied to numerical safety assessment of VVER-based reactor plants. In the semi-implicit scheme, the convective terms appearing in the phase momentum conservation equations are written in an explicit form. The mass and energy flows of the phases are represented implicitly with respect to phase velocities, and the transferred donor quantities are calculated based on the parameters taken from the previous time layer. Owing to linearization of unsteady and source terms, the linear system of finite difference equations is solved in a noniteration manner. The first of the presented algorithms compensates for numerical imbalances of phase masses and energies resulting from linearization of unsteady terms appearing in the discrete equations and ensures conservativeness of the scheme. The second algorithm introduces correction for the nonphysical redistribution of coolant mass and energy over the computational cells if a change occurs in the phase motion direction within a time step when the scheme of approximating the convective terms becomes “antidonor” with respect to flow. The numerical imbalances and refinements taking into account the correction of donor quantities are computed at each time step according to the proposed correlations and are used for making compensation at the next time step as supplementary sources in the conservation equations. Results from testing the correction algorithms in the KORSAR/GP code are presented. The effectiveness of the imbalance compensation algorithm is confirmed by solving a problem involving the natural circulation loop heating process. The adequacy of the “antidonor” scheme correction algorithm is demonstrated on problems involving steplike initial distribution of the scalar parameters of a stagnant single-phase gaseous or water coolant in a horizontal tube when the flowrate at the tube inlet has different signs at different time steps.

This paper is aimed at analyzing the consistency and convergence of a derived explicit fourth-stage fourth-order Runge-Kutta method. The analysis revealed that the method is consistent and convergent. The implementation of this method on initial-value problems was done in previous paper, and it revealed that the method compared favorably well with the existing classical fourth stage fourth order explicit Runge Kutta method.

The three-dimensional incompressible Navier-Stokes equations in conservative form are discretized with the finite volume method, and are solved in a time-accurate manner, using the pseudo-compressibility technique. A third-upwind differencing scheme, based on the Roe’s flux-splitting, is used to compute the convective term in the MUSCL scheme way which satisfies the TVD condition so as to avoid excessive numerical diffusivity. The viscous term is evaluated by a Gauss integration in the finite volume method way, resulting in the center differencing. An implicit factorization approximate method, based on the Euler implicit scheme, is used for the discretization of time derivative. In the time-accurate formulation, the time derivatives of velocity components in the momentum equations are differenced using a first-order, two-point, backward-difference formula. To satisfy the equation of continuity at each physical time step, the artificial compressibility relation is also introduced and several subiterations in pseudotime are required. The cell-centered layout is adopted, in which flow variable node is placed at the center of grid cell and the grid cell is used as control volume.

A unified gas-kinetic scheme for axisymmetric flow in all Knudsen number regimes

Variational formulations for explicit Runge-Kutta Methods

Explicit Runge–Kutta methods and implicit multistep methods utilizing a Newton–Krylov nonlinear solver are evaluated for a range of configurations of the shallow-water dynamical core of the spectral element community atmosphere model to evaluate their computational performance. These configurations are designed to explore the attributes of each method under different but relevant model usage scenarios including varied spectral order within an element, static regional refinement, and scaling to the largest problem sizes. This analysis is performed within the shallow-water dynamical core option of a full climate model code base to enable a wealth of simulations for study, with the aim of informing solver development within the more complete hydrostatic dynamical core used for climate research. The limitations and benefits to using explicit versus implicit methods, with different parameters and settings, are discussed in light of the trade-offs with Message Passing Interface (MPI) communication and memory and their inherent efficiency bottlenecks. Given the performance behavior across the configurations analyzed here, the recommendation for future work using the implicit solvers is conditional based on scale separation and the stiffness of the problem. For the regionally refined configurations, the implicit method has about the same efficiency as the explicit method, without considering efficiency gains from a preconditioner. The potential for improvement using a preconditioner is greatest for higher spectral order configurations, where more work is shifted to the linear solver. Initial simulations with OpenACC directives to utilize a Graphics Processing Unit (GPU) when performing function evaluations show improvements locally, and that overall gains are possible with adjustments to data exchanges.