Graphics Processing Units for the Real-time Linear Elastostatic Simulation of Liver

Biomedical engineering solutions like surgical simulators need High Performance Computing (HPC) to achieve real-time performance. Graphics Processing Units (GPUs) offer HPC capabilities at low cost and low power consumption. In this work, it is demonstrated that a liver which is discretized by about 2500 finite element nodes, can be graphically simulated in realtime, by making use of a GPU. Present work takes into consideration the time needed for the data transfer from CPU to GPU and back from GPU to CPU. Although behaviour of liver is very complicated, present computer simulation assumes linear elastostatics. One needs to use the commercial software ANSYS to obtain the global stiffness matrix of the liver. Results show that GPUs are useful for the real-time graphical simulation of liver, which in turn is needed in simulators that are used for training surgeons in laparoscopic surgery. Although the computer simulation should involve rendering also, neither rendering, nor the time needed for rendering and displaying the liver on a screen, is considered in the present work. The present work is just a demonstration of a concept; the concept is not really implemented and validated. Future work is to develop software which can accomplish real-time and very realistic graphical simulation of liver, with rendered image of liver on the screen changing in real-time according to the position of the surgical tool tip approximated as the mouse cursor in 3D.

New La2CuO4 derivatives La(2-2x)Sr(2x)Cu(1-x)M(x)O(4) (M=Ti, Mn, Fe, or Ru): A study of linear Cu-O-M electronic interaction in two dimensions

Oxides of the general formula La2-2xSr2xCu1-xII,M(x)(IV)O(4) (M = Ti, Mn, Fe, or Ru), crystallizing in the tetragonal K,NIF, structure, have been synthesized. For M=Ti, only the x=0,5 member could be prepared, while for M=Mn and Fe, the composition range is 0dence for an ordering of Cu(II) and M(IV) in the x = 0.5 members. While the members of the M=Ti, Mn, and Ru series are semiconducting/insulating, the members of the M=Fe series are metallic, showing a broad metal-semiconductor transition around 100 K for O Cu(III)-O-Fe(III) valence degeneracy. Increasing the strontium content at the expense of lanthanum in La2-2xSr2xCu1-xFexO4 for x less than or equal to 0.20 renders the samples metallic but not superconducting. (C) 1997 Academic Press.

Anomalous Behavior of Linear Hydrocarbon Chains in Water-DMSO Binary Mixture at Low DMSO Concentration

We find that at a mole fraction 0.05 of DMSO (x(DMSO) = 0.05) in aqueous solution, a linear hydrocarbon chain of intermediate length (n = 30-40) adopts the most stable collapsed conformation. In pure water, the same chain exhibits an intermittent oscillation between the collapsed and the extended coiled conformations. Even when the mole fraction of DMSO in the bulk is 0.05, the concentration of the same in the first hydration layer around the hydrocarbon of chain length 30 (n = 30) is as large as 17%. Formation of such hydrophobic environment around the hydrocarbon chain may be viewed as the reason for the collapsed conformation gaining additional stability. We find a second anomalous behavior to emerge near x(DMSO) = 0.15, due to a chain-like aggregation of the methyl groups of DMSO in water that lowers the relative concentration of the DMSO molecules in the hydration layer. We further find that as the concentration of DMSO is gradually increased, it progressively attains the extended coiled structure as the stable conformation. Although Flory-Huggins theory (for binary mixture solvent) fails to predict the anomaly at x(DMSO) = 0.05, it seems to capture the essence of the anomaly at 0.15.

Solving linear programming problems exactly

A linear programming problem in an inequality form having a bounded solution is solved error-free using an algorithm that sorts the inequalities, removes the redundant ones, and uses the p-adic arithmetic. (C) Elsevier Science Inc., 1997

Performance analysis of single-symbol maximum likelihood decodable linear STBCs

Performance of space-time block codes can be improved using the coordinate interleaving of the input symbols from rotated M-ary phase shift keying (MPSK) and M-ary quadrature amplitude modulation (MQAM) constellations. This paper is on the performance analysis of coordinate-interleaved space-time codes, which are a subset of single-symbol maximum likelihood decodable linear space-time block codes, for wireless multiple antenna terminals. The analytical and simulation results show that full diversity is achievable. Using the equivalent single-input single-output model, simple expressions for the average bit error rates are derived over flat uncorrelated Rayleigh fading channels. Optimum rotation angles are found by finding the minimum of the average bit error rate curves.

A DMFC with Methanol-Tolerant-Carbon-Supported-Pt-Pd-Alloy Cathode

