Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Solving Ordinary Differential Equations

SUMMARY Numerical calculations of the 2-D steady incompressible driven cavity flow are presented. The NavierStokes equations in streamfunction and vorticity formulation are solved numerically using a fine uniform grid mesh of 601 × 601. The steady driven cavity solutions are computed for Re ≤ 21,000 with a maximum absolute residuals of the governing equations that were less than 10 −10 . A new quaternary vortex at the bottom left corner and a new tertiary vortex at the top left corner of the cavity are observed in the flow field as the Reynolds number increases. Detailed results are presented and comparisons are made with benchmark solutions found in the literature.

Numerical solutions of 2‐D steady incompressible driven cavity flow at high Reynolds numbers

We consider the flapping stability and response of a thin two-dimensional flag of high extensional rigidity and low bending rigidity. The three relevant non-dimensional parameters governing the problem are the structure-to-fluid mass ratio, μ = ρ s h /(ρ f L); the Reynolds number, Re=VL/ν; and the non-dimensional bending rigidity, K B = EI / (ρfV 2 L 3 ). The soft cloth of a flag is represented by very low bending rigidity and the subsequent dominance of flow-induced tension as the main structural restoring force. We first perform linear analysis to help understand the relevant mechanisms of the problem and guide the computational investigation. To study the nonlinear stability and response, we develop a fluid-structure direct simulation (FSDS) capability, coupling a direct numerical simulation of the Navier-Stokes equations to a solver for thin-membrane dynamics of arbitrarily large motion. With the flow grid fitted to the structural boundary, external forcing to the structure is calculated from the boundary fluid dynamics. Using a systematic series of FSDS runs, we pursue a detailed analysis of the response as a function of mass ratio for the case of very low bending rigidity (K B = to-4) and relatively high Reynolds number (Re=10 3 ). We discover three distinct regimes of response as a function of mass ratio μ: (I) a small μ regime of fixed-point stability; (II) an intermediate μ regime of period-one limit-cycle flapping with amplitude increasing with increasing μ; and (III) a large μ regime of chaotic flapping. Parametric stability dependencies predicted by the linear analysis are confirmed by the nonlinear FSDS, and hysteresis in stability is explained with a nonlinear softening spring model. The chaotic flapping response shows up as a breaking of the limit cycle by inclusion of the 3/2 superharmonic. This occurs as the increased flapping amplitude yields a flapping Strouhal number (St=2Af/V) in the neighbourhood of the natural vortex wake Strouhal number, St ≃ 0.2. The limit-cycle von Karman vortex wake transitions in chaos to a wake with clusters of higher intensity vortices. For the largest mass ratios, strong vortex pairs are distributed away from the wake centreline during intermittent violent snapping events, characterized by rapid changes in tension and dynamic buckling.

https://www.cambridge.org/core/services/aop-cambridge-core/content/view/S0022112007005307

Flapping dynamics of a flag in a uniform stream

Finite differences provided the first numerical approach that permitted large-scale simulations in many applications areas, such as geophysical fluid dynamics. As accuracy and integration time requirements gradually increased, the focus shifted from finite differences to a variety of different spectral methods. During the last few years, radial basis functions, in particular in their ‘local’ RBF-FD form, have taken the major step from being mostly a curiosity approach for small-scale PDE ‘toy problems’ to becoming a major contender also for very large simulations on advanced distributed memory computer systems. Being entirely mesh-free, RBF-FD discretizations are also particularly easy to implement, even when local refinements are needed. This article gives some background to this development, and highlights some recent results.

Solving PDEs with radial basis functions

On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs

SUMMARY We note in this study that the Navier-Stokes equations, when expressed in streamfunction-vorticity fonn, can be approximated to fourth--order accuracy with stencils extending only over a 3 x 3 square of points. The key advantage of the new compact fourth-order scheme is that it allows direct iteration for low~to-mediwn Reynolds numbers. Numerical solutions are obtained for the model problem of the driven cavity and compared with solutions available in the literature. For Re $1500 point-SOR iteration is used and the convergence is fast.

/pdf/a-compact-fourth-order-finite-difference-scheme-for-the-2hdd1pul0e.pdf

A compact fourth‐order finite difference scheme for the steady incompressible Navier‐Stokes equations

In this paper, we extend a previous work on a compact scheme for the steady Navier–Stokes equations lLi, Tang, and Fornberg (1995), Int. J. Numer. Methods Fluids, 20, 1137–1151r to the unsteady case. By exploiting the coupling relation between the streamfunction and vorticity equations, the Navier–Stokes equations are discretized in space within a 3×3 stencil such that a fourth order accuracy is achieved. The time derivatives are discretized in such a way as to maintain the compactness of the stencil. We explore several known time-stepping approaches including second-order BDF method, fourth-order BDF method and the Crank–Nicolson method. Numerical solutions are obtained for the driven cavity problem and are compared with solutions available in the literature. For large values of the Reynolds number, it is found that high-order time discretizations outperform the low-order ones.

/pdf/a-compact-fourth-order-finite-difference-scheme-for-unsteady-g990xb5hia.pdf

Ming Li

Papers

A compact fourth‐order finite difference scheme for the steady incompressible Navier‐Stokes equations

A Compact Fourth-Order Finite Difference Scheme for Unsteady Viscous Incompressible Flows