BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Daniel M. Tartakovsky

Daniel M. Tartakovsky contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phcond-mat.mtrl-scicond-mat.softMachine Learningmath.APmath.DSmath.OCmath.STphysics.geo-phStatistics Theory

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2022arXiv

Machine Learning in Heterogeneous Porous Materials

The "Workshop on Machine learning in heterogeneous porous materials" brought together international scientific communities of applied mathematics, porous media, and material sciences with experts in the areas of heterogeneous materials, machine learning (ML) and applied mathematics to identify how ML can advance materials research. Within the scope of ML and materials research, the goal of the workshop was to discuss the state-of-the-art in each community, promote crosstalk and accelerate multi-disciplinary collaborative research, and identify challenges and opportunities. As the end result, four topic areas were identified: ML in predicting materials properties, and discovery and design of novel materials, ML in porous and fractured media and time-dependent phenomena, Multi-scale modeling in heterogeneous porous materials via ML, and Discovery of materials constitutive laws and new governing equations. This workshop was part of the AmeriMech Symposium series sponsored by the National Academies of Sciences, Engineering and Medicine and the U.S. National Committee on Theoretical and Applied Mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Model Reduction via Dynamic Mode Decomposition

This work proposes a new framework of model reduction for parametric complex systems. The framework employs a popular model reduction technique dynamic mode decomposition (DMD), which is capable of combining data-driven learning and physics ingredients based on the Koopman operator theory. In the offline step of the proposed framework, DMD constructs a low-rank linear surrogate model for the high dimensional quantities of interest (QoIs) derived from the (nonlinear) complex high fidelity models (HFMs) of unknown forms. Then in the online step, the resulting local reduced order bases (ROBs) and parametric reduced order models (PROMs) at the training parameter sample points are interpolated to construct a new PROM with the corresponding ROB for a new set of target/test parameter values. The interpolations need to be done on the appropriate manifolds within consistent sets of generalized coordinates. The proposed framework is illustrated by numerical examples for both linear and nonlinear problems. In particular, its advantages in computational costs and accuracy are demonstrated by the comparisons with projection-based proper orthogonal decomposition (POD)-PROM and Kriging.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Information-geometry of physics-informed statistical manifolds and its use in data assimilation

The data-aware method of distributions (DA-MD) is a low-dimension data assimilation procedure to forecast the behavior of dynamical systems described by differential equations. It combines sequential Bayesian update with the MD, such that the former utilizes available observations while the latter propagates the (joint) probability distribution of the uncertain system state(s). The core of DA-MD is the minimization of a distance between an observation and a prediction in distributional terms, with prior and posterior distributions constrained on a statistical manifold defined by the MD. We leverage the information-geometric properties of the statistical manifold to reduce predictive uncertainty via data assimilation. Specifically, we exploit the information geometric structures induced by two discrepancy metrics, the Kullback-Leibler divergence and the Wasserstein distance, which explicitly yield natural gradient descent. To further accelerate optimization, we build a deep neural network as a surrogate model for the MD that enables automatic differentiation. The manifold's geometry is quantified without sampling, yielding an accurate approximation of the gradient descent direction. Our numerical experiments demonstrate that accounting for the information-geometry of the manifold significantly reduces the computational cost of data assimilation by facilitating the calculation of gradients and by reducing the number of required iterations. Both storage needs and computational cost depend on the dimensionality of a statistical manifold, which is typically small by MD construction. When convergence is achieved, the Kullback-Leibler and $L_2$ Wasserstein metrics have similar performances, with the former being more sensitive to poor choices of the prior.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Thermal experiments for fractured rock characterization: theoretical analysis and inverse modeling

Field-scale properties of fractured rocks play crucial role in many subsurface applications, yet methodologies for identification of the statistical parameters of a discrete fracture network (DFN) are scarce. We present an inversion technique to infer two such parameters, fracture density and fractal dimension, from cross-borehole thermal experiments data. It is based on a particle-based heat-transfer model, whose evaluation is accelerated with a deep neural network (DNN) surrogate that is integrated into a grid search. The DNN is trained on a small number of heat-transfer model runs, and predicts the cumulative density function of the thermal field. The latter is used to compute fine posterior distributions of the (to-be-estimated) parameters. Our synthetic experiments reveal that fracture density is well constrained by data, while fractal dimension is harder to determine. Adding non-uniform prior information related to the DFN connectivity improves the inference of this parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Data-informed Emulators for Multi-Physics Simulations

