a new algorithm for high dimensional uncertainty quantification based on dimension adaptive sparse grid approximation and reduced basis methods in this work we develop an adaptive and reduced computational algorithm based on dimension adaptive sparse grid approximation and reduced basis methods for solving high dimensional uncertainty quantification uq problems in order to tackle the computational challenge of curse of dimensionality commonly faced by these problems we employ a dimension adaptive tensor product algorithm 16 and propose a verified version to enable effective removal of the stagnation phenomenon besides automatically detecting the importance and interaction of different dimensions to reduce the heavy computational cost of uq problems modelled by partial differential equations pde we adopt a weighted reduced basis method 7 and develop an adaptive greedy algorithm in combination with the previous verified algorithm for efficient construction of an accurate reduced basis approximation the efficiency and accuracy of the proposed algorithm are demonstrated by several numerical experiments

regularity analysis of metamaterial maxwell s equations with random coefficients and initial conditions abstract in this paper we develop and analyze the stochastic collocation method for solving the time dependent metamaterial maxwell s equations subject to random coefficients and random initial conditions we provide a rigorous regularity analysis of the solution with respect to the random variables to our best knowledge this is the first theoretical results derived for the stochastic metamaterial maxwell equations the rate of convergence is proved depending on the regularity of the solution numerical results are presented to confirm the theoretical analysis we also demonstrate that the backward wave propagation phenomenon still exists in random metamaterial

fuzzy stochastic fem based homogenization framework for materials with polymorphic uncertainties in the microstructure uncertainties in the macroscopic response of heterogeneous materials result from two sources the natural variability in the microstructure s geometry and the lack of sufficient knowledge regarding the microstructure the first type of uncertainty is denoted aleatoric uncertainty and may be characterized by a known probability density function the second type of uncertainty is denoted epistemic uncertainty this kind of uncertainty cannot be described using probabilistic methods models considering both sources of uncertainties are called polymorphic in the case of polymorphic uncertainties some combination of stochastic methods and fuzzy arithmetic should be used thus in the current work we examine a fuzzy stochastic finite element method based homogenization framework for materials with random inclusion sizes we analyze an experimental radii distribution of inclusions and develop a stochastic representative volume element the stochastic finite element method is used to obtain the material response in the case of random inclusion radii due to unavoidable noise in experimental data an insufficient number of samples and limited accuracy of the fitting procedure the radii distribution density cannot be obtained exactly thus it is described in terms of fuzzy location and scale parameters the influence of fuzzy input on the homogenized stress measures is analyzed

a stochastic galerkin method for first order quasilinear hyperbolic systems with uncertainty this paper is concerned with generalized polynomial chaos gpc approximation for first order quasilinear hyperbolic systems with uncertainty the one dimensional 1d hyperbolic system is first symmetrized with the aid of left eigenvector matrix of the jacobian matrix then the gpc stochastic galerkin method is applied to derive a provably symmetrically hyperbolic equations for the gpc expansion coefficients the resulting deterministic gpc galerkin system is discretized by a path conservative finite volume weno scheme in space and a third order total variation diminishing runge kutta method in time the method is further extended to two dimensional 2d quasilinear hyperbolic system with uncertainty where the symmetric hyperbolicity of the one dimensional gpc galerkin system is carried over via an operator splitting technique several numerical experiments are conducted to demonstrate the accuracy and effectiveness of the proposed gpc stochastic galerkin method

on the strongly damped wave equation we prove the existence of the universal attractor for the strongly damped semilinear wave equation in the presence of a quite general nonlinearity of critical growth when the nonlinearity is subcritical we prove the existence of an exponential attractor of optimal regularity having a basin of attraction coinciding with the whole phase space as a byproduct the universal attractor is regular and of finite fractal dimension moreover we carry out a detailed analysis of the asymptotic behavior of the solutions in dependence of the damping coefficient

variance reduction through robust design of boundary conditions for stochastic hyperbolic systems of equations we consider a hyperbolic system with uncertainty in the boundary and initial data our aim is to show that dierent boundary conditions gives dierent convergence rates of the variance of the solution this means that we can with the same knowledge of data get a more or less accurate description of the uncertainty in the solution a variety of boundary conditions are compared and both analytical and numerical estimates of the variance of the solution is presented as applications we study the eect of this technique on maxwell s equations as well as on a subsonic outow boundary for the euler equations

a multiscale preconditioner for stochastic mortar mixed finite elements the aim of this paper is to introduce a new approach to efficiently solve sequences of problems that typically arise when modeling flow in stochastic porous media the governing equations are based on darcy s law with a stochastic permeability field starting from a specified covariance relationship the log permeability is decomposed using a truncated karhunen loeve expansion multiscale mortar mixed finite element approximations are used in the spatial domain and a nonintrusive sampling method is used in the stochastic dimensions a multiscale mortar flux basis is computed for a single permeability called a training permeability that captures the main characteristics of the porous media and is used as a preconditioner for each stochastic realization we prove that the condition number of the preconditioned interface operator is independent of the subdomain mesh size and the mortar mesh size computational results confirm that our approach provides an efficient means to quantify the uncertainty for stochastic flow in porous media

