Lattice-free models of cell invasion : discrete simulations and travelling waves

Invasion waves of cells play an important role in development, disease and repair. Standard discrete models of such processes typically involve simulating cell motility, cell proliferation and cell-to-cell crowding effects in a lattice-based framework. The continuum-limit description is often given by a reaction–diffusion equation that is related to the Fisher–Kolmogorov equation. One of the limitations of a standard lattice-based approach is that real cells move and proliferate in continuous space and are not restricted to a predefined lattice structure. We present a lattice-free model of cell motility and proliferation, with cell-to-cell crowding effects, and we use the model to replicate invasion wave-type behaviour. The continuum-limit description of the discrete model is a reaction–diffusion equation with a proliferation term that is different from lattice-based models. Comparing lattice based and lattice-free simulations indicates that both models lead to invasion fronts that are similar at the leading edge, where the cell density is low. Conversely, the two models make different predictions in the high density region of the domain, well behind the leading edge. We analyse the continuum-limit description of the lattice based and lattice-free models to show that both give rise to invasion wave type solutions that move with the same speed but have very different shapes. We explore the significance of these differences by calibrating the parameters in the standard Fisher–Kolmogorov equation using data from the lattice-free model. We conclude that estimating parameters using this kind of standard procedure can produce misleading results.

Quantifying spatial structure in experimental observations and agent-based simulations using pair-correlation functions

We define a pair-correlation function that can be used to characterize spatiotemporal patterning in experimental images and snapshots from discrete simulations. Unlike previous pair-correlation functions, the pair-correlation functions developed here depend on the location and size of objects. The pair-correlation function can be used to indicate complete spatial randomness, aggregation or segregation over a range of length scales, and quantifies spatial structures such as the shape, size and distribution of clusters. Comparing pair-correlation data for various experimental and simulation images illustrates their potential use as a summary statistic for calibrating discrete models of various physical processes.

How to Facilitate Open Disclosure Simulations

Open Disclosure comprises two main components: Clinician Disclosure (CD), an informal process usually conducted by the treating clinician; and Formal Open Disclosure (FOD), a more structured process led by a senior clinician trained as an Open Disclosure Consultant. Training programs for both CD and FOD incorporate interactive role-play based scenarios called ‘simulations’. This section of the Open Disclosure Training Program Handbook provides guidelines and resources for facilitating the simulation components of both Clinician Disclosure and Open Disclosure Consultant training.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high-potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage. However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two dimensional (2D) numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Numerical simulations of particle deposition in metal foam heat exchangers

