A Study of Early Afterdepolarizations in a Model for Human Ventricular Tissue

Sudden cardiac death is often caused by cardiac arrhythmias. Recently, special attention has been given to a certain arrhythmogenic condition, the long-QT syndrome, which occurs as a result of genetic mutations or drug toxicity. The underlying mechanisms of arrhythmias, caused by the long-QT syndrome, are not fully understood. However, arrhythmias are often connected to special excitations of cardiac cells, called early afterdepolarizations (EADs), which are depolarizations during the repolarizing phase of the action potential. So far, EADs have been studied mainly in isolated cardiac cells. However, the question on how EADs at the single-cell level can result in fibrillation at the tissue level, especially in human cell models, has not been widely studied yet. In this paper, we study wave patterns that result from single-cell EAD dynamics in a mathematical model for human ventricular cardiac tissue. We induce EADs by modeling experimental conditions which have been shown to evoke EADs at a single-cell level: by an increase of L-type Ca currents and a decrease of the delayed rectifier potassium currents. We show that, at the tissue level and depending on these parameters, three types of abnormal wave patterns emerge. We classify them into two types of spiral fibrillation and one type of oscillatory dynamics. Moreover, we find that the emergent wave patterns can be driven by calcium or sodium currents and we find phase waves in the oscillatory excitation regime. From our simulations we predict that arrhythmias caused by EADs can occur during normal wave propagation and do not require tissue heterogeneities. Experimental verification of our results is possible for experiments at the cell-culture level, where EADs can be induced by an increase of the L-type calcium conductance and by the application of I-Kr blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium.

An enthalpy method for modeling eutectic solidification

This paper presents a new micro-scale model for solidification of eutectic alloys. The model is based on the enthalpy method and simulates the growth of adjacent alpha and beta phases from a melt of eutectic composition in a two-dimensional Eulerian framework. The evolution of the two phases is obtained from the solution of volume averaged energy and species transport equations which are formulated using the nodal enthalpy and concentration potential values. The three phases are tracked using the beta-phase fraction and the liquid fraction values in all the computational nodes. Solutal convection flow field in the domain is obtained from the solution of volume-averaged momentum and continuity equations. The governing equations are solved using a coupled explicit-implicit scheme. The model is qualitatively validated with Jackson-Hunt theory. Results show expected eutectic growth pattern and proper species transfer and diffusion field ahead of the interface. Capabilities of the model such as lamella width selection, division of lamella into thinner lamellae and the presence of solutal convection are successfully demonstrated. The present model can potentially be incorporated into the existing framework of enthalpy based micro-scale dendritic solidification models thus leading to an efficient generalized microstructure evolution model. (C) 2014 Elsevier Inc. All rights reserved.

Numerical study of the impact of inoculant and grain transport on macrosegregation and microstructure formation during solidification of an Al-22%Cu alloy

We investigate the impact of the nucleation law for nucleation on Al-Ti-B inoculant particles, of the motion of inoculant particles and of the motion of grains on the predicted macrosegregation and microstructure in a grain-refined Al-22 wt.% Cu alloy casting. We conduct the study by numerical simulations of a casting experiment in a side-cooled 76×76×254 mm sand mould. Macrosegregation and microstructure formation are studied with a volume-averaged two-phase model accounting for macroscopic heat and solute transport, melt convection, and transport of inoculant particles and equiaxed grains. On the microscopic scale it accounts for nucleation on inoculant particles with a given size distribution (and corresponding activation undercooling distribution)and for the growth of globular solid grains. The growth kinetics is described by accounting for limited solute diffusion in both liquid and solid phases and for convective effects. We show that the consideration of a size distribution of the inoculants has a strong impact on the microstructure(final grain size) prediction. The transport of inoculants significantly increases the microstructure heterogeneities and the grain motion refines the microstructure and reduces the microstructure heterogeneities.

Influence of exhaust gas heating and L/D ratios on the discharge efficiencies for an activated carbon natural gas storage system

A transient 2D axi-symmetric and lumped parameter (LP) model with constant outflow conditions have been developed to study the discharge capacity of an activated carbon bed. The predicted discharge times and variations in bed pressure and temperature are in good agreement with experimental results obtained from a 1.82 l adsorbed natural gas (ANG) storage system. Under ambient air conditions, a maximum temperature drop of 29.5 K and 45.5 K are predicted at the bed center for discharge rates of 1.0 l min(-1) and 5.0 l min(-1) respectively. The corresponding discharge efficiencies are 77% and 71.5% respectively with discharge efficiencies improving with decreasing outflow rates. Increasing the LID ratio from 1.9 to 7.8 had only a marginal increase in the discharge efficiency. Forced convection (exhaust gas) heating had a significant effect on the discharge efficiency, leading to efficiencies as high as 92.8% at a discharge of 1.0 l min(-1) and 88.7% at 5 l min(-1). Our study shows that the LP model can be reliably used to obtain discharge times due to the uniform pressure distributions in the bed. Temperature predictions with the LP model were more accurate at ambient conditions and higher discharge rates, due to greater uniformity in bed temperatures. For the low thermal conductivity carbon porous beds, our study shows that exhaust gas heating can be used as an effective and convenient strategy to improve the discharge characteristics and performance of an ANG system. (C) 2013 Elsevier Ltd. All rights reserved.