a data driven stochastic collocation approach for uncertainty quantification in mems this work presents a data driven stochastic collocation approach to include the effect of uncertain design parameters during complex multi physics simulation of micro electromechanical systems mems the proposed framework comprises of two key steps first probabilistic characterization of the input uncertain parameters based on available experimental information and second propagation of these uncertainties through the predictive model to relevant quantities of interest the uncertain input parameters are modeled as independent random variables for which the distributions are estimated based on available experimental observations using a nonparametric diffusion mixing based estimator botev nonparametric density estimation via diffusion mixing technical report 2007 the diffusion based estimator derives from the analogy between the kernel density estimation kde procedure and the heat dissipation equation and constructs density estimates that are smooth and asymptotically consistent the diffusion model allows for the incorporation of the prior density and leads to an improved density estimate in comparison with the standard kde approach as demonstrated through several numerical examples following the characterization step the uncertainties are propagated to the output variables using the stochastic collocation approach based on sparse grid interpolation smolyak soviet math dokl 1963 4 240 243 the developed framework is used to study the effect of variations in young s modulus induced as a result of variations in manufacturing process parameters or heterogeneous measurements on the performance of a mems switch copyright 2010 john wiley sons ltd

attractors for processes on time dependent spaces applications to wave equations abstract for a process u t x x t acting on a one parameter family of normed spaces we present a notion of time dependent attractor based only on the minimality with respect to the pullback attraction property such an attractor is shown to be invariant whenever u t is t closed for some t 0 a much weaker property than continuity defined in the text as a byproduct we generalize the recent theory of attractors in time dependent spaces developed in di plinio et al 2011 13 finally we exploit the new framework to study the longterm behavior of wave equations with time dependent speed of propagation

minimal multi element stochastic collocation for uncertainty quantification of discontinuous functions we propose a multi element stochastic collocation method that can be applied in high dimensional parameter space for functions with discontinuities lying along manifolds of general geometries the key feature of the method is that the parameter space is decomposed into multiple elements defined by the discontinuities and thus only the minimal number of elements are utilized on each of the resulting elements the function is smooth and can be approximated using high order methods with fast convergence properties the decomposition strategy is in direct contrast to the traditional multi element approaches which define the sub domains by repeated splitting of the axes in the parameter space such methods are more prone to the curse of dimensionality because of the fast growth of the number of elements caused by the axis based splitting the present method is a two step approach firstly a discontinuity detector is used to partition parameter space into disjoint elements in each of which the function is smooth the detector uses an efficient combination of the high order polynomial annihilation technique along with adaptive sparse grids and this allows resolution of general discontinuities with a smaller number of points when the discontinuity manifold is low dimensional after partitioning an adaptive technique based on the least orthogonal interpolant is used to construct a generalized polynomial chaos surrogate on each element the adaptive technique reuses all information from the partitioning and is variance suppressing we present numerous numerical examples that illustrate the accuracy efficiency and generality of the method when compared against standard locally adaptive sparse grid methods the present method uses many fewer number of collocation samples and is more accurate

multi fidelity non intrusive polynomial chaos based on regression abstract in this paper we present a multi fidelity mf extension of non intrusive polynomial chaos based on regression point collocation for uncertainty quantification purposes the proposed method uses the principle of a global correction function from a previous similar method that uses spectral projection to estimate the coefficients due to its usage of regression to estimate the coefficients the present method offers high flexibility in the sampling and generation of the polynomial basis the method takes advantage of a nested sampling plan to create the samples for the low fidelity lf and correction expansions where the high fidelity hf samples are a subset of the lf ones to build the polynomial basis a total order or hyperbolic truncation strategy is used with a highly flexible combination of the lf and correction polynomial expansions the method is demonstrated on some artificial test problems and aerodynamic problems of the euler flow around an airfoil and common three dimensional research models in order to derive the strategies for successful mf approximation the effect of the correlation and the errors between the lf and hf functions is also studied the results show that high correlation and moderately low errors are important to improve the mf approximation s accuracy on a common research model problem the mf approach with partially converged simulations as the lf samples can successfully reduce the computational cost to about 40 for similar accuracy compared to an approach using a single hf expansion

