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

Max Gunzburger

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

ResearcherUnclaimed source record
math.NANumerical Analysismath.APmath.STphysics.ao-phphysics.geo-phStatistics TheoryApplicationsComputationMachine Learningmath.FAmath.OCphysics.comp-ph

Catalog footprint

What is connected

23works
13topics
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

23 published item(s)

preprint2022arXiv

A Comparison of Neural Network Architectures for Data-Driven Reduced-Order Modeling

The popularity of deep convolutional autoencoders (CAEs) has engendered new and effective reduced-order models (ROMs) for the simulation of large-scale dynamical systems. Despite this, it is still unknown whether deep CAEs provide superior performance over established linear techniques or other network-based methods in all modeling scenarios. To elucidate this, the effect of autoencoder architecture on its associated ROM is studied through the comparison of deep CAEs against two alternatives: a simple fully connected autoencoder, and a novel graph convolutional autoencoder. Through benchmark experiments, it is shown that the superior autoencoder architecture for a given ROM application is highly dependent on the size of the latent space and the structure of the snapshot data, with the proposed architecture demonstrating benefits on data with irregular connectivity when the latent space is sufficiently large.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An asymptotically compatible coupling formulation for nonlocal interface problems with jumps

We introduce a mathematically rigorous formulation for a nonlocal interface problem with jumps and propose an asymptotically compatible finite element discretization for the weak form of the interface problem. After proving the well-posedness of the weak form, we demonstrate that solutions to the nonlocal interface problem converge to the corresponding local counterpart when the nonlocal data are appropriately prescribed. Several numerical tests in one and two dimensions show the applicability of our technique, its numerical convergence to exact nonlocal solutions, its convergence to the local limit when the horizons vanish, and its robustness with respect to the patch test.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-Order Multirate Explicit Time-Stepping Schemes for the Baroclinic-Barotropic Split Dynamics in Primitive Equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step \textcolor{black}{and the assistance of some intermediate approximations of the baroctropic mode obtained by substepping with the small-time step}; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode under an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Two benchmark tests from the "MPAS-Ocean" platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A cookbook for finite element methods for nonlocal problems, including quadrature rules and approximate Euclidean balls

The implementation of finite element methods (FEMs) for nonlocal models with a finite range of interaction poses challenges not faced in the partial differential equations (PDEs) setting. For example, one has to deal with weak forms involving double integrals which lead to discrete systems having higher assembly and solving costs due to possibly much lower sparsity compared to that of FEMs for PDEs. In addition, one may encounter non-smooth integrands. In many nonlocal models, nonlocal interactions are limited to bounded neighborhoods that are ubiquitously chosen to be Euclidean balls, resulting in the challenge of dealing with intersections of such balls with the finite elements. We focus on developing recipes for the efficient assembly of FEM stiffness matrices and on the choice of quadrature rules for the double integrals that contribute to the assembly efficiency and also posses sufficient accuracy. A major feature of our recipes is the use of approximate balls, e.g., several polygonal approximations of Euclidean balls, that, among other advantages, mitigate the challenge of dealing with ball-element intersections. We provide numerical illustrations of the relative accuracy and efficiency of the several approaches we develop.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A general framework for substructuring-based domain decomposition methods for models having nonlocal interactions

A rigorous mathematical framework is provided for a substructuring-based domain-decomposition approach for nonlocal problems that feature interactions between points separated by a finite distance. Here, by substructuring it is meant that a traditional geometric configuration for local partial differential equation problems is used in which a computational domain is subdivided into non-overlapping subdomains. In the nonlocal setting, this approach is substructuring-based in the sense that those subdomains interact with neighboring domains over interface regions having finite volume, in contrast to the local PDE setting in which interfaces are lower dimensional manifolds separating abutting subdomains. Key results include the equivalence between the global, single-domain nonlocal problem and its multi-domain reformulation, both at the continuous and discrete levels. These results provide the rigorous foundation necessary for the development of efficient solution strategies for nonlocal domain-decomposition methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An energy-based coupling approach to nonlocal interface problems

Nonlocal models provide accurate representations of physical phenomena ranging from fracture mechanics to complex subsurface flows, where traditional partial differential equations fail to capture effects caused by long-range forces at the microscale and mesoscale. However, the application of nonlocal models to problems involving interfaces such as multimaterial simulations and fluid-structure interaction, is hampered by the lack of a rigorous nonlocal interface theory needed to support numerical developments. In this paper, we use an energy-based approach to develop a mathematically rigorous nonlocal interface theory which provides a physically consistent extension of the classical perfect interface PDE formulation. Numerical examples validate the proposed framework and demonstrate the scope of our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Approximation of probability density functions for PDEs with random parameters using truncated series expansions

