In this paper, we investigate some micromechanical aspects of elasto-plasticity in heterogeneous geomaterials. The aim is to upscale the elasto-plastic behavior for a representative volume of the material which is indeed a very challenging task due to the irreversible deformations involved. Considering the plastic strains as eigen-strains allows us to employ the powerful tools offered by Continuum

Blades manufactured through flank and point milling will likely exhibit geometric variability. Gauging the aerodynamic repercussions of such variability, prior to manufacturing a component, is challenging enough, let alone trying to predict what the amplified impact of any in-service degradation will be. While rules of thumb that govern the tolerance band can be devised based on expected boundary layer

In this contribution we propose reduced order methods to fast and reliably solve parametrized optimal control problems governed by time dependent nonlinear partial differential equations. Our goal is to provide a general tool to deal with the time evolution of several nonlinear optimality systems in many-query context, where a system must be analysed for various physical and geometrical features. Optimal

To counter the volatile nature of renewable energy sources, gas networks take a vital role. But, to ensure fulfillment of contracts under these new circumstances, a vast number of possible scenarios, incorporating uncertain supply and demand, has to be simulated ahead of time. This many-query task can be accelerated by model order reduction, yet, large-scale, nonlinear, parametric, hyperbolic partial

Prediction of chaotic systems relies on a floating fusion of sensor data (observations) with a numerical model to decide on a good system trajectory and to compensate nonlinear feedback effects. Ensemble-based data assimilation (DA) is a major method for this concern depending on propagating an ensemble of perturbed model realizations.In this paper we develop an elastic, online, fault-tolerant and

We present a new theoretical and numerical framework for modelling mechanically-assisted corrosion in elastic-plastic solids. Both pitting and stress corrosion cracking (SCC) can be captured, as well as the pit-to-crack transition. Localised corrosion is assumed to be dissolution-driven and a formulation grounded upon the film rupture-dissolution-repassivation mechanism is presented to incorporate

This paper summarises a successful application of functional programming within a commercial environment. We report on experience at Accenture's Financial Services Solution Centre in London with simulating an object-oriented financial system in order to assist analysis and design. The work was part of a large IT project for an international investment bank and provides a pragmatic case study.

Molecular photo-switches as, e.g., azobenzene molecules allow, when embedded into a polymeric matrix, for photo-active polymer compounds responding mechanically when exposed to light of certain wavelength. Photo-mechanics, i.e. light-matter interaction in photo-active polymers holds great promise for, e.g., remote and contact-free activation of photo-driven actuators. In a series of earlier contributions

Modern financial market dynamics warrant detailed analysis due to their significant impact on the world. This, however, often proves intractable; massive numbers of agents, strategies and their change over time in reaction to each other leads to difficulties in both theoretical and simulational approaches. Notable work has been done on strategy dominance in stock markets with respect to the ratios

In this paper, we describe a stable finite element formulation for advection-diffusion-reaction problems that allows for robust automatic adaptive strategies to be easily implemented. We consider locally vanishing, heterogeneous, and anisotropic diffusivities, as well as advection-dominated diffusion problems. The general stabilized finite element framework was described and analyzed in arXiv:1907

Numerical integral operators of convolution type form the basis of most wave-equation-based methods for processing and imaging of seismic data. As several of these methods require the solution of an inverse problem, multiple forward and adjoint passes of the modelling operator must be performed to converge to a satisfactory solution. This work highlights the challenges that arise when implementing

Particle based methods such as the Discrete Element Method and the Lattice Spring Method may be used for describing the behaviour of isotropic linear elastic materials. However, the common bond models employed to describe the interaction between particles restrict the range of Poisson's ratio that can be represented. In this paper, to overcome the restriction, a modified bond model that includes the

We consider financial networks, where banks are connected by contracts such as debts or credit default swaps. We study the clearing problem in these systems: we want to know which banks end up in a default, and what portion of their liabilities can these defaulting banks fulfill. We analyze these networks in a sequential model where banks announce their default one at a time, and the system evolves

This paper presents a recently developed particle simulation code package PIFE-PIC, which is a novel three-dimensional (3-D) Parallel Immersed-Finite-Element (IFE) Particle-in-Cell (PIC) simulation model for particle simulations of plasma-material interactions. This framework is based on the recently developed non-homogeneous electrostatic IFE-PIC algorithm, which is designed to handle complex plasma-material

For multilayer structures in thin substrate systems, the interfacial failure is one of the most important reliability issues. The traction-separation relations (TSR) along fracture interface, which is often a complicated mixed-mode problem, is usually adopted as a representative of the adhesive interactions of a biomaterial system. However, the existing theoretical models lack complexity and are not

