Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNonlinear simulation of vascular tumor growth with chemotaxis and the control of necrosis
63
0
0.0
(
0
)
Added by
Min-Jhe Lu
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
In this paper, we develop a sharp interface tumor growth model to study the effect of both the intratumoral structure using a controlled necrotic core and the extratumoral nutrient supply from vasculature on tumor morphology. We first show that our model extends the benchmark results in the literature using linear stability analysis. Then we solve this generalized model numerically using a spectrally accurate boundary integral method in an evolving annular domain, not only with a Robin boundary condition on the outer boundary for the nutrient field which models tumor vasculature, but also with a static boundary condition on the inner boundary for pressure field which models the control of tumor necrosis. The discretized linear systems for both pressure and nutrient fields are shown to be well-conditioned through tracing GMRES iteration numbers. Our nonlinear simulations reveal the stabilizing effects of angiogenesis and the destabilizing ones of chemotaxis and necrosis in the development of tumor morphological instabilities if the necrotic core is fixed in a circular shape. When the necrotic core is controlled in a non-circular shape, the stabilizing effects of proliferation and the destabilizing ones of apoptosis are observed. Finally, the values of the nutrient concentration with its fluxes and the pressure level with its normal derivatives, which are solved accurately at the boundaries, help us to characterize the corresponding tumor morphology and the level of the biophysical quantities on interfaces required in keeping various shapes of the necrotic region of the tumor. Interestingly, we notice that when the necrotic region is fixed in a 3-fold non-circular shape, even if the initial shape of the tumor is circular, the tumor will evolve into a shape corresponding to the 3-fold symmetry of the shape of the fixed necrotic region.
rate research
Read More
In this paper, we develop a sharp interface tumor growth model to study the effect of the tumor microenvironment using a complex far-field geometry that mimics a heterogeneous distribution of vasculature. Together with different nutrient uptake rates inside and outside the tumor, this introduces variability in spatial diffusion gradients. Linear stability analysis suggests that the uptake rate in the tumor microenvironment, together with chemotaxis, may induce unstable growth, especially when the nutrient gradients are large. We investigate the fully nonlinear dynamics using a spectrally accurate boundary integral method. Our nonlinear simulations reveal that vascular heterogeneity plays an important role in the development of morphological instabilities that range from fingering and chain-like morphologies to compact, plate-like shapes in two-dimensions.
Herewith we discuss a network model of the ferroptosis avascular and vascular tumor growth based on our previous proposed framework. Chiefly, ferroptosis should be viewed as a first order phase transition characterized by a supercritical Andronov Hopf bifurcation, with the emergence of limit cycle. The increase of the population of the oxidized PUFA fragments, take as the control parameter, involves an inverse Feigenbaum, (a cascade of saddle foci Shilnikovs bifurcations) scenario, which results in the stabilization of the dynamics and in a decrease of complexity.
The present work was inspired by the recent developments in laboratory experiments made on chip, where culturing of multiple cell species was possible. The model is based on coupled reaction-diffusion-transport equations with chemotaxis, and takes into account the interactions among cell populations and the possibility of drug administration for drug testing effects. Our effort was devoted to the development of a simulation tool that is able to reproduce the chemotactic movement and the interactions between different cell species (immune and cancer cells) living in microfluidic chip environment. The main issues faced in this work are the introduction of mass-preserving and positivity-preserving conditions involving the balancing of incoming and outgoing fluxes passing through interfaces between 2D and 1D domains of the chip and the development of mass-preserving and positivity preserving numerical conditions at the external boundaries and at the interfaces between 2D and 1D domains.
In this article, global stabilization results for the two dimensional (2D) viscous Burgers equation, that is, convergence of unsteady solution to its constant steady state solution with any initial data, are established using a nonlinear Neumann boundary feedback control law. Then, applying $C^0$-conforming finite element method in spatial direction, optimal error estimates in $L^infty(L^2)$ and in $L^infty(H^1)$- norms for the state variable and convergence result for the boundary feedback control law are derived. All the results preserve exponential stabilization property. Finally, several numerical experiments are conducted to confirm our theoretical findings.
We introduce a simple, rigorous, and unified framework for solving nonlinear partial differential equations (PDEs), and for solving inverse problems (IPs) involving the identification of parameters in PDEs, using the framework of Gaussian processes. The proposed approach: (1) provides a natural generalization of collocation kernel methods to nonlinear PDEs and IPs; (2) has guaranteed convergence for a very general class of PDEs, and comes equipped with a path to compute error bounds for specific PDE approximations; (3) inherits the state-of-the-art computational complexity of linear solvers for dense kernel matrices. The main idea of our method is to approximate the solution of a given PDE as the maximum a posteriori (MAP) estimator of a Gaussian process conditioned on solving the PDE at a finite number of collocation points. Although this optimization problem is infinite-dimensional, it can be reduced to a finite-dimensional one by introducing additional variables corresponding to the values of the derivatives of the solution at collocation points; this generalizes the representer theorem arising in Gaussian process regression. The reduced optimization problem has the form of a quadratic objective function subject to nonlinear constraints; it is solved with a variant of the Gauss--Newton method. The resulting algorithm (a) can be interpreted as solving successive linearizations of the nonlinear PDE, and (b) in practice is found to converge in a small number of iterations (2 to 10), for a wide range of PDEs. Most traditional approaches to IPs interleave parameter updates with numerical solution of the PDE; our algorithm solves for both parameter and PDE solution simultaneously. Experiments on nonlinear elliptic PDEs, Burgers equation, a regularized Eikonal equation, and an IP for permeability identification in Darcy flow illustrate the efficacy and scope of our framework.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