The probability density function (PDF) of a random variable associated with the solution of a partial differential equation (PDE) with random parameters is approximated using a truncated series expansion. The random PDE is solved using two stochastic finite element methods, Monte Carlo sampling and the stochastic Galerkin method with global polynomials. The random variable is a functional of the solution of the random PDE, such as the average over the physical domain. The truncated series are obtained considering a finite number of terms in the Gram-Charlier or Edgeworth series expansions. These expansions approximate the PDF of a random variable in terms of another PDF, and involve coefficients that are functions of the known cumulants of the random variable. To the best of our knowledge, their use in the framework of PDEs with random parameters has not yet been explored.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Exponential Time Differencing for the Tracer Equations Appearing in Primitive Equation Ocean Models

The tracer equations are part of the primitive equations used in ocean modeling and describe the transport of tracers, such as temperature, salinity or chemicals, in the ocean. Depending on the number of tracers considered, several equations may be added to and coupled to the dynamics system. In many relevant situations, the time-step requirements of explicit methods imposed by the transport and mixing in the vertical direction are more restrictive than those for the horizontal, and this may cause the need to use very small time steps if a fully explicit method is employed. To overcome this issue, we propose an exponential time differencing (ETD) solver where the vertical terms (transport and diffusion) are treated with a matrix exponential, whereas the horizontal terms are dealt with in an explicit way. We investigate numerically the computational speed-ups that can be obtained over other semi-implicit methods, and we analyze the advantages of the method in the case of multiple tracers.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Piecewise polynomial approximation of probability density functions with application to uncertainty quantification for stochastic PDEs

The probability density function (PDF) associated with a given set of samples is approximated by a piecewise-linear polynomial constructed with respect to a binning of the sample space. The kernel functions are a compactly supported basis for the space of such polynomials, i.e. finite element hat functions, that are centered at the bin nodes rather than at the samples, as is the case for the standard kernel density estimation approach. This feature naturally provides an approximation that is scalable with respect to the sample size. On the other hand, unlike other strategies that use a finite element approach, the proposed approximation does not require the solution of a linear system. In addition, a simple rule that relates the bin size to the sample size eliminates the need for bandwidth selection procedures. The proposed density estimator has unitary integral, does not require a constraint to enforce positivity, and is consistent. The proposed approach is validated through numerical examples in which samples are drawn from known PDFs. The approach is also used to determine approximations of (unknown) PDFs associated with outputs of interest that depend on the solution of a stochastic partial differential equation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Ensemble-Proper Orthogonal Decomposition Method for the Nonstationary Navier-Stokes Equations

The definition of partial differential equation (PDE) models usually involves a set of parameters whose values may vary over a wide range. The solution of even a single set of parameter values may be quite expensive. In many cases, e.g., optimization, control, uncertainty quantification, and other settings, solutions are needed for many sets of parameter values. We consider the case of the time-dependent Navier-Stokes equations for which a recently developed ensemble-based method allows for the efficient determination of the multiple solutions corresponding to many parameter sets. The method uses the average of the multiple solutions at any time step to define a linear set of equations that determines the solutions at the next time step. To significantly further reduce the costs of determining multiple solutions of the Navier-Stokes equations, we incorporate a proper orthogonal decomposition (POD) reduced-order model into the ensemble-based method. The stability and convergence results for the ensemble-based method are extended to the ensemble-POD approach. Numerical experiments are provided that illustrate the accuracy and efficiency of computations determined using the new approach.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analysis and approximation of a fractional Laplacian-based closure model for turbulent flows and its connection to Richardson pair dispersion

We study a turbulence closure model in which the fractional Laplacian $(-Δ)^α$ of the velocity field represents the turbulence diffusivity. We investigate the energy spectrum of the model by applying Pao's energy transfer theory. For the case $α=1/3$, the corresponding power law of the energy spectrum in the inertial range has a correction exponent on the regular Kolmogorov -5/3 scaling exponent. For this case, this model represents Richardson's particle pair-distance superdiffusion of a fully developed homogeneous turbulent flow as well as Lévy jumps that lead to the superdiffusion. For other values of $α$, the power law of the energy spectrum is consistent with the regular Kolmogorov -5/3 scaling exponent. We also propose and study a modular time-stepping algorithm in semi-discretized form. The algorithm is minimally intrusive to a given legacy code for solving Navier-Stokes equations by decoupling the local part and nonlocal part of the equations for the unknowns. We prove the algorithm is unconditionally stable and unconditionally, first-order convergent. We also derive error estimates for full discretizations of the model which, in addition to the time stepping algorithm, involves a finite element spatial discretization and a domain truncation approximation to the range of the fractional Laplacian.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Identification of the diffusion parameter in nonlocal steady diffusion problems

