BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Thomas Wick

Thomas Wick contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.OCCryptography and Security

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Matrix-free Monolithic Multigrid Methods for Stokes and Generalized Stokes Problems

We consider the widely used continuous $\mathcal{Q}_{k}$-$\mathcal{Q}_{k-1}$ quadrilateral or hexahedral Taylor-Hood elements for the finite element discretization of the Stokes and generalized Stokes systems in two and three spatial dimensions. For the fast solution of the corresponding symmetric, but indefinite system of finite element equations, we propose and analyze matrix-free monolithic geometric multigrid solvers that are based on appropriately scaled Chebyshev-Jacobi smoothers. The analysis is based on results by Schöberl and Zulehner (2003). We present and discuss several numerical results for typical benchmark problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Robust preconditioning for a mixed formulation of phase-field fracture problems

In this work, we consider fracture propagation in nearly incompressible and (fully) incompressible materials using a phase-field formulation. We use a mixed form of the elasticity equation to overcome volume locking effects and develop a robust, nonlinear and linear solver scheme and preconditioner for the resulting system. The coupled variational inequality system, which is solved monolithically, consists of three unknowns: displacements, pressure, and phase-field. Nonlinearities due to coupling, constitutive laws, and crack irreversibility are solved using a combined Newton algorithm for the nonlinearities in the partial differential equation and employing a primal-dual active set strategy for the crack irreverrsibility constraint. The linear system in each Newton step is solved iteratively with a flexible generalized minimal residual method (GMRES). The key contribution of this work is the development of a problem-specific preconditioner that leverages the saddle-point structure of the displacement and pressure variable. Four numerical examples in pure solids and pressure-driven fractures are conducted on uniformly and locally refined meshes to investigate the robustness of the solver concerning the Poisson ratio as well as the discretization and regularization parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

SABMIS: Sparse approximation based blind multi-image steganography scheme

We hide grayscale secret images into a grayscale cover image, which is considered to be a challenging steganography problem. Our goal is to develop a steganography scheme with enhanced embedding capacity while preserving the visual quality of the stego-image as well as the extracted secret image, and ensuring that the stego-image is resistant to steganographic attacks. The novel embedding rule of our scheme helps to hide secret image sparse coefficients into the oversampled cover image sparse coefficients in a staggered manner. The stego-image is constructed by using ADMM to solve the LASSO formulation of the underlying minimization problem. Finally, the secret images are extracted from the constructed stego-image using the reverse of our embedding rule. Using these components together, to achieve the above mentioned competing goals, forms our most novel contribution. We term our scheme SABMIS (Sparse Approximation Blind Multi-Image Steganography). We perform extensive experiments on several standard images. By choosing the size of the secret images to be half of the of cover image, we obtain embedding capacities of 2 bpp (bits per pixel), 4 bpp, 6 bpp, and 8 bpp while embedding one, two, three, and four secret images, respectively. Our focus is on hiding multiple secret images. For the case of hiding two and three secret images, our embedding capacities are higher than all the embedding capacities obtained in the literature until now. For the case of hiding four secret images, although our capacity is slightly lower than one work, we do better on the other two goals; a) very little deterioration in the quality of the stego-images and extracted secret images, and b) inherently and designed-to-be resistant to steganographic attacks. Additionally, we demonstrate that SABMIS executes in few minutes, and show its application on two real-life problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Space-time formulation, discretization, and computational performance studies for phase-field fracture optimal control problems

The purpose of this work is the development of space-time discretization schemes for phase-field optimal control problems. Specifically in the optimal control minimization problem, a tracking-type cost functional is minimized to steer the crack via the phase-field variable into a desired pattern. To achieve such optimal solutions, Neumann type boundary conditions need to be determined. First, a time discretization of the forward problem is derived using a discontinuous Galerkin formulation. Here, a challenge is to include regularization terms and the crack irreversibility constraint. The optimal control setting is formulated by means of the Lagrangian approach from which the primal part, adjoint, tangent and adjoint Hessian are derived. Herein the overall Newton algorithm is based on a reduced approach by eliminating the state constraint, namely the displacement and phase-field unknowns, but keeping the control variable as the only unknown. From the low-order discontinuous Galerkin discretization, adjoint time-stepping schemes are finally obtained. Both our formulation and algorithmic developments are substantiated and illustrated with six numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Space-time formulation, discretization, and computations for phase-field fracture optimal control problems

