Source author record

Praveen Chandrashekar

Praveen Chandrashekar appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Geometry Aware Operator Transformer as an Efficient and Accurate Neural Surrogate for PDEs on Arbitrary Domains

The very challenging task of learning solution operators of PDEs on arbitrary domains accurately and efficiently is of vital importance to engineering and industrial simulations. Despite the existence of many operator learning algorithms to approximate such PDEs, we find that accurate models are not necessarily computationally efficient and vice versa. We address this issue by proposing a geometry aware operator transformer (GAOT) for learning PDEs on arbitrary domains. GAOT combines novel multiscale attentional graph neural operator encoders and decoders, together with geometry embeddings and (vision) transformer processors to accurately map information about the domain and the inputs into a robust approximation of the PDE solution. Multiple innovations in the implementation of GAOT also ensure computational efficiency and scalability. We demonstrate this significant gain in both accuracy and efficiency of GAOT over several baselines on a large number of learning tasks from a diverse set of PDEs, including achieving state of the art performance on three large scale three-dimensional industrial CFD datasets.

preprint2022arXiv

Exact solution for Riemann problems of the shear shallow water model

The shear shallow water model is a higher order model for shallow flows which includes some shear effects that are neglected in the classical shallow models. The model is a non-conservative hyperbolic system which can admit shocks, rarefactions, shear and contact waves. The notion of weak solution is based on a path but the choice of the correct path is not known for this problem. In this paper, we construct exact solution for the Riemann problem assuming a linear path in the space of conserved variables, which is also used in approximate Riemann solvers. We compare the exact solutions with those obtained from a path conservative finite volume scheme on some representative test cases.

preprint2022arXiv

Lax-Wendroff flux reconstruction method for hyperbolic conservation laws

The Lax-Wendroff method is a single step method for evolving time dependent solutions governed by partial differential equations, in contrast to Runge- Kutta methods that need multiple stages per time step. We develop a flux reconstruction version of the method in combination with a Jacobian-free Lax- Wendroff procedure that is applicable to general hyperbolic conservation laws. The method is of collocation type, is quadrature free and can be cast in terms of matrix and vector operations. Special attention is paid to the construction of numerical flux, including for non-linear problems, resulting in higher CFL numbers than existing methods, which is shown through Fourier analysis and yielding uniform performance at all orders. Numerical results up to fifth order of accuracy for linear and non-linear problems are given to demonstrate the performance and accuracy of the method.

preprint2020arXiv

A path conservative finite volume method for a shear shallow water model

The shear shallow water model provides an approximation for shallow water flows by including the effect of vertical shear in the model. This model can be derived from the depth averaging process by including the second order velocity fluctuations which are neglected in the classical shallow water approximation. The resulting model has a non-conservative structure which resembles the 10-moment equations from gas dynamics. This structure facilitates the development of path conservative schemes and we construct HLL, 3-wave and 5-wave HLLC-type solvers. An explicit and semi-implicit MUSCL-Hancock type second order scheme is proposed for the time integration. Several test cases including roll waves show the performance of the proposed modeling and numerical strategy.

preprint2020arXiv

Constraint preserving discontinuous Galerkin method for ideal compressible MHD on 2-D Cartesian grids

We propose a constraint preserving discontinuous Galerkin method for ideal compressible MHD in two dimensions and using Cartesian grids, which automatically maintains the global divergence-free property. The approximation of the magnetic field is achieved using Raviart-Thomas polynomials and the DG scheme is based on evolving certain moments of these polynomials which automatically guarantees divergence-free property. We also develop HLL-type multi-dimensional Riemann solvers to estimate the electric field at vertices which are consistent with the 1-D Riemann solvers. When limiters are used, the divergence-free property may be lost and it is recovered by a divergence-free reconstruction step. We show the performance of the method on a range of test cases up to fourth order of accuracy.

preprint2020arXiv

High order well-balanced finite volume methods for multi-dimensional systems of hyperbolic balance laws

We introduce a general framework for the construction of well-balanced finite volume methods for hyperbolic balance laws. We use the phrase well-balancing in a broader sense, since our proposed method can be applied to exactly follow any solution of any system of hyperbolic balance laws in multiple spatial dimensions and not only time independent solutions. The solution has to be known a priori, either as an analytical expression or as discrete data. The proposed framework modifies the standard finite volume approach such that the well-balancing property is obtained and in case the method is high order accurate, this is maintained under our modification. We present numerical tests for the compressible Euler equations with and without gravity source term and with different equations of state, and for the equations of compressible ideal magnetohydrodynamics.

preprint2016arXiv

Simulating Turbulence Using the Astrophysical Discontinuous Galerkin Code TENET