IS A DEEP ONE-CELL MERIDIONAL CIRCULATION ESSENTIAL FOR THE FLUX TRANSPORT SOLAR DYNAMO?

The solar activity cycle is successfully modeled by the flux transport dynamo, in which the meridional circulation of the Sun plays an important role. Most of the kinematic dynamo simulations assume a one-cell structure of the meridional circulation within the convection zone, with the equatorward return flow at its bottom. In view of the recent claims that the return flow occurs at a much shallower depth, we explore whether a meridional circulation with such a shallow return flow can still retain the attractive features of the flux transport dynamo (such as a proper butterfly diagram, the proper phase relation between the toroidal and poloidal fields). We consider additional cells of the meridional circulation below the shallow return flow-both the case of multiple cells radially stacked above one another and the case of more complicated cell patterns. As long as there is an equatorward flow in low latitudes at the bottom of the convection zone, we find that the solar behavior is approximately reproduced. However, if there is either no flow or a poleward flow at the bottom of the convection zone, then we cannot reproduce solar behavior. On making the turbulent diffusivity low, we still find periodic behavior, although the period of the cycle becomes unrealistically large. In addition, with a low diffusivity, we do not get the observed correlation between the polar field at the sunspot minimum and the strength of the next cycle, which is reproduced when diffusivity is high. On introducing radially downward pumping, we get a more reasonable period and more solar-like behavior even with low diffusivity.

Effect of angle of incidence on mixed convective wake dynamics and heat transfer past a square cylinder in cross flow at Re=100

In this paper, a numerical investigation is performed to study the mixed convective flow and heat transfer characteristics past a square cylinder in cross flow at incidence. Utilizing air (Pr = 0.71) as an operating fluid, computations are carried out at a representative Reynolds number (Re) of 100. Angles of incidences are varied as, 0 degrees <= alpha <= 45 degrees. Effect of superimposed positive and negative cross-flow buoyancy is brought about by varying the Richardson number (RI) in the range -1.0 <= Ri <= 1.0. The detail features of flow topology and heat transport are analyzed critically for different angles of incidences. The thermo fluidic forces acting on the cylinder during mixed convection are captured in terms of the drag (C-D), lift (C-L), and moment (C-M) coefficients. The results show that the lateral width of the cylinder wake reduces with increasing alpha and the isotherms spread out far wide. In the range 0 degrees < alpha < 45 degrees, C-D reduces with increasing Ri. The functional dependence of C-M with Ri reveals a linear relationship. Thermal boundary layer thickness reduces with increasing angle of incidences. The global rate of heat transfer from the cylinder increases with increasing alpha. (C) 2014 Elsevier Ltd. All rights reserved.

Imprints of Archean to Neoproterozoic crustal processes in the Madurai Block, Southern India

Madurai Block, the largest crustal block in the Southern Granulite Terrane (SGT) of Peninsular India, preserves the imprints of multistage tectonic evolution. Here, we present U-Pb and Hf isotope data on zircons from a charnockite-granite suite in the north-western part of this block. The oscillatory zoning, and the LREE to HREE enriched patterns of the zircons with positive Ce and negative Eu anomalies suggest that the zircon cores are of magmatic origin, with ages in the range of 2634-2435 Ma implying Neoarchean-Paleoproterozoic magmatism followed by subsequent metamorphism and protocontinent formation in the north-western part of the Madurai Block. A regional 550-500 Ma metamorphic overprint is also preserved in the zircons coinciding with the final amalgamation of the Gondwana supercontinent. The Hf isotopic data suggest that the granite and charnockite were derived from isotopically heterogeneous juvenile crustal domains and the charnockites show a significant contribution of mantle-derived components. Therefore, the Hf isotopic data reflect mixing of crustal and mantle-derived sources for the generation of Neoarchean crust in the north-western Madurai Block, possibly in a suprasubduction zone setting during continent building processes. (c) 2014 Elsevier Ltd. All rights reserved.

Local and remote impacts of direct aerosol forcing on Asian monsoon