sparse tensor discretization of elliptic spdes we propose and analyze sparse deterministic stochastic tensor galerkin finite element methods sparse sgfems for the numerical solution of elliptic partial differential equations pdes with random coefficients in a physical domain d subset mathbb r d in tensor product sgfems the variational solution to the boundary value problem is approximated in tensor product finite element spaces v gamma otimes v d where v gamma and v d denote suitable finite dimensional subspaces of the stochastic and deterministic function spaces respectively these approaches lead to sgfem algorithms of complexity o n gamma n d where n gamma dim v gamma and n d dim v d in this work we use hierarchic sequences v gamma 1 subset v gamma 2 subset ldots and v d 1 subset v d 2 subset ldots of finite dimensional spaces to approximate the law of the random solution the hierarchies of approximation spaces allow us to define sparse tensor product spaces v gamma ell hat otimes v d ell ell 1 2 dots yielding algorithms of o n gamma log n d n d log n gamma work and memory we estimate the convergence rate of sgfem for algebraic decay of the input random field karhunen loeve coefficients we give an algorithm for an input adapted a priori selection of deterministic and stochastic discretization spaces the convergence rate in terms of the total number of degrees of freedom of the proposed method is superior to monte carlo approximations numerical examples illustrate the theoretical results and demonstrate superiority of the sparse tensor product discretization proposed here versus the full tensor product approach

polynomial chaos expansions for dependent random variables abstract polynomial chaos expansions pce are well suited to quantifying uncertainty in models parameterized by independent random variables the assumption of independence leads to simple strategies for building multivariate orthonormal bases and for sampling strategies to evaluate pce coefficients in contrast the application of pce to models of dependent variables is much more challenging three approaches can be used to construct pce of models of dependent variables the first approach uses mapping methods where measure transformations such as the nataf and rosenblatt transformation can be used to map dependent random variables to independent ones however we show that this can significantly degrade performance since the jacobian of the map must be approximated a second strategy is the class of dominating support methods in these approaches a pce is built using independent random variables whose distributional support dominates the support of the true dependent joint density we provide evidence that this approach appears to produce approximations with suboptimal accuracy a third approach the novel method proposed here uses gram schmidt orthogonalization gso to numerically compute orthonormal polynomials for the dependent random variables this approach has been used successfully when solving differential equations using the intrusive stochastic galerkin method and in this paper we use gso to build pce using a non intrusive stochastic collocation method the stochastic collocation method treats the model as a black box and builds approximations of the input output map from a set of samples building pce from samples can introduce ill conditioning which does not plague stochastic galerkin methods to mitigate this ill conditioning we generate weighted leja sequences which are nested sample sets to build accurate polynomial interpolants we show that our proposed approach gso with weighted leja sequences produces pce which are orders of magnitude more accurate than pce constructed using mapping or dominating support methods

high order collocation methods for differential equations with random inputs recently there has been a growing interest in designing efficient methods for the solution of ordinary partial differential equations with random inputs to this end stochastic galerkin methods appear to be superior to other nonsampling methods and in many cases to several sampling methods however when the governing equations take complicated forms numerical implementations of stochastic galerkin methods can become nontrivial and care is needed to design robust and efficient solvers for the resulting equations on the other hand the traditional sampling methods e g monte carlo methods are straightforward to implement but they do not offer convergence as fast as stochastic galerkin methods in this paper a high order stochastic collocation approach is proposed similar to stochastic galerkin methods the collocation methods take advantage of an assumption of smoothness of the solution in random space to achieve fast convergence however the numerical implementation of stochastic collocation is trivial as it requires only repetitive runs of an existing deterministic solver similar to monte carlo methods the computational cost of the collocation methods depends on the choice of the collocation points and we present several feasible constructions one particular choice based on sparse grids depends weakly on the dimensionality of the random space and is more suitable for highly accurate computations of practical applications with large dimensional random inputs numerical examples are presented to demonstrate the accuracy and efficiency of the stochastic collocation methods

application of stochastic galerkin fem to the complete electrode model of electrical impedance tomography abstract the aim of electrical impedance tomography is to determine the internal conductivity distribution of some physical body from boundary measurements of current and voltage the most accurate forward model for impedance tomography is the complete electrode model which consists of the conductivity equation coupled with boundary conditions that take into account the electrode shapes and the contact resistances at the corresponding interfaces if the reconstruction task of impedance tomography is recast as a bayesian inference problem it is essential to be able to solve the complete electrode model forward problem with the conductivity and the contact resistances treated as a random field and random variables respectively in this work we apply a stochastic galerkin finite element method to the ensuing elliptic stochastic boundary value problem and compare the results with monte carlo simulations

parameter and state model reduction for large scale statistical inverse problems a greedy algorithm for the construction of a reduced model with reduction in both parameter and state is developed for an efficient solution of statistical inverse problems governed by partial differential equations with distributed parameters large scale models are too costly to evaluate repeatedly as is required in the statistical setting furthermore these models often have high dimensional parametric input spaces which compounds the difficulty of effectively exploring the uncertainty space we simultaneously address both challenges by constructing a projection based reduced model that accepts low dimensional parameter inputs and whose model evaluations are inexpensive the associated parameter and state bases are obtained through a greedy procedure that targets the governing equations model outputs and prior information the methodology and results are presented for groundwater inverse problems in one and two dimensions

