Researcher profile

Harbir Antil

Harbir Antil contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2025arXiv

Well-Posedness of Pseudo-Parabolic Gradient Systems with State-Dependent Dynamics

This paper develops a general mathematical framework for pseudo-parabolic gradient systems with state-dependent dynamics. The state dependence is induced by variable coefficient fields in the governing energy functional. Such coefficients arise naturally in scientific and technological models, including state-dependent mobilities in KWC-type grain boundary motion and variable orientation-adaptation operators in anisotropic image denoising. We establish two main results: the existence of energy-dissipating solutions, and the uniqueness and continuous dependence on initial data. The proposed framework yields a general well-posedness theory for a broad class of nonlinear evolutionary systems driven by state-dependent operators. As illustrative applications, we present an anisotropic image-denoising model and a new pseudo-parabolic KWC-type model for anisotropic grain boundary motion, and prove that both fit naturally within the abstract structure of $(\mathrm{S})_ν$.

preprint2022arXiv

An Optimal Time Variable Learning Framework for Deep Neural Networks

Feature propagation in Deep Neural Networks (DNNs) can be associated to nonlinear discrete dynamical systems. The novelty, in this paper, lies in letting the discretization parameter (time step-size) vary from layer to layer, which needs to be learned, in an optimization framework. The proposed framework can be applied to any of the existing networks such as ResNet, DenseNet or Fractional-DNN. This framework is shown to help overcome the vanishing and exploding gradient issues. Stability of some of the existing continuous DNNs such as Fractional-DNN is also studied. The proposed approach is applied to an ill-posed 3D-Maxwell's equation.

preprint2022arXiv

Efficient algorithms for Bayesian Inverse Problems with Whittle--Matérn Priors

This paper tackles efficient methods for Bayesian inverse problems with priors based on Whittle--Matérn Gaussian random fields. The Whittle--Matérn prior is characterized by a mean function and a covariance operator that is taken as a negative power of an elliptic differential operator. This approach is flexible in that it can incorporate a wide range of prior information including non-stationary effects, but it is currently computationally advantageous only for integer values of the exponent. In this paper, we derive an efficient method for handling all admissible noninteger values of the exponent. The method first discretizes the covariance operator using finite elements and quadrature, and uses preconditioned Krylov subspace solvers for shifted linear systems to efficiently apply the resulting covariance matrix to a vector. This approach can be used for generating samples from the distribution in two different ways: by solving a stochastic partial differential equation, and by using a truncated Karhunen-Loève expansion. We show how to incorporate this prior representation into the infinite-dimensional Bayesian formulation, and show how to efficiently compute the maximum a posteriori estimate, and approximate the posterior variance. Although the focus of this paper is on Bayesian inverse problems, the techniques developed here are applicable to solving systems with fractional Laplacians and Gaussian random fields. Numerical experiments demonstrate the performance and scalability of the solvers and their applicability to model and real-data inverse problems in tomography and a time-dependent heat equation.

preprint2022arXiv

NINNs: Nudging Induced Neural Networks

New algorithms called nudging induced neural networks (NINNs), to control and improve the accuracy of deep neural networks (DNNs), are introduced. The NINNs framework can be applied to almost all pre-existing DNNs, with forward propagation, with costs comparable to existing DNNs. NINNs work by adding a feedback control term to the forward propagation of the network. The feedback term nudges the neural network towards a desired quantity of interest. NINNs offer multiple advantages, for instance, they lead to higher accuracy when compared with existing data assimilation algorithms such as nudging. Rigorous convergence analysis is established for NINNs. The algorithmic and theoretical findings are illustrated on examples from data assimilation and chemically reacting flows.

preprint2022arXiv

Nonlocal Bounded Variations with Applications

Motivated by problems where jumps across lower dimensional subsets and sharp transitions across interfaces are of interest, this paper studies the properties of fractional bounded variation ($BV$)-type spaces. Two different natural fractional analogs of classical $BV$ are considered: $BV^α$, a space induced from the Riesz-fractional gradient that has been recently studied by Comi-Stefani; and $bv^α$, induced by the Gagliardo-type fractional gradient often used in Dirichlet forms and Peridynamics - this one is naturally related to the Caffarelli-Roquejoffre-Savin fractional perimeter. Our main theoretical result is that the latter $bv^α$ actually corresponds to the Gagliardo-Slobodeckij space $W^{α,1}$. As an application, using the properties of these spaces, novel image denoising models are introduced and their corresponding Fenchel pre-dual formulations are derived. The latter requires density of smooth functions with compact support. We establish this density property for convex domains.

preprint2022arXiv

Optimal Control of a Quasi-Variational Sweeping Process

