## Search

Now showing items 1-10 of 44

JavaScript is disabled for your browser. Some features of this site may not work without it.

AuthorTempone, Raul (32)Nobile, Fabio (8)Moraes, Alvaro (6)Vilanova, Pedro (6)Haji Ali, Abdul Lateef (5)View MoreDepartment

Applied Mathematics and Computational Science Program (44)

Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division (44)Physical Sciences and Engineering (PSE) Division (5)Mechanical Engineering Program (3)Clean Combustion Research Center (2)View MoreSubjectSampling (9)Bayesian (8)SDE (5)RDFD (2)CEM (1)View MoreTypePoster (35)Presentation (9)Year (Issue Date)2016 (44)Item AvailabilityOpen Access (44)

Now showing items 1-10 of 44

- List view
- Grid view
- Sort Options:
- Relevance
- Title Asc
- Title Desc
- Issue Date Asc
- Issue Date Desc
- Submit Date Asc
- Submit Date Desc
- Results Per Page:
- 5
- 10
- 20
- 40
- 60
- 80
- 100

New Bayesian inference method using two steps of Markov chain Monte Carlo and its application to shock tube experiment data of Furan oxidation

Kim, Daesang; El Gharamti, Iman; Bisetti, Fabrizio; Farooq, Aamir; Knio, Omar (2016-01-06) [Poster]

A new Bayesian inference method has been developed and applied to Furan shock tube experimental data for efficient statistical inferences of the Arrhenius parameters of two OH radical consumption reactions. The collected experimental data, which consist of time series signals of OH radical concentrations of 14 shock tube experiments, may require several days for MCMC computations even with the support of a fast surrogate of the combustion simulation model, while the new method reduces it to several hours by splitting the process into two steps of MCMC: the first inference of rate constants and the second inference of the Arrhenius parameters. Each step has low dimensional parameter spaces and the second step does not need the executions of the combustion simulation. Furthermore, the new approach has more flexibility in choosing the ranges of the inference parameters, and the higher speed and flexibility enable the more accurate inferences and the analyses of the propagation of errors in the measured temperatures and the alignment of the experimental time to the inference results.

Optimal design of experiments considering noisy control parameters
for the inference of Furan combustion reaction rate

Long, Quan; Kim, Daesang; Bisetti, Fabrizio; Farooq, Aamir; Tempone, Raul; Knio, Omar (2016-01-06) [Poster]

We carry out the design of experiments for the identification of the reaction
parameters in Furan combustion. The lacks of information on the true value
of the control parameters, specifically, the initial temperature and the initial TBHP concentration, are considered in the design procedure by errors-invariables models. We use two types of observables. The first is a scaler observable, i.e., half decay time of the [TBHP]. The second is the time history of the concentration.

Mathematical Model of Delamination in Composite Materials

Dia, Ben Mansour; Espath, Luis; Prudhomme, Serge; Selvakumaran, Lakshmi; Tempone, Raul (2016-01-06) [Poster]

Multi-Index Stochastic Collocation (MISC) for random elliptic PDEs

Haji Ali, Abdul Lateef; Nobile, Fabio; Tamellini, Lorenzo; Tempone, Raul (2016-01-06) [Poster]

In this work we introduce the Multi-Index Stochastic Collocation method (MISC) for computing statistics of the solution of a PDE with random data. MISC is a combination technique based on mixed differences of spatial approximations and quadratures over the space of random data. We propose an optimization procedure to select the most effective mixed differences to include in the MISC estimator: such optimization is a crucial step and allows us to build a method that, provided with sufficient solution regularity, is potentially more effective than other multi-level collocation methods already available in literature. We then provide a complexity analysis that assumes decay rates of product type for such mixed differences, showing that in the optimal case the convergence rate of MISC is only dictated by the convergence of the deterministic solver applied to a one dimensional problem. We show the effectiveness of MISC with some computational tests, comparing it with other related methods available in the literature, such as the Multi-Index and Multilevel Monte Carlo, Multilevel Stochastic Collocation, Quasi Optimal Stochastic Collocation and Sparse Composite Collocation methods.

Hybrid Chernoff Tau-Leap

Moraes, Alvaro; Tempone, Raul; Vilanova, Pedro (2016-01-06) [Poster]

Markovian pure jump processes can model many phenomena, e.g. chemical reactions at molecular level, protein transcription and translation, spread of epidemics diseases in small populations and in wireless communication networks among many others. In this work we present a novel hybrid algorithm for simulating individual trajectories which adaptively switches between the SSA and the Chernoff tauleap methods. This allows us to: (a) control the global exit probability of any simulated trajectory, (b) obtain accurate and computable estimates for the expected value of any smooth observable of the process with minimal computational work.

Multi-Index Monte Carlo (MIMC)

Haji Ali, Abdul Lateef; Nobile, Fabio; Tempone, Raul (2016-01-06) [Poster]

We propose and analyze a novel Multi-Index Monte Carlo (MIMC) method for weak approximation of stochastic models that are described in terms of differential equations either driven by random measures or with random coefficients. The MIMC method is both a stochastic version of the combination technique introduced by Zenger, Griebel and collaborators and an extension of the Multilevel Monte Carlo (MLMC) method first described by Heinrich and Giles. Inspired by Giles s seminal work, instead of using first-order differences as in MLMC, we use in MIMC high-order mixed differences to reduce the variance of the hierarchical differences dramatically. Under standard assumptions on the convergence rates of the weak error, variance and work per sample, the optimal index set turns out to be of Total Degree (TD) type. When using such sets, MIMC yields new and improved complexity results, which are natural generalizations of Giles s MLMC analysis, and which increase the domain of problem parameters for which we achieve the optimal convergence, O(TOL-2).