In astrophysics, the two main methods traditionally in use for solving the Euler equations of ideal fluid dynamics are smoothed particle hydrodynamics and finite volume discretization on a stationary mesh. However, the goal to efficiently make use of future exascale machines with their ever higher degree of parallel concurrency motivates the search for more efficient and more accurate techniques for computing hydrodynamics. Discontinuous Galerkin (DG) methods represent a promising class of methods in this regard, as they can be straightforwardly extended to arbitrarily high order while requiring only small stencils. Especially for applications involving comparatively smooth problems, higher-order approaches promise significant gains in computational speed for reaching a desired target accuracy. Here, we introduce our new astrophysical DG code TENET designed for applications in cosmology, and discuss our first results for 3D simulations of subsonic turbulence. We show that our new DG implementation provides accurate results for subsonic turbulence, at considerably reduced computational cost compared with traditional finite volume methods. In particular, we find that DG needs about 1.8 times fewer degrees of freedom to achieve the same accuracy and at the same time is more than 1.5 times faster, confirming its substantial promise for astrophysical applications.

preprint2016arXiv

Well-balanced nodal discontinuous Galerkin method for Euler equations with gravity

We present a well-balanced nodal discontinuous Galerkin (DG) scheme for compressible Euler equations with gravity. The DG scheme makes use of discontinuous Lagrange basis functions supported at Gauss-Lobatto-Legendre (GLL) nodes together with GLL quadrature using the same nodes. The well-balanced property is achieved by a specific form of source term discretization that depends on the nature of the hydrostatic solution, together with the GLL nodes for quadrature of the source term. The scheme is able to preserve isothermal and polytropic stationary solutions upto machine precision on any mesh composed of quadrilateral cells and for any gravitational potential. It is applied on several examples to demonstrate its well-balanced property and the improved resolution of small perturbations around the stationary solution.

preprint2015arXiv

Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement

Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second order finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.

preprint2012arXiv

Discontinuous Galerkin method for Navier-Stokes equations using kinetic flux vector splitting

Kinetic schemes for compressible flow of gases are constructed by exploiting the connection between Boltzmann equation and the Navier-Stokes equations. This connection allows us to construct a flux splitting for the Navier-Stokes equations based on the direction of molecular motion from which a numerical flux can be obtained. The naive use of such a numerical flux function in a discontinuous Galerkin (DG) discretization leads to an unstable scheme in the viscous dominated case. Stable schemes are constructed by adding additional terms either in a symmetric or non-symmetric manner which are motivated by the DG schemes for elliptic equations. The novelty of the present scheme is the use of kinetic fluxes to construct the stabilization terms. In the symmetric case, interior penalty terms have to be added for stability and the resulting schemes give optimal convergence rates in numerical experiments. The non-symmetric schemes lead to a cell energy/entropy inequality but exhibit sub-optimal convergence rates. These properties are studied by applying the schemes to a scalar convection-diffusion equation and the 1-D compressible Navier-Stokes equations. In the case of Navier-Stokes equations, entropy variables are used to construct stable schemes.

preprint2012arXiv

Kinetic energy preserving and entropy stable finite volume schemes for compressible Euler and Navier-Stokes equations

Centered numerical fluxes can be constructed for compressible Euler equations which preserve kinetic energy in the semi-discrete finite volume scheme. The essential feature is that the momentum flux should be of the form $f^m_\jph = \tp_\jph + \avg{u}_\jph f^ρ_\jph$ where $\avg{u}_\jph = (u_j + u_{j+1})/2$ and $\tp_\jph, f^ρ_\jph$ are {\em any} consistent approximations to the pressure and the mass flux. This scheme thus leaves most terms in the numerical flux unspecified and various authors have used simple averaging. Here we enforce approximate or exact entropy consistency which leads to a unique choice of all the terms in the numerical fluxes. As a consequence novel entropy conservative flux that also preserves kinetic energy for the semi-discrete finite volume scheme has been proposed. These fluxes are centered and some dissipation has to be added if shocks are present or if the mesh is coarse. We construct scalar artificial dissipation terms which are kinetic energy stable and satisfy approximate/exact entropy condition. Secondly, we use entropy- variable based matrix dissipation flux which leads to kinetic energy and entropy stable schemes. These schemes are shown to be free of entropy violating solutions unlike the original Roe scheme. For hypersonic flows a blended scheme is proposed which gives carbuncle free solutions for blunt body flows. Numerical results for Euler and Navier-Stokes equations are presented to demonstrate the performance of the different schemes.

preprint2011arXiv

Vertex-centroid finite volume scheme on tetrahedral grids for conservation laws

Vertex-centroid schemes are cell-centered finite volume schemes for conservation laws which make use of vertex values to construct high resolution schemes. The vertex values must be obtained through a consistent averaging (interpolation) procedure. A modified interpolation scheme is proposed which is better than existing schemes in giving positive weights in the interpolation formula. A simplified reconstruction scheme is also proposed which is also more accurate and efficient. For scalar conservation laws, we develop limited versions of the schemes which are stable in maximum norm by constructing suitable limiters. The schemes are applied to compressible flows governed by the Euler equations of inviscid gas dynamics.