BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Maria Lukacova-Medvidova appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We propose a finite volume stochastic collocation method for the random Euler system. We rigorously prove the convergence of random finite volume solutions under the assumption that the discrete differential quotients remain bounded in probability. Convergence analysis combines results on the convergence of a deterministic FV method with stochastic compactness arguments due to Skorokhod and Gyöngy-Krylov.
In this paper we study the convergence of a second order finite volume approximation of the scalar conservation law. This scheme is based on the generalized Riemann problem (GRP) solver. We firstly investigate the stability of the GRP scheme and find that it might be entropy unstable when the shock wave is generated. By adding an artificial viscosity we propose a new stabilized GRP scheme. Under the assumption that numerical solutions are uniformly bounded, we prove consistency and convergence of this new GRP method.
We prove the existence of weak solutions to a viscoelastic phase separation problem in two space dimensions. The mathematical model consists of a Cahn-Hilliard-type equation for two-phase flows and the Peterlin-Navier-Stokes equations for viscoelastic fluids. We focus on the case of a polynomial-like potential and suitably bounded coefficient functions. Using the Lagrange-Galerkin finite element method complex behavior of solution for spinodal decomposition including transient polymeric network structures is demonstrated.
The aim of this paper is to prove global in time existence of weak solutions for a viscoelastic phase separation. We consider the case with singular potentials and degenerate mobilities. Our model couples the diffusive interface model with the Peterlin-Navier-Stokes equations for viscoelastic fluids. To obtain the global in time existence of weak solutions we consider appropriate approximations by solutions of the viscoelastic phase separation with a regular potential and build on the corresponding energy and entropy estimates.
The aim of this paper is to analyze a viscoelastic phase separation model. We derive a suitable notion of the relative energy taking into account the non-convex nature of the energy law for the viscoelastic phase separation. This allows us to prove the weak-strong uniqueness principle. We will provide the estimates for the full model in two space dimensions. For a reduced model we present the estimates in three space dimensions and derive conditional relative energy estimates.
A new model for viscoelastic phase separation is proposed, based on a systematically derived conservative two-fluid model. Dissipative effects are included by phenomenological viscoelastic terms. By construction, the model is consistent with the second law of thermodynamics, and we study well-posedness of the model, i.e., existence of weak solutions, a weak-strong uniqueness principle, and stability with respect to perturbations, which are proven by means of relative energy estimates. A good qualitative agreement with mesoscopic simulations is observed in numerical tests.
In this paper we analyze the three-dimensional Peterlin viscoelastic model. By means of a mixed Galerkin and semigroup approach we prove the existence of a weak solutions. Further combining parabolic regularity with the relative energy method we derive a conditional weak-strong uniqueness result.
Based on relative energy estimates, we study the stability of solutions to the Cahn-Hilliard equation with concentration dependent mobility with respect to perturbations. As a by-product of our analysis, we obtain a weak-strong uniqueness principle on the continuous level under realistic regularity assumptions on strong solutions. We then show that the stability estimates can be further inherited almost verbatim by appropriate Galerkin approximations in space and time. This allows us to derive sharp bounds for the discretization error in terms of certain projection errors and to establish order-optimal a-priori error estimates for semi- and fully discrete approximation schemes.
We propose a multiscale model of the invasion of the extracellular matrix by two types of cancer cells, the differentiated cancer cells and the cancer stem cells. We assume that the epithelial mesenchymal-like transition between them is driven primarily by the epidermal growth factors. We moreover take into account the transidifferentiation program of the cancer stem cells and the cancer associated fibroblast cells as well as the fibroblast-driven remodelling of the extracellular matrix. The proposed haptotaxis model combines the macroscopic phenomenon of the invasion of the extracellular matrix with the microscopic dynamics of the epidermal growth factors. We analyse our model in a component-wise manner and compare our findings with the literature. We investigate pathological situations regarding the epidermal growth factors and accordingly propose "mathematical-treatment" scenarios to control the aggressiveness of the tumour.
We study convergence of a mixed finite element finite volume numerical scheme for the isentropic Navier-Stokes system under the full range of the adiabatic exponent. We establish suitable stability and consistency estimates and show that the Young measure generated by numerical solutions represents a dissipative measure-valued solutions of the limit system. In particular, using the recently established weak{strong uniqueness principle in the class of dissipative measure-valued solutions we show that the numerical solutions converge strongly to a strong solutions of the limit system as long as the latter exists.
We present a problem-suited numerical method for a particularly challenging cancer invasion model. This model is a multiscale haptotaxis advection-reaction-diffusion system that describes the macroscopic dynamics of two types of cancer cells coupled with microscopic dynamics of the cells adhesion on the extracellular matrix. The difficulties to overcome arises from the non-constant advection and diffusion coefficients, a time delay term, as well as stiff reaction terms. Our numerical method is a second order finite volume implicit-explicit scheme adjusted to include a) non-constant diffusion coefficients in the implicit part, b) an interpolation technique for the time delay, and c) a restriction on the time increment for the stiff reaction terms.
In the present work we investigate the chemotactically and proteolytically driven tissue invasion by cancer cells. The model employed is a system of advection-reaction-diffusion equations that features the role of the serine protease urokinase-type plasminogen activator. The analytical and numerical study of this system constitutes a challenge due to the merging, emerging, and travelling concentrations that the solutions exhibit. Classical numerical methods applied to this system necessitate very fine discretization grids to resolve these dynamics in an accurate way. To reduce the computational cost without sacrificing the accuracy of the solution, we apply adaptive mesh refinement techniques, in particular h-refinement. Extended numerical experiments exhibit that this approach provides with a higher order, stable, and robust numerical method for this system. We elaborate on several mesh refinement criteria and compare the results with the ones in the literature. We prove, for a simpler version of this model, $L^\infty$ bounds for the solutions, we study the stability of its conditional steady states, and conclude that it can serve as a test case for further development of mesh refinement techniques for cancer invasion simulations.
Non-uniform grids and mesh adaptation have been a growing part of numerical simulation over the past years. It has been experimentally noted that mesh adaptation leads not only to locally improved solution but also to numerical stability of the underlying method. There have been though only few results on the mathematical analysis of these schemes due to the lack of proper tools that incorporate both the time evolution and the mesh adaptation step of the overall algorithm. In this paper we provide a method to perform the analysis of the mesh adaptation method, including both the mesh reconstruction and evolution of the solution. We moreover employ this method to extract sufficient conditions -on the adaptation of the mesh- that stabilize a numerical scheme in the sense of the entropy dissipation.