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

Panagiotis Stinis

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

ResearcherUnclaimed source record
math.NAmath.AP

Catalog footprint

What is connected

10works
2topics
3close 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

10 published item(s)

preprint2014arXiv

Mesh refinement for uncertainty quantification through model reduction

We present a novel way of deciding when and where to refine a mesh in probability space in order to facilitate the uncertainty quantification in the presence of discontinuities in random space. A discontinuity in random space makes the application of generalized polynomial chaos expansion techniques prohibitively expensive. The reason is that for discontinuous problems, the expansion converges very slowly. An alternative to using higher terms in the expansion is to divide the random space in smaller elements where a lower degree polynomial is adequate to describe the randomness. In general, the partition of the random space is a dynamic process since some areas of the random space, particularly around the discontinuity, need more refinement than others as time evolves. In the current work we propose a way to decide when and where to refine the random space mesh based on the use of a reduced model. The idea is that a good reduced model can monitor accurately, within a random space element, the cascade of activity to higher degree terms in the chaos expansion. In terms, this facilitates the efficient allocation of computational sources to the areas of random space where they are more needed. For the Kraichnan-Orszag system, the prototypical system to study discontinuities in random space, we present theoretical results which show why the proposed method is sound and numerical results which corroborate the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Model reduction and mesh refinement

In recent work we have presented a novel algorithm for mesh refinement which utilizes a reduced model. In particular, the reduced model is used to monitor the transfer of activity (e.g. mass, energy) from larger to smaller scales. When the activity transfer rate exceeds a prescribed tolerance, the algorithm refines the mesh accordingly. The algorithm was applied to the inviscid Burgers and focusing Schrödinger equations to detect singularities. We offer here a simple proof of why the use of a reduced model for mesh refinement is a good idea. We do so by showing that by controlling the transfer of activity rate one controls the error caused inevitably by a finite resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mori-Zwanzig reduced models for uncertainty quantification I: Parametric uncertainty

In many time-dependent problems of practical interest the parameters entering the equations describing the evolution of the various quantities exhibit uncertainty. One way to address the problem of how this uncertainty impacts the solution is to expand the solution using polynomial chaos expansions and obtain a system of differential equations for the evolution of the expansion coefficients. We present an application of the Mori-Zwanzig formalism to the problem of constructing reduced models of such systems of differential equations. In particular, we construct reduced models for a subset of the polynomial chaos expansion coefficients that are needed for a full description of the uncertainty caused by the uncertain parameters. We also present a Markovian reformulation of the Mori-Zwanzig reduced equation which replaces the solution of the orthogonal dynamics equation with an infinite hierarchy of ordinary differential equations. The viscous Burgers equation with uncertain viscosity parameter is used to illustrate the construction. For this example we provide a way to estimate the necessary parameters that appear in the reduced model without having to solve the full system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mori-Zwanzig reduced models for uncertainty quantification II: Initial condition uncertainty

In a recent preprint (arXiv:1211.4285v1) we addressed the problem of constructing reduced models for time-dependent systems described by differential equations which involve uncertain parameters. In the current work, we focus on the construction of reduced models for systems of differential equations when the initial condition is uncertain. While for both cases the reduced models are constructed through the Mori-Zwanzig formalism, the necessary estimation of the memory parameters is quite different. For the case of uncertain initial conditions we present an algorithm which allows to estimate on the fly the parameters appearing in the reduced model. The first part of the algorithm evolves the full system until the estimation of the parameters for the reduced model has converged. At the time instant that this happens, the algorithm switches to the evolution of only the reduced model with the estimated parameter values from the first part of the algorithm. The viscous Burgers equation with uncertain initial condition is used to illustrate the construction.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A drift homotopy Monte Carlo approach to particle filtering for multi-target tracking

We present a novel approach for improving particle filters for multi-target tracking. The suggested approach is based on drift homotopy for stochastic differential equations. Drift homotopy is used to design a Markov Chain Monte Carlo step which is appended to the particle filter and aims to bring the particle filter samples closer to the observations. Also, we present a simple Metropolis Monte Carlo algorithm for tackling the target-observation association problem. We have used the proposed approach on the problem of multi-target tracking for both linear and nonlinear observation models. The numerical results show that the suggested approach can improve significantly the performance of a particle filter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Conditional path sampling of stochastic differential equations by drift relaxation

