Researcher profile

Wolfgang A. Wall

Wolfgang A. Wall contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
9topics
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)

preprint2025arXiv

Multi-Physics-Enhanced Bayesian Inverse Analysis: Information Gain from Additional Fields

Inverse analysis, such as model calibration, often suffers from a lack of informative data in complex real-world scenarios. The standard remedy, designing new experimental setups, is often costly and time-consuming, while readily available but seemingly useless data are ignored. This work proposes incorporating such data from additional physical fields into the inverse analysis, even when the forward model solves a single-physics problem. A Bayesian framework easily incorporates the additional data and quantifies the resulting uncertainty reduction. We formally introduce the proposed method, which we denote as multi-physics-enhanced Bayesian inverse analysis. Moreover, this work is the first to quantify the reduction in parameter uncertainty by comparing the information gain from the prior to the posterior when using single-physics versus multi-physics data. We demonstrate the potential of the proposed method in two exemplary applications. Our results show that even a few or noisy data points from an additional physical field can considerably increase the information gain, even when the physical field is only weakly or one-way coupled. Overall, this work proposes and promotes the future use of multi-physics-enhanced Bayesian inverse analysis as a cost- and time-saving game-changer across various fields of science and industry, particularly in medicine.

preprint2022arXiv

A Generalized Probabilistic Learning Approach for Multi-Fidelity Uncertainty Propagation in Complex Physical Simulations

Two of the most significant challenges in uncertainty quantification pertain to the high computational cost for simulating complex physical models and the high dimension of the random inputs. In applications of practical interest, both of these problems are encountered, and standard methods either fail or are not feasible. To overcome the current limitations, we present a generalized formulation of a Bayesian multi-fidelity Monte-Carlo (BMFMC) framework that can exploit lower-fidelity model versions in a small data regime. The goal of our analysis is an efficient and accurate estimation of the complete probabilistic response for high-fidelity models. BMFMC circumvents the curse of dimensionality by learning the relationship between the outputs of a reference high-fidelity model and potentially several lower-fidelity models. While the continuous formulation is mathematically exact and independent of the low-fidelity model's accuracy, we address challenges associated with the small data regime (i.e., only a small number of 50 to 300 high-fidelity model runs can be performed). Specifically, we complement the formulation with a set of informative input features at no extra cost. Despite the inaccurate and noisy information that some low-fidelity models provide, we demonstrate that accurate and certifiable estimates for the quantities of interest can be obtained for uncertainty quantification problems in high stochastic dimensions, with significantly fewer high-fidelity model runs than state-of-the-art methods for uncertainty quantification. We illustrate our approach by applying it to challenging numerical examples such as Navier-Stokes flow simulations and fluid-structure interaction problems.

preprint2022arXiv

A Mortar Finite Element Formulation for Large Deformation Lubricated Contact Problems with Smooth Transition Between Mixed, Elasto-Hydrodynamic and Full Hydrodynamic Lubrication

This work proposes a novel model and numerical formulation for lubricated contact problems describing the mutual interaction between two deformable 3D solid bodies and an interposed fluid film. The solid bodies are consistently described based on nonlinear continuum mechanics allowing for finite deformations and arbitrary constitutive laws. The fluid film is modelled as a quasi-2D flow problem on the interface between the solids governed by the averaged Reynolds equation. The averaged Reynolds equation accounts for surface roughness utilizing spatially homogenized, effective fluid parameters and for cavitation through a positivity constraint imposed on the pressure field. In contrast to existing approaches, the proposed model accounts for the co-existence of frictional contact tractions and hydrodynamic fluid tractions at every local point on the contact surface of the interacting bodies and covers the entire range from boundary lubrication to mixed, elastohydrodynamic, and eventually to full film hydrodynamic lubrication in one unified modelling framework with smooth transition between these different regimes. Critically, the model relies on a recently proposed regularization scheme for the mechanical contact constraint combining the advantages of classical penalty and Lagrange multiplier approaches by expressing the mechanical contact pressure as a function of the effective gap between the solid bodies while at the same time limiting the minimal gap value occurring at the (theoretical) limit of infinitely high contact pressures. From a physical point of view, this approach can be considered as a model for the elastic deformation of surface asperities, with a bounded magnitude depending on the interacting solids' surface roughness. A consistent and accurate model behavior is demonstrated and validated by employing several challenging and practically relevant benchmark test cases.

