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

Bedřich Sousedík

Bedřich Sousedík appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.PRNumerical Analysisphysics.ao-ph

Catalog footprint

What is connected

20works
4topics
4close 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

20 published item(s)

preprint2022arXiv

A stochastic Galerkin method with adaptive time-stepping for the Navier-Stokes equations

We study the time-dependent Navier-Stokes equations in the context of stochastic finite element discretizations. Specifically, we assume that the viscosity is a random field given in the form of a generalized polynomial chaos expansion, and we use the stochastic Galerkin method to extend the methodology from [D. A. Kay et al., \textit{SIAM J. Sci. Comput.} 32(1), pp. 111--128, 2010] into this framework. For the resulting stochastic problem, we explore the properties of the resulting stochastic solutions, and we also compare the results with that of Monte Carlo and stochastic collocation. Since the time-stepping scheme is fully implicit, we also propose strategies for efficient solution of the stochastic Galerkin linear systems using a preconditioned Krylov subspace method. The effectiveness of the stochastic Galerkin method is illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Inexact and primal multilevel FETI-DP methods: a multilevel extension and interplay with BDDC

We study a framework that allows to solve the coarse problem in the FETI-DP method approximately. It is based on the saddle-point formulation of the FETI-DP system with a block-triangular preconditioner. One of the blocks approximates the coarse problem, for which we use the multilevel BDDC method as the main tool. This strategy then naturally leads to a version of multilevel FETI-DP method, and we show that the spectra of the multilevel FETI-DP and BDDC preconditioned operators are essentially the same. The theory is illustrated by a set of numerical experiments, and we also present a few experiments when the coarse solve is approximated by algebraic multigrid.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On Surrogate Learning for Linear Stability Assessment of Navier-Stokes Equations with Stochastic Viscosity

We study linear stability of solutions to the Navier\textendash Stokes equations with stochastic viscosity. Specifically, we assume that the viscosity is given in the form of a~stochastic expansion. Stability analysis requires a solution of the steady-state Navier-Stokes equation and then leads to a generalized eigenvalue problem, from which we wish to characterize the real part of the rightmost eigenvalue. While this can be achieved by Monte Carlo simulation, due to its computational cost we study three surrogates based on generalized polynomial chaos, Gaussian process regression and a shallow neural network. The results of linear stability analysis assessment obtained by the surrogates are compared to that of Monte Carlo simulation using a set of numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stochastic Galerkin finite element method for nonlinear elasticity and application to reinforced concrete members

We develop a stochastic Galerkin finite element method for nonlinear elasticity and apply it to reinforced concrete members with random material properties. The strategy is based on the modified Newton-Raphson method, which consists of an incremental loading process and a linearization scheme applied at each load increment. We consider that the material properties are given by a stochastic expansion in the so-called generalized polynomial chaos (gPC) framework. We search the gPC expansion of the displacement, which is then used to update the gPC expansions of the stress, strain and internal forces. The proposed method is applied to a reinforced concrete beam with uncertain initial concrete modulus of elasticity and a shear wall with uncertain maximum compressive stress of concrete, and the results are compared to those of stochastic collocation and Monte Carlo methods. Since the systems of equations obtained in the linearization scheme using the stochastic Galerkin method are very large, and there are typically many load increments, we also studied iterative solution using preconditioned conjugate gradients. The efficiency of the proposed method is illustrated by a set of numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stochastic Galerkin methods for linear stability analysis of systems with parametric uncertainty

We present a method for linear stability analysis of systems with parametric uncertainty formulated in the stochastic Galerkin framework. Specifically, we assume that for a model partial differential equation, the parameter is given in the form of generalized polynomial chaos expansion. The stability analysis leads to the solution of a stochastic eigenvalue problem, and we wish to characterize the rightmost eigenvalue. We focus, in particular, on problems with nonsymmetric matrix operators, for which the eigenvalue of interest may be a complex conjugate pair, and we develop methods for their efficient solution. These methods are based on inexact, line-search Newton iteration, which entails use of preconditioned GMRES. The method is applied to linear stability analysis of Navier-Stokes equation with stochastic viscosity, its accuracy is compared to that of Monte Carlo and stochastic collocation, and the efficiency is illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic Galerkin methods for the steady-state Navier-Stokes equations

We study the steady-state Navier-Stokes equations in the context of stochastic finite element discretizations. Specifically, we assume that the viscosity is a random field given in the form of a generalized polynomial chaos expansion. For the resulting stochastic problem, we formulate the model and linearization schemes using Picard and Newton iterations in the framework of the stochastic Galerkin method, and we explore properties of the resulting stochastic solutions. We also propose a preconditioner for solving the linear systems of equations arising at each step of the stochastic (Galerkin) nonlinear iteration and demonstrate its effectiveness for solving a set of benchmark problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

