Source author record

Sashikumaar Ganesan

Sashikumaar Ganesan 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

10works
8topics
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

10 published item(s)

preprint2022arXiv

An Operator-Splitting Finite Element Method for the Numerical Solution of Radiative Transfer Equation

An operator-splitting finite element scheme for the time-dependent, high-dimensional radiative transfer equation is presented in this paper. The streamline upwind Petrov-Galerkin finite element method and discontinuous Galerkin finite element method are used for the spatial-angular discretization of the radiative transfer equation, whereas the implicit backward Euler scheme is used for temporal discretization. Error analysis of the proposed numerical scheme for the fully discrete radiative transfer equation is presented. The stability and convergence estimates for the fully discrete problem are derived. Moreover, an operator-splitting algorithm for numerical simulation of high-dimensional equations is also presented. The validation of the derived estimates and implementation is demonstrated with appropriate numerical experiments.

preprint2020arXiv

An overlapping local projection stabilization for Galerkin approximations of Darcy and Stokes problems

A priori analysis for a generalized local projection stabilized conforming finite element approximation of Darcy flow and Stokes problems is presented in this paper. A first-order conforming P1 finite element space is used to approximate both the velocity and the pressure. It is shown that the stabilized discrete bilinear form satisfy the inf-sup condition with respect to a generalized local projection norm. Moreover, a priori error estimates are derived for both problems. Finally, the validation of the proposed stabilization scheme is demonstrated with appropriate numerical examples.

preprint2020arXiv

Generalized local projection stabilized finite element method for advection-reaction problems

A priori analysis for a generalized local projection stabilized finite element approximations for the solution of an advection-reaction equation is presented in this article. The stability and a priori error estimates are established for both the conforming and the nonconforming (Crouzeix-Raviart) approximations with respect to the local projection streamline derivative norm. Finally, the validation of the proposed stabilization scheme and verification of the the derived estimates are presented with appropriate numerical experiments.

preprint2020arXiv

SParSH-AMG: A library for hybrid CPU-GPU algebraic multigrid and preconditioned iterative methods

Hybrid CPU-GPU algorithms for Algebraic Multigrid methods (AMG) to efficiently utilize both CPU and GPU resources are presented. In particular, hybrid AMG framework focusing on minimal utilization of GPU memory with performance on par with GPU-only implementations is developed. The hybrid AMG framework can be tuned to operate at a significantly lower GPU-memory, consequently, enables to solve large algebraic systems. Combining the hybrid AMG framework as a preconditioner with Krylov Subspace solvers like Conjugate Gradient, BiCG methods provides a solver stack to solve a large class of problems. The performance of the proposed hybrid AMG framework is analysed for an array of matrices with different properties and size. Further, the performance of CPU-GPU algorithms are compared with the GPU-only implementations to illustrate the significantly lower memory requirements.

preprint2020arXiv

Spatio-temporal predictive modeling framework for infectious disease spread

A novel predictive modeling framework for the spread of infectious diseases using high dimensional partial differential equations is developed and implemented. A scalar function representing the infected population is defined on a high-dimensional space and its evolution over all directions is described by a population balance equation (PBE). New infections are introduced among the susceptible population from non-quarantined infected population based on their interaction, adherence to distancing norms, hygiene levels and any other societal interventions. Moreover, recovery, death, immunity and all aforementioned parameters are modeled on the high-dimensional space. To epitomize the capabilities and features of the above framework, prognostic estimates of Covid-19 spread using a six-dimensional (time, 2D space, infection severity, duration of infection, and population age) PBE is presented. Further, scenario analysis for different policy interventions and population behavior is presented, throwing more insights into the spatio-temporal spread of infections across disease age, intensity and age of population. These insights could be used for science-informed policy planning.

preprint2016arXiv

ALE-SUPG finite element method for convection-diffusion problems in time-dependent domains: Non-conservative form

Stability estimates for Streamline Upwind Petrov-Galerkin (SUPG) finite element method with different time integration schemes for the solution of a scalar transient convection-diffusion-reaction equation in a time-dependent domain are derived. The deformation of the domain is handled with the arbitrary Lagrangian-Eulerian (ALE) approach. In particular, the non-conservative form of the ALE scheme is considered. The implicit Euler, the Crank-Nicolson, and the backward-difference~(BDF-2) methods are used for temporal discretization. It is shown that the stability of the semi-discrete (continuous in time) ALE-SUPG equation is independent of the mesh velocity, whereas the stability of the fully discrete problem is only conditionally stable. The theoretical considerations are illustrated by a numerical example. Further, the dependence of numerical solution on the choice of stabilization parameter $δ_k$ is also presented.

