A semi-alternating direction method for a 2-D fractional FitzHugh–Nagumo monodomain model on an approximate irregular domain

A FitzHugh-Nagumo monodomain model has been used to describe the propagation of the electrical potential in heterogeneous cardiac tissue. In this paper, we consider a two-dimensional fractional FitzHugh-Nagumo monodomain model on an irregular domain. The model consists of a coupled Riesz space fractional nonlinear reaction-diffusion model and an ordinary differential equation, describing the ionic fluxes as a function of the membrane potential. Secondly, we use a decoupling technique and focus on solving the Riesz space fractional nonlinear reaction-diffusion model. A novel spatially second-order accurate semi-implicit alternating direction method (SIADM) for this model on an approximate irregular domain is proposed. Thirdly, stability and convergence of the SIADM are proved. Finally, some numerical examples are given to support our theoretical analysis and these numerical techniques are employed to simulate a two-dimensional fractional Fitzhugh-Nagumo model on both an approximate circular and an approximate irregular domain.

Leaving family behind: Understanding the irregular migration of unaccompanied asylum-seeking minors - A qualitative study

This qualitative study investigated the drivers and determinants of irregular maritime migration among 17 protection visa holders who arrived in Australia as unaccompanied asylum-seeking minors. Semi-structured interviews were also conducted with eight non-government service providers working with unaccompanied minors in the Greater Brisbane area.

Transition to an unsteady flow induced by a fin on the sidewall of a differentially heated air-filled square cavity and heat transfer

The transition from a steady to an unsteady flow induced by an adiabatic fin on the sidewall of a differentially heated air-filled cavity is numerically investigated. Numerical simulations have been performed over the range of Rayleigh numbers from Ra = 105–109. The temporal development and spatial structures of natural convection flows in the cavity with a fin are described. It has been demonstrated that the fin may induce the transition to an unsteady flow and the critical Rayleigh number for the occurrence of the transition is between 3.72 × 106 and 3.73 × 106. Furthermore, the peak frequencies of the oscillations triggered by different mechanisms are obtained through spectral analysis. It has been found that the flow rate through the cavity with a fin is larger than that without a fin under the unsteady flow, indicating that the fin may improve the unsteady flow in the cavity.

Double diffusive Marangoni convection flow of electrically conducting fluid in a square cavity with chemical reaction

Double diffusive Marangoni convection flow of viscous incompressible electrically conducting fluid in a square cavity is studied in this paper by taking into consideration of the effect of applied magnetic field in arbitrary direction and the chemical reaction. The governing equations are solved numerically by using alternate direct implicit (ADI) method together with the successive over relaxation (SOR) technique. The flow pattern with the effect of governing parameters, namely the buoyancy ratio W, diffusocapillary ratio w, and the Hartmann number Ha, is investigated. It is revealed from the numerical simulations that the average Nusselt number decreases; whereas the average Sherwood number increases as the orientation of magnetic field is shifted from horizontal to vertical. Moreover, the effect of buoyancy due to species concentration on the flow is stronger than the one due to thermal buoyancy. The increase in diffusocapillary parameter, w caus

Comparison of oropharyngeal and oral cavity squamous cell cancer incidence and trends in New Zealand and Queensland, Australia

Background Increases in the incidence of squamous cell oropharyngeal cancer (OPC) have been reported from some countries, but have not been assessed in Australia or New Zealand. This study examines trends for squamous cell OPC and squamous cell oral cavity cancer (OCC) in two similarly sized populations, New Zealand and Queensland, Australia. Methods Incidence data for 1982–2010 were obtained from the respective population-based cancer registries for squamous cell OPC and OCC, by subsite, sex, and age. Time trends and annual percentage changes (APCs) were assessed by joinpoint regression. Results The incidence rates of squamous cell OPC in males in New Zealand since 2005 and Queensland since 2006 have increased rapidly, with APCs of 11.9% and 10.6% respectively. The trends were greatest at ages 50–69 and followed more gradual increases previously. In females, rates increased by 2.1% per year in New Zealand from 1982, but by only 0.9% (not significant) in Queensland. In contrast, incidence rates for OCC decreased by 1.2% per year in males in Queensland since 1982, but remained stable for females in Queensland and for both sexes in New Zealand. Overall, incidence rates for both OCC and OPC were substantially higher in Queensland than in New Zealand. In males in both areas, OPC incidence is now higher than that of OCC. Conclusions Incidence rates of squamous cell OPC have increased rapidly in men, while rates of OCC have been stable or reducing, showing distinct etiologies. This has both clinical and public health importance, including implications for the extension of human papilloma virus (HPV) vaccination to males.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