Computation of High-Frequency Waves with Random Uncertainty

Malenova, Gabriela; Motamed, Mohammad; Runborg, Olof; Tempone, Raul (2016-01-06) [Poster]

We consider the forward propagation of uncertainty in high-frequency waves, described by the second order wave equation with highly oscillatory initial data. The main sources of uncertainty are the wave speed and/or the initial phase and amplitude, described by a finite number of random variables with known joint probability distribution. We propose a stochastic spectral asymptotic method [1] for computing the statistics of uncertain output quantities of interest (QoIs), which are often linear or nonlinear functionals of the wave solution and its spatial/temporal derivatives. The numerical scheme combines two techniques: a high-frequency method based on Gaussian beams [2, 3], a sparse stochastic collocation method [4]. The fast spectral convergence of the proposed method depends crucially on the presence of high stochastic regularity of the QoI independent of the wave frequency. In general, the high-frequency wave solutions to parametric hyperbolic equations are highly oscillatory and non-smooth in both physical and stochastic spaces. Consequently, the stochastic regularity of the QoI, which is a functional of the wave solution, may in principle below and depend on frequency. In the present work, we provide theoretical arguments and numerical evidence that physically motivated QoIs based on local averages of |uE|2 are smooth, with derivatives in the stochastic space uniformly bounded in E, where uE and E denote the highly oscillatory wave solution and the short wavelength, respectively. This observable related regularity makes the proposed approach more efficient than current asymptotic approaches based on Monte Carlo sampling techniques.

Tight Error Bounds for Fourier Methods for Option Pricing for Exponential Levy Processes

Flores, Fabian Crocce; Häppölä, Juho; Keissling, Jonas; Tempone, Raul (2016-01-06) [Poster]

Prices of European options whose underlying asset is driven by the L´evy process are solutions to partial integrodifferential Equations (PIDEs) that generalise the Black-Scholes equation by incorporating a non-local integral term to account for the discontinuities in the asset price. The Levy -Khintchine formula provides an explicit representation of the characteristic function of a L´evy process (cf, [6]): One can derive an exact expression for the Fourier transform of the solution of the relevant PIDE. The rapid rate of convergence of the trapezoid quadrature and the speedup provide efficient methods for evaluationg option prices, possibly for a range of parameter configurations simultaneously. A couple of works have been devoted to the error analysis and parameter selection for these transform-based methods. In [5] several payoff functions are considered for a rather general set of models, whose characteristic function is assumed to be known. [4] presents the framework and theoretical approach for the error analysis, and establishes polynomial convergence rates for approximations of the option prices. [1] presents FT-related methods with curved integration contour. The classical flat FT-methods have been, on the other hand, extended for option pricing problems beyond the European framework [3]. We present a methodology for studying and bounding the error committed when using FT methods to compute option prices. We also provide a systematic way of choosing the parameters of the numerical method, minimising the error bound and guaranteeing adherence to a pre-described error tolerance. We focus on exponential L´evy processes that may be of either diffusive or pure jump in type. Our contribution is to derive a tight error bound for a Fourier transform method when pricing options under risk-neutral Levy dynamics. We present a simplified bound that separates the contributions of the payoff and of the process in an easily processed and extensible product form that is independent of the asymptotic behaviour of the option price at extreme prices and at strike parameters. We also provide a proof for the existence of optimal parameters of the numerical computation that minimise the presented error bound.

A Stochastic Multiscale Method for the Elastic Wave Equations Arising from Fiber Composites

Babuska, Ivo; Motamed, Mohammad; Tempone, Raul (2016-01-06) [Poster]

We present a stochastic multilevel global-local algorithm [1] for computing elastic waves propagating in fiber-reinforced polymer composites, where the material properties and the size and distribution of fibers in the polymer matrix may be random. The method aims at approximating statistical moments of some given quantities of interest, such as stresses, in regions of relatively small size, e.g. hot spots or zones that are deemed vulnerable to failure. For a fiber-reinforced cross-plied laminate, we introduce three problems: 1) macro; 2) meso; and 3) micro problems, corresponding to the three natural length scales: 1) the sizes of plate; 2) the tickles of plies; and 3) and the diameter of fibers. The algorithm uses a homogenized global solution to construct a local approximation that captures the microscale features of the problem. We perform numerical experiments to show the applicability and efficiency of the method.

Some Numerical Aspects on Crowd Motion - The Hughes Model

Gomes, Diogo A.; Machado Velho, Roberto (2016-01-06) [Poster]

Here, we study a crowd model proposed by R. Hughes in [5] and we describe a numerical approach to solve it. This model comprises a Fokker-Planck equation coupled with an Eikonal equation with Dirichlet or Neumann data. First, we establish a priori estimates for the solution. Second, we study radial solutions and identify a shock formation mechanism. Third, we illustrate the existence of congestion, the breakdown of the model, and the trend to the equilibrium. Finally, we propose a new numerical method and consider two numerical examples.

The export option will allow you to export the current search results of the entered query to a file. Different formats are available for download. To export the items, click on the button corresponding with the preferred download format.

By default, clicking on the export buttons will result in a download of the allowed maximum amount of items. For anonymous users the allowed maximum amount is 50 search results.

To select a subset of the search results, click "Selective Export" button and make a selection of the items you want to export. The amount of items that can be exported at once is similarly restricted as the full export.

After making a selection, click one of the export format buttons. The amount of items that will be exported is indicated in the bubble next to export format.