preprint2016arXiv

An object oriented parallel finite element scheme for computations of PDEs: Design and implementation

Parallel finite element algorithms based on object-oriented concepts are presented. Moreover, the design and implementation of a data structure proposed are utilized in realizing a parallel geometric multigrid method. The ParFEMapper and the ParFECommunicator are the key components of the data structure in the proposed parallel scheme. These classes are constructed based on the type of finite elements (continuous or nonconforming or discontinuous) used. The proposed solver is compared with the open source direct solvers, MUMPS and PasTiX. Further, the performance of the parallel multigrid solver is analyzed up to 1080 processors. The solver shows a very good speedup up to 960 processors and the problem size has to be increased in order to maintain the good speedup when the number of processors are increased further. As a result, the parallel solver is able to handle large scale problems on massively parallel supercomputers. The proposed parallel finite element algorithms and multigrid solver are implemented in our in-house package ParMooN.

preprint2015arXiv

ALE-SUPG finite element method for convection-diffusion problems in time-dependent domains: Conservative form

A Streamline Upwind Petrov-Galerkin (SUPG) finite element method for transient convection-diffusion-reaction equation in time-dependent domains is proposed. In particular, a convection dominated transient scalar problem is considered. The time-dependent domain is handled by the arbitrary Lagrangian-Eulerian (ALE) approach, whereas the SUPG finite element method is used for the spatial discretization. Further, the first order backward Euler and the second order Crank-Nicolson methods are used for the temporal discretization. It is shown that the stability of the semidiscrete (continuous in time) conservative ALE-SUPG equation is independent of the mesh velocity, whereas the stability of the fully discrete problem is unconditionally stable for implicit Euler method and is only conditionally stable for Crank-Nicolson time discretization. Numerical results are presented to show the influence of the SUPG stabilization parameter in a time-dependent domain. Further, the proposed numerical scheme is applied to a boundary/layer problem in a time-dependent domain.

preprint2015arXiv

Computation of impinging droplet with soluble surfactants and dynamic contact angle

An arbitrary Lagrangian--Eulerian (ALE) finite element scheme for computations of soluble surfactant droplet impingement on a horizontal surface is presented. The numerical scheme solves the time-dependent Navier--Stokes equations for the fluid flow, scalar convection-diffusion equation for the surfactant transport in the bulk phase, and simultaneously, surface evolution equations for the surfactants on the free surface and on the liquid-solid interface. The effects of surfactants on the flow dynamics are included into the model through the surfactant-dependent surface tension and dynamic contact angle. In particular, the dynamic contact angle $(θ_d)$ of the droplet is defined in terms of the surfactant concentration at the contact line and the equilibrium contact angle $(θ_e^0)$ of the clean surface using the nonlinear equation of state for surface tension. Further, the surface forces are included in the model as the surface divergence of the surface stress tensor that allows to incorporate the Marangoni effects without calculating the surface gradient of the surfactant concentration on the free surface. In addition to a mesh convergence study and validation of the numerical results with experiments, the effects of adsorption and desorption surfactant coefficients on the flow dynamics in wetting, partially wetting and non-wetting droplets are studied in detail. It is observed that the effect of surfactants are more in wetting droplets than in the non-wetting droplets. Further, the presence of surfactants at the contact line reduces the equilibrium contact angle further when $ θ_e^0 $ is less than $90^\circ$, and increases it further when $ θ_e^0 $ is greater than $90^\circ$. The presence of surfactants has no effect on the contact angle, when $θ_e^0 =90^\circ$. The numerical study clearly demonstrates that the surfactant-dependent contact angle ...

preprint2015arXiv

On the Navier-slip boundary condition for computations of impinging droplets

A mesh-dependent relation for the slip number in the Navier-slip with friction boundary condition for computations of impinging droplets with sharp interface methods is proposed. The relation is obtained as a function of Reynolds number, Weber number and the mesh size. The proposed relation is validated for several test cases by comparing the numerically obtained wetting diameter with the experimental results. Further, the computationally obtained maximum wetting diameter using the proposed slip relation is verified with the theoretical predictions. The relative error between the computationally obtained maximum wetting diameter and the theoretical predictions is less than 10 % for impinging droplet on a hydrophilic surface, and the error increases in the case of hydrophobic surface.