We present an algorithm for the efficient sampling of conditional paths of stochastic differential equations (SDEs). While unconditional path sampling of SDEs is straightforward, albeit expensive for high dimensional systems of SDEs, conditional path sampling can be difficult even for low dimensional systems. This is because we need to produce sample paths of the SDE which respect both the dynamics of the SDE and the initial and endpoint conditions. The dynamics of a SDE are governed by the deterministic term (drift) and the stochastic term (noise). Instead of producing conditional paths directly from the original SDE, one can consider a sequence of SDEs with modified drifts. The modified drifts should be chosen so that it is easier to produce sample paths which satisfy the initial and endpoint conditions. Also, the sequence of modified drifts converges to the drift of the original SDE. We construct a simple Markov Chain Monte Carlo (MCMC) algorithm which samples, in sequence, conditional paths from the modified SDEs, by taking the last sampled path at each level of the sequence as an initial condition for the sampling at the next level in the sequence. The algorithm can be thought of as a stochastic analog of deterministic homotopy methods for solving nonlinear algebraic equations or as a SDE generalization of simulated annealing. The algorithm is particularly suited for filtering/smoothing applications. We show how it can be used to improve the performance of particle filters. Numerical results for filtering of a stochastic differential equation are included.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Path sampling for particle filters with application to multi-target tracking

In recent work (arXiv:1006.3100v1), we have presented a novel approach for improving particle filters for multi-target tracking. The suggested approach was based on drift homotopy for stochastic differential equations. Drift homotopy was used to design a Markov Chain Monte Carlo step which is appended to the particle filter and aims to bring the particle filter samples closer to the observations. In the current work, we present an alternative way to append a Markov Chain Monte Carlo step to a particle filter to bring the particle filter samples closer to the observations. Both current and previous approaches stem from the general formulation of the filtering problem. We have used the currently proposed approach on the problem of multi-target tracking for both linear and nonlinear observation models. The numerical results show that the suggested approach can improve significantly the performance of a particle filter.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Renormalized reduced models for singular PDEs

We present a novel way of constructing reduced models for systems of ordinary differential equations. The reduced models we construct depend on coefficients which measure the importance of the different terms appearing in the model and need to be estimated. The proposed approach allows the estimation of these coefficients on the fly by enforcing the equality of integral quantities of the solution as computed from the original system and the reduced model. In particular, the approach combines the concepts of renormalization and effective field theory developed in the context of high energy physics and the Mori-Zwanzig formalism of irreversible statistical mechanics. It allows to construct stable reduced models of higher order than was previously possible. The method is applied to the problem of computing reduced models for ordinary differential equation systems resulting from Fourier expansions of singular (or near-singular) time-dependent partial differential equations. Results for the 1D Burgers and the 3D incompressible Euler equations are used to illustrate the construction. We also present, for the 1D Burgers and the 3D Euler equations, a simple and efficient recursive algorithm for calculating the higher order terms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical computation of solutions of the critical nonlinear Schrodinger equation after the singularity

We present numerical results for the solution of the 1D critical nonlinear Schrodinger with periodic boundary conditions and initial data that give rise to a finite time singularity. We construct, through the Mori-Zwanzig formalism, a reduced model which allows us to follow the solution after the formation of the singularity. The computed post-singularity solution exhibits the same characteristics as the post-singularity solutions constructed recently by Terence Tao.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Application of the t-model of optimal prediction to the estimation of the rate of decay of solutions of the Euler equations in two and three dimensions

The "t-model" for dimensional reduction is applied to the estimation of the rate of decay of solutions of the Burgers equation and of the Euler equations in two and three space dimensions. The model was first derived in a statistical mechanics context, but here we analyze it purely as a numerical tool and prove its convergence. In the Burgers case the model captures the rate of decay exactly, as was already previously shown. For the Euler equations in two space dimensions, the model preserves energy as it should. In three dimensions, we find a power law decay in time and observe a temporal intermittency.

0 citations0 reviewsNo code linkedOpen access