The impact of heating by black carbon aerosols on Indian summer monsoon has remained inconclusive. Some investigators have predicted that black carbon aerosols reduce monsoon rainfall while others have argued that it will increase monsoon rainfall. These conclusions have been based on local influence of aerosols on the radiative fluxes. The impact of aerosol-like heating in one region on the rainfall in a remote region has not been examined in detail. Here, using an atmospheric general circulation model, it has been shown that remote influence of aerosol-like heating can be as important as local influence on Indian summer monsoon. Precipitation in northern Arabian Sea and north-west Indian region increased by 16% in June to July when aerosol-like heating were present globally. The corresponding increase in precipitation due to presence of aerosol-like heating only over South Asia (local impact) and East Asia (remote impact) were 28 and 13%, respectively. This enhancement in precipitation was due to destabilization of the atmosphere in pre-monsoon season that affected subsequent convection. Moreover, pre-monsoon heating of the lower troposphere changed the circulation substantially that enabled influx of more moisture over certain regions and reduced the moist static stability of the atmosphere. It has been shown that regional aerosol heating can have large impact on the phase of upper tropospheric Rossby wave in pre-monsoon season, which acts as a primary mechanism behind teleconnection and leads to the change in precipitation during monsoon season. These results demonstrate that changes in aerosol in one region can influence the precipitation in a remote region through changes in circulation.

A DYNAMO MODEL OF MAGNETIC ACTIVITY IN SOLAR-LIKE STARS WITH DIFFERENT ROTATIONAL VELOCITIES

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of 1M(circle dot) stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude f(m) of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of f(m). Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of the magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes rise radially from the bottom of the convection zone. Taking into account the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

Effect of Rayleigh numbers on the evolution of double-diffusive salt fingers

This is a transient two-dimensional numerical study of double-diffusive salt fingers in a two-layer heat-salt system for a wide range of initial density stability ratio (R-rho 0) and thermal Rayleigh numbers (Ra-T similar to 10(3) - 10(11)). Salt fingers have been studied for several decades now, but several perplexing features of this rich and complex system remain unexplained. The work in question studies this problem and shows the morphological variation in fingers from low to high thermal Rayleigh numbers, which have been missed by the previous investigators. Considerable variations in convective structures and evolution pattern were observed in the range of Ra-T used in the simulation. Evolution of salt fingers was studied by monitoring the finger structures, kinetic energy, vertical profiles, velocity fields, and transient variation of R-rho(t). The results show that large scale convection that limits the finger length was observed only at high Rayleigh numbers. The transition from nonlinear to linear convection occurs at about Ra-T similar to 10(8). Contrary to the popular notion, R-rho(t) first decrease during diffusion before the onset time and then increase when convection begins at the interface. Decrease in R-rho(t) is substantial at low Ra-T and it decreases even below unity resulting in overturning of the system. Interestingly, all the finger system passes through the same state before the onset of convection irrespective of Rayleigh number and density stability ratio of the system. (C) 2014 AIP Publishing LLC.

POLAR NETWORK INDEX AS A MAGNETIC PROXY FOR THE SOLAR CYCLE STUDIES

The Sun has a polar magnetic field which oscillates with the 11 yr sunspot cycle. This polar magnetic field is an important component of the dynamo process which operates in the solar convection zone and produces the sunspot cycle. We have direct systematic measurements of the Sun's polar magnetic field only from about the mid-1970s. There are, however, indirect proxies which give us information about this field at earlier times. The Ca-K spectroheliograms taken at the Kodaikanal Solar Observatory during 1904-2007 have now been digitized with 4k x 4k CCD and have higher resolution (similar to 0.86 arcsec) than the other available historical data sets. From these Ca-K spectroheliograms, we have developed a completely new proxy (polar network index, hereafter PNI) for the Sun's polar magnetic field. We calculate PNI from the digitized images using an automated algorithm and calibrate our measured PNI against the polar field as measured by the Wilcox Solar Observatory for the period 1976-1990. This calibration allows us to estimate the polar fields for the earlier period up to 1904. The dynamo calculations performed with this proxy as input data reproduce reasonably well the Sun's magnetic behavior for the past century.

Time-Domain Approach to Estimate Series Capacitance of an Isolated Phase Winding of a Transformer

Information available in frequency response data is equivalently available in the time domain as a response due to an impulse excitation. The idea to pursue this equivalence to estimate series capacitance is linked to the well-known fact that under impulse excitation, the line/neutral current in a transformer has three distinct components, of which, the initial capacitive component is the first to manifest, followed by the oscillatory and inductive components. Of these, the capacitive component is temporally well separated from the rest-a crucial feature permitting its direct access and analysis. Further, the winding initially behaves as a pure capacitive network, so the initial component must obviously originate from only the (series and shunt) capacitances. With this logic, it should therefore be possible to estimate series capacitance, just by measuring the initial capacitive component of line current and the total shunt capacitance. The principle of the method and details of its implementation on two actual isolated transformerwindings (uniformly wound) are presented. For implementation, a low-voltage recurrent surge generator, a current probe, and a digital oscilloscope are all that is needed. The method is simple and requires no programming and needs least user intervention, thus paving the way for its widespread use.

