Search
Results
-
Existence of solutions to a model of two-phase flow in porous media
(2011)SAM Research ReportWe consider the flow of two-phases in a porous medium and propose a modi ed version of the fractional flow model for incompressible, two-phase flow based on a Helmholtz regularization of the Darcy phase velocities. We show the existence of global-in-time entropy solutions for this model with suitable assumptions on the boundary conditions. Numerical experiments demonstrating the approximation of the classical two-phase flow equations with ...Report -
Arbitrarily high order accurate entropy stable essentially non-oscillatory schemes for systems of conservation laws
(2011)SAM Research ReportWe design arbitrarily high-order accurate entropy stable schemes for systems of conservation laws. The schemes, termed TeCNO schemes, are based on two main ingredients: (i) high-order accurate entropy conservative uxes, and (ii) suitable numerical di usion operators involving ENO reconstructed cell-interface values of scaled entropy variables. Numerical experiments in one and two space dimensions are presented to illustrate the robust ...Report -
Accurate numerical schemes for approximating intitial-boundary value problems for systems of conservation law
(2011)SAM Research ReportSolutions of initial-boundary value problems for systems of conservation laws depend on the underlying viscous mechanism, namely different viscosity operators lead to different limit solutions. Standard numerical schemes for approximating conservation laws do not take into account this fact and converge to solutions that are not necessarily physically relevant. We design numerical schemes that incorporate explicit information about the ...Report -
Graph-Coupled Oscillator Networks
(2022)SAM Research ReportWe propose Graph-Coupled Oscillator Networks (GraphCON), a novel framework for deep learning on graphs. It is based on discretizations of a second-order system of ordinary differential equations (ODEs), which model a network of nonlinear forced and damped oscillators, coupled via the adjacency structure of the underlying graph. The flexibility of our framework permits any basic GNN layer (e.g. convolutional or attentional) as the coupling ...Report -
Scaling pattern with size by a morphogen directed cell division rule
(2011)CMA Collected Preprint seriesReport -
A Multi-level procedure for enhancing accuracy of machine learning algorithms
(2019)SAM Research ReportWe propose a multi-level method to increase the accuracy of machine learning algorithms for approximating observables in scientific computing, particularly those that arise in systems modeled by differential equations. The algorithm relies on judiciously combining a large number of computationally cheap training data on coarse resolutions with a few expensive training samples on fine grid resolutions. Theoretical arguments for lowering ...Report -
On the approximation of functions by tanh neural networks
(2021)SAM Research ReportWe derive bounds on the error, in high-order Sobolev norms, incurred in the approximation of Sobolev-regular as well as analytic functions by neural networks with the hyperbolic tangent activation function. These bounds provide explicit estimates on the approximation error with respect to the size of the neural networks. We show that tanh neural networks with only two hidden layers suffice to approximate functions at comparable or better ...Report -
Numerical methods for conservation laws with discontinuous coefficients
(2016)Research reports / Seminar for Applied MathematicsReport -
An operator preconditioning perspective on training in physics-informed machine learning
(2023)SAM Research ReportIn this paper, we investigate the behavior of gradient descent algorithms in physics-informed machine learning methods like PINNs, which minimize resid uals connected to partial differential equations (PDEs). Our key result is that the difficulty in training these models is closely related to the conditioning of a specific differential operator. This operator, in turn, is associated to the Hermi tian square of the differential operator ...Report -
Error analysis for physics informed neural networks (PINNs) approximating Kolmogorov PDEs
(2021)SAM Research ReportPhysics informed neural networks approximate solutions of PDEs by minimizing pointwise residuals. We derive rigorous bounds on the error, incurred by PINNs in approximating the solutions of a large class of linear parabolic PDEs, namely Kolmogorov equations that include the heat equation and Black-Scholes equation of option pricing, as examples. We construct neural networks, whose PINN residual (generalization error) can be made as small ...Report