The problem of identifying the diffusion parameter appearing in a nonlocal steady diffusion equation is considered. The identification problem is formulated as an optimal control problem having a matching functional as the objective of the control and the parameter function as the control variable. The analysis makes use of a nonlocal vector calculus that allows one to define a variational formulation of the nonlocal problem. In a manner analogous to the local partial differential equations counterpart, we demonstrate, for certain kernel functions, the existence of at least one optimal solution in the space of admissible parameters. We introduce a Galerkin finite element discretization of the optimal control problem and derive a priori error estimates for the approximate state and control variables. Using one-dimensional numerical experiments, we illustrate the theoretical results and show that by using nonlocal models it is possible to estimate non-smooth and discontinuous diffusion parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Methods for a Class of Nonlocal Diffusion Problems with the Use of Backward SDEs

We propose a novel numerical approach for nonlocal diffusion equations [8] with integrable kernels, based on the relationship between the backward Kolmogorov equation and backward stochastic differential equations (BSDEs) driven by Lèvy processes with jumps. The nonlocal diffusion problem under consideration is converted to a BSDE,for which numerical schemes are developed and applied directly. As a stochastic approach, the proposed method does not require the solution of linear systems, which allows for embarrassingly parallel implementations and also enables adaptive approximation techniques to be incorporated in a straightforward fashion. Moreover, our method is more accurate than classic stochastic approaches due to the use of high-order temporal and spatial discretization schemes. In addition, our approach can handle a broad class of problems with general nonlinear forcing terms as long as they are globally Lipchitz continuous. Rigorous error analysis of the new method is provided as several numerical examples that illustrate the effectiveness and efficiency of the proposed approach.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Generalized Nonlocal Calculus with Application to the Peridynamics Model for Solid Mechanics

A nonlocal vector calculus was introduced in [2] that has proved useful for the analysis of the peridynamics model of nonlocal mechanics and nonlocal diffusion models. A generalization is developed that provides a more general setting for the nonlocal vector calculus that is independent of particular nonlocal models. It is shown that general nonlocal calculus operators are integral operators with specific integral kernels. General nonlocal calculus properties are developed, including nonlocal integration by parts formula and Green's identities. The nonlocal vector calculus introduced in [2] is shown to be recoverable from the general formulation as a special example. This special nonlocal vector calculus is used to reformulate the peridynamics equation of motion in terms of the nonlocal gradient operator and its adjoint. A new example of nonlocal vector calculus operators is introduced, which shows the potential use of the general formulation for general nonlocal models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Multilevel Stochastic Collocation Method for Partial Differential Equations with Random Input Data

Stochastic collocation methods for approximating the solution of partial differential equations with random input data (e.g., coefficients and forcing terms) suffer from the curse of dimensionality whereby increases in the stochastic dimension cause an explosion of the computational effort. We propose and analyze a multilevel version of the stochastic collocation method that, as is the case for multilevel Monte Carlo (MLMC) methods, uses hierarchies of spatial approximations to reduce the overall computational complexity. In addition, our proposed approach utilizes, for approximation in stochastic space, a sequence of multi-dimensional interpolants of increasing fidelity which can then be used for approximating statistics of the solution as well as for building high-order surrogates featuring faster convergence rates. A rigorous convergence and computational cost analysis of the new multilevel stochastic collocation method is provided, demonstrating its advantages compared to standard single-level stochastic collocation approximations as well as MLMC methods. Numerical results are provided that illustrate the theory and the effectiveness of the new multilevel method.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Optimal Point Sets for Total Degree Polynomial Interpolation in Moderate Dimensions

This paper is concerned with Lagrange interpolation by total degree polynomials in moderate dimensions. In particular, we are interested in characterising the optimal choice of points for the interpolation problem, where we define the optimal interpolation points as those which minimise the Lebesgue constant. We give a novel algorithm for numerically computing the location of the optimal points, which is independent of the shape of the domain and does not require computations with Vandermonde matrices. We perform a numerical study of the growth of the minimal Lebesgue constant with respect to the degree of the polynomials and the dimension, and report the lowest values known as yet of the Lebesgue constant in the unit cube and the unit ball in up to 10 dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Peridynamics and Material Interfaces