Assimilation of continuously streamed monitored data is an essential component of a digital twin; the assimilated data are used to ensure the digital twin is a true representation of the monitored system. One way this is achieved is by calibration of simulation models, whether data-derived or physics-based, or a combination of both. Traditional manual calibration is not possible in this context hence

In this paper, we present a case study demonstrating how dynamic and uncertain criteria can be incorporated into a multicriteria analysis with the help of discrete event simulation. The simulation guided multicriteria analysis can include both monetary and non-monetary criteria that are static or dynamic, whereas standard multi criteria analysis only deals with static criteria and cost benefit analysis

We evaluate a number of different finite element approaches for fluid-structure (contact) interaction problems against data from physical experiments. For this we take the data from experiments by Hagemeier [Mendeley Data, doi: 10.17632/mf27c92nc3.1]. This consists of trajectories of single particles falling through a highly viscous fluid and rebounding off the bottom fluid tank wall. The resulting

A new predictor-corrector type incremental algorithm is proposed for the exact construction of weighted straight skeletons of 2D general planar polygons of arbitrary complexity based on the notion of deforming polygon. In the proposed algorithm, the raw input provided by the polygon itself is enough to resolve edge collapse and edge split events. Neither the construction of a kinetic triangulation

This work proposes a novel, general and robust method of determining bond micromoduli for anisotropic linear elastic bond-based peridynamics. The problem of finding a discrete distribution of bond micromoduli that reproduces an anisotropic peridynamic stiffness tensor is cast as a least-squares problem. The proposed numerical method is able to find a distribution of bond micromoduli that is able to

We introduce a new preconditioned Newton-GMRES method for solving the nonlinear systems arising from incompressible Navier-Stokes equations. When the Reynolds number is relatively high, these systems often involve millions of degrees of freedom (DOFs), and the nonlinear systems are difficult to converge, partially due to their saddle-point structure. In this work, we propose to alleviate these issues

Ionic models, the set of ordinary differential equations (ODEs) describing the time evolution of the state of excitable cells, are the cornerstone of modeling in neuro- and cardiac electrophysiology. Modern ionic models can have tens of state variables and hundreds of tunable parameters. Fitting ionic models to experimental data, which usually covers only a limited subset of state variables, remains

The phase field fracture method has emerged as a promising computational tool for modelling a variety of problems including, since recently, hydrogen embrittlement and stress corrosion cracking. In this work, we demonstrate the potential of phase field-based multi-physics models in transforming the engineering assessment and design of structural components in hydrogen-containing environments. First

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary

Battery Energy Storage Systems (BESS) are more and more competitive due to their increasing performances and decreasing costs. Although certain battery storage technologies may be mature and reliable from a technological perspective, with further cost reductions expected, the economic concern of battery systems is still a major barrier to be overcome before BESS can be fully utilized as a mainstream

Classical approaches and procedures for testing of automated vehicles of SAE levels 1 and 2 were based on defined scenarios with specific maneuvers, depending on the function under test. For automated driving systems (ADS) of SAE level 3+, the scenario space is infinite and calling for virtual testing and verification. However, even in simulation, the generation of safety-relevant scenarios for ADS

Propagation of a fluid-driven crack in an impermeable linear elastic medium under axis-symmetric conditions is investigated in the present work. The fluid exerting the pressure inside the crack is an incompressible Newtonian one and its front is allowed to lag behind the propagating fracture tip. The tip cavity is considered as filled by fluid vapors under constant pressure having a negligible value

The continuously growing amount of monitored data in the Industry 4.0 context requires strong and reliable anomaly detection techniques. The advancement of Digital Twin technologies allows for realistic simulations of complex machinery, therefore, it is ideally suited to generate synthetic datasets for the use in anomaly detection approaches when compared to actual measurement data. In this paper,

The sequential fully implicit (SFI) scheme was introduced (Jenny et al. 2006) for solving coupled flow and transport problems. Each time step for SFI consists of an outer loop, in which there are inner Newton loops to implicitly and sequentially solve the pressure and transport sub-problems. In standard SFI, the sub-problems are usually solved with tight tolerances at every outer iteration. This can

A method (Ember) for non-stationary spatial modelling with multiple secondary variables by combining Geostatistics with Random Forests is applied to a three-dimensional Reservoir Model. It extends the Random Forest method to an interpolation algorithm retaining similar consistency properties to both Geostatistical algorithms and Random Forests. It allows embedding of simpler interpolation algorithms

