Researcher profile

Ilya Peshkov

Ilya Peshkov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

11 published item(s)

preprint2026arXiv

Modeling phononic band gap in microstructured solids using the Riemann-Cartan geometric framework

This paper discusses the modeling of acoustic wave fields in microstructured elastic solids within the framework of Riemann-Cartan geometry. We consider a scenario in which microstructural deformations occur significantly faster than those of the bulk material. This time-scale separation creates apparent geometric incompatibilities at the macroscopic level, even in the absence of permanent inelastic deformation or damage. We formalize this phenomenon by using a non-holonomic frame field to represent macroscopic elastic deformations and an associated torsion field to characterize the resulting geometric incompatibilities. The spatial components of the torsion tensor quantify the instantaneous geometric incompatibility of the macroscopic deformations, while its temporal components capture the inertial effects arising from the reversible energy exchange between the micro- and macro-scales. A key finding is that the model's dispersion relation predicts a complete frequency band gap. Furthermore, the governing equations exhibit a mathematical analogy to Maxwell's equations, potentially bridging the modeling of phononic and photonic metamaterials.

preprint2023arXiv

On thermodynamically compatible finite volume schemes for continuum mechanics

In this paper we present a new family of semi-discrete and fully-discrete finite volume schemes for overdetermined, hyperbolic and thermodynamically compatible PDE systems. In the following we will denote these methods as HTC schemes. In particular, we consider the Euler equations of compressible gasdynamics, as well as the more complex Godunov-Peshkov-Romenski (GPR) model of continuum mechanics, which, at the aid of suitable relaxation source terms, is able to describe nonlinear elasto-plastic solids at large deformations as well as viscous fluids as two special cases of a more general first order hyperbolic model of continuum mechanics. The main novelty of the schemes presented in this paper lies in the fact that we solve the \textit{entropy inequality} as a primary evolution equation rather than the usual total energy conservation law. Instead, total energy conservation is achieved as a mere consequence of a thermodynamically compatible discretization of all the other equations. For this, we first construct a discrete framework for the compressible Euler equations that mimics the continuous framework of Godunov's seminal paper \textit{An interesting class of quasilinear systems} of 1961 \textit{exactly} at the discrete level. All other terms in the governing equations of the more general GPR model, including non-conservative products, are judiciously discretized in order to achieve discrete thermodynamic compatibility, with the exact conservation of total energy density as a direct consequence of all the other equations. As a result, the HTC schemes proposed in this paper are provably marginally stable in the energy norm and satisfy a discrete entropy inequality by construction. We show some computational results obtained with HTC schemes in one and two space dimensions, considering both the fluid limit as well as the solid limit of the governing partial differential equations.

preprint2022arXiv

Comparison of the Symmetric Hyperbolic Thermodynamically Compatible framework with Hamiltonian mechanics of binary mixtures

How to properly describe continuum thermodynamics of binary mixtures where each constituent has its own momentum? The Symmetric Hyperbolic Thermodynamically Consistent (SHTC) framework and Hamiltonian mechanics in the form of the General Equation for Non-Equilibrium Reversible-Irreversible Coupling (GENERIC) provide two answers, which are similar but not identical, and are compared in this article. They are compared both analytically and numerically on several levels of description, varying in the amount of detail. The GENERIC equations, stemming from the Liouville equation, contain terms expressing self-advection of the relative velocity by itself, which lead to a vorticity-dependent diffusion matrix after a reduction. The SHTC equations, on the other hand, do not contain such terms. We also show how to formulate a theory of mixtures with two momenta and only one temperature that is compatible with the Liouville equation and possesses the Hamiltonian structure, including Jacobi identity.

preprint2021arXiv

A cell-centered implicit-explicit Lagrangian scheme for a unified model of nonlinear continuum mechanics on unstructured meshes