The convergence of a peridynamic model for solid mechanics inside heterogeneous media in the limit of vanishing nonlocality is analyzed. It is shown that the operator of linear peridynamics for an isotropic heterogeneous medium converges to the corresponding operator of linear elasticity when the material properties are sufficiently regular. On the other hand, when the material properties are discontinuous, i.e., when material interfaces are present, it is shown that the operator of linear peridynamics diverges, in the limit of vanishing nonlocality, at material interfaces. Nonlocal interface conditions, whose local limit implies the classical interface conditions of elasticity, are then developed and discussed. A peridynamics material interface model is introduced which generalizes the classical interface model of elasticity. The model consists of a new peridynamics operator along with nonlocal interface conditions. The new peridynamics interface model converges to the classical interface model of linear elasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Power-Law Noises over General Spatial Domains and on Non-Standard Meshes

Power-law noises abound in nature and have been observed extensively in both time series and spatially varying environmental parameters. Although, recent years have seen the extension of traditional stochastic partial differential equations to include systems driven by fractional Brownian motion, spatially distributed scale-invariance has received comparatively little attention, especially for parameters defined over non-standard spatial domains. This paper discusses the generalization of power-law noises to general spatial domains by outlining their theoretical underpinnings as well as addressing their numerical simulation on arbitrary meshes. Three computational algorithms are presented for efficiently generating their sample paths, accompanied by numerous numerical illustrations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The fractional Laplacian operator on bounded domains as a special case of the nonlocal diffusion operator

We analyze a nonlocal diffusion operator having as special cases the fractional Laplacian and fractional differential operators that arise in several applications. In our analysis, a nonlocal vector calculus is exploited to define a weak formulation of the nonlocal problem. We demonstrate that, when sufficient conditions on certain kernel functions hold, the solution of the nonlocal equation converges to the solution of the fractional Laplacian equation on bounded domains as the nonlocal interactions become infinite. We also introduce a continuous Galerkin finite element discretization of the nonlocal weak formulation and we derive a priori error estimates. Through several numerical examples we illustrate the theoretical results and we show that by solving the nonlocal problem it is possible to obtain accurate approximations of the solutions of fractional differential equations circumventing the problem of treating infinite-volume constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Efficient and Long-Time Accurate Second-Order Methods for Stokes-Darcy System

We propose and study two second-order in time implicit-explicit (IMEX) methods for the coupled Stokes-Darcy system that governs flows in karst aquifers. The first is a combination of a second-order backward differentiation formula and the second-order Gear's extrapolation approach. The second is a combination of the second-order Adams-Moulton and second-order Adams-Bashforth methods. Both algorithms only require the solution of two decoupled problems at each time step, one Stokes and the other Darcy. Hence, these schemes are very efficient and can be easily implemented using legacy codes. We establish the unconditional and uniform in time stability for both schemes. The uniform in time stability leads to uniform in time control of the error which is highly desirable for modeling physical processes, e.g., contaminant sequestration and release, that occur over very long time scales. Error estimates for fully-discretized schemes using finite element spatial discretizations are derived. Numerical examples are provided that illustrate the accuracy, efficiency, and long-time stability of the two schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A scale-invariant formulation of the anticipated potential vorticity method

The long-term success of climate models that operate on multi-resolution grids will depend on access to subgrid parameterizations that act appropriately across a wide range of spatial and temporal scales. As the first step in a series of efforts to obtain such scale-aware subgrid parametrizations, in this article we focus on the anticipated potential vorticity method (APVM) on quasi-uniform grids with varying resolutions. By a scale analysis technique and the phenomenological theories for two-dimensional turbulent flows, we derive a new formulation of the APVM, which depends on a single parameter that is, in principle, invariant with respect to the time step size, the grid resolution and the flow itself. Results of numerical experiments with this new formulation demonstrate that the optimal parameter of the APVM is invariant with respect to the time step size, insensitive to the flows, and is only weakly dependent on the grid resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Goal-oriented A Posteriori Error Estimation for Finite Volume Methods

A general framework for goal-oriented a posteriori error estimation for finite volume methods is presented. The framework does not rely on recasting finite volume methods as special cases of finite element methods, but instead directly determines error estimators from the discretized finite volume equations. Thus, the framework can be ap- plied to arbitrary finite volume methods. It also provides the proper functional settings to address well-posedness issues for the primal and adjoint problems. Numerical results are presented to illustrate the validity and effectiveness of the a posteriori error estimates and their applicability to adaptive mesh refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Partial and spectral-viscosity models for geophysical flows

Two models based on the hydrostatic primitive equa- tions are proposed. The first model is the primitive equations with partial viscosity only, and is oriented towards large-scale wave structures in the ocean and atmosphere. The second model is the viscous primitive equations with spectral eddy viscosity, and is oriented towards turbulent geophysical flows. For both models, the existence and uniqueness of global strong solutions is estab- lished. For the second model, the convergence of the solutions to the solutions of the classical primitive equations as eddy viscosity parameters tend to zero is also established.

0 citations0 reviewsNo code linkedOpen access