optimal control with stochastic pde constraints and uncertain controls abstract the optimal control of problems that are constrained by partial differential equations with uncertainties and with uncertain controls is addressed the lagrangian that defines the problem is postulated in terms of stochastic functions with the control function possibly decomposed into an unknown deterministic component and a known zero mean stochastic component the extra freedom provided by the stochastic dimension in defining cost functionals is explored demonstrating the scope for controlling statistical aspects of the system response one shot stochastic finite element methods are used to find approximate solutions to control problems it is shown that applying the stochastic collocation finite element method to the formulated problem leads to a coupling between stochastic collocation points when a deterministic optimal control is considered or when moments are included in the cost functional thereby forgoing the primary advantage of the collocation method over the stochastic galerkin method for the considered problem the application of the presented methods is demonstrated through a number of numerical examples the presented framework is sufficiently general to also consider a class of inverse problems and numerical examples of this type are also presented

probabilistic collocation for period 1 limit cycle oscillations in this paper probabilistic collocation for limit cycle oscillations pclco is proposed probabilistic collocation pc is a non intrusive approach to compute the polynomial chaos description of uncertainty numerically polynomial chaos can require impractical high orders to approximate long term time integration problems due to the fast increase of required polynomial chaos order with time pclco is a pc formulation for modeling the long term stochastic behavior of dynamical systems exhibiting a periodic response i e a limit cycle oscillation lco in the pc method deterministic time series are computed at collocation points in probability space in pclco pc is applied to a time independent parametrization of the periodic response of the deterministic solves instead of to the time dependent functions themselves due to the time independent parametrization the accuracy of pclco is independent of time the approach is applied to period 1 oscillations with one main frequency subject to a random parameter numerical results are presented for the harmonic oscillator a two dof airfoil flutter model and the fluid structure interaction of an elastically mounted cylinder

three dimensional stochastic finite element method for elasto plastic bodies a new stochastic finite element method sfem is formulated for three dimensional softening elasto plastic bodies with random material properties the method is based on the karhunen loeve and polynomial chaos expansions and able to efficiently estimate complete probabilistic characteristics of the response such as moments or pdfs to reduce the computational complexity in the three dimensional setting two alterations are made with respect to the two dimensional sfem proposed earlier by the authors first a variability preserving modification of the karhunen loeve expansion is rigorously derived and applied in the stochastic discretization of random fields representing material properties second an efficient algorithm for parallel processing is developed with time consumption being the same order as for an ordinary fem rendering the proposed sfem an effective alternative to monte carlo simulation the applicability of the proposed method to stochastic analysis of strain localization is examined using monte carlo simulation then it is applied to a fault formation problem which is a recent concern of earthquake engineering ground surface layers are modelled by a softening elasto plastic body and the evolution of probabilistic characteristics of the rupture process is analysed in detail some practical observations are made regarding the nature of the fault formation from the stochastic viewpoint copyright 2001 john wiley sons ltd

on stability and monotonicity requirements of finite difference approximations of stochastic conservation laws with random viscosity on stability and monotonicity requirements of finite difference approximations of stochastic conservation laws with random viscosity

a dynamically bi orthogonal method for time dependent stochastic partial differential equations ii adaptivity and generalizations this is part ii of our paper in which we propose and develop a dynamically bi orthogonal method dybo to study a class of time dependent stochastic partial differential equations spdes whose solutions enjoy a low dimensional structure in part i of our paper 9 we derived the dybo formulation and proposed numerical algorithms based on this formulation some important theoretical results regarding consistency and bi orthogonality preservation were also established in the first part along with a range of numerical examples to illustrate the effectiveness of the dybo method in this paper we focus on the computational complexity analysis and develop an effective adaptivity strategy to add or remove modes dynamically our complexity analysis shows that the ratio of computational complexities between the dybo method and a generalized polynomial chaos method gpc is roughly of order o m n p 3 for a quadratic nonlinear spde where m is the number of mode pairs used in the dybo method and n p is the number of elements in the polynomial basis in gpc the effective dimensions of the stochastic solutions have been found to be small in many applications so we can expect m is much smaller than n p and computational savings of our dybo method against gpc are dramatic the adaptive strategy plays an essential role for the dybo method to be effective in solving some challenging problems another important contribution of this paper is the generalization of the dybo formulation for a system of time dependent spdes several numerical examples are provided to demonstrate the effectiveness of our method including the navier stokes equations and the boussinesq approximation with brownian forcing

a fast collocation method for solving stochastic integral equations based on sparse grid multiscale piecewise polynomial bases we develop a fast collocation method for solving fredholm integral equations of the second kind with stochastic loading terms it is proved that the proposed method preserves the optimal rate of convergence and has linear up to a logarithmic factor computational complexity the reduction in computational costs come from three aspects the adoption of the collocation principle the use of the sparse grid for the solution space and the employment of a truncation strategy for the coefficient matrix of the resulting linear system numerical experiments confirm the theoretical results on approximation accuracy and computational efficiency of the proposed method

