BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Harbir Antil

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

ResearcherUnclaimed source record
math.NAmath.OCNumerical Analysismath.APMachine LearningBiological Physicsphysics.flu-dynphysics.soc-phComputational Engineering, Finance, and Scienceeess.IVgr-qcmath.CANeural and Evolutionary Computingphysics.geo-phphysics.med-ph

Catalog footprint

What is connected

31works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 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})_ν$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A note on semilinear fractional elliptic equation: analysis and discretization

In this paper we study existence, regularity, and approximation of solution to a fractional semilinear elliptic equation of order $s \in (0,1)$. We identify minimal conditions on the nonlinear term and the source which leads to existence of weak solutions and uniform $L^\infty$-bound on the solutions. Next we realize the fractional Laplacian as a Dirichlet-to-Neumann map via the Caffarelli-Silvestre extension. We introduce a first-degree tensor product finite elements space to approximate the truncated problem. We derive a priori error estimates and conclude with an illustrative numerical example.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An a posteriori error analysis for an optimal control problem involving the fractional Laplacian

In a previous work, we introduced a discretization scheme for a constrained optimal control problem involving the fractional Laplacian. For such a control problem, we derived optimal a priori error estimates that demand the convexity of the domain and some compatibility conditions on the data. To relax such restrictions, in this paper, we introduce and analyze an efficient and, under certain assumptions, reliable a posteriori error estimator. We realize the fractional Laplacian as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on a semi--infinite cylinder in one more spatial dimension. This extra dimension further motivates the design of an posteriori error indicator. The latter is defined as the sum of three contributions, which come from the discretization of the state and adjoint equations and the control variable. The indicator for the state and adjoint equations relies on an anisotropic error estimator in Muckenhoupt weighted Sobolev spaces. The analysis is valid in any dimension. On the basis of the devised a posteriori error estimator, we design a simple adaptive strategy that exhibits optimal experimental rates of convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Galerkin v. least-squares Petrov--Galerkin projection in nonlinear model reduction

