A one-equation turbulence model for aerodynamic flows

Preface 1. Introduction 2. Conservation laws and differential equations 3. Characteristics and Riemann problems for linear hyperbolic equations 4. Finite-volume methods 5. Introduction to the CLAWPACK software 6. High resolution methods 7. Boundary conditions and ghost cells 8. Convergence, accuracy, and stability 9. Variable-coefficient linear equations 10. Other approaches to high resolution 11. Nonlinear scalar conservation laws 12. Finite-volume methods for nonlinear scalar conservation laws 13. Nonlinear systems of conservation laws 14. Gas dynamics and the Euler equations 15. Finite-volume methods for nonlinear systems 16. Some nonclassical hyperbolic problems 17. Source terms and balance laws 18. Multidimensional hyperbolic problems 19. Multidimensional numerical methods 20. Multidimensional scalar equations 21. Multidimensional systems 22. Elastic waves 23. Finite-volume methods on quadrilateral grids Bibliography Index.

https://depts.washington.edu/clawpack/book2/sample.pdf

Finite Volume Methods for Hyperbolic Problems

A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.

Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes

This Key Note presents a summary of the development of the Finite Element Method in the field of Electromagnet ic Engineering, together with a description of several contributions of the authors to the Finite Element Method and its application to the solution of electromagnetic problems. First, a self-adaptive mesh scheme is presented in the context of the quasi-static and full-wave analysis of general anisotropic multiconductor arbitrary shaped waveguiding structures. A comparison between two a posteriori error estimates is done. The first one is based on the complete residual of the differential equations defining the problem. The second one is based on a recovery or smoothing technique of the electromagnetic field. Next, an implementation of the first family of Nedelec's curl-conforming elements done by the authors is outlined. Its features are highlighted and compared with other curl-conforming elements. A presentation of an iterative procedure using a numerically exact radiation condition for the analysis of open (scattering and radiation) problems follows. Other contributions of the authors, like the use of wavelet like basis functions and an implementation of a Time Domain Finite Element Method in the context of two-dimensional scattering problems are only mentioned due to the lack of space.

The finite element method in electromagnetics

A new set of benchmarks has been developed for the performance evaluation of highly parallel supercom puters. These consist of five "parallel kernel" bench marks and three "simulated application" benchmarks. Together they mimic the computation and data move ment characteristics of large-scale computational fluid dynamics applications. The principal distinguishing feature of these benchmarks is their "pencil and paper" specification-all details of these benchmarks are specified only algorithmically. In this way many of the difficulties associated with conventional bench- marking approaches on highly parallel systems are avoided.

/pdf/the-nas-parallel-benchmarks-1v48ldcocu.pdf

The Nas Parallel Benchmarks

Preconditioners for hyperbolic systems are numerical artifices intended to accelerate the convergence path to a steady state. In addition, in some cases, the preconditioner can also be included in the artificial viscosity/upwinding terms in order to improve the accuracy of the steady state solution. For time dependent problems we use a dual time stepping approach; and therefore the preconditioner affects the convergence rate and accuracy of the subiterations at each physical time step. We consider two types of local preconditioners that couple the governing equations at a node point: Jacobi and low-speed preconditioning. We consider their effectiveness for both steady state and time dependent problems with regard to the convergence rate and the numerical accuracy. We also consider the effect of the far field boundary conditions on both steady state and time dependent problems.

Assessment of local preconditioners for steady state and time dependent flows

Artificial viscosity is added either implicity or explicitly in practically every numerical scheme for suppressing spurious oscillations in the solution of fluid-dynamics equations. In the present central-difference scheme, artificial viscosity is added explicitly for suppressing high-frequency oscillations and achieving good convergence properties. The amount of artificial viscosity added is controlled through the use of preselected coefficients. In the standard scheme, scalar coefficients based on the spectral radii of the Jacobian of the convective fluxes are used. However, this can add too much viscosity to the slower waves. Hence, the use of matrix-valued coefficients, which give appropriate viscosity for each wave component, is suggested. With the matrix-valued coefficients, the central-difference scheme produces more accurate solutions on a given grid, particularly in the vicinity of shocks and boundary layers, while still maintaining good convergence properties.

Effect of artificial viscosity on three-dimensional flow solutions

Acceleration methods for multi-physics compressible flow

Most surviving biblical period Hebrew inscriptions are ostraca—ink-on-clay texts. They are poorly preserved and once unearthed, fade rapidly. Therefore, proper and timely documentation of ostraca is essential. Here we show a striking example of a hitherto invisible text on the back side of an ostracon revealed via multispectral imaging. This ostracon, found at the desert fortress of Arad and dated to ca. 600 BCE (the eve of Judah’s destruction by Nebuchadnezzar), has been on display for half a century. Its front side has been thoroughly studied, while its back side was considered blank. Our research revealed three lines of text on the supposedly blank side and four "new" lines on the front side. Our results demonstrate the need for multispectral image acquisition for both sides of all ancient ink ostraca. Moreover, in certain cases we recommend employing multispectral techniques for screening newly unearthed ceramic potsherds prior to disposal.

/pdf/multispectral-imaging-reveals-biblical-period-inscription-b2rg7s7m8r.pdf

Eli Turkel

Papers

Assessment of local preconditioners for steady state and time dependent flows

Effect of artificial viscosity on three-dimensional flow solutions

Acceleration methods for multi-physics compressible flow

Multispectral imaging reveals biblical-period inscription unnoticed for half a century

Numerical Solution of the Shallow-Water Equations.