This article was accepted into the corpus but its outbound wikilinks were never NER-processed — typical at the deepest BFS hop or when the run's entity cap was reached. No expansion funnel to show.
| HLLC solver | |
|---|---|
| Name | HLLC solver |
| Developer | Einfeldt et al. |
| First appeared | 1990s |
| Type | approximate Riemann solver |
| Applications | Computational fluid dynamics, aerodynamics, astrophysics |
HLLC solver
The HLLC solver is an approximate Riemann solver widely used in Computational fluid dynamics and aerodynamics for hyperbolic systems of conservation laws, notable for restoring the contact and shear waves omitted by its predecessor, the HLL method. It was developed to combine robustness and simplicity, and it has been adopted in codes for high-speed flows, magnetohydrodynamics, and reactive flows across institutions such as NASA, ETH Zurich, Imperial College London, Princeton University, and Stanford University. The solver plays a central role in methods used alongside schemes by Godunov, Roe, Lax-Friedrichs, and MUSCL, and informs modern codes like those from OpenFOAM, FLASH, PLUTO, and Athena.
The HLLC solver originates from modifications to ideas introduced by P. D. Lax and R. Courant in early finite-volume methods and by the approximate Riemann framework of Harten, Lax, and van Leer; its specific contact-restoring extension is attributed to researchers including Toro and Einfeldt. It occupies a place in numerical methods alongside solvers such as Roe solver, HLL solver, AUSM, and AUSM+-up, and is often compared with exact Riemann solvers used in codes developed at Los Alamos National Laboratory and Lawrence Livermore National Laboratory. The technique is used in research groups at MIT, Caltech, University of Cambridge, University of Tokyo, and Tsinghua University.
The HLLC solver is cast within the conservation-law framework formalized by Jean Leray and later applied in computational contexts by Lax and Godunov. For a hyperbolic system with left and right states informed by reconstruction techniques from van Leer and Monotone Upstream-centered Schemes for Conservation Laws (MUSCL), HLLC estimates wave speeds similarly to approaches used by Davis and Einfeldt. The formulation introduces three-wave structure inspired by discontinuous solutions studied by Riemann and extended in works by Rankine and Hugoniot; it enforces conservation across approximate intermediate states analogous to constraints used in the Rankine–Hugoniot condition literature. Flux evaluation employs signal velocities related to eigenvalue estimates used in analyses by Liu and Toro, ensuring coupling with limiters developed by Sweby and Van Leer.
Practical implementation follows finite-volume paradigms popularized at Los Alamos National Laboratory and in software like OpenFOAM and FLASH. Typical codebases interoperate with reconstruction modules by ENO and WENO families attributed to Harten and Shu, and time integration schemes such as Runge–Kutta variants studied by Butcher. Implementation choices often reflect best practices from workshops at ICCS and AIAA conferences, and adopt slope limiters inspired by Minmod and Superbee formulations discussed in texts by Toro and LeVeque. Boundary condition treatments borrow methods developed at CEA and CNRS research groups, while parallelization strategies reference paradigms from MPI and OpenMP users at Argonne National Laboratory.
Extensions adapt the basic HLLC concept to systems beyond the Euler equations, paralleling developments in ideal magnetohydrodynamics by researchers at Princeton University and NASA Goddard Space Flight Center. Variants include entropy-fix modifications similar to those proposed by Harten and Hyman, positivity-preserving adaptations investigated at Stanford University and University of Michigan, and multidimensional reconstructions used in codes from Duke University and University of California, Berkeley. Specializations address reactive flows studied at Sandia National Laboratories and Lawrence Berkeley National Laboratory, shallow-water adaptations linked to NOAA and USGS modeling, and relativistic extensions connected to groups at Max Planck Institute for Astrophysics and Kavli Institute.
The solver is applied in aerospace design workflows at Boeing and Airbus, atmospheric reentry simulations at ESA and JAXA, and astrophysical modeling at Harvard-Smithsonian Center for Astrophysics and Harvard University. It supports multiphysics simulations for turbomachinery by teams at GE Aviation and Siemens and is used in environmental modeling projects coordinated by EPA and European Environment Agency. In academia, HLLC-based schemes support studies at University of Oxford, Yale University, Columbia University, University of Toronto, and National University of Singapore.
Performance assessments often compare HLLC with Roe solver, HLL solver, AUSM+, and exact solvers employed in benchmark suites from CFD General Notation System and competitions organized by AIAA. Studies at Imperial College London and ETH Zurich show HLLC provides improved contact resolution over HLL and greater robustness than Roe in strong-shock regimes examined in wind-tunnel tests at NASA Langley Research Center and DNW. Comparisons in magnetohydrodynamic contexts reference work at Los Alamos National Laboratory and Princeton University, while studies in reactive flows cite research from Sandia National Laboratories and Stanford University.
Challenges include difficulties in complex geometries addressed in collaborations with CERN and European XFEL, ensuring positivity under extreme conditions studied at Lawrence Livermore National Laboratory, and tuning for low-Mach-number flows investigated at University of Cambridge and Caltech. Numerical pathologies analogous to carbuncle phenomena analyzed by Quirk and mitigation strategies proposed by Liou remain active research topics at University of Illinois Urbana-Champaign and University of Maryland. Scaling to exascale computing environments involves efforts led by Argonne National Laboratory and Oak Ridge National Laboratory.
Category:Numerical methods