The purpose of this work is the development of space-time discretization schemes for phase-field optimal control problems. First, a time discretization of the forward problem is derived using a discontinuous Galerkin formulation. Here, a challenge is to include regularization termsand the crack irreversibility constraint. The optimal control setting is formulated by means of the Lagrangian approach from which the primal part, adjoint, tangent and adjoint Hessian are derived. Herein the overall Newton algorithm is based on a reduced approach by eliminating the state constraint. From the low-order discontinuous Galerkin discretization, adjoint time-stepping schemes are finally obtained. Our algorithmic developments are substantiated and illustrated with some numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Neural network guided adjoint computations in dual weighted residual error estimation

In this work, we are concerned with neural network guided goal-oriented a posteriori error estimation and adaptivity using the dual weighted residual method. The primal problem is solved using classical Galerkin finite elements. The adjoint problem is solved in strong form with a feedforward neural network using two or three hidden layers. The main objective of our approach is to explore alternatives for solving the adjoint problem with greater potential of a numerical cost reduction. The proposed algorithm is based on the general goal-oriented error estimation theorem including both linear and nonlinear stationary partial differential equations and goal functionals. Our developments are substantiated with some numerical experiments that include comparisons of neural network computed adjoints and classical finite element solutions of the adjoints. In the programming software, the open-source library deal.II is successfully coupled with LibTorch, the PyTorch C++ application programming interface.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Bayesian estimation method for variational phase-field fracture problems

In this work, we propose a parameter estimation framework for fracture propagation problems. The fracture problem is described by a phase-field method. Parameter estimation is realized with a Bayesian framework. Here, the focus is on uncertainties arising in the solid material parameters and the critical energy release rate. A reference value (obtained on a sufficiently small mesh) as the replacement of measurement will be chosen, and their posterior distribution is obtained. Due to time- and mesh dependency of the problem, the computational costs can be high. Using Bayesian inversion, we solve the problem on a relatively coarse mesh and fit the parameters. The obtained load-displacement curve that is usually the target function is matched with the reference values. Finally, our algorithmic approach is substantiated with several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A global-local approach for hydraulic phase-field fracture in poroelastic media

In this work, phase-field modeling of hydraulic fractures in porous media is extended towards a global-local approach. Therein, the failure behavior is solely analyzed in a (small) local domain. In the surrounding medium, a simplified and linearized system of equations is solved. Both domains are coupled by Robin-type interface conditions. The fracture(s) inside the local domain are allowed to propagate and consequently both subdomains change within time. Here, a predictor-corrector strategy is adopted in which the local domain is dynamically adjusted to the current fracture pattern. The resulting framework is algorithmically described in detail and substantiated with some numerical tests.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Parallel matrix-free higher-order finite element solvers for phase-field fracture problems

Phase-field fracture models lead to variational problems that can be written as a coupled variational equality and inequality system. Numerically, such problems can be treated with Galerkin finite elements and primal-dual active set methods. Specifically, low-order and high-order finite elements may be employed, where, for the latter, only few studies exist to date. The most time-consuming part in the discrete version of the primal-dual active set (semi-smooth Newton) algorithm consists in the solutions of changing linear systems arising at each semi-smooth Newton step. We propose a new parallel matrix-free monolithic multigrid preconditioner for these systems. We provide two numerical tests, and discuss the performance of the parallel solver proposed in the paper. Furthermore, we compare our new preconditioner with a block-AMG preconditioner available in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reliability and efficiency of DWR-type a posteriori error estimates with smart sensitivity weight recovering

