Bio-Mechanical Model of Osteosarcoma Tumor Microenvironment: A Porous Media Approach.

Osteosarcoma is the most common malignant bone tumor in children and adolescents with a poor prognosis. To describe the progression of osteosarcoma, we expanded a system of data-driven ODE from a previous study into a system of Reaction-Diffusion-Advection (RDA) equations and coupled it with Biot equations of poroelasticity to form a bio-mechanical model. The RDA system includes the spatio-temporal information of the key components of the tumor microenvironment. The Biot equations are comprised of an equation for the solid phase, which governs the movement of the solid tumor, and an equation for the fluid phase, which relates to the motion of cells. The model predicts the total number of cells and cytokines of the tumor microenvironment and simulates the tumor's size growth. We simulated different scenarios using this model to investigate the impact of several biomedical settings on tumors' growth. The results indicate the importance of macrophages in tumors' growth. Particularly, we have observed a high co-localization of macrophages and cancer cells, and the concentration of tumor cells increases as the number of macrophages increases.

Explicit meshfree u - p w solution of the dynamic Biot formulation at large strain.

In this paper, an efficient and robust methodology to simulate saturated soils subjected to low-medium frequency dynamic loadings under large deformation regime is presented. The coupling between solid and fluid phases is solved through the dynamic reduced formulation u - p w (solid displacement - pore water pressure) of the Biot's equations. The additional novelty lies in the employment of an explicit two-steps Newmark predictor-corrector time integration scheme that enables accurate solutions of related geomechanical problems at large strain without the usually high computational cost associated with the implicit counterparts. Shape functions based on the elegant Local Maximum Entropy approach, through the Optimal Transportation Meshfree framework, are considered to solve numerically different dynamic problems in fluid saturated porous media.

Parameter-robust multiphysics algorithms for Biot model with application in brain edema simulation.

In this paper, we develop parameter-robust numerical algorithms for Biot model and apply the algorithms in brain edema simulations. By introducing an intermediate variable, we derive a multiphysics reformulation of the Biot model. Based on the reformulation, the Biot model is viewed as a generalized Stokes subproblem combining with a reaction-diffusion subproblem. Solving the two subproblems together or separately leads to a coupled or a decoupled algorithm. We conduct extensive numerical experiments to show that the two algorithms are robust with respect to the key physical parameters. The algorithms are applied to study the brain swelling caused by abnormal accumulation of cerebrospinal fluid in injured areas. The effects of the key physical parameters on brain swelling are carefully investigated. It is observed that the permeability has the biggest influence on intracranial pressure (ICP) and tissue deformation; the Young's modulus and the Poisson ratio do not affect the maximum value of ICP too much but have big influence on the tissue deformation and the developing speed of brain swelling.

Finite volume hydromechanical simulation in porous media.

Cell-centered finite volume methods are prevailing in numerical simulation of flow in porous media. However, due to the lack of cell-centered finite volume methods for mechanics, coupled flow and deformation is usually treated either by coupled finite-volume-finite element discretizations, or within a finite element setting. The former approach is unfavorable as it introduces two separate grid structures, while the latter approach loses the advantages of finite volume methods for the flow equation. Recently, we proposed a cell-centered finite volume method for elasticity. Herein, we explore the applicability of this novel method to provide a compatible finite volume discretization for coupled hydromechanic flows in porous media. We detail in particular the issue of coupling terms, and show how this is naturally handled. Furthermore, we observe how the cell-centered finite volume framework naturally allows for modeling fractured and fracturing porous media through internal boundary conditions. We support the discussion with a set of numerical examples: the convergence properties of the coupled scheme are first investigated; second, we illustrate the practical applicability of the method both for fractured and heterogeneous media.