Recent work on material interface reconstruction
- Los Alamos National Lab., NM (United States)
For the last 15 years, many Eulerian codes have relied on a series of piecewise linear interface reconstruction algorithms developed by David Youngs. In a typical Youngs` method, the material interfaces were reconstructed based upon nearly cell values of volume fractions of each material. The interfaces were locally represented by linear segments in two dimensions and by pieces of planes in three dimensions. The first step in such reconstruction was to locally approximate an interface normal. In Youngs` 3D method, a local gradient of a cell-volume-fraction function was estimated and taken to be the local interface normal. A linear interface was moved perpendicular to the now known normal until the mass behind it matched the material volume fraction for the cell in question. But for distorted or nonorthogonal meshes, the gradient normal estimate didn`t accurately match that of linear material interfaces. Moreover, curved material interfaces were also poorly represented. The authors will present some recent work in the computation of more accurate interface normals, without necessarily increasing stencil size. Their estimate of the normal is made using an iterative process that, given mass fractions for nearby cells of known but arbitrary variable density, converges in 3 or 4 passes in practice (and quadratically--like Newton`s method--in principle). The method reproduces a linear interface in both orthogonal and nonorthogonal meshes. The local linear approximation is generally 2nd-order accurate, with a 1st-order accurate normal for curved interfaces in both two and three dimensional polyhedral meshes. Recent work demonstrating the interface reconstruction for curved surfaces will /be discussed.
- Research Organization:
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- OSTI ID:
- 332744
- Report Number(s):
- SAND-98-1591; CONF-9709141-PROC.; ON: DE99000778; TRN: IM9916%%46
- Resource Relation:
- Conference: 5. joint Russian-American computational mathematics conference, Albuquerque, NM (United States), 2-5 Sep 1997; Other Information: PBD: [1997]; Related Information: Is Part Of Proceedings of the 5. joint Russian-American computational mathematics conference; PB: 312 p.
- Country of Publication:
- United States
- Language:
- English
Similar Records
Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method
Direction-aware slope limiter for three-dimensional cubic grids with adaptive mesh refinement