preprint2022arXiv

A Versatile SPH Modeling Framework for Coupled Microfluid-Powder Dynamics in Additive Manufacturing: Binder Jetting, Material Jetting, Directed Energy Deposition and Powder Bed Fusion

Many additive manufacturing (AM) technologies rely on powder feedstock, which is fused to form the final part either by melting or by chemical binding with subsequent sintering. In both cases, process stability and resulting part quality depend on dynamic interactions between powder particles and a fluid phase, i.e., molten metal or liquid binder. The present work proposes a versatile computational modeling framework for simulating such coupled microfluid-powder dynamics problems involving thermo-capillary flow and reversible phase transitions. In particular, a liquid and a gas phase are interacting with a solid phase that consists of a substrate and mobile powder particles while simultaneously considering temperature-dependent surface tension and wetting effects. In case of laser-metal interactions, the effect of rapid evaporation is incorporated through additional mechanical and thermal interface fluxes. All phase domains are spatially discretized using smoothed particle hydrodynamics. The method's Lagrangian nature is beneficial in the context of dynamically changing interface topologies. Special care is taken in the formulation of phase transitions, which is crucial for the robustness of the computational scheme. While the underlying model equations are of a very general nature, the proposed framework is especially suitable for the mesoscale modeling of various AM processes. To this end, the generality and robustness of the computational modeling framework is demonstrated by several application-motivated examples representing the specific AM processes binder jetting, material jetting, directed energy deposition, and powder bed fusion. Among others, it is shown how the dynamic impact of droplets in binder jetting or the evaporation-induced recoil pressure in powder bed fusion leads to powder motion, distortion of the powder packing structure, and powder particle ejection.

preprint2022arXiv

Analytical disk-cylinder interaction potential laws for the computational modeling of adhesive, deformable (nano)fibers

The analysis of complex fibrous systems or materials on the micro- and nanoscale, which have a high practical relevance for many technical or biological systems, requires accurate analytical descriptions of the adhesive and repulsive forces acting on the fiber surfaces. While such analytical expressions are generally needed both for theoretical studies and for computer-based simulations, the latter motivates us here to derive disk-cylinder interaction potential laws that are valid for arbitrary mutual orientations in the decisive regime of small surface separations. The chosen type of fundamental point-pair interaction follows the simple Lennard-Jones model with inverse power laws for both the adhesive van der Waals part and the steric, repulsive part. We present three different solutions, ranging from highest accuracy to the best trade-off between simplicity of the expression and sufficient accuracy for our intended use. The validity of simplifying approximations and the accuracy of the derived potential laws is thoroughly analyzed, using both numerical and analytical reference solutions for specific interaction cases. Most importantly, the correct asymptotic scaling behavior in the decisive regime of small separations is achieved, and also the theoretically predicted $(1\!/\!\sin\!α)$-angle dependence (for non-parallel cylinders) is obtained by the proposed analytical solutions. As we show in the outlook to our current research, the derived analytical disk-cylinder interaction potential laws may be used to formulate highly efficient computational models for the interaction of arbitrarily curved fibers, such that the disk represents the cross-section of the first and the cylinder a local approximation to the shape of the second fiber.

preprint2022arXiv

Asymptotically consistent and computationally efficient modeling of short-ranged molecular interactions between curved slender fibers undergoing large 3D deformations

