Source author record

Dustin D. Keck

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

ResearcherUnclaimed source record

math.NA Quantitative Methods

Catalog footprint

What is connected

2works

2topics

1close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2014arXiv

Generalized sensitivity functions for size-structured population models

Size-structured population models provide a popular means to mathematically describe phenomena such as bacterial aggregation, schooling fish, and planetesimal evolution. For parameter estimation, generalized sensitivity functions (GSFs) provide a tool that quantifies the impact of data from specific regions of the experimental domain. These functions help identify the most relevant data subdomains, which enhances the optimization of experimental design. To our knowledge, GSFs have not been used in the partial differential equation (PDE) realm, so we provide a novel PDE extension of the discrete and continuous ordinary differential equation (ODE) concepts of Thomaseth and Cobelli and Banks et al. respectively. We analyze the GSFs in the context of size-structured population models, and specifically analyze the Smoluchowski coagulation equation to determine the most relevant time and volume domains for three, distinct aggregation kernels. Finally, we provide evidence that parameter estimation for the Smoluchowski coagulation equation does not require post-gelation data.

preprint2013arXiv

Numerical simulation of solutions and moments of the smoluchowski coagulation equation

Researchers have employed variations of the Smoluchowski coagulation equation to model a wide variety of both organic and inorganic phenomena and with relatively few known analytical solutions, numerical solutions play an important role in studying this equation. In this article, we consider numerical approximations, focusing on how different discretization schemes impact the accuracy of approximate solution moments. Pursuing the eventual goal of comparing simulated solutions to experimental data, we must carefully choose the numerical method most appropriate to the type of data we attain. Within this context, we compare and contrast the accuracy and computational cost of a finite element approach and a finite volume-based scheme. Our study provides theoretical and numerical evidence that the finite element approach achieves much more accuracy when the system aggregates slowly, and it does so with much less computation cost. Conversely, the finite volume method is slightly more accurate approximating the zeroth moment when the system aggregates quickly and is much more accurate approximating the first moment in general. Lastly, our study also provides numerical evidence that the finite element method (conventionally considered first order) actually belongs to a class of discontinuous Galerkin methods that exhibit superconvergence, or second order in our case.