Machine learning techniques are powerful tools for construction of emulators for complex systems. We explore different machine learning methods and conceptual methodologies, ranging from functional approximations to dynamical approximations, to build such emulators for coupled thermal, hydrological, mechanical and chemical processes that occur near an engineered barrier system in the nuclear waste repository. Two nonlinear approximators, random forests and neural networks, are deployed to capture the complexity of the physics-based model and to identify its most significant hydrological and geochemical parameters. Our emulators capture the temporal evolution of the Uranium distribution coefficient of the clay buffer, and identify its functional dependence on these key parameters. The emulators' accuracy is further enhanced by assimilating relevant simulated predictors and clustering strategy. The relative performance of random forests and neural networks shows the advantage of ensemble learning in random forests algorithm, especially for highly nonlinear problems with limited data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamic Mode Decomposition for Construction of Reduced-Order Models of Hyperbolic Problems with Shocks

Construction of reduced-order models (ROMs) for hyperbolic conservation laws is notoriously challenging mainly due to the translational property and nonlinearity of the governing equations. While the Lagrangian framework for ROM construction resolves the translational issue, it is valid only before a shock forms. Once that occurs, characteristic lines cross each other and projection from a high-fidelity model space onto a ROM space distorts a moving grid, resulting in numerical instabilities. We address this grid distortion issue by developing a physics-aware dynamic mode decomposition (DMD) method based on hodograph transformation. The latter provides a map between the original nonlinear system and its linear counterpart, which coincides with the Koopman operator. This strategy is consistent with the spirit of physics-aware DMDs in that it retains information about shock dynamics. Several numerical examples are presented to validate the proposed physics-aware DMD approach to constructing accurate ROMs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of Data-driven Ambiguity Sets for Hyperbolic Conservation Laws with Uncertain Inputs

Ambiguity sets of probability distributions are used to hedge against uncertainty about the true probabilities of random quantities of interest (QoIs). When available, these ambiguity sets are constructed from both data (collected at the initial time and along the boundaries of the physical domain) and concentration-of-measure results on the Wasserstein metric. To propagate the ambiguity sets into the future, we use a physics-dependent equation governing the evolution of cumulative distribution functions (CDF) obtained through the method of distributions. This study focuses on the latter step by investigating the spatio-temporal evolution of data-driven ambiguity sets and their associated guarantees when the random QoIs they describe obey hyperbolic partial-differential equations with random inputs. For general nonlinear hyperbolic equations with smooth solutions, the CDF equation is used to propagate the upper and lower envelopes of pointwise ambiguity bands. For linear dynamics, the CDF equation allows us to construct an evolution equation for tighter ambiguity balls. We demonstrate that, in both cases, the ambiguity sets are guaranteed to contain the true (unknown) distributions within a prescribed confidence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic Pore Collapse Models in Granular Materials

Stochastic models for pore collapse in granular materials are developed. First, a general fluctuating stress-strain relation for a plastic flow rule is derived. The fluctuations account for non-associativity in plastic deformations typically observed in heterogeneous materials. Second, an axisymmetric spherical shell compaction model is extended to account for fluctuations in the material microstructure due to granular interactions at the pore scale. This changes the stress-strain constitutive equation determining the dynamics of pore collapse. Results show that stochastic differential equations can account for multiscale interactions in a statistical sense.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Estimation of distributions via multilevel Monte Carlo with stratified sampling

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the cumulative distribution function of a quantity of interest in problems with random input parameters or initial conditions. Our approach combines a standard MLMC method with stratified sampling by replacing standard Monte Carlo at each level with stratified Monte Carlo with proportional allocation. We show that the resulting stratified MLMC algorithm is more efficient than its standard MLMC counterpart, due to the reduction in variance at each level provided by the stratification of the random parameter's domain. A smoothing approximation for the indicator function based on kernel density estimation yields a more efficient algorithm compared to the typically used polynomial smoothing. The difference in computational cost between the smoothing methods depends on the required error tolerance.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Lagrangian Dynamic Mode Decomposition for Construction of Reduced-Order Models of Advection-Dominated Phenomena

Proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are two complementary singular-value decomposition (SVD) techniques that are widely used to construct reduced-order models (ROMs) in a variety of fields of science and engineering. Despite their popularity, both DMD and POD struggle to formulate accurate ROMs for advection-dominated problems because of the nature of SVD-based methods. We investigate this shortcoming of conventional POD and DMD methods formulated within the Eulerian framework. Then we propose a Lagrangian-based DMD method to overcome this so-called translational issues. Our approach is consistent with the spirit of physics-aware DMD since it accounts for the evolution of characteristic lines. Several numerical tests are presented to demonstrate the accuracy and efficiency of the proposed Lagrangian DMD method.

0 citations0 reviewsNo code linkedOpen access