Shock-resolved Navier–Stokes simulation of the Richtmyer–Meshkov instability start-up at a light–heavy interface

The single-mode Richtmyer–Meshkov instability is investigated using a first-order perturbation of the two-dimensional Navier–Stokes equations about a one-dimensional unsteady shock-resolved base flow. A feature-tracking local refinement scheme is used to fully resolve the viscous internal structure of the shock. This method captures perturbations on the shocks and their influence on the interface growth throughout the simulation, to accurately examine the start-up and early linear growth phases of the instability. Results are compared to analytic models of the instability, showing some agreement with predicted asymptotic growth rates towards the inviscid limit, but significant discrepancies are noted in the transient growth phase. Viscous effects are found to be inadequately predicted by existing models

An exploration of the IGA method for efficient reservoir simulation

Novel numerical methods present exciting opportunities to improve the efficiency of reservoir simulators. Because potentially significant gains to computational speed and accuracy may be obtained, it is worthwhile explore alternative computational algorithms for both general and case-by-case application to the discretization of the equations of porous media flow, fluid-structure interaction, and/or production. In the present work, the fairly new concept of isogeometric analysis (IGA) is evaluated for its suitability to reservoir simulation via direct comparison with the industry standard finite difference (FD) method and 1st order standard finite element method (SFEM). To this end, two main studies are carried out to observe IGA’s performance with regards to geometrical modeling and ability to capture steep saturation fronts. The first study explores IGA’s ability to model complex reservoir geometries, observing L2 error convergence rates under a variety of refinement schemes. The numerical experimental setup includes an 'S' shaped line sink of varying curvature from which water is produced in a 2D homogenous domain. The accompanying study simplifies the domain to 1D, but adds in multiphase physics that traditionally introduce difficulties associated with modeling of a moving saturation front. Results overall demonstrate promise for the IGA method to be a particularly effective tool in handling geometrically difficult features while also managing typically challenging numerical phenomena.Petroleum and Geosystems Engineerin

Status of research at the Institute for Computer Applications in Science and Engineering (ICASE)

Research conducted at the Institute for Computer Applications in Science and Engineering in applied mathematics, numerical analysis and computer science is summarized

Characteristic Evolution and Matching

I review the development of numerical evolution codes for general relativity based upon the characteristic initial value problem. Progress in characteristic evolution is traced from the early stage of 1D feasibility studies to 2D axisymmetric codes that accurately simulate the oscillations and gravitational collapse of relativistic stars and to current 3D codes that provide pieces of a binary black hole spacetime. Cauchy codes have now been successful at simulating all aspects of the binary black hole problem inside an artificially constructed outer boundary. A prime application of characteristic evolution is to extend such simulations to null infinity where the waveform from the binary inspiral and merger can be unambiguously computed. This has now been accomplished by Cauchy-characteristic extraction, where data for the characteristic evolution is supplied by Cauchy data on an extraction worldtube inside the artificial outer boundary. The ultimate application of characteristic evolution is to eliminate the role of this outer boundary by constructing a global solution via Cauchy-characteristic matching. Progress in this direction is discussed.Comment: New version to appear in Living Reviews 2012. arXiv admin note: updated version of arXiv:gr-qc/050809

Institute for Computational Mechanics in Propulsion (ICOMP) fourth annual review, 1989

The Institute for Computational Mechanics in Propulsion (ICOMP) is operated jointly by Case Western Reserve University and the NASA Lewis Research Center. The purpose of ICOMP is to develop techniques to improve problem solving capabilities in all aspects of computational mechanics related to propulsion. The activities at ICOMP during 1989 are described

Cumulative reports and publications through December 31, 1990

This document contains a complete list of ICASE reports. Since ICASE reports are intended to be preprints of articles that will appear in journals or conference proceedings, the published reference is included when it is available

Spectral methods in fluid dynamics

Fundamental aspects of spectral methods are introduced. Recent developments in spectral methods are reviewed with an emphasis on collocation techniques. Their applications to both compressible and incompressible flows, to viscous as well as inviscid flows, and also to chemically reacting flows are surveyed. The key role that these methods play in the simulation of stability, transition, and turbulence is brought out. A perspective is provided on some of the obstacles that prohibit a wider use of these methods, and how these obstacles are being overcome