risk averse structural topology optimization under random fields using stochastic expansion methods abstract this work proposes a level set based approach for solving risk averse structural topology optimization problems considering random field loading and material uncertainty the use of random fields increases the dimensionality of the stochastic domain which poses several computational challenges related to the minimization of the excess probability as a measure of risk awareness this problem is addressed both from the theoretical and numerical viewpoints first an existence result under a typical geometrical constraint on the set of admissible shapes is proved second a level set continuous approach to find the numerical solution of the problem is proposed since the considered cost functional has a discontinuous integrand the numerical approximation of the functional and its sensitivity combine an adaptive anisotropic polynomial chaos pc approach with a monte carlo mc sampling method for uncertainty propagation furthermore to address the increment of dimensionality induced by the random field an anisotropic sparse grid stochastic collocation method is used for the efficient computation of the pc coefficients a key point is that the non intrusive nature of such an approach facilitates the use of high performance computing hpc to alleviate the computational burden of the problem several numerical experiments including random field loading and material uncertainty are presented to show the feasibility of the proposal

time reversal frameness and superselection we show that appropriate superpositions of motional states are a reference frame resource that enables breaking of time reversal superselection so that two parties lacking knowledge about the other s direction of time can still communicate we identify the time reversal reference frame resource states and determine the corresponding frameness monotone which connects time reversal frameness to entanglement in contradistinction to other studies of reference frame quantum resources this is the first analysis that involves an antiunitary rather than unitary representation

robust shape optimization of continuous structures via the level set method abstract this work proposes a stochastic shape optimization method for continuous structures using the level set method such a method aims to minimize the expected compliance and its variance as measures of the structural robustness the behavior of continuous structures is modeled by linear elasticity equations with uncertain loading and material this uncertainty can be modeled using random variables with different probability distributions as well as random fields the proper problem formulation is ensured by the proof of the existence colorrev of solution under certain geometrical constraints on the set of admissible shapes the proposed method addresses the stochastic linear elasticity problem in its weak form obtaining the explicit expressions for the continuous shape derivatives some numerical examples are presented to show the effectiveness of the proposed approach

a robust bi orthogonal dynamically orthogonal method using the covariance pseudo inverse with application to stochastic flow problems abstract we develop a new robust methodology for the stochastic navier stokes equations based on the dynamically orthogonal do and bi orthogonal bo methods 1 2 3 both approaches are variants of a generalized karhunen loeve kl expansion in which both the stochastic coefficients and the spatial basis evolve according to system dynamics hence capturing the low dimensional structure of the solution the do and bo formulations are mathematically equivalent 3 but they exhibit computationally complimentary properties specifically the bo formulation may fail due to crossing of the eigenvalues of the covariance matrix while both bo and do become unstable when there is a high condition number of the covariance matrix or zero eigenvalues to this end we combine the two methods into a robust hybrid framework and in addition we employ a pseudo inverse technique to invert the covariance matrix the robustness of the proposed method stems from addressing the following issues in the do bo formulation i eigenvalue crossing we resolve the issue of eigenvalue crossing in the bo formulation by switching to the do near eigenvalue crossing using the equivalence theorem and switching back to bo when the distance between eigenvalues is larger than a threshold value ii ill conditioned covariance matrix we utilize a pseudo inverse strategy to invert the covariance matrix iii adaptivity we utilize an adaptive strategy to add remove modes to resolve the covariance matrix up to a threshold value in particular we introduce a soft threshold criterion to allow the system to adapt to the newly added removed mode and therefore avoid repetitive and unnecessary mode addition removal when the total variance approaches zero we show that the do bo formulation becomes equivalent to the evolution equation of the optimally time dependent modes 4 we demonstrate the capability of the proposed methodology with several numerical examples namely i stochastic burgers equation we analyze the performance of the method in the presence of eigenvalue crossing and zero eigenvalues ii stochastic kovasznay flow we examine the method in the presence of a singular covariance matrix and iii we examine the adaptivity of the method for an incompressible flow over a cylinder where for large stochastic forcing thirteen do bo modes are active

a variational inequality approach to stochastic boundary value problems with inequality constraints and its application to contact and elastoplasticity summary this paper is concerned with stochastic boundary value problems sbvps whose formulation involves inequality constraints a class of stochastic variational inequalities svis is defined which is well adapted to characterize the solution of specified inequality constrained sbvps a methodology for solving such svis is proposed which involves their discretization by projection onto polynomial chaos and collocation of the inequality constraints followed by the solution of a finite dimensional constrained optimization problem simulation studies in contact and elastoplasticity are provided to demonstrate the proposed framework copyright 2011 john wiley sons ltd

complex symmetric operators and applications ii a bounded linear operator t on a complex hilbert space h is called complex symmetric if t ct c where c is a conjugation an isometric antilinear involution of h we prove that t cj t where j is an auxiliary conjugation commuting with t t t we consider numerous examples including the poincare neumann singular integral bounded operator and the jordan model operator compressed shift the decomposition t cj t also extends to the class of unbounded c selfadjoint operators originally introduced by glazman in this context it provides a method for estimating the norms of the resolvents of certain unbounded operators

