Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Tracking the states of a nonlinear and nonstationary system in the weight-space of artificial neural networks

We propose a novel interpretation and usage of Neural Network (NN) in modeling physiological signals, which are allowed to be nonlinear and/or nonstationary. The method consists of training a NN for the k-step prediction of a physiological signal, and then examining the connection-weight-space (CWS) of the NN to extract information about the signal generator mechanism. We de. ne a novel feature, Normalized Vector Separation (gamma(ij)), to measure the separation of two arbitrary states i and j in the CWS and use it to track the state changes of the generating system. The performance of the method is examined via synthetic signals and clinical EEG. Synthetic data indicates that gamma(ij) can track the system down to a SNR of 3.5 dB. Clinical data obtained from three patients undergoing carotid endarterectomy of the brain showed that EEG could be modeled (within a root-means-squared-error of 0.01) by the proposed method, and the blood perfusion state of the brain could be monitored via gamma(ij), with small NNs having no more than 21 connection weight altogether.

Transient temperature rise in a mouse due to low-frequency regional hyperthermia

A refined nonlinear heat transfer model of a mouse has been developed to simulate the transient temperature rise in a neoplastic tumour and neighbouring tissue during regional hyperthermia using a 150 kHz inductive coil. In this study, we incorporate various bio-energetic enhancements to the heat transfer equation and numerical validations based on experimental findings for the mouse, in terms of nonlinear metabolic heat production, homeothermy, blood perfusion parameters, thermoregulation, psychological and physiological effects. The discretized bio-heat transfer equation has been validated with the commercial software FEMLAB on a canonical multi-sphere object before applying the scheme to the inhomogeneous mouse voxel phantom. The time-dependent numerical results of regional hyperthermia of mouse thigh have been compared with the available experimental temperature results with only a few small disparities. During the first 20 min of local unfocused heating, the temperature in the tumour and the surrounding tissue increased by around 7.5 degrees C. The objective of this preliminary study was to develop a validated electrothermal numerical scheme for inductive hyperthermia of a small mammal with the intention of expanding the model into a complete numerical solution involving ferromagnetic nanoparticles for targeted heating of tumours at low frequencies. In addition, the numerical scheme herein could assist in optimizing and tailoring of focused electromagnetic fields for hyperthermia.

Yield criterion of porous materials subjected to complex stress states

Finite-element simulations are used to obtain many thousands of yield points for porous materials with arbitrary void-volume fractions with spherical voids arranged in simple cubic, body-centred cubic and face-centred cubic three-dimensional arrays. Multi-axial stress states are explored. We show that the data may be fitted by a yield function which is similar to the Gurson-Tvergaard-Needleman (GTN) form, but which also depends on the determinant of the stress tensor, and all additional parameters may be expressed in terms of standard GTN-like parameters. The dependence of these parameters on the void-volume fraction is found. (c) 2006 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Yield functions for porous materials with cubic symmetry using different definitions of yield

Plastic yield criteria for porous ductile materials are explored numerically using the finite-element technique. The cases of spherical voids arranged in simple cubic, body-centred cubic and face-centred cubic arrays are investigated with void volume fractions ranging from 2 % through to the percolation limit (over 90 %). Arbitrary triaxial macroscopic stress states and two definitions of yield are explored. The numerical data demonstrates that the yield criteria depend linearly on the determinant of the macroscopic stress tensor for the case of simple-cubic and body-centred cubic arrays - in contrast to the famous Gurson-Tvergaard-Needleman (GTN) formula - while there is no such dependence for face-centred cubic arrays within the accuracy of the finite-element discretisation. The data are well fit by a simple extension of the GTN formula which is valid for all void volume fractions, with yield-function convexity constraining the form of the extension in terms of parameters in the original formula. Simple cubic structures are more resistant to shear, while body-centred and face-centred structures are more resistant to hydrostatic pressure. The two yield surfaces corresponding to the two definitions of yield are not related by a simple scaling.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.

Performance of a scramjet combustor with combined normal and tangential fuel injection