High dimensional integrals are abundant in many fields of research including quantum physics. The aim of this paper is to develop efficient recursive strategies to tackle a class of high dimensional integrals having a special product structure with low order couplings, motivated by models in lattice gauge theory from quantum field theory. A novel element of this work is the potential benefit in using

This paper presents nonlinear iterative methods for the fundamental thermal radiative transfer (TRT) model defined by the time-dependent multifrequency radiative transfer (RT) equation and the material energy balance (MEB) equation. The iterative methods are based on the nonlinear projection approach and use multiple grids in photon frequency. They are formulated by the high-order RT equation on a

The ductile fracture process in porous metals due to growth and coalescence of micron scale voids is not only affected by the imposed stress state but also by the distribution of the voids and the material size effect. The objective of this work is to understand the interaction of the inter-void spacing (or ligaments) and the resultant gradient induced material size effect on void coalescence for a

This paper presents a tractable model of non-linear dynamics of market returns using a Langevin approach.Due to non-linearity of an interaction potential, the model admits regimes of both small and large return fluctuations. Langevin dynamics are mapped onto an equivalent quantum mechanical (QM) system. Borrowing ideas from supersymmetric quantum mechanics (SUSY QM), we use a parameterized ground state

The monodomain model is widely used in in-silico cardiology to describe excitation propagation in the myocardium. Frequently, operator splitting is used to decouple the stiff reaction term and the diffusion term in the monodomain model so that they can be solved separately. Commonly, the diffusion term is solved implicitly with a large time step while the reaction term is solved by using an explicit

Decentralized electricity markets are often dominated by a small set of generator companies who control the majority of the capacity. In this paper, we explore the effect of the total controlled electricity capacity by a single, or group, of generator companies can have on the average electricity price. We demonstrate this through the use of ElecSim, a simulation of a country-wide energy market. We

This manuscript presents the formulation, implementation, calibration and application of a multiscale reduced-order model to simulate a titanium panel structure subjected to thermo-mechanical loading associated with high-speed flight. The formulation is based on the eigenstrain-based reduced-order homogenization model (EHM) and further considers thermal strain as well as temperature dependent material

A transient magneto-quasistatic vector potential formulation involving nonlinear material is spatially discretized using the finite element method of first and second polynomial order. By applying a generalized Schur complement the resulting system of differential algebraic equations is reformulated into a system of ordinary differential equations (ODE). The ODE system is integrated in time using the

While non-parametric models, such as neural networks, are sufficient in the load forecasting, separate estimates of fixed and shiftable loads are beneficial to a wide range of applications such as distribution system operational planning, load scheduling, energy trading, and utility demand response programs. A semi-parametric estimation model is usually required, where cost sensitivities of demands

The propagation of sound waves in a horizontally stratified environment, a classical problem in ocean acoustics, has traditionally been calculated using normal modes. Most programs based on the normal mode model are discretized using the finite difference method (FDM). This paper proposes a Chebyshev collocation method (CCM) based on domain decomposition to solve this problem. A set of collocation

This paper proposes a hierarchical adaptive sampling scheme for passivity characterization of large-scale linear lumped macromodels. Here, large-scale is intended both in terms of dynamic order and especially number of input/output ports. Standard passivity characterization approaches based on spectral properties of associated Hamiltonian matrices are either inefficient or non-applicable for large-scale

The numerical simulation of nonlinear ultrasound is important in the treatment planning for high-intensity focused ultrasound (HIFU) therapies in the abdomen. However, the large domain sizes and generation of higher harmonics at the focus make these problems extremely computationally demanding. Numerical methods typically employ a uniform mesh fine enough the resolve the highest harmonic present in

We present a fast spectral Galerkin scheme for the discretization of boundary integral equations arising from two-dimensional Helmholtz transmission problems in multi-layered periodic structures or gratings. Employing suitably parametrized Fourier basis and excluding Rayleigh-Wood anomalies, we rigorously establish the well-posedness of both continuous and discrete problems, and prove super-algebraic

ESCAPE is a free python package and framework for creating applications for simulating and fitting of X-ray and neutron scattering data with current support for specular reflectivity, polarized neutron reflectometry, high resolution X-ray diffraction, small angle scattering with future support for off-specular scattering from structured samples with complicated morphology. Dedicated to professionals

We present a computational framework to explore the effect of microstructure and constituent properties upon the fracture toughness of fibre-reinforced polymer composites. To capture microscopic matrix cracking and fibre-matrix debonding, the framework couples the phase field fracture method and a cohesive zone model in the context of the finite element method. Virtual single-notched three point bending

