Related papers: Subdifferential-based implicit return-mapping oper…
This paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a…
We study an elastic Calderon-type inverse problem: recover the mass density $\rho(x)$ in a bounded domain $\Omega\subset\mathbb{R}^3$ from the Neumann-to-Dirichlet map associated with the isotropic Lam\'e system…
In this paper, we study the preferential stiffness and the crack-tip fields for an elastic porous solid of which material properties are dependent upon the density. Such a description is necessary to describe the failure that can be caused…
We consider the problem of estimating an additive regression function in an inverse regres- sion model with a convolution type operator. A smooth backfitting procedure is developed and asymptotic normality of the resulting estimator is…
Numerical simulation of fracture contact poromechanics is essential for various applications, including CO2 sequestration, geothermal energy production and underground gas storage. Modeling this problem accurately presents significant…
We herein propose a variant of the projected inexact Levenberg--Marquardt method (ILMM) for solving constrained nonsmooth equations. Since the orthogonal projection onto the feasible set may be computationally expensive, we propose a local…
In this work, a consistent viscoplasticity formulation is derived from thermodynamical principles and employing the concept of continuum elastic corrector rate. The proposed model is developed based on the principle of maximum viscoplastic…
A multiscale (micro-to-macro) analysis is proposed for the prediction of the finite strain behavior of composites with hyperelastic constituents and embedded localized damage. The composites are assumed to possess periodic microstructure…
This paper aims to develop a Newton-type method to solve a class of nonconvex composite programs. In particular, the nonsmooth part is possibly nonconvex. To tackle the nonconvexity, we develop a notion of strong prox-regularity which is…
Current methods of estimating the change in stress caused by injecting fluid into subsurface formations require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the…
We present a practical algorithm for partially relaxing multiwell energy densities such as pertain to materials undergoing martensitic phase transitions. The algorithm is based on sequential lamination, but the evolution of the…
We propose an optimization proxy in terms of iterative implicit gradient methods for solving constrained optimization problems with nonconvex loss functions. This framework can be applied to a broad range of machine learning settings,…
High-order implicit shock tracking (fitting) is a class of high-order, optimization-based numerical methods to approximate solutions of conservation laws with non-smooth features by aligning elements of the computational mesh with…
The transport of charged particles, which can be described by the Maxwell-Ampere Nernst-Planck (MANP) framework, is essential in various applications including ion channels and semiconductors. We propose a decoupled structure-preserving…
In this paper, we consider a class of systems of nonlinear equations, which arise in discretized mixed formulations of problems in solid mechanics by $hp$-finite elements. We introduce a semismooth Newton solver for this specific class and…
The derivation of the time-dependent variational equations of the Multi-Configuration Time-Dependent Hartree (MCTDH) method for high-dimensional quantum propagation is revisited from the perspective of tangent space projection methods. In…
We deal with the problem of determining the shape of an inclusion embedded in a homogenous background medium. The multifre-quency electrical impedance tomography is used to image the inclusion. For different frequencies, a current is…
The constitutive relation of the quasi-static deformation on two dimensional packed samples of polygons is calculated using molecular dynamic simulations. The stress values at which the system remains stable are bounded by a failure…
We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying…
We consider the compressible Euler system with anelastic scaling, modeling isentropic flows under the influence of gravity. In the zero-Mach-number limit, the solution of the compressible Euler system converges to a variable density…