This article proposes a novel computational modeling approach for short-ranged molecular interactions between curved slender fibers undergoing large 3D deformations, and gives a detailed overview how it fits into the framework of existing fiber or beam interaction models, either considering microscale molecular or macroscale contact effects. The direct evaluation of a molecular interaction potential between two general bodies in 3D space would require to integrate molecule densities over two 3D volumes, leading to a sixfold integral to be solved numerically. By exploiting the short-range nature of the considered class of interaction potentials as well as the fundamental kinematic assumption of undeformable fiber cross-sections, as typically applied in mechanical beam theories, a recently derived, closed-form analytical solution is applied for the interaction potential between a given section of the first fiber (slave beam) and the entire second fiber (master beam). This novel approach based on a pre-defined section-beam interaction potential (SBIP) requires only one single integration step along the slave beam length to be performed numerically. In terms of accuracy, the total beam-beam interaction potential resulting from this approach is shown to exhibit an asymptotically consistent angular and distance scaling behavior. In addition to elementary two-fiber systems, carefully chosen to verify accuracy and asymptotic consistence of the proposed SBIP approach, a potential practical application in form of adhesive nanofiber-grafted surfaces is studied. Involving a large number of helicoidal fibers undergoing large 3D deformations, arbitrary mutual fiber orientations as well as frequent local fiber pull-off and snap-into-contact events, this example demonstrates the robustness and computational efficiency of the new approach.

preprint2022arXiv

Generalized Section-Section Interaction Potentials in the Geometrically Exact Beam Theory: Modeling of Intermolecular Forces, Asymptotic Limit as Strain-Energy Function, and Formulation of Rotational Constraints

The present contribution proposes a universal framework to formulate generalized section-section interaction potentials (SSIP) within the geometrically exact beam theory. By exploiting the fundamental kinematic assumption of undeformable cross-sections, an objective (i.e., frame-invariant) description of SSIPs via a minimal set of six (translational and rotational) relative coordinates, either in spatial or in material form, is proposed. Based on work-pairing, work-conjugated section-section interaction forces and moments, either in spatial or in material form, are identified that can be consistently derived from a variational principle. Interestingly, it is shown that hyperelastic stored-energy functions relating the deformation measures and stress-resultants of the well-known geometrically exact Simo-Reissner beam theory can also be identified as SSIPs when considering the asymptotic limit of small relative distances and rotations between the interacting cross-sections. Moreover, the proposed variational problem formulation is demonstrated to be of a very general nature, thus allowing for the formulation of translational and rotational constraints between arbitrarily oriented cross-sections based on either a penalty or a Lagrange multiplier potential. Possible applications include fiber-based structures and materials in technical and biological systems, where the proposed approach allows to model short- or long-ranged inter-molecular (e.g., electrostatic, van der Waals or repulsive steric) interactions between fibers in geometrically complex arrangements and to formulate translational and rotational coupling constraints between different fibers (e.g., cross-linked polymer chains) or between fibers and a matrix phase (e.g., fiber-reinforced composites).

preprint2022arXiv

Global sensitivity analysis based on Gaussian-process metamodelling for complex biomechanical problems

Biomechanical models often need to describe very complex systems, organs or diseases, and hence also include a large number of parameters. One of the attractive features of physics-based models is that in those models (most) parameters have a clear physical meaning. Nevertheless, the determination of these parameters is often very elaborate and costly and shows a large scatter within the population. Hence, it is essential to identify the most important parameter for a particular problem at hand. In order to distinguish parameters which have a significant influence on a specific model output from non-influential parameters, we use sensitivity analysis, in particular the Sobol method as a global variance-based method. However, the Sobol method requires a large number of model evaluations, which is prohibitive for computationally expensive models. We therefore employ Gaussian processes as a metamodel for the underlying full model. Metamodelling introduces further uncertainty, which we also quantify. We demonstrate the approach by applying it to two different problems: nanoparticle-mediated drug delivery in a multiphase tumour-growth model, and arterial growth and remodelling. Even relatively small numbers of evaluations of the full model suffice to identify the influential parameters in both cases and to separate them from non-influential parameters. The approach also allows the quantification of higher-order interaction effects. We thus show that a variance-based global sensitivity analysis is feasible for computationally expensive biomechanical models. Different aspects of sensitivity analysis are covered including a transparent declaration of the uncertainties involved in the estimation process. Such a global sensitivity analysis not only helps to massively reduce costs for experimental determination of parameters but is also highly beneficial for inverse analysis of such complex models.

