Search
Results
-
Numerical approximation of statistical solutions of incompressible flow
(2015)Research ReportWe present a finite difference-(Multi-level) Monte Carlo algorithm to efficiently compute statistical solutions of the two dimensional Navier-Stokes equations, with periodic bound- ary conditions and for arbitrarily high Reynolds number. We propose a reformulation of statistical solutions in the vorticity-stream function form. The vorticity-stream function for- mulation is discretized with a finite difference scheme. We obtain a convergence ...Report -
Efficient computation of all speed flows using an entropy stable shock-capturing space-time discontinuous Galerkin method
(2014)Research ReportWe present a shock-capturing space-time Discontinuous Galerkin method to approximate all speed flows modeled by systems of conservation laws with multiple time scales. The method provides a very general and computationally efficient framework for approximating such systems on account of its ability to incorporate large time steps. Numerical examples ranging from computing the incompressible limit (robustness with respect to Mach number) ...Report -
Multi-level Monte Carlo finite volume method for shallow water equations with uncertain parameters applied to landslides-generated tsunamis
(2014)SAM Research ReportTwo layer Savage-Hutter type shallow water PDEs model flows such as tsunamis generated by rockslides. On account of heterogeneities in the composition of the granular matter, these models contain uncertain parameters like the ratio of densities of layers, Coulomb and interlayer friction. These parameters are modeled statistically and quantifying the resulting solution uncertainty (UQ) is a crucial task in geophysics. We propose a novel ...Report -
ENO reconstruction and ENO interpolation are stable
(2011)SAM Research ReportWe prove stability estimates for the ENO reconstruction and ENO interpolation procedures. In particular, we show that the jump of the reconstructed ENO pointvalues at each cell interface has the same sign as the jump of the underlying cell averages across that interface. We also prove that the jump of the reconstructed values can be upper-bounded in terms of the jump of the underlying cell averages. Similar sign properties hold for the ...Report -
Entropy stable schemes for initial-boundary-value conservation laws
(2010)SAM Research ReportWe consider initial boundary value problems for systems of conservation laws and design entropy stable finite difference schemes to approximate them. The schemes are shown to be entropy stable for a large class of systems that are equipped with a symmetric splitting, derived from the entropy formulation. Numerical examples for the Euler equations of gas dynamics are presented to illustrate the robust performance of the proposed method.Report -
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 -
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