The dependence of computed tomography number to relative electron density conversion on phantom geometry and its impact on planned dose

A computed tomography number to relative electron density (CT-RED) calibration is performed when commissioning a radiotherapy CT scanner by imaging a calibration phantom with inserts of specified RED and recording the CT number displayed. In this work, CT-RED calibrations were generated using several commercially available phantoms to observe the effect of phantom geometry on conversion to electron density and, ultimately, the dose calculation in a treatment planning system. Using an anthropomorphic phantom as a gold standard, the CT number of a material was found to depend strongly on the amount and type of scattering material surrounding the volume of interest, with the largest variation observed for the highest density material tested, cortical bone. Cortical bone gave a maximum CT number difference of 1,110 when a cylindrical insert of diameter 28 mm scanned free in air was compared to that in the form of a 30 × 30 cm2 slab. The effect of using each CT-RED calibration on planned dose to a patient was quantified using a commercially available treatment planning system. When all calibrations were compared to the anthropomorphic calibration, the largest percentage dose difference was 4.2 % which occurred when the CT-RED calibration curve was acquired with heterogeneity inserts removed from the phantom and scanned free in air. The maximum dose difference observed between two dedicated CT-RED phantoms was ±2.1 %. A phantom that is to be used for CT-RED calibrations must have sufficient water equivalent scattering material surrounding the heterogeneous objects that are to be used for calibration.

Void Geometry Driven Spin Crossover in Zeolite-Encapsulated Cobalt Tris(bipyridyl) Complex Ion

The cobalt(II) tris(bipyridyl) complex ion encapsulated in zeolite-Y supercages exhibits a thermally driven interconversion between a low-spin and a high-spin state-a phenomenon not observed for this ion either in solid state or in solution. From a comparative study of the magnetism and optical spectroscopy of the encapsulated and unencapsulated complex ion, supported by molecular modeling, such spin behavior is shown to be intramolecular in origin. In the unencapsulated or free state, the [Co(bipy)(3)](2+) ion exhibits a marked trigonal prismatic distortion, but on encapsulation, the topology of the supercage forces it to adopt a near-octahedral geometry. An analysis using the angular overlap ligand field model with spectroscopically derived parameters shows that the geometry does indeed give rise to a low-spin ground state, and suggests a possible scenario for the spin state interconversion.

The equilibrium geometry and electronic structure of Bi nanolines on clean and hydrogenated Si(001) surfaces

The equilibrium geometry, electronic structure and energetic stability of Bi nanolines on clean and hydrogenated Si(001) surfaces have been examined by means of ab initio total energy calculations and scanning tunnelling microscopy. For the Bi nanolines on a clean Si surface the two most plausible structural models, the Miki or M model (Miki et al 1999 Phys. Rev. B 59 14868) and the Haiku or H model (Owen et al 2002 Phys. Rev. Lett. 88 226104), have been examined in detail. The results of the total energy calculations support the stability of the H model over the M model, in agreement with previous theoretical results. For Bi nanolines on the hydrogenated Si(001) surface, we find that an atomic configuration derived from the H model is also more stable than an atomic configuration derived from the M model. However, the energetically less stable (M) model exhibits better agreement with experimental measurements for equilibrium geometry. The electronic structures of the H and M models are very similar. Both models exhibit a semiconducting character, with the highest occupied Bi-derived bands lying at ~0.5 eV below the valence band maximum. Simulated and experimental STM images confirm that at a low negative bias the Bi lines exhibit an 'antiwire' property for both structural models.

The geometry of Bi nanolines on Si(0 0 1)

A study of the Bi nanoline geometry on Si(0 0 1) has been performed using a combination of ab initio theoretical technique and scanning tunnelling microscopy (STM). Our calculations demonstrate decisively that the recently proposed Haiku geometry is a lower energy configuration than any of the previously proposed line geometries. Furthermore, we have made comparisons between STM constant-current topographs of the lines and Tersoff–Haman STM simulations. Although the Haiku and the Miki geometries both reproduce the main features of the constant-current topographs, the simulated STM images of the Miki geometry have a dark stripe between the dimer rows that does not correspond well with experiment.

