-
A Sequential Global Programming Approach for Two-scale Optimization of Homogenized Multiphysics Problems with Application to Biot Porous Media
Authors:
Bich Ngoc Vu,
Vladimir Lukeš,
Michael Stingl,
Eduard Rohan
Abstract:
We present a new approach and an algorithm for optimizing the material configuration and behaviour of a fluid saturated porous medium in a two-scale setting. The state problem is governed by the Biot model describing the fluid-structure interaction in homogenized poroelastic structures. However, the approach is widely applicable to multiphysics problems involving several macroscopic fields where h…
▽ More
We present a new approach and an algorithm for optimizing the material configuration and behaviour of a fluid saturated porous medium in a two-scale setting. The state problem is governed by the Biot model describing the fluid-structure interaction in homogenized poroelastic structures. However, the approach is widely applicable to multiphysics problems involving several macroscopic fields where homogenization provides the relationship between the microconfigurations and the macroscopic mathematical model. The optimization variables describe the local microstructure design by virtue of the pore shape which determines the effective medium properties - the material coefficients - computed by the homogenization method. The main idea of the numerical optimization strategy consists in a) employing a precomputed database of the material coefficients associated to the geometric parameters and b) applying the sequential global programming (SGP) method for solving the problem of macroscopically optimized distribution of material coefficients. Although there are similarities with the free material optimization (FMO) approach, only effective material coefficients are considered admissible, for which a well-defined set of corresponding configurable microstructures exist. Due to the flexibility of the SGP approach, different types of microstructures with fully independent parametrizations can easily be handled. The efficiency of the concept is demonstrated by a series of numerical experiments. We show that the SGP method can handle simultaneously multiple types of microstructures with nontrivial parametrizations using a considerably low and stable number of state problems to be solved.
△ Less
Submitted 27 January, 2023;
originally announced January 2023.
-
Homogenization of the vibro-acoustic transmission on periodically perforated elastic plates with arrays of resonators
Authors:
Eduard Rohan,
Vladimír Lukeš
Abstract:
Based on our previous work, we propose a homogenized model of acoustic waves propagating through periodically perforated elastic plates with metamaterial properties due to embedded arrays of soft elastic inclusions serving for resonators. Such structures enable to suppress the acoustic transmission for selected frequency bands. Homogenization of the vibro-acoustic fluid-structure interaction probl…
▽ More
Based on our previous work, we propose a homogenized model of acoustic waves propagating through periodically perforated elastic plates with metamaterial properties due to embedded arrays of soft elastic inclusions serving for resonators. Such structures enable to suppress the acoustic transmission for selected frequency bands. Homogenization of the vibro-acoustic fluid-structure interaction problem in a 3D complex geometry of the transmission layer leads to effective transmission conditions prescribed on the acoustic meta-surface associated with the mid-plane of the Reissner-Mindlin plate. Asymptotic analysis with respect to the layer thickness, proportional to the plate thickness and to the perforation period, yields an implicit Dirichlet-to-Neumann operator defined on the homogenized metasurface. An efficient method is proposed for computing frequency-dependent effective parameters involved in the homogenized model of the layer. These can change their signs, thus modifying the acoustic impedance and the effective mass of the metasurface. The global problem of the acoustic wave propagation in a waveguide fitted with the plate is solved using the finite element method. The homogenized interface allows for a significant reduction of the computational model. Numerical illustrations are presented.
△ Less
Submitted 3 April, 2021;
originally announced April 2021.
-
Homogenization of the vibro-acoustic transmission on perforated plates
Authors:
Eduard Rohan,
Vladimír Lukeš
Abstract:
The paper deals with modelling of acoustic waves which propagate in inviscid fluids interacting with perforated elastic plates. The plate can be replaced by an interface on which transmission conditions are derived by homogenization of a problem describing vibroacoustic fluid-structure interactions in a transmission layer in which the plate is embedded. The Reissner-Mindlin theory of plates is ado…
▽ More
The paper deals with modelling of acoustic waves which propagate in inviscid fluids interacting with perforated elastic plates. The plate can be replaced by an interface on which transmission conditions are derived by homogenization of a problem describing vibroacoustic fluid-structure interactions in a transmission layer in which the plate is embedded. The Reissner-Mindlin theory of plates is adopted for periodic perforations designed by arbitrary cylindrical holes with axes orthogonal to the plate midplane. The homogenized model of the vibroacoustic transmission is obtained using the two-scale asymptotic analysis with respect to the layer thickness which is proportional to the plate thickness and to the perforation period. The nonlocal, implicit transmission conditions involve a jump in the acoustic potential and its normal one-side derivatives across the interface which represents the plate with a given thickness. The homogenized model was implemented using the finite element method and validated using direct numerical simulations of the non-homogenized problem. Numerical illustrations of the vibroacoustic transmission are presented.
△ Less
Submitted 1 January, 2019;
originally announced January 2019.
-
Multiscale finite element calculations in Python using SfePy
Authors:
Robert Cimrman,
Vladimír Lukeš,
Eduard Rohan
Abstract:
SfePy (Simple finite elements in Python) is a software for solving various kinds of problems described by partial differential equations in one, two or three spatial dimensions by the finite element method. Its source code is mostly (85\%) Python and relies on fast vectorized operations provided by the NumPy package. For a particular problem two interfaces can be used: a declarative application pr…
▽ More
SfePy (Simple finite elements in Python) is a software for solving various kinds of problems described by partial differential equations in one, two or three spatial dimensions by the finite element method. Its source code is mostly (85\%) Python and relies on fast vectorized operations provided by the NumPy package. For a particular problem two interfaces can be used: a declarative application programming interface (API), where problem description/definition files (Python modules) are used to define a calculation, and an imperative API, that can be used for interactive commands, or in scripts and libraries. After outlining the SfePy package development, the paper introduces its implementation, structure and general features. The components for defining a partial differential equation are described using an example of a simple heat conduction problem. Specifically, the declarative API of SfePy is presented in the example. To illustrate one of SfePy's main assets, the framework for implementing complex multiscale models based on the theory of homogenization, an example of a two-scale piezoelastic model is presented, showing both the mathematical description of the problem and the corresponding code.
△ Less
Submitted 1 October, 2018;
originally announced October 2018.
-
Homogenization of the fluid-saturated piezoelectric porous media
Authors:
Eduard Rohan,
Vladimír Lukeš
Abstract:
The paper is devoted to the homogenization of porous piezoelectric materials saturated by electrically inert fluid. The solid part of a representative volume element consists of the piezoelectric skeleton with embedded conductors. The pore fluid in the periodic structure can constitute a single connected domain, or an array of inclusions. Also the conducting parts are represented by several mutual…
▽ More
The paper is devoted to the homogenization of porous piezoelectric materials saturated by electrically inert fluid. The solid part of a representative volume element consists of the piezoelectric skeleton with embedded conductors. The pore fluid in the periodic structure can constitute a single connected domain, or an array of inclusions. Also the conducting parts are represented by several mutually separated connected domains, or by inclusions. Two of four possible arrangements are considered for upscaling by the homogenization method. The macroscopic model of the first type involves coefficients responsible for interactions between the electric field and the pore pressure, or the pore volume. For the second type, the electrodes can be used for controlling the electric field at the pore level, so that the deformation and the pore volume can be influenced locally. Effective constitutive coefficients are computed using characteristic responses of the microstructure. The two-scale modelling procedure is implemented numerically using the finite element method. The macroscopic strain and electric fields are used to reconstruct the corresponding local responses at the pore level. For validation of the models, these are compared with results obtained by direct numerical simulations of the heterogeneous structure; a good agreement is demonstrated, showing relevance of the two-scale numerical modelling approach.
△ Less
Submitted 16 January, 2018;
originally announced January 2018.
-
Modeling of the contrast-enhanced perfusion test in liver based on the multi-compartment flow in porous media
Authors:
Eduard Rohan,
Vladimír Lukeš,
Alena Jonášová
Abstract:
The paper deals with modeling the liver perfusion intended to improve quantitative analysis of the tissue scans provided by the contrast-enhanced computed tomography (CT). For this purpose, we developed a model of dynamic transport of the contrast fluid through the hierarchies of the perfusion trees. Conceptually, computed time-space distributions of the so-called tissue density can be compared wi…
▽ More
The paper deals with modeling the liver perfusion intended to improve quantitative analysis of the tissue scans provided by the contrast-enhanced computed tomography (CT). For this purpose, we developed a model of dynamic transport of the contrast fluid through the hierarchies of the perfusion trees. Conceptually, computed time-space distributions of the so-called tissue density can be compared with the measured data obtained from CT; such a modeling feedback can be used for model parameter identification. The blood flow is characterized at several scales for which different models are used. Flows in upper hierarchies represented by larger branching vessels are described using simple 1D models based on the Bernoulli equation extended by correction terms to respect the local pressure losses. To describe flows in smaller vessels and in the tissue parenchyma, we propose a 3D continuum model of porous medium defined in terms of hierarchically matched compartments characterized by hydraulic permeabilities. The 1D models corresponding to the portal and hepatic veins are coupled with the 3D model through point sources, or sinks. The contrast fluid saturation is governed by transport equations adapted for the 1D and 3D flow models. The complex perfusion model has been implemented using the finite element and finite volume methods. We report numerical examples computed for anatomically relevant geometries of the liver organ and of the principal vascular trees. The simulated tissue density corresponding to the CT examination output reflects a pathology modeled as a localized permeability deficiency.
△ Less
Submitted 30 May, 2016;
originally announced May 2016.
-
Numerical simulation of liver perfusion: from CT scans to FE model
Authors:
Vladimír Lukeš,
Miroslav Jiřík,
Alena Jonášová,
Eduard Rohan,
Ondřej Bublík,
Robert Cimrman
Abstract:
We use a collection of Python programs for numerical simulation of liver perfusion. We have an application for semi-automatic generation of a finite element mesh of the human liver from computed tomography scans and for reconstruction of the liver vascular structure. When the real vascular trees can not be obtained from the CT data we generate artificial trees using the constructive optimization m…
▽ More
We use a collection of Python programs for numerical simulation of liver perfusion. We have an application for semi-automatic generation of a finite element mesh of the human liver from computed tomography scans and for reconstruction of the liver vascular structure. When the real vascular trees can not be obtained from the CT data we generate artificial trees using the constructive optimization method. The generated FE mesh and vascular trees are imported into SfePy (Simple Finite Elements in Python) and numerical simulations are performed in order to get the pressure distribution and perfusion flows in the liver tissue. In the post-processing steps we calculate transport of a contrast fluid through the liver parenchyma.
△ Less
Submitted 19 December, 2014;
originally announced December 2014.