a scalable framework for the solution of stochastic inverse problems using a sparse grid collocation approach experimental evidence suggests that the dynamics of many physical phenomena are significantly affected by the underlying uncertainties associated with variations in properties and fluctuations in operating conditions recent developments in stochastic analysis have opened the possibility of realistic modeling of such systems in the presence of multiple sources of uncertainties these advances raise the possibility of solving the corresponding stochastic inverse problem the problem of designing estimating the evolution of a system in the presence of multiple sources of uncertainty given limited information a scalable parallel methodology for stochastic inverse design problems is developed in this article the representation of the underlying uncertainties and the resultant stochastic dependant variables is performed using a sparse grid collocation methodology a novel stochastic sensitivity method is introduced based on multiple solutions to deterministic sensitivity problems the stochastic inverse design problem is transformed to a deterministic optimization problem in a larger dimensional space that is subsequently solved using deterministic optimization algorithms the design framework relies entirely on deterministic direct and sensitivity analysis of the continuum systems thereby significantly enhancing the range of applicability of the framework for the design in the presence of uncertainty of many other systems usually analyzed with legacy codes various illustrative examples with multiple sources of uncertainty including inverse heat conduction problems in random heterogeneous media are provided to showcase the developed framework

a stochastic mixed finite element heterogeneous multiscale method for flow in porous media a computational methodology is developed to efficiently perform uncertainty quantification for fluid transport in porous media in the presence of both stochastic permeability and multiple scales in order to capture the small scale heterogeneity a new mixed multiscale finite element method is developed within the framework of the heterogeneous multiscale method hmm in the spatial domain this new method ensures both local and global mass conservation starting from a specified covariance function the stochastic log permeability is discretized in the stochastic space using a truncated karhunen loeve expansion with several random variables due to the small correlation length of the covariance function this often results in a high stochastic dimensionality therefore a newly developed adaptive high dimensional stochastic model representation technique hdmr is used in the stochastic space this results in a set of low stochastic dimensional subproblems which are efficiently solved using the adaptive sparse grid collocation method asgc numerical examples are presented for both deterministic and stochastic permeability to show the accuracy and efficiency of the developed stochastic multiscale method

a two stage adaptive stochastic collocation method on nested sparse grids for multiphase flow in randomly heterogeneous porous media a new computational method is proposed for efficient uncertainty quantification of multiphase flow in porous media with stochastic permeability for pressure estimation it combines the dimension adaptive stochastic collocation method on smolyak sparse grids and the kronrod patterson hermite nested quadrature formulas for saturation estimation an additional stage is developed in which the pressure and velocity samples are first generated by the sparse grid interpolation and then substituted into the transport equation to solve for the saturation samples to address the low regularity problem of the saturation numerical examples are presented for multiphase flow with stochastic permeability fields to demonstrate accuracy and efficiency of the proposed two stage adaptive stochastic collocation method on nested sparse grids

entanglement of pure states for a single copy an optimal local conversion strategy between any two pure states of a bipartite system is presented it is optimal in that the probability of success is the largest achievable if the parties which share the system and which can communicate classically are only allowed to act locally on it the study of optimal local conversions sheds some light on the entanglement of a single copy of a pure state we propose a quantification of such an entanglement by means of a finite minimal set of new measures from which the optimal probability of conversion follows

efficient posterior exploration of a high dimensional groundwater model from two stage markov chain monte carlo simulation and polynomial chaos expansion 1 this study reports on two strategies for accelerating posterior inference of a highly parameterized and cpu demanding groundwater flow model our method builds on previous stochastic collocation approaches e g marzouk and xiu 2009 and marzouk and najm 2009 and uses generalized polynomial chaos gpc theory and dimensionality reduction to emulate the output of a large scale groundwater flow model the resulting surrogate model is cpu efficient and serves to explore the posterior distribution at a much lower computational cost using two stage mcmc simulation the case study reported in this paper demonstrates a two to five times speed up in sampling efficiency

accurate and efficient evaluation of failure probability for partial different equations with random input data several computational challenges arise when evaluating the failure probability of a given system in the context of risk prediction or reliability analysis when the dimension of the uncertainties becomes high well established direct numerical methods can not be employed because of the curse of dimensionality many surrogate models have been proposed with the aim of reducing computational effort however most of them fail in computing an accurate failure probability when the limit state surface defined by the failure event in the probability space lacks smoothness in addition for a stochastic system modeled by partial differential equations pdes with random input only a limited number of the underlying pdes order of a few tens are affordable to solve in practice due to the considerable computational effort therefore preventing the application of many numerical methods especially for high dimensional random inputs in this work we develop hybrid and goal oriented adaptive reduced basis methods to tackle these challenges by accurately and efficiently computing the failure probability of a stochastic pde the curse of dimensionality is significantly alleviated by reduced basis approximation whose bases are constructed by goal oriented adaptation moreover an accurate evaluation of the failure probability for pde system with solution of low regularity in probability space is guaranteed by the certified a posteriori error bound for the output approximation error at the end of this paper we suitably extend our proposed method to deal with more general pde models finally we perform several numerical experiments to illustrate its computational accuracy and efficiency c 2013 elsevier b v all rights reserved

