Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
August 26, 2021
### Abstract

A variational formulation is analysed for the Oseen equations written in terms of vorticity and Bernoulli pressure. The velocity is fully decoupled using the momentum balance equation, and it is later recovered by a post-process. A finite element method is also proposed, consisting in equal-order Nédélec finite elements and piecewise continuous polynomials for the vorticity and the Bernoulli pressure, respectively. The a priori error analysis is carried out in the L 2 -norm for vorticity, pressure, and velocity; under a smallness assumption either on the convecting velocity, or on the mesh parameter. Furthermore, an a posteriori error estimator is designed and its robustness and efficiency are studied using weighted norms. Finally, a set of numerical examples in 2D and 3D is given, where the error indicator serves to guide adaptive mesh refinement. These tests illustrate the behaviour of the new formulation in typical flow conditions, and also confirm the theoretical findings.

Unable to retrieve citations for this document

Retrieving citations for document...

Requires Authentication
Accessible
August 1, 2021
### Abstract

We present simplified formulæ for the analytic integration of the Newton potential of polynomials over boxes in two- and three-dimensional space. These are implemented in an easy-to-use C++ library that allows computations in arbitrary precision arithmetic which is also documented here. We describe how these results can be combined with fast multipole methods to evaluate the Newton potential of more general, non-polynomial densities.

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
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.