Note on Logarithmic Switchback Terms in Regular and Singular Perturbation Expansions

The occurrence of logarithmic switchback is studied for ordinary differential equations containing a parameter k which is allowed to take any value in a continuum of real numbers and with boundary conditions imposed at x = ε and x = ∞. Classical theory tells us that if the equation has a regular singular point at the origin there is a family of solutions which varies continuously with k, and the expansion around the origin has log x terms for a discrete set of values of k. It is shown here how nonlinearity enlarges this set so that it may even be dense in some interval of the real numbers. A log x term in the expansion in x leads to expansion coefficients containing log ε (switchback) in the perturbation expansion. If for a given value of k logarithmic terms in x and ε occur they may be obtained by continuity from neighboring values of k. Switchback terms occurred conspicuously in singular-perturbation solutions of problems posed for semi-infinite domain x ≥ ε. This connection is historical rather than logical. In particular we study here switchback terms for a specific example using methods of both singular and regular perturbations

Spectral stability of noncharacteristic isentropic Navier-Stokes boundary layers

Building on work of Barker, Humpherys, Lafitte, Rudd, and Zumbrun in the shock wave case, we study stability of compressive, or "shock-like", boundary layers of the isentropic compressible Navier-Stokes equations with gamma-law pressure by a combination of asymptotic ODE estimates and numerical Evans function computations. Our results indicate stability for gamma in the interval [1, 3] for all compressive boundary-layers, independent of amplitude, save for inflow layers in the characteristic limit (not treated). Expansive inflow boundary-layers have been shown to be stable for all amplitudes by Matsumura and Nishihara using energy estimates. Besides the parameter of amplitude appearing in the shock case, the boundary-layer case features an additional parameter measuring displacement of the background profile, which greatly complicates the resulting case structure. Moreover, inflow boundary layers turn out to have quite delicate stability in both large-displacement and large-amplitude limits, necessitating the additional use of a mod-two stability index studied earlier by Serre and Zumbrun in order to decide stability

Nonexistence of Chaotic Solutions of Nonlinear Differential Equations

We discuss some important issues arising from computational efforts in dynamical systems and fluid dynamics. Various individuals have misunderstood these issues since the onset of these problem areas; indeed, they have been routinely misinterpreted, and even viewed as "laws" by some. This paper hopes to stimulate appropriate corrections and to realign thinking, with the overall goal being sound future progress in dynamical systems and fluid dynamics.Comment: 12 page

Some analytical results for an algebraic flux correction scheme for a steady convection-diffusion equation in 1D

Algebraic flux correction schemes are nonlinear discretizations of convection dominated problems. In this work, a scheme from this class is studied for a steady-state convection--diffusion equation in one dimension. It is proved that this scheme satisfies the discrete maximum principle. Also, as it is a nonlinear scheme, the solvability of the linear subproblems arising in a Picard iteration is studied, where positive and negative results are proved. Furthermore, the non-existence of solutions for the nonlinear scheme is proved by means of counterexamples. Therefore, a modification of the method, which ensures the existence of a solution, is proposed. A weak version of the discrete maximum principle is proved for this modified method

Global Asymptotic Behavior of Iterative Implicit Schemes

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing three models of 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations (ODEs) is analyzed using the theory of dynamical systems. The iterative procedures include simple iteration and full and modified Newton iterations. The results are compared with standard Runge-Kutta explicit methods, a noniterative implicit procedure, and the Newton method of solving the steady part of the ODEs. Studies showed that aside from exhibiting spurious asymptotes, all of the four implicit LMMs can change the type and stability of the steady states of the differential equations (DEs). They also exhibit a drastic distortion but less shrinkage of the basin of attraction of the true solution than standard nonLMM explicit methods. The simple iteration procedure exhibits behavior which is similar to standard nonLMM explicit methods except that spurious steady-state numerical solutions cannot occur. The numerical basins of attraction of the noniterative implicit procedure mimic more closely the basins of attraction of the DEs and are more efficient than the three iterative implicit procedures for the four implicit LMMs. Contrary to popular belief, the initial data using the Newton method of solving the steady part of the DEs may not have to be close to the exact steady state for convergence. These results can be used as an explanation for possible causes and cures of slow convergence and nonconvergence of steady-state numerical solutions when using an implicit LMM time-dependent approach in computational fluid dynamics