We derive efficient and reliable goal-oriented error estimations, and devise adaptive mesh procedures for the finite element method that are based on the localization of a posteriori estimates. In our previous work [SIAM J. Sci. Comput., 42(1), A371--A394, 2020], we showed efficiency and reliability for error estimators based on enriched finite element spaces. However, the solution of problems on a enriched finite element space is expensive. In the literature, it is well known that one can use some higher-order interpolation to overcome this bottleneck. Using a saturation assumption, we extend the proofs of efficiency and reliability to such higher-order interpolations. The results can be used to create a new family of algorithms, where one of them is tested on three numerical examples (Poisson problem, p-Laplace equation, Navier-Stokes benschmark), and is compared to our previous algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An Adaptive Global-Local Approach for Phase-Field Modeling of Anisotropic Brittle Fracture

This work addresses an efficient Global-Local approach supplemented with predictor-corrector adaptivity applied to anisotropic phase-field brittle fracture. The phase-field formulation is used to resolve the sharp crack surface topology on the anisotropic/non-uniform local state in the regularized concept. To resolve the crack phase-field by a given single preferred direction, second-order structural tensors are imposed to both the bulk and crack surface density functions. Accordingly, a split in tension and compression modes in anisotropic materials is considered. A Global-Local formulation is proposed, in which the full displacement/phase-field problem is solved on a lower (local) scale, while dealing with a purely linear elastic problem on an upper (global) scale. Robin-type boundary conditions are introduced to relax the stiff local response at the global scale and enhancing its stabilization. Another important aspect of this contribution is the development of an adaptive Global-Local approach, where a predictor-corrector scheme is designed in which the local domains are dynamically updated during the computation. To cope with different finite element discretizations at the interface between the two nested scales, a non-matching dual mortar method is formulated. Hence, more regularity is achieved on the interface. Several numerical results substantiate our developments.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An iterative staggered scheme for phase field brittle fracture propagation with stabilizing parameters

This paper concerns the analysis and implementation of a novel iterative staggered scheme for quasi-static brittle fracture propagation models, where the fracture evolution is tracked by a phase field variable. The model we consider is a two-field variational inequality system, with the phase field function and the elastic displacements of the solid material as independent variables. Using a penalization strategy, this variational inequality system is transformed into a variational equality system, which is the formulation we take as the starting point for our algorithmic developments. The proposed scheme involves a partitioning of this model into two subproblems; phase field and mechanics, with added stabilization terms to both subproblems for improved efficiency and robustness. We analyze the convergence of the proposed scheme using a fixed point argument, and find that under a natural condition, the elastic mechanical energy remains bounded, and, if the diffusive zone around crack surfaces is sufficiently thick, monotonic convergence is achieved. Finally, the proposed scheme is validated numerically with several bench-mark problems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Multigoal-oriented optimal control problems with nonlinear PDE constraints

In this work, we consider an optimal control problem subject to a nonlinear PDE constraint and apply it to the regularized $p$-Laplace equation. To this end, a reduced unconstrained optimization problem in terms of the control variable is formulated. Based on the reduced approach, we then derive an a posteriori error representation and mesh adaptivity for multiple quantities of interest. All quantities are combined to one, and then the dual-weighted residual (DWR) method is applied to this combined functional. Furthermore, the estimator allows for balancing the discretization error and the nonlinear iteration error. These developments allow us to formulate an adaptive solution strategy, which is finally substantiated via several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Optimization with nonstationary, nonlinear monolithic fluid-structure interaction

Within this work, we consider optimization settings for nonlinear, nonstationary fluid-structure interaction. The problem is formulated in a monolithic fashion using the arbitrary Lagrangian-Eulerian framework to set-up the fluid-structure forward problem. In the optimization approach, either optimal control or parameter estimation problems are treated. In the latter, the stiffness of the solid is estimated from given reference values. In the numerical solution, the optimization problem is solved with a gradient-based solution algorithm. The nonlinear subproblems of the FSI forward problem are solved with a Newton method including line search. Specifically, we will formally provide the backward-in-time running adjoint state used for gradient computations. Our algorithmic developments are demonstrated with some numerical examples as for instance extensions of the well-known fluid-structure benchmark settings and a flapping membrane test in a channel flow with elastic walls.

0 citations0 reviewsNo code linkedOpen access