Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
July 3, 2021
### Abstract

We introduce and analyze a stress-based formulation for Zener’s model in linear viscoelasticity. The method is aimed to tackle efficiently heterogeneous materials that admit purely elastic and viscoelastic parts in their composition. We write the mixed variational formulation of the problem in terms of a class of tensorial wave equation and obtain an energy estimate that guaranties the well-posedness of the problem through a standard Galerkin procedure. We propose and analyze mixed continuous and discontinuous Galerkin space discretizations of the problem and derive optimal error bounds for each semidiscrete solution in the corresponding energy norm. Finally, we discuss full discretization strategies for both Galerkin methods.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
March 9, 2021
### Abstract

We present a novel numerical scheme to approximate the solution map s ↦ u( s ) := 𝓛 − s f to fractional PDEs involving elliptic operators. Reinterpreting 𝓛 − s as an interpolation operator allows us to write u ( s ) as an integral including solutions to a parametrized family of local PDEs. We propose a reduced basis strategy on top of a finite element method to approximate its integrand. Unlike prior works, we deduce the choice of snapshots for the reduced basis procedure analytically. The integral is interpreted in a spectral setting to evaluate the surrogate directly. Its computation boils down to a matrix approximation L of the operator whose inverse is projected to the s -independent reduced space, where explicit diagonalization is feasible. Exponential convergence rates are proven rigorously. A second algorithm is presented to avoid inversion of L . Instead, we directly project the matrix to the subspace, where its negative fractional power is evaluated. A numerical comparison with the predecessor highlights its competitive performance.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
January 16, 2021
### Abstract

Bounds on the spectrum of Schur complements of subdomain stiffness matrices with respect to the interior variables are key ingredients of the convergence analysis of FETI (finite element tearing and interconnecting) based domain decomposition methods. Here we give bounds on the regular condition number of Schur complements of “floating” clusters arising from the discretization of 3D Laplacian on a cube decomposed into cube subdomains. The results show that the condition number of the cluster defined on a fixed domain decomposed into m × m × m cube subdomains connected by face and optionally edge averages increases proportionally to m . The estimates support scalability of unpreconditioned H-FETI-DP (hybrid FETI dual-primal) method. Though the research is most important for the solution of variational inequalities, the results of numerical experiments indicate that unpreconditioned H-FETI-DP with large clusters can be useful also for the solution of huge linear problems.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
January 16, 2021
### Abstract

Convergence of a finite element method with mass-lumping and flux upwinding is formulated for solving the immiscible two-phase flow problem in porous media. The method approximates directly the wetting phase pressure and saturation, which are the primary unknowns. Well-posedness is obtained in [7]. Theoretical convergence is proved via a compactness argument. The numerical phase saturation converges strongly to a weak solution in L 2 in space and in time whereas the numerical phase pressures converge strongly to weak solutions in L 2 in space almost everywhere in time. The proof is not straightforward because of the degeneracy of the phase mobilities and the unboundedness of the derivative of the capillary pressure.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
November 26, 2020
### Abstract

We explore algebraic strategies for numerically solving linear elliptic partial differential equations in polygonal domains. To discretize the polygon by means of structured meshes, we employ Schwarz-Christoffel conformal mappings, leading to a multiterm linear equation possibly including Hadamard products of some of the terms. This new algebraic formulation allows us to clearly distinguish between the role of the discretized operators and that of the domain meshing. Various algebraic strategies are discussed for the solution of the resulting matrix equation.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
November 14, 2020
### Abstract

This paper develops an efficient and robust solution technique for the steady Boussinesq model of non-isothermal flow using Anderson acceleration applied to a Picard iteration. After analyzing the fixed point operator associated with the nonlinear iteration to prove that certain stability and regularity properties hold, we apply the authors’ recently constructed theory for Anderson acceleration, which yields a convergence result for the Anderson accelerated Picard iteration for the Boussinesq system. The result shows that the leading term in the residual is improved by the gain in the optimization problem, but at the cost of additional higher order terms that can be significant when the residual is large. We perform numerical tests that illustrate the theory, and show that a 2-stage choice of Anderson depth can be advantageous. We also consider Anderson acceleration applied to the Newton iteration for the Boussinesq equations, and observe that the acceleration allows the Newton iteration to converge for significantly higher Rayleigh numbers that it could without acceleration, even with a standard line search.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
November 3, 2020
### Abstract

The methodology proposed in this paper bridges the gap between entropy stable and positivity-preserving discontinuous Galerkin (DG) methods for nonlinear hyperbolic problems. The entropy stability property and, optionally, preservation of local bounds for cell averages are enforced using flux limiters based on entropy conditions and discrete maximum principles, respectively. Entropy production by the (limited) gradients of the piecewise-linear DG approximation is constrained using Rusanov-type entropy viscosity. The Taylor basis representation of the entropy stabilization term reveals that it penalizes the solution gradients in a manner similar to slope limiting and requires implicit treatment to avoid severe time step restrictions. The optional application of a vertex-based slope limiter constrains the DG solution to be bounded by local maxima and minima of the cell averages. Numerical studies are performed for two scalar two-dimensional test problems with nonlinear and nonconvex flux functions.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
October 29, 2020
### Abstract

With the climate change processes over times, all professions and trades in Three Gorges Reservoir Area will be influenced. One of the biggest challenges is the risk of rising sea level. In this situation, a large number of uncertainties for climate changes will be faced in Three Gorges Reservoir Area. Therefore, it is of importance to investigate the complexity of decision making on investing in the long term rising sea level risk related projects in Three Gorges Reservoir Area. This paper investigates the sea level and the temperature as the underlying assets in Three Gorges Reservoir Area. A real option model is constructed to evaluate potential sea level rising risk. We formulate European and American real option models into a linear parabolic variational inequalities and propose a power penalty approach to solve it. Then we obtain a nonlinear parabolic equation. It shows that the nonlinear parabolic equation is unique and solvable. Also, the solutions of the nonlinear parabolic equation converge to the solutions of the parabolic variational inequalities at the rate of order O(λ − k2).$\begin{array}{} \displaystyle O(\lambda^{-\frac{k}{2}}). \end{array}$ Since the analytic solution of nonlinear parabolic equation is difficult to obtain, a fitted finite volume method is developed to solve it in case of European and American options, and the convergence of the nonlinear parabolic equation is obtained. An empirical analysis is presented to illustrate our theoretical results.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
October 4, 2020
### Abstract

The paper studies a higher order unfitted finite element method for the Stokes system posed on a surface in ℝ 3 . The method employs parametric P k - P k −1 finite element pairs on tetrahedral bulk mesh to discretize the Stokes system on embedded surface. Stability and optimal order convergence results are proved. The proofs include a complete quantification of geometric errors stemming from approximate parametric representation of the surface. Numerical experiments include formal convergence studies and an example of the Kelvin--Helmholtz instability problem on the unit sphere.