931 research outputs found
Dynamic model adaptation for multiscale simulation of hyperbolic systems with relaxation
International audienceIn numerous industrial CFD applications, it is usual to use two (or more)different codes to solve a physical phenomenon: where the flow is a priori assumed to have a simple behavior, a code based on a coarse model is applied, while a code based on a fine model is used elsewhere. This leads to a complex coupling problem with fixed interfaces. The aim of the present work is to provide a numerical indicator to optimize to position of these coupling interfaces. In other words, thanks to this numerical indicator, one could verify if the use of the coarser model and of the resulting coupling does not introduce spurious effects. In order to validate this indicator, we use it in a dynamical multiscale method with moving coupling interfaces. The principle of this method is to use as much as possible a coarse model instead of the fine model in the computational domain, in order to obtain an accuracy which is comparable with the one provided by the fine model. We focus here on general hyperbolic systems with stiff relaxation source terms together with the corresponding hyperbolic equilibrium systems. Using a numerical Chapman-Enskog expansion and the distance to the equilibrium manifold, we construct the numerical indicator. Based on several works on the coupling of different hyperbolic models, an original numerical method of dynamic model adaptation is proposed. We prove that this multiscale method preserves invariant domains and that the entropy of the numerical solution decreases with respect to time. The reliability of the adaptation procedure is assessed on various 1D and 2D test cases coming from two-phase flow modeling
Taxis Equations for Amoeboid Cells
The classical macroscopic chemotaxis equations have previously been derived
from an individual-based description of the tactic response of cells that use a
"run-and-tumble" strategy in response to environmental cues. Here we derive
macroscopic equations for the more complex type of behavioral response
characteristic of crawling cells, which detect a signal, extract directional
information from a scalar concentration field, and change their motile behavior
accordingly. We present several models of increasing complexity for which the
derivation of population-level equations is possible, and we show how
experimentally-measured statistics can be obtained from the transport equation
formalism. We also show that amoeboid cells that do not adapt to constant
signals can still aggregate in steady gradients, but not in response to
periodic waves. This is in contrast to the case of cells that use a
"run-and-tumble" strategy, where adaptation is essential.Comment: 35 pages, submitted to the Journal of Mathematical Biolog
A new numerical strategy with space-time adaptivity and error control for multi-scale streamer discharge simulations
This paper presents a new resolution strategy for multi-scale streamer
discharge simulations based on a second order time adaptive integration and
space adaptive multiresolution. A classical fluid model is used to describe
plasma discharges, considering drift-diffusion equations and the computation of
electric field. The proposed numerical method provides a time-space accuracy
control of the solution, and thus, an effective accurate resolution independent
of the fastest physical time scale. An important improvement of the
computational efficiency is achieved whenever the required time steps go beyond
standard stability constraints associated with mesh size or source time scales
for the resolution of the drift-diffusion equations, whereas the stability
constraint related to the dielectric relaxation time scale is respected but
with a second order precision. Numerical illustrations show that the strategy
can be efficiently applied to simulate the propagation of highly nonlinear
ionizing waves as streamer discharges, as well as highly multi-scale nanosecond
repetitively pulsed discharges, describing consistently a broad spectrum of
space and time scales as well as different physical scenarios for consecutive
discharge/post-discharge phases, out of reach of standard non-adaptive methods.Comment: Support of Ecole Centrale Paris is gratefully acknowledged for
several month stay of Z. Bonaventura at Laboratory EM2C as visiting
Professor. Authors express special thanks to Christian Tenaud (LIMSI-CNRS)
for providing the basis of the multiresolution kernel of MR CHORUS, code
developed for compressible Navier-Stokes equations (D\'eclaration d'Invention
DI 03760-01). Accepted for publication; Journal of Computational Physics
(2011) 1-2
Comparison of Subgrid-scale Viscosity Models and Selective Filtering Strategy for Large-eddy Simulations
Explicitly filtered large-eddy simulations (LES), combining high-accuracy schemes with the use of a selective filtering without adding an explicit subgrid-scales (SGS) model, are carried out for the Taylor-Green-vortex and the supersonic-boundary-layer cases. First, the present approach is validated against direct numerical simulation (DNS) results. Subsequently, several SGS models are implemented in order to investigate if they can improve the initial filter-based methodology. It is shown that the most accurate results are obtained when the filtering is used alone as an implicit model, and for a minimal cost. Moreover, the tests for the Taylor-Green vortex indicate that the discretization error from the numerical methods, notably the dissipation error from the high-order filtering, can have a greater influence than the SGS models
Error analysis of a dynamic model adaptation procedure for nonlinear hyperbolic equations
International audienceWe propose a dynamic model adaptation method for a nonlinear conservation law coupled with an ordinary differential equation. This model, called the ''fine model", involves a small time scale and setting this time scale to 0 leads to a classical conservation law, called the ''coarse model", with a flux which depends on the unknown and on space and time. The dynamic model adaptation consists in detecting the regions where the fine model can be replaced by the coarse one in an automatic way, without deteriorating the accuracy of the result. To do so, we provide an error estimate between the solution of the fine model and the solution of the adaptive method, enabling a sharp control of the different parameters. This estimate rests upon stability results for conservation laws with respect to the flux function. Numerical results are presented at the end and show that our estimate is optimal
SOLID-SHELL FINITE ELEMENT MODELS FOR EXPLICIT SIMULATIONS OF CRACK PROPAGATION IN THIN STRUCTURES
Crack propagation in thin shell structures due to cutting is conveniently simulated
using explicit finite element approaches, in view of the high nonlinearity of the problem. Solidshell
elements are usually preferred for the discretization in the presence of complex material
behavior and degradation phenomena such as delamination, since they allow for a correct
representation of the thickness geometry. However, in solid-shell elements the small thickness
leads to a very high maximum eigenfrequency, which imply very small stable time-steps. A new
selective mass scaling technique is proposed to increase the time-step size without affecting
accuracy. New ”directional” cohesive interface elements are used in conjunction with selective
mass scaling to account for the interaction with a sharp blade in cutting processes of thin ductile
shells
OSAMOAL: optimized simulations by adapted models using asymptotic limits
We propose in this work to address the problem of model adaptation, dedicated to hyper- bolic models with relaxation and to their parabolic limit. The goal is to replace a hyperbolic system of balance laws (the so-called fine model) by its parabolic limit (the so-called coarse model), in delimited parts of the computational domain. Our method is based on the construction of asymptotic preserving schemes and on interfacial coupling methods between hyperbolic and parabolic models. We study in parallel the cases of the Goldstein-Taylor model and of the p-system with friction
Fully adaptive multiresolution schemes for strongly degenerate parabolic equations with discontinuous flux
A fully adaptive finite volume multiresolution scheme for one-dimensional
strongly degenerate parabolic equations with discontinuous flux is presented.
The numerical scheme is based on a finite volume discretization using the
Engquist--Osher approximation for the flux and explicit time--stepping. An
adaptivemultiresolution scheme with cell averages is then used to speed up CPU
time and meet memory requirements. A particular feature of our scheme is the
storage of the multiresolution representation of the solution in a dynamic
graded tree, for the sake of data compression and to facilitate navigation.
Applications to traffic flow with driver reaction and a clarifier--thickener
model illustrate the efficiency of this method
An equation-free computational approach for extracting population-level behavior from individual-based models of biological dispersal
The movement of many organisms can be described as a random walk at either or
both the individual and population level. The rules for this random walk are
based on complex biological processes and it may be difficult to develop a
tractable, quantitatively-accurate, individual-level model. However, important
problems in areas ranging from ecology to medicine involve large collections of
individuals, and a further intellectual challenge is to model population-level
behavior based on a detailed individual-level model. Because of the large
number of interacting individuals and because the individual-level model is
complex, classical direct Monte Carlo simulations can be very slow, and often
of little practical use. In this case, an equation-free approach may provide
effective methods for the analysis and simulation of individual-based models.
In this paper we analyze equation-free coarse projective integration. For
analytical purposes, we start with known partial differential equations
describing biological random walks and we study the projective integration of
these equations. In particular, we illustrate how to accelerate explicit
numerical methods for solving these equations. Then we present illustrative
kinetic Monte Carlo simulations of these random walks and show a decrease in
computational time by as much as a factor of a thousand can be obtained by
exploiting the ideas developed by analysis of the closed form PDEs. The
illustrative biological example here is chemotaxis, but it could be any random
walker which biases its movement in response to environmental cues.Comment: 30 pages, submitted to Physica
Multiphase modelling of vascular tumour growth in two spatial dimensions
In this paper we present a continuum mathematical model of vascular tumour growth which is based on a multiphase framework in which the tissue is decomposed into four distinct phases and the principles of conservation of mass and momentum are applied to the normal/healthy cells, tumour cells, blood vessels and extracellular material. The inclusion of a diffusible nutrient, supplied by the blood vessels, allows the vasculature to have a nonlocal influence on the other phases. Two-dimensional computational simulations are carried out on unstructured, triangular meshes to allow a natural treatment of irregular geometries, and the tumour boundary is captured as a diffuse interface on this mesh, thereby obviating the need to explicitly track the (potentially highly irregular and ill-defined) tumour boundary. A hybrid finite volume/finite element algorithm is used to discretise the continuum model: the application of a conservative, upwind, finite volume scheme to the hyperbolic mass balance equations and a finite element scheme with a stable element pair to the generalised Stokes equations derived from momentum balance, leads to a robust algorithm which does not use any form of artificial stabilisation. The use of a matrix-free Newton iteration with a finite element scheme for the nutrient reaction-diffusion equations allows full nonlinearity in the source terms of the mathematical model. Numerical simulations reveal that this four-phase model reproduces the characteristic pattern of tumour growth in which a necrotic core forms behind an expanding rim of well-vascularised proliferating tumour cells. The simulations consistently predict linear tumour growth rates. The dependence of both the speed with which the tumour grows and the irregularity of the invading tumour front on the model parameters are investigated
- …