A cell-centered implicit-explicit updated Lagrangian finite volume scheme on unstructured grids is proposed for a unified first order hyperbolic formulation of continuum fluid and solid mechanics. The scheme provably respects the stiff relaxation limits of the continuous model at the fully discrete level, thus it is asymptotic preserving. Furthermore, the GCL is satisfied by a compatible discretization that makes use of a nodal solver to compute vertex-based fluxes that are used both for the motion of the computational mesh as well as for the time evolution of the governing PDEs. Second-order accuracy in space is achieved using a TVD piecewise linear reconstruction, while an implicit-explicit (IMEX) Runge-Kutta time discretization allows the scheme to obtain higher accuracy also in time. Particular care is devoted to the design of a stiff ODE solver, based on approximate analytical solutions of the governing equations, that plays a crucial role when the visco-plastic limit of the model is approached. We demonstrate the accuracy and robustness of the scheme on a wide spectrum of material responses covered by the unified continuum model that includes inviscid hydrodynamics, viscous heat conducting fluids, elastic and elasto-plastic solids in multidimensional settings.

preprint2021arXiv

Two-phase hyperbolic model for porous media saturated with a viscous fluid and its application to wavefields simulation

We derive and study a new hyperbolic two-phase model of a porous deformable medium saturated by a viscous fluid. The governing equations of the model are derived in the framework of Symmetric Hyperbolic Thermodynamically Compatible (SHTC) systems and by generalizing the unified hyperbolic model of continuum fluid and solid mechanics. Similarly to the unified model, the presented model takes into account the viscosity of the saturating fluid through a hyperbolic reformulation. The model accounts for such dissipative mechanisms as interfacial friction and viscous dissipation of the saturated fluid. Using the presented nonlinear finite-strain SHTC model, the governing equations for the propagation of small-amplitude waves in a porous medium saturated with a viscous fluid are derived. As in the conventional Biot theory of porous media, three types of waves can be found: fast and slow compression waves and shear waves. It turns out that the shear wave attenuates rapidly due to the viscosity of the saturating fluid, and this wave is difficult to see in typical test cases. However, some test cases are presented in which shear waves can be observed in the vicinity of interfaces between regions with different porosity.

preprint2020arXiv

A new continuum model for general relativistic viscous heat-conducting media

The lack of formulation of macroscopic equations for irreversible dynamics of viscous heat-conducting media compatible with the causality principle of Einstein's Special Relativity and the Euler-Lagrange structure of General Relativity is a long-lasting problem. In this paper, we propose a possible solution to this problem in the framework of SHTC equations. The approach does not rely on postulates of equilibrium irreversible thermodynamics but treats irreversible processes from the non-equilibrium point of view. Thus, each transfer process is characterized by a characteristic velocity of perturbation propagation in the non-equilibrium state, as well as by an intrinsic time/length scale of the dissipative dynamics. The resulting system of governing equations is formulated as a first-order system of hyperbolic equations with relaxation-type irreversible terms. Via a formal asymptotic analysis, we demonstrate that classical transport coefficients such as the viscosity and heat conductivity are recovered in leading terms of our theory as effective transport coefficients. Some numerical examples are presented in order to demonstrate the viability of the approach.

preprint2020arXiv

High order ADER schemes and GLM curl cleaning for a first order hyperbolic formulation of compressible flow with surface tension

In this work, we introduce two novel reformulations of a recent weakly hyperbolic model for two-phase flow with surface tension. In the model, the tracking of phase boundaries is achieved by using a vector interface field, rather than a scalar tracer, so that the surface-force stress tensor can be expressed as an algebraic function of the state variables, without requiring the computation of gradients of the tracer. An interesting and important feature of the model is that this interface field obeys a curl involution constraint, that is, the vector field is required to be curl-free at all times. The proposed modifications are intended to restore the strong hyperbolicity of the model, and are closely related to divergence-preserving numerical approaches developed in the field of numerical magnetohydrodynamics (MHD). The first strategy is based on the theory of Symmetric Hyperbolic and Thermodynamically Compatible (SHTC) systems forwarded by Godunov in the 60s and 70s and yields a modified system of governing equations which includes some symmetrisation terms, in analogy to the approach adopted later by Powell et al for the ideal MHD equations. The second technique is an extension of the hyperbolic Generalized Lagrangian Multiplier (GLM) divergence cleaning approach, forwarded by Munz et al in applications to the Maxwell and MHD equations. We solve the resulting nonconservative hyperbolic PDE systems with high order ADER Discontinuous Galerkin (DG) methods with a posteriori Finite Volume subcell limiting and carry out a set of numerical tests concerning flows dominated by surface tension as well as shock-driven flows. We also provide a new exact solution to the equations, show convergence of the schemes for orders of accuracy up to ten in space and time, and investigate the role of hyperbolicity and of curl constraints in the long-term stability of the computations.