Australia is a high potential country for geothermal power with reserves currently estimated in the tens of millions of petajoules, enough to power the nation for at least 1000 years at current usage.However, these resources are mainly located in isolated arid regions where water is scarce. Therefore, wet cooling systems for geothermal plants in Australia are the least attractive solution and thus air-cooled heat exchangers are preferred. In order to increase the efficiency of such heat exchangers, metal foams have been used. One issue raised by this solution is the fouling caused by dust deposition. In this case, the heat transfer characteristics of the metal foam heat exchanger can dramatically deteriorate. Exploring the particle deposition property in the metal foam exchanger becomes crucial. This paper is a numerical investigation aimed to address this issue. Two-dimensional(2D numerical simulations of a standard one-row tube bundle wrapped with metal foam in cross-flow are performed and highlight preferential particle deposition areas.

Simulations of the spontaneous emission of a quantum dot near a gap plasmon waveguide

In this paper we modeled a quantum dot at near proximity to a gap plasmon waveguide to study the quantum dot-plasmon interactions. Assuming that the waveguide is single mode, this paper is concerned about the dependence of spontaneous emission rate of the quantum dot on waveguide dimensions such as width and height. We compare coupling efficiency of a gap waveguide with symmetric configuration and asymmetric configuration illustrating that symmetric waveguide has a better coupling efficiency to the quantum dot. We also demonstrate that optimally placed quantum dot near a symmetric waveguide with 50 nm x 50 nm cross section can capture 80% of the spontaneous emission into a guided plasmon mode.

Using narrow beam profiles to quantify focal spot size, for accurate Monte Carlo simulations of SRS/SRT systems

This study investigates the variation of photon field penumbra shape with initial electron beam diameter, for very narrow beams. A Varian Millenium MLC (Varian Medical Systems, Palo Alto, USA) and a Brainlab m3 microMLC (Brainlab AB. Feldkirchen, Germany) were used, with one Varian iX linear accelerator, to produce fields that were (nominally) 0.20 cm across. Dose profiles for these fields were measured using radiochromic film and compared with the results of simulations completed using BEAMnrc and DOSXYZnrc, where the initial electron beam was set to FWHM = 0.02, 0.10, 0.12, 0.15, 0.20 and 0.50 cm. Increasing the electron-beam FWHM produced increasing occlusion of the photon source by the closely spaced collimator leaves and resulted in blurring of the simulated profile widths from 0.26 to 0.64 cm, for the MLC, from 0.12 to 0.43 cm, for the microMLC. Comparison with measurement data suggested that the electron spot size in the clinical linear accelerator was between FWHM = 0.10 and 0.15 cm, encompassing the result of our previous output-factor based work, which identified a FWHM of 0.12. Investigation of narrow-beam penumbra variation has been found to be a useful procedure, with results varying noticeably with linear accelerator spot size and allowing FWHM estimates obtained using other methods to be verified.

Monte Carlo simulations of small field output factors in solid phantoms

Introduction Total scatter factor (or output factor) in megavoltage photon dosimetry is a measure of relative dose relating a certain field size to a reference field size. The use of solid phantoms has been well established for output factor measurements, however to date these phantoms have not been tested with small fields. In this work, we evaluate the water equivalency of a number of solid phantoms for small field output factor measurements using the EGSnrc Monte Carlo code. Methods The following small square field sizes were simulated using BEAMnrc: 5, 6, 7, 8, 10 and 30 mm. Each simulated phantom geometry was created in DOSXYZnrc and consisted of a silicon diode (of length and width 1.5 mm and depth 0.5 mm) submersed in the phantom at a depth of 5 g/cm2. The source-to-detector distance was 100 cm for all simulations. The dose was scored in a single voxel at the location of the diode. Interaction probabilities and radiation transport parameters for each material were created using custom PEGS4 files. Results A comparison of the resultant output factors in the solid phantoms, compared to the same factors in a water phantom are shown in Fig. 1. The statistical uncertainty in each point was less than or equal to 0.4 %. The results in Fig. 1 show that the density of the phantoms affected the output factor results, with higher density materials (such as PMMA) resulting in higher output factors. Additionally, it was also calculated that scaling the depth for equivalent path length had negligible effect on the output factor results at these field sizes. Discussion and conclusions Electron stopping power and photon mass energy absorption change minimally with small field size [1]. Also, it can be seen from Fig. 1 that the difference from water decreases with increasing field size. Therefore, the most likely cause for the observed discrepancies in output factors is differing electron disequilibrium as a function of phantom density. When measuring small field output factors in a solid phantom, it is important that the density is very close to that of water.

Short-term changes in axial length during simulations of typical far, intermediate and near tasks

Purpose: To investigate the changes in axial length with the combined effect of accommodation and angle of gaze (convergence and downward gaze) over 5 minutes in groups of myopes and emmetropes. Methods: A total of 31 subjects (nine emmetropes, 10 low myopes, and 12 moderate to high myopes) aged from 18 to 31 years were recruited. To measure ocular biometrics in inferonasal gaze with accommodation, an optical biometer (Lenstar LS900) was inclined on a tilt and height adjustable stage, with the subject’s chinrest mounted on a rotary stage to induce various levels of convergence by rotation of the subject’s head in primary or downward gaze. Initially, the subjects performed a distance viewing task in primary gaze for 10 minutes to provide a ‘wash-out’ period for prior visual tasks, and then the subject’s axial length and ocular biometrics were measured in nine different combinations of gaze/accommodation over 5 minutes. These nine sessions for all gaze measurements (i.e. three levels of accommodation 9 three levels of convergence) were completed across 3 days of testing (one accommodation condition on each day).The nine combinations of gaze/accommodation were based on those required to view the centre, right and left edges of a distant TV at 6 m in primary gaze, an intermediate task (i.e. computer at 50 cm in 10° downward gaze) and a near task (i.e. reading A4 page at 20 cm in 20° downward gaze). Subjects were wearing a custom built three-axes head tracker throughout the experiment that monitored subjects’ relative head movements (roll, pitch and yaw) during measurements. Results: A significant increase in axial length occurred with the combined effect of accommodation, convergence and downward gaze (repeated measures ANOVA, p < 0.001), with the greatest axial elongation during the near task in downward gaze with convergence (i.e. downward 20°/inward 33°, with 5 D accommodation) (mean change 33 ± 13 lm, after 5 minutes task) followed by the intermediate task (i.e. downward 10°/inward 25°, with 2 D accommodation) (mean change 14 ± 11 lm, after 5 minutes task).Changes in axial length for the distance task (i.e. primary gaze/9° convergence, with 0.16 D accommodation) were not statistically significant (mean change 4 ± 8 lm, after 5 minutes task, p > 0.05). Moderate to high myopes had a greater change in the axial length (mean change 40 ± 11 lm after 5 minutes of near task) than that of emmetropes (mean change 29 ± 15 lm after 5 minutes of near task) and low myopes (mean change 29 ± 16 lm after 5 minutes of near task) associated with time (p = 0.02) and accommodation by time (p = 0.03). Conclusions: The combination of accommodation, convergence and downward angle has a significant short term effect on axial length over time. The near task in downward gaze with convergence caused a greater change in axial length than the intermediate and distant visual tasks. The greater axial elongation measured in the infero-nasal direction with accommodation is most likely associated with a combination of biomechanical factors such as, extraocular muscle forces and ciliary muscle contraction.

Controlling current and voltage type interfaces in power-hardware-in-the-loop simulations

The usual practice to study a large power system is through digital computer simulation. However, the impact of large scale use of small distributed generators on a power network cannot be evaluated strictly by simulation since many of these components cannot be accurately modelled. Moreover, the network complexity makes the task of practical testing on a physical network nearly impossible. This study discusses the paradigm of interfacing a real-time simulation of a power system to real-life hardware devices. This type of splitting a network into two parts and running a real-time simulation with a physical system in parallel is usually termed as power-hardware-in-the-loop (PHIL) simulation. The hardware part is driven by a voltage source converter that amplifies the signals of the simulator. In this paper, the effects of suitable control strategy on the performance of PHIL and the associated stability aspects are analysed in detail. The analyses are validated through several experimental tests using an real-time digital simulator.

CFD simulations of flow and heat transfer through the porous interface of a metal foam heat exchanger

This paper offers numerical modelling of a waste heat recovery system. A thin layer of metal foam is attached to a cold plate to absorb heat from hot gases leaving the system. The heat transferred from the exhaust gas is then transferred to a cold liquid flowing in a secondary loop. Two different foam PPI (Pores Per Inch) values are examined over a range of fluid velocities. Numerical results are then compared to both experimental data and theoretical results available in the literature. Challenges in getting the simulation results to match those of the experiments are addressed and discussed in detail. In particular, interface boundary conditions specified between a porous layer and a fluid layer are investigated. While physically one expects much lower fluid velocity in the pores compared to that of free flow, capturing this sharp gradient at the interface can add to the difficulties of numerical simulation. The existing models in the literature are modified by considering the pressure gradient inside and outside the foam. Comparisons against the numerical modelling are presented. Finally, based on experimentally-validated numerical results, thermo-hydraulic performance of foam heat exchangers as waste heat recovery units is discussed with the main goal of reducing the excess pressure drop and maximising the amount of heat that can be recovered from the hot gas stream.

3D CFD simulations of a candidate R143a radial-inflow turbine for geothermal power applications

Optimisation of Organic Rankine Cycles (ORCs) for binary cycle applications could play a major role in determining the competitiveness of low to moderate renewable sources. An important aspect of the optimisation is to maximise the turbine output power for a given resource. This requires careful attention to the turbine design notably through numerical simulations. Challenges in the numerical modelling of radial-inflow turbines using high-density working fluids still need to be addressed in order to improve the turbine design and better optimise ORCs. This paper presents preliminary 3D numerical simulations of a radial-inflow turbine working with high-density fluids in realistic geothermal ORCs. Following extensive investigation of the operating conditions and thermodynamic cycle analysis, the refrigerant R143a is chosen as the high-density working fluid. The 1D design of the candidate radial-inflow turbine is presented in details. Furthermore, commercially-available software Ansys-CFX is used to perform the 3D CFD simulations for a number of operating conditions including off-design conditions. The real-gas properties are obtained using the Peng-Robinson equations of state. The preliminary design created using dedicated radial-inflow turbine software Concepts-Rital is discussed and the 3D CFD results are presented and compared against the meanline analysis.

Prediction of thermal expansion properties of carbon nanotubes using molecular dynamics simulations

The axial coefficients of thermal expansion (CTE) of various carbon nanotubes (CNTs), i.e., single-wall carbon nanotubes (SWCNTs), and some multi-wall carbon nanotubes (MWCNTs), were predicted using molecular dynamics (MDs) simulations. The effects of two parameters, i.e., temperature and the CNT diameter, on CTE were investigated extensively. For all SWCNTs and MWCNTs, the obtained results clearly revealed that within a wide low temperature range, their axial CTEs are negative. As the diameter of CNTs decreases, this temperature range for negative axial CTEs becomes narrow, and positive axial CTEs appear in high temperature range. It was found that the axial CTEs vary nonlinearly with the temperature, however, they decrease linearly as the CNT diameter increases. Moreover, within a wide temperature range, a set of empirical formulations was proposed for evaluating the axial CTEs of armchair and zigzag SWCNTs using the above two parameters. Finally, it was found that the absolute value of the negative axial CTE of any MWCNT is much smaller than those of its constituent SWCNTs, and the average value of the CTEs of its constituent SWCNTs. The present fundamental study is very important for understanding the thermal behaviors of CNTs in such as nanocomposite temperature sensors, or nanoelectronics devices using CNTs.

Supramolecular organization in fluorene/indenofluorene-oligothiophene alternating conjugated copolymers

A series of conjugated copolymers containing fluorene or indenofluorene units alternating with oligothiophene segments, with potential interest for use as the active layer in field-effect transistors, is investigated. Atomic force microscopy analysis of the morphology of thin deposits shows either the formation of fibrillar structures, which are the signature of long-range π stacking, or the presence of untextured aggregates, resulting from disordered assembly. These morphologies are interpreted in terms of the supramolecular organization of the conjugated chains. Molecular modeling simulations indicate that the commensurability between the lengths of the monomer units and the presence of alkyl side groups are the two key structural factors governing the chain organization into highly ordered assemblies. The most favorable structures are those combining fluorene (indenofluorene) units with unsubstituted bithiophene (terthiophene) segments.

Bridging experiments, models and simulations : an integrative approach to validation in computational cardiac electrophysiology

Computational models in physiology often integrate functional and structural information from a large range of spatio-temporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and scepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace and refine animal experiments. A fundamental requirement to fulfil these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations between experiments, models and simulations in cardiac electrophysiology. We describe the processes, data and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. Validation must therefore take into account the complex interplay between models, simulations and experiments. Key points for developing strategies for validation are: 1) understanding sources of bio-variability is crucial to the comparison between simulation and experimental results; 2) robustness of techniques and tools is a pre-requisite to conducting physiological investigations using the MSE system; 3) definition and adoption of standards facilitates interoperability of experiments, models and simulations; 4) physiological validation must be understood as an iterative process that defines the specific aspects of electrophysiology the MSE system targets, and is driven by advancements in experimental and computational methods and the combination of both.