A multiscale model for real-time simulation of terrain dynamics is explored. To represent the dynamics on different scales the model combines the description of soil as a continuous solid, as distinct particles and as rigid multibodies. The models are dynamically coupled to each other and to the earthmoving equipment. Agitated soil is represented by a hybrid of contacting particles and continuum solid

Machine learning with missing data has been approached in two different ways, including feature imputation where missing feature values are estimated based on observed values, and label prediction where downstream labels are learned directly from incomplete data. However, existing imputation models tend to have strong prior assumptions and cannot learn from downstream tasks, while models targeting

In electrocardiography, the "classic" inverse problem consists of finding electric potentials on a surface enclosing the heart from remote recordings on the body surface and an accurate description of the anatomy. The latter being affected by noise and obtained with limited resolution due to clinical constraints, a possibly large uncertainty may be perpetuated in the inverse reconstruction. The purpose

Data-driven quantitative defect reconstructions using ultrasonic guided waves has recently demonstrated great potential in the area of non-destructive testing. In this paper, we develop an efficient deep learning-based defect reconstruction framework, called NetInv, which recasts the inverse guided wave scattering problem as a data-driven supervised learning progress that realizes a mapping between

A physics-informed neural network is presented for poroelastic problems with coupled flow and deformation processes. The governing equilibrium and mass balance equations are discussed and specific derivations for two-dimensional cases are presented. A fully-connected deep neural network is used for training. Barry and Mercer's source problem with time-dependent fluid injection/extraction in an idealized

3D free-form surfaces with high aesthetical and functional requirements are widely used in automotive and aerospace industry. Geometric and dimensional variations of these free-form surfaces caused by inevitable uncertainties in the manufacturing process often leads to product quality issues. Failing to model the effect of non-ideal parts, i.e., parts with geometric and dimensional errors, during design

The present work proposes an extension of the third medium contact method for solving structural topology optimization problems that involve and exploit self-contact. A new regularization of the void region, which acts as the contact medium, makes the method suitable for cases with very large deformations. The proposed contact method is implemented in a second order topology optimization framework

We consider phylogeny estimation under a two-state model of sequence evolution by site substitution on a tree. In the asymptotic regime where the sequence lengths tend to infinity, we show that for any fixed $k$ no statistically consistent phylogeny estimation is possible from $k$-mer counts of the leaf sequences alone. Formally, we establish that the joint leaf distributions of $k$-mer counts on two

An isogeometric finite element formulation for geometrically and materially nonlinear Timoshenko beams is presented, which incorporates in-plane deformation of the cross-section described by two extensible director vectors. Since those directors belong to the space ${\Bbb R}^3$, a configuration can be additively updated. The developed formulation allows direct application of nonlinear three-dimensional

Machine learning with application to questions in the physical sciences has become a widely used tool, successfully applied to classification, regression and optimization tasks in many areas. Research focus mostly lies in improving the accuracy of the machine learning models in numerical predictions, while scientific understanding is still almost exclusively generated by human researchers analysing

A novel finite element framework is proposed for the numerical simulation of two phase flows with surface tension. The Level-Set (LS) method with piece-wise quadratic (P2) interpolation for the liquid-gas interface is used in order to reach higher-order convergence rates in regions with smooth interface. A balanced-force implementation of the continuum surface force model is used to take into account

This paper considers the approximation of partial differential equations with a point collocation framework based on high-order local maximum-entropy schemes (HOLMES). In this approach, smooth basis functions are computed through an optimization procedure and the strong form of the problem is directly imposed at the collocation points, reducing significantly the computational times with respect to

We employ textbook multigrid efficiency (TME), as introduced by Achi Brandt, to construct an asymptotically optimal monolithic multigrid solver for the Stokes system. The geometric multigrid solver builds upon the concept of hierarchical hybrid grids (HHG), which is extended to higher-order finite-element discretizations, and a corresponding matrix-free implementation. The computational cost of the

We propose a new approach for data-driven automated discovery of hyperelastic constitutive laws. The approach is unsupervised, i.e., it requires no stress data but only displacement and global force data, which are realistically available through mechanical testing and digital image correlation techniques; it delivers interpretable models, i.e., models that are embodied by parsimonious mathematical

The exposure of a human by magneto-quasistatic fields from wireless charging systems is to be determined from magnetic field measurements in near real-time. This requires a fast linear equations solver for the discrete Poisson system of the Co-Simulation Scalar Potential Finite Difference (Co-Sim. SPFD) scheme. Here, the use of the AmgX library on NVIDIA GPUs is presented for this task. It enables