preprint2020arXiv

High order ADER schemes for continuum mechanics

In this paper we first review the development of high order ADER finite volume and ADER discontinuous Galerkin schemes on fixed and moving meshes, since their introduction in 1999 by Toro et al. We show the modern variant of ADER based on a space-time predictor-corrector formulation in the context of ADER discontinuous Galerkin schemes with a posteriori subcell finite volume limiter on fixed and moving grids, as well as on space-time adaptive Cartesian AMR meshes. We then present and discuss the unified symmetric hyperbolic and thermodynamically compatible (SHTC) formulation of continuum mechanics developed by Godunov, Peshkov and Romenski (GPR model), which allows to describe fluid and solid mechanics in one single and unified first order hyperbolic system. In order to deal with free surface and moving boundary problems, a simple diffuse interface approach is employed, which is compatible with Eulerian schemes on fixed grids as well as direct Arbitrary-Lagrangian-Eulerian methods on moving meshes. We show some examples of moving boundary problems in fluid and solid mechanics.

preprint2020arXiv

Modeling wavefields in saturated elastic porous media based on thermodynamically compatible system theory for multiphase mixtures

A two-phase model and its application to wavefields numerical simulation are discussed in the context of modeling of compressible fluid flows in elastic porous media. The derivation of the model is based on a theory of thermodynamically compatible systems and on a model of nonlinear elastoplasticity combined with a two-phase compressible fluid flow model. The governing equations of the model include phase mass conservation laws, a total momentum conservation law, an equation for the relative velocities of the phases, an equation for mixture distortion, and a balance equation for porosity. They form a hyperbolic system of conservation equations that satisfy the fundamental laws of thermodynamics. Two types of phase interaction are introduced in the model: phase pressure relaxation to a common value and interfacial friction. Inelastic deformations also can be accounted for by source terms in the equation for distortion. The thus formulated model can be used for studying general compressible fluid flows in a deformable elastoplastic porous medium, and for modeling wave propagation in a saturated porous medium. Governing equations for small-amplitude wave propagation in a uniform porous medium saturated with a single fluid are derived. They form a first-order hyperbolic PDE system written in terms of stress and velocities and, like in Biot's model, predict three types of waves existing in real fluid-saturated porous media: fast and slow longitudinal waves and shear waves. For the numerical solution of these equations, an efficient numerical method based on a staggered-grid finite difference scheme is used. The results of solving some numerical test problems are presented and discussed

preprint2020arXiv

On Hamiltonian continuum mechanics

Continuum mechanics can be formulated in the Lagrangian frame (addressing motion of individual continuum particles) or in the Eulerian frame (addressing evolution of fields in an inertial frame). There is a canonical Hamiltonian structure in the Lagrangian frame. By transformation to the Eulerian frame we find the Poisson bracket for Eulerian continuum mechanics with deformation gradient (or the related distortion matrix). Both Lagrangian and Eulerian Hamiltonian structures are then discussed from the perspective of space-time variational formulation and by means of semidirect products and Lie algebras. Finally, we discuss the importance of the Jacobi identity in continuum mechanics and approaches to prove hyperbolicity of the evolution equations and their gauge invariance.

preprint2020arXiv

On numerical methods for hyperbolic PDE with curl involutions

In this paper we present three different numerical approaches to account for curl-type involution constraints in hyperbolic partial differential equations for continuum physics. All approaches have a direct analogy to existing and well-known divergence-preserving schemes for the Maxwell and MHD equations. The first method consists in a generalization of the Godunov-Powell terms, which means adding suitable multiples of the involution constraints to the PDE system in order to achieve the symmetric Godunov form. The second method is an extension of the generalized Lagrangian multiplier (GLM) approach of Munz et al., where the numerical errors in the involution constraint are propagated away via an augmented PDE system. The last method is an exactly involution preserving discretization, similar to the exactly divergence-free schemes for the Maxwell and MHD equations, making use of appropriately staggered meshes. We present some numerical results that allow to compare all three approaches with each other.