Methanol-tolerant Pt-Pd alloy catalysts supported on to carbon with varying Pt:Pd atomic ratios of 1:1, 2:1 and 3:1 are prepared by a novel wet-chemical method and characterized using powder XRD, XPS, FESEM, EDAX and TEM techniques. The optimum atomic weight ratio for Pt to Pd in the carbon-supported alloy catalyst as established by linear-sweep voltammetry (LSV) and cell polarization studies is found to be 2:1. A direct methanol fuel cell (DMFC) employing carbon-supported Pt-Pd (2:1) alloy (Pt-Pd/C) catalyst as the cathode catalyst delivers a peak-power density of 115 mW/cm(2) at 70 degrees C as compared to peak-power density of 60 mW/cm(2) obtained with the DMFC employing carbon-supported Pt (Pt/C) catalyst operating under similar conditions. In the literature, DMFCs operating with Pt-TiO2 (2:1)/C and Pt-Au (2:1)/C methanol-tolerant cathodes are reported to exhibit maximum ORR activity among the group of these methanol-tolerant cathodes with varying catalysts compositions. Accordingly, the present study also provides an effective route to design methanol-tolerant-oxygen-reduction catalysts for DMFCs. (C) 2011 The Electrochemical Society. DOI: 10.1149/1.3596542] All rights reserved.

Weakly non-linear stability of a hydrodynamic accretion disc