Collision-Geometry-Based Pulsed Guidance Law for Exoatmospheric Interception

No abstract is available.

Ultrafast Solvation Dynamics of an Ion in the gamma-Cyclodextrin Cavity: The Role of Restricted Environment

A detailed study of the solvation dynamics of a charged coumarin dye molecule in gamma-cyclodextrin/water has been carried out by using two different theoretical approaches. The first approach is based on a multishell continuum model (MSCM). This model predicts the time scales of the dynamics rather well, provided an accurate description of the frequency-dependent dielectric function is supplied. The reason for this rather surprising agreement is 2-fold. First, there is a cancellation of errors, second, the two-zone model mimics the heterogeneous microenvironment surrounding the ion rather well. The second approach is based on the molecular hydrodynamics theory (MI-IT). In this molecular approach, the solvation dynamics has been studied by restricting the translational motion of the solvent molecules enclosed within the cavity. The results from the molecular theory are also in good agreement with the experimental results. Our study indicates that, in the present case, the restricted environment affects only the long time decay of the solvation time correlation function. The short time dynamics is still governed by the librational (and/or vibrational) modes present in bulk water.

Unsteady Natural Convection Flow in a Square Cavity Filled with a Porous Medium Due to Impulsive Change in Wall Temperature

Unsteady natural convection flow in a two- dimensional square cavity filled with a porous material has been studied. The flow is initially steady where the left- hand vertical wall has temperature T-h and the right- hand vertical wall is maintained at temperature T-c ( T-h > T-c) and the horizontal walls are insulated. At time t > 0, the left- hand vertical wall temperature is suddenly raised to (T-h) over bar ((T-h) over bar > T-h) which introduces unsteadiness in the flow field. The partial differential equations governing the unsteady natural convection flow have been solved numerically using a finite control volume method. The computation has been carried out until the final steady state is reached. It is found that the average Nusselt number attains a minimum during the transient period and that the time required to reach the final steady state is longer for low Rayleigh number and shorter for high Rayleigh number.

A follow up MRI-based geometry and computational fluid dynamics study of carotid bifurcation

Cardiovascular disease is the leading causes of death in the developed world. Wall shear stress (WSS) is associated with the initiation and progression of atherogenesis. This study combined the recent advances in MR imaging and computational fluid dynamics (CFD) and evaluated the patient-specific carotid bifurcation. The patient was followed up for 3 years. The geometry changes (tortuosity, curvature, ICA/CCA area ratios, central to the cross-sectional curvature, maximum stenosis) and the CFD factors (Velocity distribute, Wall Shear Stress (WSS) and Oscillatory Shear Index (OSI)) were compared at different time points.The carotid stenosis was a slight increase in the central to the cross-sectional curvature, and it was minor and variable curvature changes for carotid centerline. The OSI distribution presents ahigh-values in the same region where carotid stenosis and normal border, indicating complex flow and recirculation.The significant geometric changes observed during the follow-up may also cause significant changes in bifurcation hemodynamics.

Study of tracheal collapsibility, compliance and stress by considering its asymmetric geometry

The shape of tracheal cartilage has been widely treated as symmetric in analytical and numerical models. However, according to both histological images and in vivo medical image, tracheal cartilage is of highly asymmetric shape. Taking the cartilage as symmetric structure will induce bias in calculation of the collapse behavior, as well as compliance and muscular stress. However, this has been rarely discussed. In this paper, tracheal collapse is represented by considering its asymmetric shape. For comparison, the symmetric shape, which is reconstructed by half of the cartilage, is also presented. A comparison of cross-sectional area, compliance of airway and stress in the muscular membrane, determined by asymmetric shape and symmetric shape is made. The result indicates that the symmetric assumption brings a small error, around 5% in predicting the cross-sectional area under loading conditions. The relative error of compliance is more than 10%. Particularly when the pressure is close to zero, the error could be more than 50%. The model considering the symmetric shape results in a significant difference in predicting stress in muscular membrane by either under- or over-estimating it. In conclusion, tracheal cartilage should not be treated as a symmetric structure. The results obtained in this study are helpful in evaluating the error induced by the assumption in geometry.