on the time dependent cattaneo law in space dimension one we consider the one dimensional wave equation e u tt u xx 1 e f u u t f u h where e e t is a decreasing function vanishing at infinity providing a model for heat conduction of cattaneo type with thermal resistance decreasing in time within the theory of processes on time dependent spaces we prove the existence of an invariant time dependent attractor which converges in a suitable sense to the attractor of the classical fourier equation u t u xx f u h formally arising in the limit t

a domain decomposition method of stochastic pdes recent advances in high performance computing systems and sensing technologies motivate computational simulations with extremely high resolution models with capabilities to quantify uncertainties for credible numerical predictions a two level domain decomposition method is reported in this investigation to devise a linear solver for the large scale system in the galerkin spectral stochastic finite element method ssfem in particular a two level scalable preconditioner is introduced in order to iteratively solve the large scale linear system in the intrusive ssfem using an iterative substructuring based domain decomposition solver the implementation of the algorithm involves solving a local problem on each subdomain that constructs the local part of the preconditioner and a coarse problem that propagates information globally among the subdomains the numerical and parallel scalabilities of the two level preconditioner are contrasted with the previously developed one level preconditioner for two dimensional flow through porous media and elasticity problems with spatially varying non gaussian material properties a distributed implementation of the parallel algorithm is carried out using mpi and petsc parallel libraries the scalabilities of the algorithm are investigated in a linux cluster

stochastic elastic plastic finite elements abstract a computational framework has been developed for simulations of the behavior of solids and structures made of stochastic elastic plastic materials uncertain elastic plastic material properties are modeled as random fields which appear as the coefficient term in the governing partial differential equation of mechanics a spectral stochastic elastic plastic finite element method with fokker planck kolmogorov equation based probabilistic constitutive integrator is proposed for solution of this non linear elastic plastic partial differential equation with stochastic coefficient to this end the random field material properties are discretized in both spatial and stochastic dimension into finite numbers of independent basic random variables using karhunen loeve expansion those random variables are then propagated through the elastic plastic constitutive rate equation using fokker planck kolmogorov equation approach to obtain the evolutionary material properties as the material plastifies the unknown displacement solution random field is then assembled using polynomial chaos as a function of known basic random variables and unknown deterministic coefficients which are obtained by minimizing the error of finite representation by galerkin technique

acceleration of the spectral stochastic fem using pod and element based discrete empirical approximation for a micromechanical model of heterogeneous materials with random geometry abstract the spectral stochastic fem with local basis functions in the stochastic domain sl fem is one of the most flexible and accurate stochastic methods however also the most computationally expensive these expenses are traditionally associated with the extra large tangent stiffness matrix and a huge number of elements which need to be re integrated in every iteration in this work we incorporate the proper orthogonal decomposition pod into the sl fem thus performing a drastic reduction of the stiffness matrix in order to reduce the integration costs by hyperreduction a novel element based modification of the discrete empirical interpolation the so called element based empirical approximation method edeam is developed and combined with the pod particular advantages of the sl fem for order reduction and hyperreduction compared to other stochastic techniques are discussed the new reduced order sl fem is applied to the computational homogenization of materials with random geometry of the microstructure i e to a general class of problems exhibiting strongly nonlinear non smooth and sometimes discontinuous dependency of the solution on some random parameters the reduced order sl fem demonstrates a high accuracy and a high solution speed whereby the solution time for the reduced order sl fem is comparable to the solution time of only one single monte carlo sample

limiting behavior for strongly damped nonlinear wave equations in this paper we investigate both the existence and the limiting behavior for the equation utt aut au t u ut where a is a sectorial operator is periodic in t and satisfies certain regularity and growth assumptions in most results on limiting behavior we will assume a has compact resolvent we consider the equation as an abstract ode defined on a paired space x x 0 1 with regard to the limiting behavior one of our principal results will be to show that if there is a bounded set in one of the spaces considered for which all points or trajectories enter into and remain then there is a set j consisting of very smooth functions defined on all of the spaces considered which is the maximum compact invariant set uniformly asymptotically stable connected and having very strong attractivity properties in all these spaces we will often show it attracts all points in a bounded set uniformly we will give a few sharper results for the case where a the work is motivated by recent papers of webb and fitzgibbon and applies techniques found in recent papers by the author