Flexibility in flapping foil suppresses meandering of induced jet in absence of free stream

Thrust-generating flapping foils are known to produce jets inclined to the free stream at high Strouhal numbers St = fA/U-infinity, where f is the frequency and A is the amplitude of flapping and U-infinity is the free-stream velocity. Our experiments, in the limiting case of St —> infinity (zero free-stream speed), show that a purely oscillatory pitching motion of a chordwise flexible foil produces a coherent jet composed of a reverse Benard-Karman vortex street along the centreline, albeit over a specific range of effective flap stiffnesses. We obtain flexibility by attaching a thin flap to the trailing edge of a rigid NACA0015 foil; length of flap is 0.79 c where c is rigid foil chord length. It is the time-varying deflections of the flexible flap that suppress the meandering found in the jets produced by a pitching rigid foil for zero free-stream condition. Recent experiments (Marais et al., J. Fluid Mech., vol. 710, 2012, p. 659) have also shown that the flexibility increases the St at which non-deflected jets are obtained. Analysing the near-wake vortex dynamics from flow visualization and particle image velocimetry (PIV) measurements, we identify the mechanisms by which flexibility suppresses jet deflection and meandering. A convenient characterization of flap deformation, caused by fluid-flap interaction, is through a non-dimensional effective stiffness', EI* = 8 EI/(rho V-TEmax(2) s(f) c(f)(3)/2), representing the inverse of the flap deflection due to the fluid-dynamic loading; here, EI is the bending stiffness of flap, rho is fluid density, V-TEmax is the maximum velocity of rigid foil trailing edge, s(f) is span and c(f) is chord length of the flexible flap. By varying the amplitude and frequency of pitching, we obtain a variation in EI* over nearly two orders of magnitude and show that only moderate EI*. (0.1 less than or similar to EI * less than or similar to 1 generates a sustained, coherent, orderly jet. Relatively `stiff' flaps (EI* greater than or similar to 1), including the extreme case of no flap, produce meandering jets, whereas highly `flexible' flaps (EI* less than or similar to 0.1) produce spread-out jets. Obtained from the measured mean velocity fields, we present values of thrust coefficients for the cases for which orderly jets are observed.

The role of buoyancy in polarity reversals of the geodynamo

We investigate polarity reversals in the geodynamo using a rotating, convection-driven dynamo model. As the flow in rapidly rotating convection is dominated by columns aligned with the axis of rotation, the focus is on the dynamics of columnar vortices. By studying the growth of a seed magnetic field to a stable axial dipole field, we show that the magnetic field acts in ways that significantly enhance the relative helicity between cyclonic and anticyclonic vortices. This flow asymmetry is the hallmark of a dipolar dynamo. Strong buoyancy, on the other hand, offsets the effect of the magnetic field, establishing parity between positive and negative vortices. As the dipole field is deprived of the helicity required to support itself, the dynamo is pushed into a reversing state. This is a likely regime for polarity reversals in the Earth's core. The integral lengthscale at which buoyancy injects energy is not significantly different from the convective flow lengthscale, which implies that buoyancy does not feed vortices at the small scales where non-linear inertia is present. The lengthscale at which the Lorentz force acts in the reversing dynamo is small, which may allow the passive presence of non-linear inertia in the small scales.

Influence of melt pouring temperature and plate inclination on solidification and microstructure of A356 aluminum alloy produced using oblique plate

The preparation of semisolid slurry of A356 aluminum alloy using an oblique plate was investigated. A356 alloy melt undergoes partial solidification when it flows down on an oblique plate cooled from underneath by counter flowing water. It results in continuous formation of columnar dendrites on plate wall. Due to forced convection, these dendrites are sheared off into equiaxed/fragmented grains and then washed away continuously to produce semisolid slurry at plate exit. Melt pouring temperature provides required condition of solidification whereas plate inclination enables necessary shear for producing semisolid slurry of desired quality. Slurry obtained was solidified in metal mould to produce semisolid-cast billets of desired microstructure. Furthermore, semisolid-cast billets were heat treated to improve surface quality. Microstructures of both semisolid-cast and heat-treated billets were analyzed. Effects of melt pouring temperature and plate inclination on solidification and microstructure of billets produced using oblique plate were described. The investigations involved four different melt pouring temperatures (620, 625, 630 and 635 degrees C) associated with four different plate inclinations (30 degrees, 45 degrees, 60 degrees and 75 degrees). Melt pouring temperature of 625 degrees C with plate inclination of 60 degrees shows fine and globular microstructures and it is the optimum.