preprint2022arXiv

Quantitative analysis of thin metal powder layers via transmission X-ray imaging and discrete element simulation: Blade-based spreading approaches

Spreading uniform and dense layers is of paramount importance to creating high-quality components using powder bed additive manufacturing (AM). Blade-like tools are often employed for spreading powder metal feedstocks, especially in laser powder bed fusion and electron beam melting, where powders are characterized by a D50 of 30 microns or greater. Along with variations in boundary conditions introduced by the layer-wise geometry and surface topography of the printed component, stochastic interactions between the spreading tool and powder result in spatial variations of layer quality that are still not well understood. Here, to study powder spreading under conditions representative of powder bed AM, we employ a modular, mechanized apparatus to create powder layers from moderately and highly cohesive powders with a selection of blade-like spreading tools. Powder layer effective depth is spatially mapped using transmission X-ray imaging, and uniformity is quantified via a statistical approach. We first compare layer density, or the effective depth of powder layer, and show that blade geometries with a curved profile lead to increased material deposition. Second, this approach enables quantification of local fluctuations, or layer defect severity. For example, we observe that the primary benefit of a V-shaped rubber blade, as compared to a 45 degree rigid blade, lies in enabling local deflection of the blade edge to eliminate streaking from large particles, while also increasing deposition. Additionally, we employ a custom DEM simulation to elucidate the opposing roles of particle density and surface energy with a pseudo-material approach, where the balance of inertial and cohesive forces determine macro-scale powder flowability. For specific alloy densities, we find a critical surface energy beyond which layer density is greatly impaired when powder spreading is performed using a blade.

preprint2022arXiv

Quantitative analysis of thin metal powder layers via transmission X-ray imaging and discrete element simulation: Roller-based spreading approaches

A variety of tools can be used for spreading metal, ceramic, and polymer feedstocks in powder bed additive manufacturing methods. Rollers are often employed when spreading powders with limited flowability, as arises in powders comprising fine particle sizes or high surface energy materials. Here, we study roller-based powder spreading for powder bed AM using the unique combination of a purpose-built powder spreading testbed with a proven method for X-ray mapping of powder layer depth. We focus on the density and uniformity of nominally 100 micrometer thick layers of roller-spread Ti-6Al-4V and Al-10Si-Mg powders. Our results indicate that when rotation is too rapid, roller-applied shear force impedes the creation of dense and uniform layers from powders of high innate flowability, or where inertial forces driven by particle density dominate cohesive forces. Roller counter-rotation augments the uniformity of cohesive powder layers, primarily though reducing the influence of particle clusters in the flowing powder, which are otherwise shown to cause deep, trench-like streaks. Companion discrete element method (DEM) simulations further contextualize the experiments through isolation of the effects of cohesion on layer attributes. Results suggest that roller motion parameters could apply a strategic level of additional shear force to the flowing powder, thereby mitigating the clumping behavior characteristic of highly cohesive feedstocks while maintaining high layer uniformity.

preprint2020arXiv

A weakly compressible hybridizable discontinuous Galerkin formulation for fluid-structure interaction problems

A scheme for the solution of fluid-structure interaction (FSI) problems with weakly compressible flows is proposed in this work. A novel hybridizable discontinuous Galerkin (HDG) method is derived for the discretization of the fluid equations, while the standard continuous Galerkin (CG) approach is adopted for the structural problem. The chosen HDG solver combines robustness of discontinuous Galerkin (DG) approaches in advection-dominated flows with higher order accuracy and efficient implementations. Two coupling strategies are examined in this contribution, namely a partitioned Dirichlet-Neumann scheme in the context of hybrid HDG-CG discretizations and a monolithic approach based on Nitsche's method, exploiting the definition of the numerical flux and the trace of the solution to impose the coupling conditions. Numerical experiments show optimal convergence of the HDG and CG primal and mixed variables and superconvergence of the postprocessed fluid velocity. The robustness and the efficiency of the proposed weakly compressible formulation, in comparison to a fully incompressible one, are also highlighted on a selection of two and three dimensional FSI benchmark problems.

