We consider a three-dimensional composite material made of small inclusions periodically embedded in an elastic matrix, the whole structure presents strong heterogeneities between its different components. In the general framework of linearized elasticity we show that, when the size of the microstructures tends to zero, the limit homogeneous structure presents, for some wavelengths, a negative mass density tensor. Hence we are able to rigorously justify the existence of forbidden bands, i.e., intervals of frequencies in which there is no propagation of elastic waves. In particular, we show how to compute these band gaps and we illustrate the theoretical results with some numerical simulations.
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.
This paper presents a theoretical investigation of the multiphysical phenomena that govern cortical bone behaviour. Taking into account the piezoelectricity of the collagen-apatite matrix and the electrokinetics governing the interstitial fluid movement, we adopt a multiscale approach to derive a coupled poroelastic model of cortical tissue. Following how the phenomena propagate from the microscale to the tissue scale, we are able to determine the nature of macroscopically observed electric phenomena in bone.
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.
We consider periodically heterogeneous fluid-saturated poroelastic media described by the Biot model with inertia effects. The weak and semistrong formulations for displacement, seepage and pressure fields involve three equations expressing the momentum and mass balance and the Darcy law. Using the two-scale homogenization method, we obtain the limit two-scale problem and prove the existence and uniqueness of its weak solutions. The Laplace transformation in time is used to decouple the macroscopic and microscopic scales. It is shown that the seepage velocity is eliminated from the macroscopic equations involving strain and pressure fields only. The plane harmonic wave propagation is studied using an example of layered medium. Illustrations show some influence of the orthotropy on the dispersion phenomena.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.