BDDC for Mixed-Hybrid Formulation of Flow in Porous Media with Combined Mesh Dimensions

We extend the Balancing Domain Decomposition by Constraints (BDDC) method to flows in porous media discretised by mixed-hybrid finite elements with combined mesh dimensions. Such discretisations appear when major geological fractures are modelled by 1D or 2D elements inside three-dimensional domains. In this set-up, the global problem as well as the substructure problems have a symmetric saddle-point structure, containing a `penalty' block due to the combination of meshes. We show that the problem can be reduced by means of iterative substructuring to an interface problem, which is symmetric and positive definite. The interface problem can thus be solved by conjugate gradients with the BDDC method as a preconditioner. A parallel implementation of this algorithm is incorporated into an existing software package for subsurface flow simulations. We study the performance of the iterative solver on several academic and real-world problems. Numerical experiments illustrate its efficiency and scalability.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse subspace iteration for spectral stochastic finite element methods

We study random eigenvalue problems in the context of spectral stochastic finite elements. In particular, given a parameter-dependent, symmetric positive-definite matrix operator, we explore the performance of algorithms for computing its eigenvalues and eigenvectors represented using polynomial chaos expansions. We formulate a version of stochastic inverse subspace iteration, which is based on the stochastic Galerkin finite element method, and we compare its accuracy with that of Monte Carlo and stochastic collocation methods. The coefficients of the eigenvalue expansions are computed from a stochastic Rayleigh quotient. Our approach allows the computation of interior eigenvalues by deflation methods, and we can also compute the coefficients of multiple eigenvectors using a stochastic variant of the modified Gram-Schmidt process. The effectiveness of the methods is illustrated by numerical experiments on benchmark problems arising from vibration analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adaptive-Multilevel BDDC and its parallel implementation

We combine the adaptive and multilevel approaches to the BDDC and formulate a method which allows an adaptive selection of constraints on each decomposition level. We also present a strategy for the solution of local eigenvalue problems in the adaptive algorithm using the LOBPCG method with a preconditioner based on standard components of the BDDC. The effectiveness of the method is illustrated on several engineering problems. It appears that the Adaptive-Multilevel BDDC algorithm is able to effectively detect troublesome parts on each decomposition level and improve convergence of the method. The developed open-source parallel implementation shows a good scalability as well as applicability to very large problems and core counts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hierarchical Schur complement preconditioner for the stochastic Galerkin finite element methods

Use of the stochastic Galerkin finite element methods leads to large systems of linear equations obtained by the discretization of tensor product solution spaces along their spatial and stochastic dimensions. These systems are typically solved iteratively by a Krylov subspace method. We propose a preconditioner which takes an advantage of the recursive hierarchy in the structure of the global matrices. In particular, the matrices posses a recursive hierarchical two-by-two structure, with one of the submatrices block diagonal. Each one of the diagonal blocks in this submatrix is closely related to the deterministic mean-value problem, and the action of its inverse is in the implementation approximated by inner loops of Krylov iterations. Thus our hierarchical Schur complement preconditioner combines, on each level in the approximation of the hierarchical structure of the global matrix, the idea of Schur complement with loops for a number of mutually independent inner Krylov iterations, and several matrix-vector multiplications for the off-diagonal blocks. Neither the global matrix, nor the matrix of the preconditioner need to be formed explicitly. The ingredients include only the number of stiffness matrices from the truncated Karhunen-Loève expansion and a good preconditioned for the mean-value deterministic problem. We provide a condition number bound for a model elliptic problem and the performance of the method is illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nested BDDC for a saddle-point problem

We propose a Nested BDDC for a class of saddle-point problems. The method solves for both flux and pressure variables. The fluxes are resolved in three-steps: the coarse solve is followed by subdomain solves, and last we look for a divergence-free flux correction and pressure variables using conjugate gradients with a Multilevel BDDC preconditioner. Because the coarse solve in the first step has the same structure as the original problem, we can use this procedure recursively and solve (a hierarchy of) coarse problems only approximately, utilizing the coarse problems known from the BDDC. The resulting algorithm thus first performs several upscaling steps, and then solves a hierarchy of problems that have the same structure but increase in size while sweeping down the levels, using the same components in the first and in the third step on each level, and also reusing the components from the higher levels. Because the coarsening can be quite aggressive, the number of levels can be kept small and the additional computational cost is significantly reduced due to the reuse of the components. We also provide the condition number bound and numerical experiments confirming the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Truncated hierarchical preconditioning for the stochastic Galerkin FEM

Stochastic Galerkin finite element discretizations of partial differential equations with coefficients characterized by arbitrary distributions lead, in general, to fully block dense linear systems. We propose two novel strategies for constructing preconditioners for these systems to be used with Krylov subspace iterative solvers. In particular, we present a variation on of the hierarchical Schur complement preconditioner, developed recently by the authors, and an adaptation of the symmetric block Gauss-Seidel method. Both preconditioners take advantage of the hierarchical structure of global stochastic Galerkin matrices, and also, when applicable, of the decay of the norms of the stiffness matrices obtained from the polynomial chaos expansion of the coefficients. This decay allows to truncate the matrix-vector multiplications in the action of the preconditioners. Also, throughout the global matrix hierarchy, we approximate solves with certain submatrices by the associated diagonal block solves. The preconditioners thus require only a limited number of stiffness matrices obtained from the polynomial chaos expansion of the coefficients, and a preconditioner for the diagonal blocks of the global matrix. The performance is illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parallel implementation of Multilevel BDDC

In application of the Balancing Domain Decomposition by Constraints (BDDC) to a case with many substructures, solving the coarse problem exactly becomes the bottleneck which spoils scalability of the solver. However, it is straightforward for BDDC to substitute the exact solution of the coarse problem by another step of BDDC method with subdomains playing the role of elements. In this way, the algorithm of three-level BDDC method is obtained. If this approach is applied recursively, multilevel BDDC method is derived. We present a detailed description of a recently developed parallel implementation of this algorithm. The implementation is applied to an engineering problem of linear elasticity and a benchmark problem of Stokes flow in a cavity. Results by the multilevel approach are compared to those by the standard (two-level) BDDC method.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Adaptive BDDC in Three Dimensions

The adaptive BDDC method is extended to the selection of face constraints in three dimensions. A new implementation of the BDDC method is presented based on a global formulation without an explicit coarse problem, with massive parallelism provided by a multifrontal solver. Constraints are implemented by a projection and sparsity of the projected operator is preserved by a generalized change of variables. The effectiveness of the method is illustrated on several engineering problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Application of the parallel BDDC preconditioner to the Stokes flow

A parallel implementation of the Balancing Domain Decomposition by Constraints (BDDC) method is described. It is based on formulation of BDDC with global matrices without explicit coarse problem. The implementation is based on the MUMPS parallel solver for computing the approximate inverse used for preconditioning. It is successfully applied to several problems of Stokes flow discretized by Taylor-Hood finite elements and BDDC is shown to be a promising method also for this class of problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulation of the 2009 Harmanli fire (Bulgaria)

We use a coupled atmosphere-fire model to simulate a fire that occurred on August 14--17, 2009, in the Harmanli region, Bulgaria. Data was obtained from GIS and satellites imagery, and from standard atmospheric data sources. Fuel data was classified in the 13 Anderson categories. For correct fire behavior, the spatial resolution of the models needed to be fine enough to resolve the essential micrometeorological effects. The simulation results are compared to available incident data. The code runs faster than real time on a cluster. The model is available from openwfm.org and it extends WRF-Fire from WRF 3.3 release.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Multispace and Multilevel BDDC

BDDC method is the most advanced method from the Balancing family of iterative substructuring methods for the solution of large systems of linear algebraic equations arising from discretization of elliptic boundary value problems. In the case of many substructures, solving the coarse problem exactly becomes a bottleneck. Since the coarse problem in BDDC has the same structure as the original problem, it is straightforward to apply the BDDC method recursively to solve the coarse problem only approximately. In this paper, we formulate a new family of abstract Multispace BDDC methods and give condition number bounds from the abstract additive Schwarz preconditioning theory. The Multilevel BDDC is then treated as a special case of the Multispace BDDC and abstract multilevel condition number bounds are given. The abstract bounds yield polylogarithmic condition number bounds for an arbitrary fixed number of levels and scalar elliptic problems discretized by finite elements in two and three spatial dimensions. Numerical experiments confirm the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

BDDC and FETI-DP under Minimalist Assumptions

The FETI-DP, BDDC and P-FETI-DP preconditioners are derived in a particulary simple abstract form. It is shown that their properties can be obtained from only on a very small set of algebraic assumptions. The presentation is purely algebraic and it does not use any particular definition of method components, such as substructures and coarse degrees of freedom. It is then shown that P-FETI-DP and BDDC are in fact the same. The FETI-DP and the BDDC preconditioned operators are of the same algebraic form, and the standard condition number bound carries over to arbitrary abstract operators of this form. The equality of eigenvalues of BDDC and FETI-DP also holds in the minimalist abstract setting. The abstract framework is explained on a standard substructuring example.

0 citations0 reviewsNo code linkedOpen access