24,876 research outputs found
An efficient parallel immersed boundary algorithm using a pseudo-compressible fluid solver
We propose an efficient algorithm for the immersed boundary method on
distributed-memory architectures, with the computational complexity of a
completely explicit method and excellent parallel scaling. The algorithm
utilizes the pseudo-compressibility method recently proposed by Guermond and
Minev [Comptes Rendus Mathematique, 348:581-585, 2010] that uses a directional
splitting strategy to discretize the incompressible Navier-Stokes equations,
thereby reducing the linear systems to a series of one-dimensional tridiagonal
systems. We perform numerical simulations of several fluid-structure
interaction problems in two and three dimensions and study the accuracy and
convergence rates of the proposed algorithm. For these problems, we compare the
proposed algorithm against other second-order projection-based fluid solvers.
Lastly, the strong and weak scaling properties of the proposed algorithm are
investigated
A computational method for the coupled solution of reaction–diffusion equations on evolving domains and manifolds: application to a model of cell migration and chemotaxis
In this paper, we devise a moving mesh finite element method for the approximate solution of coupled bulk–surface reaction–diffusion equations on an evolving two dimensional domain. Fundamental to the success of the method is the robust generation of bulk and surface meshes. For this purpose, we use a novel moving mesh partial differential equation (MMPDE) approach. The developed method is applied to model problems with known analytical solutions; these experiments indicate second-order spatial and temporal accuracy. Coupled bulk–surface problems occur frequently in many areas; in particular, in the modelling of eukaryotic cell migration and chemotaxis. We apply the method to a model of the two-way interaction of a migrating cell in a chemotactic field, where the bulk region corresponds to the extracellular region and the surface to the cell membrane
A Moving Frame Algorithm for High Mach Number Hydrodynamics
We present a new approach to Eulerian computational fluid dynamics that is
designed to work at high Mach numbers encountered in astrophysical hydrodynamic
simulations. The Eulerian fluid conservation equations are solved in an
adaptive frame moving with the fluid where Mach numbers are minimized. The
moving frame approach uses a velocity decomposition technique to define local
kinetic variables while storing the bulk kinetic components in a smoothed
background velocity field that is associated with the grid velocity.
Gravitationally induced accelerations are added to the grid, thereby minimizing
the spurious heating problem encountered in cold gas flows. Separately tracking
local and bulk flow components allows thermodynamic variables to be accurately
calculated in both subsonic and supersonic regions. A main feature of the
algorithm, that is not possible in previous Eulerian implementations, is the
ability to resolve shocks and prevent spurious heating where both the preshock
and postshock Mach numbers are high. The hybrid algorithm combines the high
resolution shock capturing ability of the second-order accurate Eulerian TVD
scheme with a low-diffusion Lagrangian advection scheme. We have implemented a
cosmological code where the hydrodynamic evolution of the baryons is captured
using the moving frame algorithm while the gravitational evolution of the
collisionless dark matter is tracked using a particle-mesh N-body algorithm.
The MACH code is highly suited for simulating the evolution of the IGM where
accurate thermodynamic evolution is needed for studies of the Lyman alpha
forest, the Sunyaev-Zeldovich effect, and the X-ray background. Hydrodynamic
and cosmological tests are described and results presented. The current code is
fast, memory-friendly, and parallelized for shared-memory machines.Comment: 19 pages, 5 figure
Exponential integrators for a Markov chain model of the fast sodium channel of cardiomyocytes
The modern Markov chain models of ionic channels in excitable membranes are
numerically stiff. The popular numerical methods for these models require very
small time steps to ensure stability. Our objective is to formulate and test
two methods addressing this issue, so that the timestep can be chosen based on
accuracy rather than stability.
Both proposed methods extend Rush-Larsen technique, which was originally
developed to Hogdkin-Huxley type gate models. One method, "Matrix Rush-Larsen"
(MRL) uses a matrix reformulation of the Rush-Larsen scheme, where the matrix
exponentials are calculated using precomputed tables of eigenvalues and
eigenvectors. The other, "hybrid operator splitting" (HOS) method exploits
asymptotic properties of a particular Markov chain model, allowing explicit
analytical expressions for the substeps.
We test both methods on the Clancy and Rudy (2002) INa Markov chain model.
With precomputed tables for functions of the transmembrane voltage, both
methods are comparable to the forward Euler method in accuracy and
computational cost, but allow longer time steps without numerical instability.
We conclude that both methods are of practical interest. MRL requires more
computations than HOS, but is formulated in general terms which can be readily
extended to other Markov Chain channel models, whereas the utility of HOS
depends on the asymptotic properties of a particular model.
The significance of the methods is that they allow a considerable speed-up of
large-scale computations of cardiac excitation models by increasing the time
step, while maintaining acceptable accuracy and preserving numerical stability.Comment: 9 pages, 5 figures main text + 14 pages, 1 figure appendix, as
submitted in final form to IEEE TBME 2014/11/11. Copyright IEEE (2014
Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids
In this review, we describe and analyze a mesoscale simulation method for
fluid flow, which was introduced by Malevanets and Kapral in 1999, and is now
called multi-particle collision dynamics (MPC) or stochastic rotation dynamics
(SRD). The method consists of alternating streaming and collision steps in an
ensemble of point particles. The multi-particle collisions are performed by
grouping particles in collision cells, and mass, momentum, and energy are
locally conserved. This simulation technique captures both full hydrodynamic
interactions and thermal fluctuations. The first part of the review begins with
a description of several widely used MPC algorithms and then discusses
important features of the original SRD algorithm and frequently used
variations. Two complementary approaches for deriving the hydrodynamic
equations and evaluating the transport coefficients are reviewed. It is then
shown how MPC algorithms can be generalized to model non-ideal fluids, and
binary mixtures with a consolute point. The importance of angular-momentum
conservation for systems like phase-separated liquids with different
viscosities is discussed. The second part of the review describes a number of
recent applications of MPC algorithms to study colloid and polymer dynamics,
the behavior of vesicles and cells in hydrodynamic flows, and the dynamics of
viscoelastic fluids
Model-based analysis for the thermal management of open-cathode proton exchange membrane fuel cell systems concerning efficiency and stability
In this work we present a dynamic, control-oriented, concentrated parameter model of an open-cathode proton exchange membrane fuel cell system for the study of stability and efficiency improvement with respect to thermal management. The system model consists of two dynamic states which are the fuel cell temperature and the liquid water saturation in the cathode catalyst layer. The control action of the system is the inlet air velocity of the cathode air flow manifold, set by the cooling fan, and the system output is the stack voltage. From the model we derive the equilibrium points and eigenvalues within a set of operating conditions and subsequently discuss stability and the possibility of efficiency improvement. The model confirms the existence of a temperature-dependent maximum power in the moderate temperature region. The stability analysis shows that the maximum power line decomposes the phase plane in two parts, namely stable and unstable equilibrium points. The model is capable of predicting the temperature of a stable steady-state voltage maximum and the simulation results serve for the design of optimal thermal management strategies.Postprint (author's final draft
- …