Least-squares Petrov--Galerkin (LSPG) model-reduction techniques such as the Gauss--Newton with Approximated Tensors (GNAT) method have shown promise, as they have generated stable, accurate solutions for large-scale turbulent, compressible flow problems where standard Galerkin techniques have failed. However, there has been limited comparative analysis of the two approaches. This is due in part to difficulties arising from the fact that Galerkin techniques perform optimal projection associated with residual minimization at the time-continuous level, while LSPG techniques do so at the time-discrete level. This work provides a detailed theoretical and computational comparison of the two techniques for two common classes of time integrators: linear multistep schemes and Runge--Kutta schemes. We present a number of new findings, including conditions under which the LSPG ROM has a time-continuous representation, conditions under which the two techniques are equivalent, and time-discrete error bounds for the two approaches. Perhaps most surprisingly, we demonstrate both theoretically and computationally that decreasing the time step does not necessarily decrease the error for the LSPG ROM; instead, the time step should be `matched' to the spectral content of the reduced basis. In numerical experiments carried out on a turbulent compressible-flow problem with over one million unknowns, we show that increasing the time step to an intermediate value decreases both the error and the simulation time of the LSPG reduced-order model by an order of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimization with respect to order in a fractional diffusion model: analysis, approximation and algorithmic aspects

We consider an identification problem, where the state $u$ is governed by a fractional elliptic equation and the unknown variable corresponds to the order $s \in (0,1)$ of the underlying operator. We study the existence of an optimal pair $(\bar s, \bar u)$ and provide sufficient conditions for its local uniqueness. We develop semi-discrete and fully discrete algorithms to approximate the solutions to our identification problem and provide a convergence analysis. We present numerical illustrations that confirm and extend our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimizing the Kelvin force in a moving target subdomain

In order to generate a desired Kelvin (magnetic) force in a target subdomain moving along a prescribed trajectory, we propose a minimization problem with a tracking type cost functional. We use the so-called dipole approximation to realize the magnetic field, where the location and the direction of the magnetic sources are assumed to be fixed. The magnetic field intensity acts as the control and exhibits limiting pointwise constraints. We address two specific problems: the first one corresponds to a fixed final time whereas the second one deals with an unknown force to minimize the final time. We prove existence of solutions and deduce local uniqueness provided that a second order sufficient condition is valid. We use the classical backward Euler scheme for time discretization. For both problems we prove the $H^1$-weak convergence of this semi-discrete numerical scheme. This result is motivated by $Γ$-convergence and does not require second order sufficient condition. If the latter holds then we prove $H^1$-strong local convergence. We report computational results to assess the performance of the numerical methods. As an application, we study the control of magnetic nanoparticles as those used in magnetic drug delivery, where the optimized Kelvin force is used to transport the drug to a desired location.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Some applications of weighted norm inequalities to the error analysis of PDE constrained optimization problems

The purpose of this work is to illustrate how the theory of Muckenhoupt weights, Muckenhoupt weighted Sobolev spaces and the corresponding weighted norm inequalities can be used in the analysis and discretization of PDE constrained optimization problems. We consider: a linear quadratic constrained optimization problem where the state solves a nonuniformly elliptic equation; a problem where the cost involves pointwise observations of the state and one where the state has singular sources, e.g. point masses. For all three examples we propose and analyze numerical schemes and provide error estimates in two and three dimensions. While some of these problems might have been considered before in the literature, our approach allows for a simpler, Hilbert space-based, analysis and discretization and further generalizations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A FEM for an optimal control problem of fractional powers of elliptic operators

We study solution techniques for a linear-quadratic optimal control problem involving fractional powers of elliptic operators. These fractional operators can be realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on a semi-infinite cylinder in one more spatial dimension. Thus, we consider an equivalent formulation with a nonuniformly elliptic operator as state equation. The rapid decay of the solution to this problem suggests a truncation that is suitable for numerical approximation. We discretize the proposed truncated state equation using first degree tensor product finite elements on anisotropic meshes. For the control problem we analyze two approaches: one that is semi-discrete based on the so-called variational approach, where the control is not discretized, and the other one is fully discrete via the discretization of the control by piecewise constant functions. For both approaches, we derive a priori error estimates with respect to the degrees of freedom. Numerical experiments validate the derived error estimates and reveal a competitive performance of anisotropic over quasi-uniform refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A fractional space-time optimal control problem: analysis and discretization\

We study a linear-quadratic optimal control problem involving a parabolic equation with fractional diffusion and Caputo fractional time derivative of orders $s \in (0,1)$ and $γ\in (0,1]$, respectively. The spatial fractional diffusion is realized as the Dirichlet-to-Neumann map for a nonuniformly elliptic operator. Thus, we consider an equivalent formulation with a quasi-stationary elliptic problem with a dynamic boundary condition as state equation. The rapid decay of the solution to this problem suggests a truncation that is suitable for numerical approximation. We consider a fully-discrete scheme: piecewise constant functions for the control and, for the state, first-degree tensor product finite elements in space and a finite difference discretization in time. We show convergence of this scheme and, for $s \in (0,1)$ and $γ= 1$, we derive a priori error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optimal Control of a Free Boundary Problem with Surface Tension Effects: A Priori Error Analysis

We present a finite element method along with its analysis for the optimal control of a model free boundary problem with surface tension effects, formulated and studied in \cite{HAntil_RHNochetto_PSodre_2014a}. The state system couples the Laplace equation in the bulk with the Young-Laplace equation on the free boundary to account for surface tension. We first prove that the state and adjoint system have the requisite regularity for the error analysis (strong solutions). We discretize the state, adjoint and control variables via piecewise linear finite elements and show optimal $O(h)$ error estimates for all variables, including the control. This entails using the second order sufficient optimality conditions of \cite{HAntil_RHNochetto_PSodre_2014a}, and the first order necessary optimality conditions for both the continuous and discrete systems. We conclude with two numerical examples which examine the various error estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Stokes problem with Navier slip boundary condition: Minimal fractional Sobolev regularity of the domain

We prove well-posedness in reflexive Sobolev spaces of weak solutions to the stationary Stokes problem with Navier slip boundary condition over bounded domains $Ω$ of $\mathbb{R}^n$ of class $W^{2-1/s}_s$, $s>n$. Since such domains are of class $C^{1,1-n/s}$, our result improves upon the recent one by Amrouche-Seloula, who assume $Ω$ to be of class $C^{1,1}$. We deal with the slip boundary condition directly via a new localization technique, which features domain, space and operator decompositions. To flatten the boundary of $Ω$ locally, we construct a novel $W^2_s$ diffeomorphism for $Ω$ of class $W^{2-1/s}_s$. The fractional regularity gain, from $2-1/s$ to $2$, guarantees that the Piola transform is of class $W^1_s$. This allows us to transform $W^1_r$ vector fields without changing their regularity, provided $r\le s$, and preserve the unit normal which is Hölder. It is in this sense that the boundary regularity $W^{2-1/s}_s$ seems to be minimal.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal Control of a Free Boundary Problem: Analysis with Second Order Sufficient Conditions

We consider a PDE-constrained optimization problem governed by a free boundary problem. The state system is based on coupling the Laplace equation in the bulk with a Young-Laplace equation on the free boundary to account for surface tension, as proposed by P.\ Saavedra and L.\ R.\ Scott \cite{PSaavedra_RScott_1991}. This amounts to solving a second order system both in the bulk and on the interface. Our analysis hinges on a convex control constraint such that the state constraints are always satisfied. Using only first order regularity we show that the control to state operator is twice continuously Fréchet differentiable. We improve slightly the regularity of the state variables and exploit it to show existence of a control together with second order sufficient optimality conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal Control Of Surface Shape

Controlling the shapes of surfaces provides a novel way to direct self-assembly of colloidal particles on those surfaces and may be useful for material design. This motivates the investigation of an optimal control problem for surface shape in this paper. Specifically, we consider an objective (tracking) functional for surface shape with the prescribed mean curvature equation in graph form as a state constraint. The control variable is the prescribed curvature. We prove existence of an optimal control, and using improved regularity estimates, we show sufficient differentiability to make sense of the first order optimality conditions. This allows us to rigorously compute the gradient of the objective functional for both the continuous and discrete (finite element) formulations of the problem. Moreover, we provide error estimates for the state variable and adjoint state. Numerical results are shown to illustrate the minimizers and optimal controls on different domains.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two-step greedy algorithm for reduced order quadratures

We present an algorithm to generate application-specific, global reduced order quadratures (ROQ) for multiple fast evaluations of weighted inner products between parameterized functions. If a reduced basis (RB) or any other projection-based model reduction technique is applied, the dimensionality of integrands is reduced dramatically; however, the cost of approximating the integrands by projection still scales as the size of the original problem. In contrast, using discrete empirical interpolation (DEIM) points as ROQ nodes leads to a computational cost which depends linearly on the dimension of the reduced space. Generation of a reduced basis via a greedy procedure requires a training set, which for products of functions can be very large. Since this direct approach can be impractical in many applications, we propose instead a two-step greedy targeted towards approximation of such products. We present numerical experiments demonstrating the accuracy and the efficiency of the two-step approach. The presented ROQ are expected to display very fast convergence whenever there is regularity with respect to parameter variation. We find that for the particular application here considered, one driven by gravitational wave physics, the two-step approach speeds up the offline computations to build the ROQ by more than two orders of magnitude. Furthermore, the resulting ROQ rule is found to converge exponentially with the number of nodes, and a factor of ~50 savings, without loss of accuracy, is observed in evaluations of inner products when ROQ are used as a downsampling strategy for equidistant samples using the trapezoidal rule. While the primary focus of this paper is on quadrature rules for inner products of parameterized functions, our method can be easily adapted to integrations of single parameterized functions, and some examples of this type are considered.

0 citations0 reviewsNo code linkedOpen access