preprint2020arXiv

Directed force propagation in semiflexible networks

We consider the propagation of tension along specific filament of a semiflexible filament network in response to the application of a point force using a combination of numerical simulations and analytic theory. We find the distribution of force within the network is highly heterogenous, with a small number of fibers supporting a significant fraction of the applied load over distances of multiple mesh sizes surrounding the point of force application. We suggest that these structures may be thought of as tensile force chains, whose structure we explore via simulation. We develop self-consistent calculations of the point-force response function and introduce a transfer matrix approach to explore the decay of tension (into bending) energy and the branching of tensile force chains in the network.

preprint2019arXiv

Hybrid multigrid methods for high-order discontinuous Galerkin discretizations

The present work develops hybrid multigrid methods for high-order discontinuous Galerkin discretizations of elliptic problems. Fast matrix-free operator evaluation on tensor product elements is used to devise a computationally efficient PDE solver. The multigrid hierarchy exploits all possibilities of geometric, polynomial, and algebraic coarsening, targeting engineering applications on complex geometries. Additionally, a transfer from discontinuous to continuous function spaces is performed within the multigrid hierarchy. This does not only further reduce the problem size of the coarse-grid problem, but also leads to a discretization most suitable for state-of-the-art algebraic multigrid methods applied as coarse-grid solver. The relevant design choices regarding the selection of optimal multigrid coarsening strategies among the various possibilities are discussed with the metric of computational costs as the driving force for algorithmic selections. We find that a transfer to a continuous function space at highest polynomial degree (or on the finest mesh), followed by polynomial and geometric coarsening, shows the best overall performance. The success of this particular multigrid strategy is due to a significant reduction in iteration counts as compared to a transfer from discontinuous to continuous function spaces at lowest polynomial degree (or on the coarsest mesh). The coarsening strategy with transfer to a continuous function space on the finest level leads to a multigrid algorithm that is robust with respect to the penalty parameter of the SIPG method. Detailed numerical investigations are conducted for a series of examples ranging from academic test cases to more complex, practically relevant geometries. Performance comparisons to state-of-the-art methods from the literature demonstrate the versatility and computational efficiency of the proposed multigrid algorithms.

preprint2018arXiv

Parametric model order reduction and its application to inverse analysis of large nonlinear coupled cardiac problems

Predictive high-fidelity finite element simulations of human cardiac mechanics co\-mmon\-ly require a large number of structural degrees of freedom. Additionally, these models are often coupled with lumped-parameter models of hemodynamics. High computational demands, however, slow down model calibration and therefore limit the use of cardiac simulations in clinical practice. As cardiac models rely on several patient-specific parameters, just one solution corresponding to one specific parameter set does not at all meet clinical demands. Moreover, while solving the nonlinear problem, 90\% of the computation time is spent solving linear systems of equations. We propose a novel approach to reduce only the structural dimension of the monolithically coupled structure-windkessel system by projection onto a lower-dimensional subspace. We obtain a good approximation of the displacement field as well as of key scalar cardiac outputs even with very few reduced degrees of freedom while achieving considerable speedups. For subspace generation, we use proper orthogonal decomposition of displacement snapshots. To incorporate changes in the parameter set into our reduced order model, we provide a comparison of subspace interpolation methods. We further show how projection-based model order reduction can be easily integrated into a gradient-based optimization and demonstrate its performance in a real-world multivariate inverse analysis scenario. Using the presented projection-based model order reduction approach can significantly speed up model personalization and could be used for many-query tasks in a clinical setting.