The performance of a scramjet combustor with combined normal and tangential injection was experimentally investigated. Experiments were performed on a 500-mm cylindrical scramjet combustor at a freestream Mach number of 4.5, a nozzle supply pressure of 35.8 MPa, and a nozzle supply enthalpy of 5.8 MJ/kg. Hydrogen fuel was injected normally through portholes to promote combustion and tangentially through a slot to reduce viscous drag. A series of fuel injectors were used to vary the proportion of tangential to normal fuel between 45 and 100%. Reductions in the viscous drag of up to 25% were observed with the greatest reductions occurring at the lowest total equivalence ratio tested for each injector. However, the average pressure produced by combustion with combined normal and tangential injection was approximately 50% less than that produced by normal injection alone. An analysis of the change in specific impulse of the scramjet combustor indicated that the best overall performance was produced by 100% normal injection.

Getting under the skin

The needle and syringe may be on their way out as materials scientists begin to find better ways of delivering drugs through the skin.

Assessment of epidermal cell viability by near infrared multi-photon microscopy following ballistic delivery of gold micro-particles

The use of gene guns in ballistically delivering DNA vaccine coated gold micro-particles to skin can potentially damage targeted cells, therefore influencing transfection efficiencies. In this paper, we assess cell death in the viable epidermis by non-invasive near infrared two-photon microscopy following micro-particle bombardment of murine skin. We show that the ballistic delivery of micro-particles to the viable epidermis can result in localised cell death. Furthermore, experimental results show the degree of cell death is dependant on the number of micro-particles delivered per unit of tissue surface area. Micro-particles densities of 0.16 +/- 0.27 (mean +/- S.D.), 1.35 +/- 0.285 and 2.72 +/- 0.47 per 1000 mu m(2) resulted in percent deaths of 3.96 +/- 5.22, 45.91 +/- 10.89, 90.52 +/- 12.28, respectively. These results suggest that optimization of transfection by genes administered with gene guns is - among other effects - a compromise of micro-particle payload and cell death. (c) 2005 Elsevier Ltd. All rights reserved.

Evaluation of the indirect measures of rock brittleness and fracture toughness in rock cutting

Performance prediction models for partial face mechanical excavators, when developed in laboratory conditions, depend on relating the results of a set of rock property tests and indices to specific cutting energy (SE) for various rock types. There exist some studies in the literature aiming to correlate the geotechnical properties of intact rocks with the SE, especially for massive and widely jointed rock environments. However, those including direct and/or indirect measures of rock fracture parameters such as rock brittleness and fracture toughness, along with the other rock parameters expressing different aspects of rock behavior under drag tools (picks), are rather limited. With this study, it was aimed to investigate the relationships between the indirect measures of rock brittleness and fracture toughness and the SE depending on the results of a new and two previous linear rock cutting programmes. Relationships between the SE, rock strength parameters, and the rock index tests have also been investigated in this study. Sandstone samples taken from the different fields around Ankara, Turkey were used in the new testing programme. Detailed mineralogical analyses, petrographic studies, and rock mechanics and rock cutting tests were performed on these selected sandstone specimens. The assessment of rock cuttability was based on the SE. Three different brittleness indices (B1, B2, and B4) were calculated for sandstones samples, whereas a toughness index (T-i), being developed by Atkinson et al.(1), was employed to represent the indirect rock fracture toughness. The relationships between the SE and the large amounts of new data obtained from the mineralogical analyses, petrographic studies, rock mechanics, and linear rock cutting tests were evaluated by using bivariate correlation and curve fitting techniques, variance analysis, and Student's t-test. Rock cutting and rock property testing data that came from well-known studies of McFeat-Smith and Fowell(2) and Roxborough and Philips(3) have also been employed in statistical analyses together with the new data. Laboratory tests and subsequent analyses revealed that there were close correlations between the SE and B4 whereas no statistically significant correlation has been found between the SE and T-i. Uniaxial compressive and Brazilian tensile strengths and Shore scleroscope hardness of sandstones also exhibited strong relationships with the SE. NCB cone indenter test had the greatest influence on the SE among the other engineering properties of rocks, confirming the previous studies in rock cutting and mechanical excavation. Therefore, it was recommended to employ easy-to-use index tests of NCB cone indenter and Shore scleroscope in the estimation of laboratory SE of sandstones ranging from very low to high strengths in the absence of a rock cutting rig to measure it until the easy-to-use universal measures of the rock brittleness and especially the rock fracture toughness, being an intrinsic rock property, are developed.