efficient solvers for a linear stochastic galerkin mixed formulation of diffusion problems with random data we introduce a stochastic galerkin mixed formulation of the steady state diffusion equation and focus on the efficient iterative solution of the saddle point systems obtained by combining standard finite element discretizations with two distinct types of stochastic basis functions so called mean based preconditioners based on fast solvers for scalar diffusion problems are introduced for use with the minimum residual method we derive eigenvalue bounds for the preconditioned system matrices and report on the efficiency of the chosen preconditioning schemes with respect to all the discretization parameters

a stochastic collocation method for elliptic partial differential equations with random input data in this paper we propose and analyze a stochastic collocation method to solve elliptic partial differential equations with random coefficients and forcing terms input data of the model the input data are assumed to depend on a finite number of random variables the method consists in a galerkin approximation in space and a collocation in the zeros of suitable tensor product orthogonal polynomials gauss points in the probability space and naturally leads to the solution of uncoupled deterministic problems as in the monte carlo approach it can be seen as a generalization of the stochastic galerkin method proposed in i babuscka r tempone and g e zouraris siam j numer anal 42 2004 pp 800 825 and allows one to treat easily a wider range of situations such as input data that depend nonlinearly on the random variables diffusivity coefficients with unbounded second moments and random variables that are correlated or even unbounded we provide a rigorous convergence analysis and demonstrate exponential convergence of the probability error with respect to the number of gauss points in each direction in the probability space under some regularity assumptions on the random input data numerical examples show the effectiveness of the method

a note on stochastic elliptic models there exist two types of commonly studied stochastic elliptic models in literature i a x u x f x and ii a x u x f x where indicates randomness the wick product and a x is a positive random process model i is widely used in engineering and physical applications while model ii is usually studied from the mathematical point of view the difference between the above two stochastic elliptic models has not been fully clarified in this work we discuss the difference between models i and ii when a x is a log normal random process we show that the difference between models i and ii is mainly characterized by a scaling factor which is an exponential function of the degree of perturbation of a x we then construct a new stochastic elliptic model iii a 1 1 u x f x which has the same scaling factor as model i after removing the divergence from the scaling factor models i and iii can be highly comparable for many cases we demonstrate this by a numerical study for a one dimensional problem

a feti preconditioned conjugate gradient method for large scale stochastic finite element problems in the spectral stochastic finite element method for analyzing an uncertain system the uncertainty is represented by a set of random variables and a quantity of interest such as the system response is considered as a function of these random variables consequently the underlying galerkin projection yields a block system of deterministic equations where the blocks are sparse but coupled the solution of this algebraic system of equations becomes rapidly challenging when the size of the physical system and or the level of uncertainty is increased this paper addresses this challenge by presenting a preconditioned conjugate gradient method for such block systems where the preconditioning step is based on the dual primal finite element tearing and interconnecting method equipped with a krylov subspace reusage technique for accelerating the iterative solution of systems with multiple and repeated right hand sides preliminary performance results on a linux cluster suggest that the proposed solution method is numerically scalable and demonstrate its potential for making the uncertainty quantification of realistic systems tractable

finite element error analysis of elliptic pdes with random coefficients and its application to multilevel monte carlo methods we consider a finite element approximation of elliptic partial differential equations with random coefficients such equations arise for example in uncertainty quantification in subsurface flow modeling models for random coefficients frequently used in these applications such as log normal random fields with exponential covariance have only very limited spatial regularity and lead to variational problems that lack uniform coercivity and boundedness with respect to the random parameter in our analysis we overcome these challenges by a careful treatment of the model problem almost surely in the random parameter which then enables us to prove uniform bounds on the finite element error in standard bochner spaces these new bounds can then be used to perform a rigorous analysis of the multilevel monte carlo method for these elliptic problems that lack full regularity and uniform coercivity and boundedness to conclude we give some numerical results that confirm the new bounds

tensor train decomposition a simple nonrecursive form of the tensor decomposition in d dimensions is presented it does not inherently suffer from the curse of dimensionality it has asymptotically the same number of parameters as the canonical decomposition but it is stable and its computation is based on low rank approximation of auxiliary unfolding matrices the new form gives a clear and convenient way to implement all basic operations efficiently a fast rounding procedure is presented as well as basic linear algebra operations examples showing the benefits of the decomposition are given and the efficiency is demonstrated by the computation of the smallest eigenvalue of a 19 dimensional operator

propagation of mechanical waves through a stochastic medium with spherical symmetry abstract we theoretically analyze the propagation of outgoing mechanical waves through an infinite isotropic elastic medium possessing spherical symmetry whose lame coefficients and density are spatial random functions characterized by well defined statistical parameters we derive the differential equation that governs the average displacement for a system whose properties depend on the radial coordinate we show that such an equation is an extended version of the well known bessel differential equation whose perturbative additional terms contain coefficients that depend directly on the squared noise intensities and the autocorrelation lengths in an exponential decay fashion we numerically solve the second order differential equation for several values of noise intensities and autocorrelation lengths and compare the corresponding displacement profiles with that of the exact analytic solution for the case of absent inhomogeneities