When the cold accretion disc coupling between neutral gas and a magnetic field is so weak that the magnetorotational instability is less effective or even stops working, it is of prime interest to investigate the pure hydrodynamic origin of turbulence and transport phenomena. As the Reynolds number increases, the relative importance of the non-linear term in the hydrodynamic equation increases. In an accretion disc where the molecular viscosity is too small, the Reynolds number is large enough for the non-linear term to have new effects. We investigate the scenario of the `weakly non-linear' evolution of the amplitude of the linear mode when the flow is bounded by two parallel walls. The unperturbed flow is similar to the plane Couette flow, but with the Coriolis force included in the hydrodynamic equation. Although there is no exponentially growing eigenmode, because of the self-interaction, the least stable eigenmode will grow in an intermediate phase. Later, this will lead to higher-order non-linearity and plausible turbulence. Although the non-linear term in the hydrodynamic equation is energy-conserving, within the weakly non-linear analysis it is possible to define a lower bound of the energy (alpha A(c)(2), where A(c) is the threshold amplitude) needed for the flow to transform to the turbulent phase. Such an unstable phase is possible only if the Reynolds number >= 10(3-4). The numerical difficulties in obtaining such a large Reynolds number might be the reason for the negative result of numerical simulations on a pure hydrodynamic Keplerian accretion disc.

Organically templated linear and layered cadmium sulfates

By reacting cadmium salts with H2SO4 in the presence of organic amines or directly with amine sulfates under hydrothermal conditions, it has been possible to prepare three linear cadmium sulfates of linarite topology, with the compositions [H3N(CH2)(2)NH3](2)[CdCl2(SO4)][SO4].H2O, I, [HN(CH2)(6)NH][CdBr2(SO4)], II, [HN(CH2)(6)NH][CdCl2-(SO4)], III. A layered cadmium sulfate of composition [H3N(CH2)(3)NH3][Cd-2(H2O)(2)(SO4)(3)], IV, has also been obtained. These sulfates are the first examples of a family of organically templated metal sulfates with interesting structural features. In the linarite chains, the CdX4O2 (X = Cl, Br) octahedron shares two trans-edges to form an [Mphi(4)] (phi = anionic ligand) chain decorated by the SO4 tetrahedron that adopts a staggered arrangement on either side of the chain. IV is constructed by the fusion of four-membered ring ladders involving edge sharing between the sulfate tetrahedron and metal octahedron. IV appears to be the first member of a family of organically templated metal sulfates containing an octahedral-tetrahedral 2D net wherein the sulfate tetrahedron is connected at all four corners.

Phase diagram of a two-species lattice model with a linear instability

We discuss the properties of a one-dimensional lattice model of a driven system with two species of particles in which the mobility of one species depends on the density of the other. This model was introduced by Lahiri and Ramaswamy (Phys. Rev. Lett., 79, 1150 (1997)) in the context of sedimenting colloidal crystals, and its continuum version was shown to exhibit an instability arising from linear gradient couplings. In this paper we review recent progress in understanding the full phase diagram of the model. There are three phases. In the first, the steady state can be determined exactly along a representative locus using the condition of detailed balance. The system shows phase separation of an exceptionally robust sort, termed strong phase separation, which survives at all temperatures. The second phase arises in the threshold case where the first species evolves independently of the second, but the fluctuations of the first influence the evolution of the second, as in the passive scalar problem. The second species then shows phase separation of a delicate sort, in which long-range order coexists with fluctuations which do not damp down in the large-size limit. This fluctuation-dominated phase ordering is associated with power law decays in cluster size distributions and a breakdown of the Porod law. The third phase is one with a uniform overall density, and along a representative locus the steady state is shown to have product measure form. Density fluctuations are transported by two kinematic waves, each involving both species and coupled at the nonlinear level. Their dissipation properties are governed by the symmetries of these couplings, which depend on the overall densities. In the most interesting case,, the dissipation of the two modes is characterized by different critical exponents, despite the nonlinear coupling.

One-dimensional zinc phosphates with linear chain structure

Three one-dimensional zinc phosphates, [C5N2H14][Zn(HPO4)2], I, [C10N4H26][Zn(HPO4)2].2H2O II, and [C4N2H6]2[Zn(HPO4)], III, have been prepared employing hydro/solvothermal methods in the presence of organic amines. While I and II consist of linear chains of corner-shared four-membered rings, III is a polymeric wire where the amine molecule is directly bonded to the metal center. The wire, as well as the chain in these structures, are held together by hydrogen bond interactions involving the amine and the framework oxygens. The polymeric zinc phosphate with wire-like architecture, III, is only the second example of such architecture. Crystal data: I, monoclinic, P21/c (no. 14), a=8.603(2), b=13.529(2), c=10.880(1) Å, β=94.9(1)°, V=1261.6(1) Å3, Z=4, ρcalc.=1.893 gcm−3, μ(MoKα)=2.234 mm−1, R1=0.032, wR2=0.086, [1532 observed reflections with I>2σ(I)], II, orthorhombic, Pbca (no. 61), a=8.393(1), b=15.286(1), c=22.659(1) Å, V=2906.9(2) Å3, Z=8, ρcalc.=1.794 gcm−3, μ(MoKα)=1.957 mm−1, R1=0.055, wR2=0.11, [1565 observed reflections with I>2σ(I) and III, monoclinic, P21/c (no. 14), a=8.241(1), b=13.750(2), c=10.572(1) Å, β=90.9(1)°, V=1197.7(2) Å3, Z=4, ρcalc.=1.805 gcm−3, μ(MoKα)=2.197 mm−1, R1=0.036, wR2=0.10, [1423 observed reflections with I>2σ(I)].

Quadratic Lyapunov functions for linear systems Journal

The paper deals with the existence of a quadratic Lyapunov function V = x′P(t)x for an exponentially stable linear system with varying coefficients described by the vector differential equation S0305004100044777_inline1 The derivative dV/dt is allowed to be strictly semi-(F) and the locus dV/dt = 0 does not contain any arc of the system trajectory. It is then shown that the coefficient matrix A(t) of the exponentially stable sy

Effect of screening of intermicellar interactions on the linear and nonlinear rheology of a viscoelastic gel

We report our studies of the linear and nonlinear rheology of aqueous solutions of the surfactant cetyl trimethylammonium tosylate (CTAT) with varying amounts of sodium chloride (NaCl). The CTAT concentration is fixed at 42 mM, and the salt concentration is varied between 0 and 120 mM. On increasing the salt (NaCl) concentration, we see three distinct regimes in the zero-shear viscosity and the high-frequency plateau modulus data. In regime 1, the zero-shear viscosity shows a weak increase with salt concentration due to enhanced micellar growth. The decrease in the zero-shear viscosities with salt concentration in regimes II and III can be explained in terms of intermicellar branching. The most intriguing feature of our data, however, is the anomalous behavior of the high-frequency plateau modulus in regime II (0.12 less than or equal to [NaCl]/[CTAT] less than or equal to 1.42). In this regime, the plateau modulus increases with an increase in NaCl concentration. This is highly interesting, since the correlation length of concentration fluctuations and hence the plateau modulus G(0) are not expected to change appreciably in the semidilute regime. We propose to explain the changes in regime II in terms of a possible unbinding of the organic counterions (tosylate) from the CTA(+) surfaces on the addition of NaCl. In the nonlinear flow curves of the samples with high salt content, significant deviations from the predictions of the Giesekus model for entangled micelles are observed.

Determination of solid fraction evolution during solidification of aluminum alloy in presence of linear electromagnetic stirring

In many industrial casting processes, knowledge of the solid fraction evolution during the solidification process is a key factor in determining the process parameters such as cooling rate, stirring intensity and in estimating the total solidification time. In the present work, a new method of estimating solid fraction is presented, which is based on calorimetric principles. In this method, the cooling curve data at each point in the melt, along with the thermal boundary conditions, are used to perform energy balance in the mould, from which solid fraction generation during any time interval can be estimated. This method is applied to the case of a rheocasting process, in which Al-Si alloy (A356 alloy) is solidified by stirring in a cylindrical mould placed in the annulus of a linear electromagnetic stirrer. The metal in the mould is simultaneously cooled and stirred to produce a cylindrical billet with non-dendritic globular microstructure. Temperature is measured at key locations in the mould to assess the various heat exchange processes prevalent in the mould and to monitor the solidification rate. The results obtained by energy balance method are compared with those by the conventional procedure of calculating solid fraction using the Schiel equation.

Qualitative theorem proving in linear constraints

We know, from the classical work of Tarski on real closed fields, that elimination is, in principle, a fundamental engine for mechanized deduction. But, in practice, the high complexity of elimination algorithms has limited their use in the realization of mechanical theorem proving. We advocate qualitative theorem proving, where elimination is attractive since most processes of reasoning take place through the elimination of middle terms, and because the computational complexity of the proof is not an issue. Indeed what we need is the existence of the proof and not its mechanization. In this paper, we treat the linear case and illustrate the power of this paradigm by giving extremely simple proofs of two central theorems in the complexity and geometry of linear programming.