The paper addresses the study of a class of evolutionary quasi-variational inequalities of the parabolic type arising in the formation and growth models of granular and cohensionless materials. Such models and their mathematical descriptions are highly challenging and require powerful tools of their analysis and implementation. We formulate a space-time continuous optimal control problem for a basic model of this type, develop several regularization and approximation procedures, and establish the existence of optimal solutions for the time-continuous and space-discrete problem. Viewing a version of this problem as a controlled quasi-variational sweeping process leads us to deriving necessary optimality conditions for the fully discrete problem by using the advanced machinery of variational analysis and generalized differentiation.

preprint2021arXiv

Novel Deep neural networks for solving Bayesian statistical inverse

We consider the simulation of Bayesian statistical inverse problems governed by large-scale linear and nonlinear partial differential equations (PDEs). Markov chain Monte Carlo (MCMC) algorithms are standard techniques to solve such problems. However, MCMC techniques are computationally challenging as they require several thousands of forward PDE solves. The goal of this paper is to introduce a fractional deep neural network based approach for the forward solves within an MCMC routine. Moreover, we discuss some approximation error estimates and illustrate the efficiency of our approach via several numerical examples.

preprint2020arXiv

Detailed Simulation of Viral Propagation In The Built Environment

A summary is given of the mechanical characteristics of virus contaminants and the transmission via droplets and aerosols. The ordinary and partial differential equations describing the physics of these processes with high fidelity are presented, as well as appropriate numerical schemes to solve them. Several examples taken from recent evaluations of the built environment are shown, as well as the optimal placement of sensors.

preprint2020arXiv

Fractional Deep Neural Network via Constrained Optimization

This paper introduces a novel algorithmic framework for a deep neural network (DNN), which in a mathematically rigorous manner, allows us to incorporate history (or memory) into the network -- it ensures all layers are connected to one another. This DNN, called Fractional-DNN, can be viewed as a time-discretization of a fractional in time nonlinear ordinary differential equation (ODE). The learning problem then is a minimization problem subject to that fractional ODE as constraints. We emphasize that an analogy between the existing DNN and ODEs, with standard time derivative, is well-known by now. The focus of our work is the Fractional-DNN. Using the Lagrangian approach, we provide a derivation of the backward propagation and the design equations. We test our network on several datasets for classification problems. Fractional-DNN offers various advantages over the existing DNN. The key benefits are a significant improvement to the vanishing gradient issue due to the memory effect, and better handling of nonsmooth data due to the network's ability to approximate non-smooth functions.

preprint2020arXiv

Fractional Diffusion Maps

In this paper, we extend the diffusion maps algorithm on a family of heat kernels that are either local (having exponential decay) or nonlocal (having polynomial decay), arising in various applications. For example, these kernels have been used as a regularizer in various supervised learning tasks for denoising images. Importantly, these heat kernels give rise to operators that include (but are not restricted to) the generators of the classical Laplacian associated to Brownian processes as well as the fractional Laplacian associated with $β$-stable Lévy processes. For local kernels, while the method is a version of the diffusion maps algorithm, we show that the applications with non-Gaussian local heat kernels approximate temporally rescaled Laplace-Beltrami operators. For the non-local heat kernels, we modify the diffusion maps algorithm to estimate fractional Laplacian operators. Here, the graph distance is used to approximate the geodesic distance with appropriate error bounds. While this approximation becomes numerically expensive as the number of data points increases, it produces an accurate operator estimation that is robust to the choice of the kernel bandwidth parameter value. In contrast, the local kernels are numerically more efficient but more sensitive to the choice of kernel bandwidth parameter value. In an application to estimate non-smooth regression functions, we find that using the nonlocal kernel as a regularizer produces a more robust and accurate estimate than using local kernels. For manifolds with boundary, we find that the proposed fractional diffusion maps framework implemented with non-local kernels approximates the regional fractional Laplacian.

preprint2020arXiv

High Fidelity Modeling of Aerosol Pathogen Propagation in Built Environments with Moving Pedestrians

A high fidelity model for the propagation of pathogens via aerosols in the presence of moving pedestrians is proposed. The key idea is the tight coupling of computational fluid dynamics and computational crowd dynamics in order to capture the emission, transport and inhalation of pathogen loads in space and time. An example simulating pathogen propagation in a narrow corridor with moving pedestrians clearly shows the considerable effect that pedestrian motion has on airflow, and hence on pathogen propagation and potential infectivity.

preprint2020arXiv

Moreau-Yosida regularization for optimal control of fractional PDEs with state constraints: parabolic case