On the multicomponent polynomial solution models

The present study addresses the problem of predicting the properties of multicomponent systems from those of corresponding binary systems. Two types of multicomponent polynomial models have been analysed. A probabilistic interpretation of the parameters of the Polynomial model, which explicitly relates them with the Gibbs free energies of the generalised quasichemical reactions, is proposed. The presented treatment provides a theoretical justification for such parameters. A methodology of estimating the ternary interaction parameter from the binary ones is presented. The methodology provides a way in which the power series multicomponent models, where no projection is required, could be incorporated into the Calphad approach.

Beach face and berm morphodynamics fronting a coastal lagoon

This study documents two different modes of berm development: (I) vertical growth at spring tides or following significant beach cut due to substantial swash overtopping, and (2) horizontal progradation at neap tides through the formation of a proto-berm located lower and further seaward of the principal berm. Concurrent high-frequency measurements of bed elevation and the associated wave runup distribution reveal the details of each of these berm growth modes. In mode I sediment is eroded from the inner surf and lower swash zone where swash interactions are prevalent. The net transport of this sediment is landward only, resulting in accretion onto the upper beach face and over the berm crest. The final outcome is a steepening of the beach face gradient, a change in the profile shape towards concave and rapid vertical and horizontal growth of the berm. In mode 2 sediment is eroded from the lower two-thirds of the active swash zone during the rising tide and is transported both landward and seaward. On the falling tide sediment is eroded from the inner surf and transported landward to backfill the zone eroded on the rising tide. The net result is relatively slow steepening of the beach face, a change of the profile shape towards convex, and horizontal progradation through the formation of a neap berm. The primary factor determining which mode of berm growth occurs is the presence or absence of swash overtopping at the time of sediment accumulation on the beach face. This depends on the current phase of the spring-neap tide cycle, the wave runup height (and indirectly offshore wave conditions) and the height of the pre-existing berm. A conceptual model for berm morphodynamics is presented, based on sediment transport shape functions measured during the two modes of berm growth. (c) 2006 Elsevier B.V. All rights reserved.

Swash-aquifer interaction in the vicinity of the water table exit point on a sandy beach

The coupling of sandy beach aquifers with the swash zone in the vicinity of the water table exit point is investigated through simultaneous measurements of the instantaneous shoreline (swash front) location, pore pressures and the water table exit point. The field observations reveal new insights into swash-aquifer coupling not previously gleaned from measurements of pore pressure only. In particular, for the case where the exit point is seaward of the observation point, the pore pressure response is correlated with the distance between the exit point and the shoreline in that when the distance is large the rate of pressure drop is fast and when the distance is small the rate decreases. The observations expose limitations in a simple model describing exit point dynamics which is based only on the force balance on a particle of water at the sand surface and neglects subsurface pressures. A new modified form of the model is shown to significantly improve the model-data comparison through a parameterization of the effects of capillarity into the aquifer storage coefficient. The model enables sufficiently accurate predictions of the exit point to determine when the swash uprush propagates over a saturated or a partially saturated sand surface, potentially an important factor in the morphological evolution of the beach face. Observations of the shoreward propagation of the swash-induced pore pressure waves ahead of the runup limit shows that the magnitude of the pressure fluctuation decays exponentially and that there is a linear increase in time lags, behavior similar to that of tidally induced water table waves. The location of the exit point and the intermittency of wave runup events is also shown to be significant in terms of the shore-normal energy distribution. Seaward of the mean exit point location, peak energies are small because of the saturated sand surface within the seepage face acting as a "rigid lid'' and limiting pressure fluctuations. Landward of the mean exit point the peak energies grow before decreasing landward of the maximum shoreline position.