20,948 research outputs found
A comparative study fourth order runge kutta-tvd Scheme and fluent software case of inlet flow problems
Inlet as part of aircraft engine plays important role in controlling the rate of airflow
entering to the engine. The shape of inlet has to be designed in such away to make the
rate of airflow does not change too much with angle of attack and also not much
pressure losses at the time, the airflow entering to the compressor section. It is therefore
understanding on the flow pattern inside the inlet is important. The present work
presents on the use of the Fourth Order Runge Kutta – Harten Yee TVD scheme
for
the flow analysis inside inlet. The flow is assumed as an inviscid quasi two dimensional
compressible flow. As an initial stage of computer code development, here uses three
generic inlet models. The first inlet model to allow the problem in hand solved as the
case of inlet with expansion wave case. The second inlet model will relate to the case of
expansion compression wave. The last inlet model concerns with the inlet which
produce series of weak shock wave and end up with a normal shock wave. The
comparison result for the same test case with Fluent Software
[1, 2]
indicates that the
developed computer code based on the Fourth Order Runge Kutta – Harten – Yee TVD
scheme are very close to each other. However for complex inlet geometry, the problem
is in the way how to provide an appropriate mesh model
Cancer modelling: Getting to the heart of the problem
Paradoxically, improvements in healthcare that have enhanced the life expectancy of humans in the Western world have, indirectly, increased the prevalence of certain types of cancer such as prostate and breast. It remains unclear whether this phenomenon should be attributed to the ageing process itself or the cumulative effect of prolonged exposure to harmful environmental stimuli such as ultraviolet light, radiation and carcinogens (Franks and Teich, 1988). Equally, there is also compelling evidence that certain genetic abnormalities can predispose individuals to specific cancers (Ilyas et al., 1999). The variety of factors that have been implicated in the development of solid tumours stems, to a large extent, from the fact that ‘cancer’ is a generic term, often used to characterize a series of disorders that share common features. At this generic level of description, cancer may be viewed as a cellular disease in which controls that usually regulate growth and maintain homeostasis are disrupted. Cancer is typically initiated by genetic mutations that lead to enhanced mitosis of a cell lineage and the formation of an avascular tumour. Since it receives nutrients by diffusion from the surrounding tissue, the size of an avascular tumour is limited to several millimeters in diameter. Further growth relies on the tumour acquiring the ability to stimulate the ingrowth of a new, circulating blood supply from the host vasculature via a process termed angiogenesis (Folkman, 1974). Once vascularised, the tumour has access to a vast nutrient source and rapid growth ensues. Further, tumour fragments that break away from the primary tumour, on entering the vasculature, may be transported to other organs in which they may establish secondary tumours or metastases that further compromise the host. Invasion is another key feature of solid tumours whereby contact with the tissue stimulates the production of enzymes that digest the tissue, liberating space into which the tumour cells migrate. Thus, cancer is a complex, multiscale process. The spatial scales of interest range from the subcellular level, to the cellular and macroscopic (or tissue) levels while the timescales may vary from seconds (or less) for signal transduction pathways to months for tumour doubling times The variety of phenomena involved, the range of spatial and temporal scales over which they act and the complex way in which they are inter-related mean that the development of realistic theoretical models of solid tumour growth is extremely challenging. While there is now a large literature focused on modelling solid tumour growth (for a review, see, for example, Preziosi, 2003), existing models typically focus on a single spatial scale and, as a result, are unable to address the fundamental problem of how phenomena at different scales are coupled or to combine, in a systematic manner, data from the various scales. In this article, a theoretical framework will be presented that is capable of integrating a hierarchy of processes occurring at different scales into a detailed model of solid tumour growth (Alarcon et al., 2004). The model is formulated as a hybrid cellular automaton and contains interlinked elements that describe processes at each spatial scale: progress through the cell cycle and the production of proteins that stimulate angiogenesis are accounted for at the subcellular level; cell-cell interactions are treated at the cellular level; and, at the tissue scale, attention focuses on the vascular network whose structure adapts in response to blood flow and angiogenic factors produced at the subcellular level. Further coupling between the different spatial scales arises from the transport of blood-borne oxygen into the tissue and its uptake at the cellular level. Model simulations will be presented to illustrate the effect that spatial heterogeneity induced by blood flow through the vascular network has on the tumour’s growth dynamics and explain how the model may be used to compare the efficacy of different anti-cancer treatment protocols
Simulation of networks of spiking neurons: A review of tools and strategies
We review different aspects of the simulation of spiking neural networks. We
start by reviewing the different types of simulation strategies and algorithms
that are currently implemented. We next review the precision of those
simulation strategies, in particular in cases where plasticity depends on the
exact timing of the spikes. We overview different simulators and simulation
environments presently available (restricted to those freely available, open
source and documented). For each simulation tool, its advantages and pitfalls
are reviewed, with an aim to allow the reader to identify which simulator is
appropriate for a given task. Finally, we provide a series of benchmark
simulations of different types of networks of spiking neurons, including
Hodgkin-Huxley type, integrate-and-fire models, interacting with current-based
or conductance-based synapses, using clock-driven or event-driven integration
strategies. The same set of models are implemented on the different simulators,
and the codes are made available. The ultimate goal of this review is to
provide a resource to facilitate identifying the appropriate integration
strategy and simulation tool to use for a given modeling problem related to
spiking neural networks.Comment: 49 pages, 24 figures, 1 table; review article, Journal of
Computational Neuroscience, in press (2007
Mean Field description of and propagation of chaos in recurrent multipopulation networks of Hodgkin-Huxley and Fitzhugh-Nagumo neurons
We derive the mean-field equations arising as the limit of a network of
interacting spiking neurons, as the number of neurons goes to infinity. The
neurons belong to a fixed number of populations and are represented either by
the Hodgkin-Huxley model or by one of its simplified version, the
Fitzhugh-Nagumo model. The synapses between neurons are either electrical or
chemical. The network is assumed to be fully connected. The maximum
conductances vary randomly. Under the condition that all neurons initial
conditions are drawn independently from the same law that depends only on the
population they belong to, we prove that a propagation of chaos phenomenon
takes places, namely that in the mean-field limit, any finite number of neurons
become independent and, within each population, have the same probability
distribution. This probability distribution is solution of a set of implicit
equations, either nonlinear stochastic differential equations resembling the
McKean-Vlasov equations, or non-local partial differential equations resembling
the McKean-Vlasov-Fokker- Planck equations. We prove the well-posedness of
these equations, i.e. the existence and uniqueness of a solution. We also show
the results of some preliminary numerical experiments that indicate that the
mean-field equations are a good representation of the mean activity of a finite
size network, even for modest sizes. These experiment also indicate that the
McKean-Vlasov-Fokker- Planck equations may be a good way to understand the
mean-field dynamics through, e.g., a bifurcation analysis.Comment: 55 pages, 9 figure
The what and where of adding channel noise to the Hodgkin-Huxley equations
One of the most celebrated successes in computational biology is the
Hodgkin-Huxley framework for modeling electrically active cells. This
framework, expressed through a set of differential equations, synthesizes the
impact of ionic currents on a cell's voltage -- and the highly nonlinear impact
of that voltage back on the currents themselves -- into the rapid push and pull
of the action potential. Latter studies confirmed that these cellular dynamics
are orchestrated by individual ion channels, whose conformational changes
regulate the conductance of each ionic current. Thus, kinetic equations
familiar from physical chemistry are the natural setting for describing
conductances; for small-to-moderate numbers of channels, these will predict
fluctuations in conductances and stochasticity in the resulting action
potentials. At first glance, the kinetic equations provide a far more complex
(and higher-dimensional) description than the original Hodgkin-Huxley
equations. This has prompted more than a decade of efforts to capture channel
fluctuations with noise terms added to the Hodgkin-Huxley equations. Many of
these approaches, while intuitively appealing, produce quantitative errors when
compared to kinetic equations; others, as only very recently demonstrated, are
both accurate and relatively simple. We review what works, what doesn't, and
why, seeking to build a bridge to well-established results for the
deterministic Hodgkin-Huxley equations. As such, we hope that this review will
speed emerging studies of how channel noise modulates electrophysiological
dynamics and function. We supply user-friendly Matlab simulation code of these
stochastic versions of the Hodgkin-Huxley equations on the ModelDB website
(accession number 138950) and
http://www.amath.washington.edu/~etsb/tutorials.html.Comment: 14 pages, 3 figures, review articl
- …