This paper considers optimal control of fractional parabolic PDEs with both state and control constraints. The key challenge is how to handle the state constraints. Similarly, to the elliptic case, in this paper, we establish several new mathematical tools in the parabolic setting that are of wider interest. For example, existence of solution to the fractional parabolic equation with measure data on the right-hand-side. We employ the Moreau-Yosida regularization to handle the state constraints. We establish convergence, with rate, of the regularized optimal control problem to the original one. Numerical experiments confirm what we have proven theoretically.

preprint2020arXiv

Optimal Control of Parameterized Maxwell's System: Reduced Basis, Convergence Analysis, and A Posteriori Error Estimates

We consider control constrained optimal control problems governed by parameterized stationary Maxwell's system with the Gauss's law. The parameters enter through dielectric, magnetic permeability, and charge density. Moreover, the parameter set is assumed to be compact. We discretize the electric field by a finite element method and use variational discretization concept to discretize the control. We create a reduced basis method for the optimal control problem and establish uniform convergence of the reduced order solutions to that of the original high dimensional problem provided that the snapshot parameter sample is dense in the parameter set, with an appropriate parameter separability rule. Finally, we establish the absolute a posteriori error estimator for the reduced order solutions and the corresponding cost functions in terms of the state and adjoint residuals.

preprint2020arXiv

Optimal Control Problems Governed by 1-D Kobayashi-Warren-Carter Type Systems

This paper is devoted to the study of a class of optimal control problems governed by 1-D Kobayashi-Warren-Carter type systems, which are based on a phase-field model of grain boundary motion, proposed by [Kobayashi et al, Physica D, 140, 141-150, 2000]. The class consists of an optimal control problem for a physically realistic state-system of Kobayashi-Warren-Carter type, and its regularized approximating problems. The results of this paper are stated in three Main Theorems 1-3. The first Main Theorem 1 is concerned with the solvability and continuous dependence for the state-systems. Meanwhile, the second Main Theorem 2 is concerned with the solvability of optimal control problems, and some semi-continuous association in the class of our optimal control problems. Finally, in the third Main Theorem 3, we derive the first order necessary optimality conditions for optimal controls of the regularized approximating problems. By taking the approximating limit, we also derive the optimality conditions for the optimal controls for the physically realistic problem.

preprint2019arXiv

Bilevel Optimization, Deep Learning and Fractional Laplacian Regularization with Applications in Tomography

In this work we consider a generalized bilevel optimization framework for solving inverse problems. We introduce fractional Laplacian as a regularizer to improve the reconstruction quality, and compare it with the total variation regularization. We emphasize that the key advantage of using fractional Laplacian as a regularizer is that it leads to a linear operator, as opposed to the total variation regularization which results in a nonlinear degenerate operator. Inspired by residual neural networks, to learn the optimal strength of regularization and the exponent of fractional Laplacian, we develop a dedicated bilevel optimization neural network with a variable depth for a general regularized inverse problem. We also draw some parallels between an activation function in a neural network and regularization. We illustrate how to incorporate various regularizer choices into our proposed network. As an example, we consider tomographic reconstruction as a model problem and show an improvement in reconstruction quality, especially for limited data, via fractional Laplacian regularization. We successfully learn the regularization strength and the fractional exponent via our proposed bilevel optimization neural network. We observe that the fractional Laplacian regularization outperforms total variation regularization. This is specially encouraging, and important, in the case of limited and noisy data.

preprint2019arXiv

Fractional Operators Applied to Geophysical Electromagnetics

A growing body of applied mathematics literature in recent years has focussed on the application of fractional calculus to problems of anomalous transport. In these analyses, the anomalous transport (of charge, tracers, fluid, etc.) is presumed attributable to long-range correlations of material properties within an inherently complex, and in some cases self-similar, conducting medium. Rather than considering an exquisitely discretized (and computationally intractable) representation of the medium, the complex and spatially correlated heterogeneity is represented through reformulation of the PDE governing the relevant transport physics such that its coefficients are, instead, smooth but paired with fractional-order space derivatives. Here we apply these concepts to the scalar Helmholtz equation and its use in electromagnetic interrogation of Earth's interior through the magnetotelluric method. We outline a practical algorithm for solving the Helmholtz equation using spectral methods coupled with finite element discretization. Execution of this algorithm for the magnetotelluric problem reveals several interesting features observable in field data: long--range correlation of the predicted electromagnetic fields; a power-law relationship between the squared impedance amplitude and squared wavenumber whose slope is a function of the fractional exponent within the governing Helmholtz equation; and, a non-constant apparent resistivity spectrum whose variability arises solely from the fractional exponent. In geologic settings characterized by self--similarity (e.g. fracture systems; thick and richly-textured sedimentary sequences, etc.) we posit that diagnostics are useful for geologic characterization of features far below the typical resolution limit of electromagnetic methods in geophysics.