To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
In this work, we address the numerical solution of fluid-structure interaction problems. This issue is particularly difficulty to tackle when the fluid and the solid densities are of the same order, for instance as it happens in hemodynamic applications, since fully implicit coupling schemes are required to ensure stability of the resulting method. Thus, at each time step, we have to solve a highly non-linear coupled system, since the fluid domain depends on the unknown displacement of the structure. Standard strategies for solving this non-linear problems, are fixed point based methods such as Block-Gauss-Seidel (BGS) iterations. Unfortunately, these methods are very CPU time consuming and usually show slow convergence. We propose a modified fixed-point algorithm which combines the standard BGS iterations with a transpiration formulation. Numerical experiments show the great improvement in computing time with respect to the standard BGS method.
We consider a simple model for the immune systemin which virus are able to undergo mutations and are in competitionwith leukocytes. These mutations are related to several other concepts which havebeen proposed in the literature like those of shape or ofvirulence – a continuous notion. For a given species, the system admits aglobally attractive critical point. We prove that mutations do not affect thispicture for small perturbations and under strong structural assumptions.Based on numerical and theoretical arguments, we also examine how, releasing these assumptions, the system can blow-up.
We propose a quasi-Newton algorithm for solving fluid-structure interaction problems. The basic idea of the method is to build an approximate tangent operator which is cost effective and which takes into account the so-called added mass effect. Various test cases show that the method allows a significant reduction of the computational effort compared to relaxed fixed point algorithms. We present 2D and 3D fluid-structure simulations performed either with a simple 1D structure model or with shells in large displacements.
Our purpose is to estimate numerically the influence of particles on the global viscosity of fluid–particle mixtures. Particles are supposed to rigid, and the surrounding fluid is newtonian. The motion of the mixture is computed directly, i.e. all the particle motions are computed explicitly. Apparent viscosity, based on the force exerted by the fluid on the sliding walls, is computed at each time step of the simulation.In order to perform long–time simulations and still control the solid fraction, we assume periodicity of the flow in the shear direction.Direct simulations are based on the so–called Arbitrary Lagrangian Eulerian approach, which we adapted to make it suitable to periodic domains.As a first step toward modelling of interacting red cells in the blood, we propose a simple model of circular particles submitted to an attractive force which tends to form aggregates.
In this paper we outline the hyperbolic system of governing equationsdescribing one-dimensional blood flow in arterial networks. Thissystem is numerically discretised using a discontinuous Galerkinformulation with a spectral/hp element spatial approximation. Weapply the numerical model to arterial networks in theplacenta. Starting with a single placenta we investigate the velocity waveformin the umbilical artery and its relationship with the distalbifurcation geometry and the terminal resistance. We then present resultsfor the waveform patterns and the volume fluxes throughout a simplifiedmodel of the arterial placental network in a monochorionic twin pregnancy with an arterio-arterial anastomosis and an arterio-venous anastomosis. Theeffects of varying the time period of the two fetus' heart beats, increasing the input flux of one fetus and the role ofterminal resistance in the network are investigated and discussed. The results show that the main features of the in vivo, physiological waves are captured by thecomputational model and demonstrate the applicability of themethods to the simulation of flows in arterial networks.
In Parts I and II of this book, we studied in some detail the Richardson extrapolation and its generalizations and various important sequence transformations. We also mentioned several applications of them. Actually, we discussed in detail the Romberg integration of infinite-range integrals of regular integrands, numerical differentiation, and the computation of infinite-range integrals by the D-transformation. We discussed the application of the various generalizations of the D-transformation to the computation of oscillatory infinite-range integrals, including some important integral transforms. We also treated in detail the acceleration of convergence of infinite series, including power series and Fourier series and their generalizations, by the d-transformation and other methods, such as the Shanks transformation, the θ-algorithm, the Baker-Gammel approximants and their extensions, and so on. In connection with acceleration of convergence of power series, we also discussed in some detail the subject of prediction via the d-transformation and mentioned that the approach presented can be used with any sequence transformation. In this chapter, we add further applications of special interest.
We would like to note that extensive surveys and bibliographies covering the application of extrapolation methods to numerical integration can be found in Joyce [145], Davis and Rabinowitz [63], and Rabinowitz [234].
In Section 4.4, we gave a brief convergence study of GREP(m) for both Process I and Process II. In this study, we treated the cases in which GREP(m) was stable. In addition, we made some practical remarks on stability of GREP(m) in Section 4.5. The aim of the study was to justify the preference given to Process I and Process II as the relevant limiting processes to be used for approximating A, the limit or antilimit of A(y) as y → 0+. We also mentioned that stability was not necessary for convergence and that convergence could be proved at least in some cases in which the extrapolation process is clearly unstable.
In this chapter as well as the next, we would like to make more refined statements about the convergence and stability properties of GREP(1), the simplest form and prototype of GREP, as it is being applied to functions A(y) ∈ F(1).
Before going on, we mention that this chapter is an almost exact reproduction of the recent paper Sidi [306].
As we will be using the notation and results of Section 7.2 on the W-algorithm, we believe a review of this material is advisable at this point.