A study on the usefulness of Support Vector Machines for the realtime computational simulation of soft biological organs

Realistic and realtime computational simulation of soft biological organs (e.g., liver, kidney) is necessary when one tries to build a quality surgical simulator that can simulate surgical procedures involving these organs. Since the realistic simulation of these soft biological organs should account for both nonlinear material behavior and large deformation, achieving realistic simulations in realtime using continuum mechanics based numerical techniques necessitates the use of a supercomputer or a high end computer cluster which are costly. Hence there is a need to employ soft computing techniques like Support Vector Machines (SVMs) which can do function approximation, and hence could achieve physically realistic simulations in realtime by making use of just a desktop computer. Present work tries to simulate a pig liver in realtime. Liver is assumed to be homogeneous, isotropic, and hyperelastic. Hyperelastic material constants are taken from the literature. An SVM is employed to achieve realistic simulations in realtime, using just a desktop computer. The code for the SVM is obtained from [1]. The SVM is trained using the dataset generated by performing hyperelastic analyses on the liver geometry, using the commercial finite element software package ANSYS. The methodology followed in the present work closely follows the one followed in [2] except that [2] uses Artificial Neural Networks (ANNs) while the present work uses SVMs to achieve realistic simulations in realtime. Results indicate the speed and accuracy that is obtained by employing the SVM for the targeted realistic and realtime simulation of the liver.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

An integrated model for optimal reservoir operation for irrigation of multiple crops

An integratedm odel is developed,b asedo n seasonailn puts of reservoiri nflow and rainfall in the irrigated area, to determine the optimal reservoir release policies and irrigation allocationst o multiple crops.T he model is conceptuallym ade up of two modules. Module 1 is an intraseasonal allocation model to maximize the sum of relative yieldso f all crops,f or a givens tateo f the systemu, singl inear programming(L P). The module takes into account reservoir storage continuity, soil moisture balance, and crop root growthw ith time. Module 2 is a seasonaal llocationm odel to derive the steadys tate reservoiro peratingp olicyu sings tochastidc ynamicp rogramming(S DP). Reservoir storage, seasonal inflow, and seasonal rainfall are the state variables in the SDP. The objective in SDP is to maximize the expected sum of relative yields of all crops in a year.The resultso f module 1 and the transitionp robabilitieso f seasonailn flow and rainfall form the input for module 2. The use of seasonailn puts coupledw ith the LP-SDP solution strategy in the present formulation facilitates in relaxing the limitations of an earlier study,w hile affectinga dditionali mprovementsT. he model is applied to an existing reservoir in Karnataka State, India.

MODELING AND SIMULATION OF HYDRODYNAMIC INTERACTION OF DNA IN A MICRO-FLUIDIC CHANNEL

The motion of DNA (in the bulk solution) and the non-Newtonian effective fluid behavior are considered separately and self-consistently with the fluid motion satisfying the no-slip boundary condition on the surface of the confining geometry in the presence of channel pressure gradients. A different approach has been developed to model DNA in the micro-channel. In this study the DNA is assumed as an elastic chain with its characteristic Young's modulus, Poisson's ratio and density. The force which results from the fluid dynamic pressure, viscous forces and electromotive forces is applied to the elastic chain in a coupled manner. The velocity fields in the micro-channel are influenced by the transport properties. Simulations are carried out for the DNAs attached to the micro-fluidic wall. Numerical solutions based on a coupled multiphysics finite element scheme are presented. The modeling scheme is derived based on mass conservation including biomolecular mass, momentum balance including stress due to Coulomb force field and DNA-fluid interaction, and charge transport associated to DNA and other ionic complexes in the fluid. Variation in the velocity field for the non-Newtonian flow and the deformation of the DNA strand which results from the fluid-structure interaction are first studied considering a single DNA strand. Motion of the effective center of mass is analyzed considering various straight and coil geometries. Effects of DNA statistical parameters (geometry and spatial distribution of DNAs along the channel) on the effective flow behavior are analyzed. In particular, the dynamics of different DNA physical properties such as radius of gyration, end-to-end length etc. which are obtained from various different models (Kratky-Porod, Gaussian bead-spring etc.) are correlated to the nature of interaction and physical properties under the same background fluid environment.

Layerwise decomposition of water dynamics in reverse micelles: A simulation study of two-dimensional infrared spectrum

We present computer simulation study of two-dimensional infrared spectroscopy (2D-IR) of water confined in reverse micelles (RMs) of various sizes. The present study is motivated by the need to understand the altered dynamics of confined water by performing layerwise decomposition of water, with an aim to quantify the relative contributions of different layers water molecules to the calculated 2D-IR spectrum. The 0-1 transition spectra clearly show substantial elongation, due to in-homogeneous broadening and incomplete spectral diffusion, along the diagonal in the surface water layer of different sized RMs. Fitting of the frequency fluctuation correlation functions reveal that the motion of the surface water molecules is sub-diffusive and indicate the constrained nature of their dynamics. This is further supported by two peak nature of the angular analogue of van Hove correlation function. With increasing system size, the water molecules become more diffusive in nature and spectral diffusion almost completes in the central layer of the larger size RMs. Comparisons between experiments and simulations establish the correspondence between the spectral decomposition available in experiments with the spatial decomposition available in simulations. Simulations also allow a quantitative exploration of the relative role of water, sodium ions, and sulfonate head groups in vibrational dephasing. Interestingly, the negative cross correlation between force on oxygen and hydrogen of O-H bond in bulk water significantly decreases in the surface layer of each RM. This negative cross correlation gradually increases in the central water pool with increasing RMs size and this is found to be partly responsible for the faster relaxation rate of water in the central pool. (C) 2013 AIP Publishing LLC.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact' solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.

Role of the cloud adjustment time scale in simulation of the interannual variability of Indian summer monsoon

The simulation of precipitation in a general circulation model relying on relaxed mass flux cumulus parameterization scheme is sensitive to cloud adjustment time scale (CATS). In this study, the frequency of the dominant intra-seasonal mode and interannual variability of Indian summer monsoon rainfall (ISMR) simulated by an atmospheric general circulation model is shown to be sensitive to the CATS. It has been shown that a longer CATS of about 5 h simulates the spatial distribution of the ISMR better. El Nio Southern Oscillation-ISMR relationship is also sensitive to CATS. The equatorial Indian Ocean rainfall and ISMR coupling is sensitive to CATS. Our study suggests that a careful choice of CATS is necessary for adequate simulation of spatial pattern as well as interannual variation of Indian summer monsoon precipitation.

Monte Carlo simulation of light transport in tissue for optimizing light delivery in photoacoustic imaging of the sentinel lymph node

Noninvasive or minimally invasive identification of sentinel lymph node (SLN) is essential to reduce the surgical effects of SLN biopsy. Photoacoustic (PA) imaging of SLN in animal models has shown its promise for clinical use in the future. Here, we present a Monte Carlo simulation for light transport in the SLN for various light delivery configurations with a clinical ultrasound probe. Our simulation assumes a realistic tissue layer model and also can handle the transmission/reflectance at SLN-tissue boundary due to the mismatch of refractive index. Various light incidence angles show that for deeply situated SLNs the maximum absorption of light in the SLN is for normal incidence. We also show that if a part of the diffused reflected photons is reflected back into the skin using a reflector, the absorption of light in the SLN can be increased significantly to enhance the PA signal. (C) 2013 Society of Photo-Optical Instrumentation Engineers (SPIE)

System-level SoC near-field (NF) emissions: Simulation to measurement correlation

As System-on-Chip (SoC) designs migrate to 28nm process node and beyond, the electromagnetic (EM) co-interactions of the Chip-Package-Printed Circuit Board (PCB) becomes critical and require accurate and efficient characterization and verification. In this paper a fast, scalable, and parallelized boundary element based integral EM solutions to Maxwell equations is presented. The accuracy of the full-wave formulation, for complete EM characterization, has been validated on both canonical structures and real-world 3-D system (viz. Chip + Package + PCB). Good correlation between numerical simulation and measurement has been achieved. A few examples of the applicability of the formulation to high speed digital and analog serial interfaces on a 45nm SoC are also presented.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated, with specific reference to the speed and the accuracy offered by BEM. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. A pig liver is the biological organ considered. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

Computationally Efficient Model for Simulation of Boost Converter

This paper presents a computationally efficient model for a dc-dc boost converter, which is valid for continuous and discontinuous conduction modes; the model also incorporates significant non-idealities of the converter. Simulation of the dc-dc boost converter using an average model provides practically all the details, which are available from the simulation using the switching (instantaneous) model, except for the quantum of ripple in currents and voltages. A harmonic model of the converter can be used to evaluate the ripple quantities. This paper proposes a combined (average-cum-harmonic) model of the boost converter. The accuracy of the combined model is validated through extensive simulations and experiments. A quantitative comparison of the computation times of the average, combined and switching models are presented. The combined model is shown to be more computationally efficient than the switching model for simulation of transient and steady-state responses of the converter under various conditions.

Comparative Evaluation of Switching and Average Models of a DC-DC Boost Converter for Real-Time Simulation

This paper presents a comparative evaluation of the average and switching models of a dc-dc boost converter from the point of view of real-time simulation. Both the models are used to simulate the converter in real-time on a Field Programmable Gate Array (FPGA) platform. The converter is considered to function over a wide range of operating conditions, and could do transition between continuous conduction mode (CCM) and discontinuous conduction mode (DCM). While the average model is known to be computationally efficient from the perspective of off-line simulation, the same is shown here to consume more logical resources than the switching model for real-time simulation of the dc-dc converter. Further, evaluation of the boundary condition between CCM and DCM is found to be the main reason for the increased consumption of resources by the average model.

A Quick-Simulation Tool for Induction Motor Drives Controlled Using Advanced Space-Vector-Based PWM Techniques

Space-vector-based pulse width modulation (PWM) for a voltage source inverter (VSI) offers flexibility in terms of different switching sequences. Numerical simulation is helpful to assess the performance of a PWM method before actual implementation. A quick-simulation tool to simulate a variety of space-vector-based PWM strategies for a two-level VSI-fed squirrel cage induction motor drive is presented. The simulator is developed using C and Python programming languages, and has a graphical user interface (GUI) also. The prime focus being PWM strategies, the simulator developed is 40 times faster than MATLAB in terms of the actual time taken for a simulation. Simulation and experimental results are presented on a 5-hp ac motor drive.

A Study of Speed of the Boundary Element Method as applied to the Realtime Computational Simulation of Biological Organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated. Biological organs are assumed to follow linear elastostatic material behavior, and constant boundary element is the element type used. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. Results indicate that BEM is fast enough to provide for realtime graphics if biological organs are assumed to follow linear elastostatic material behavior. Although the present work does not conduct any simulation using nonlinear material models, results from using the linear elastostatic material model imply that it would be difficult to obtain realtime performance if highly nonlinear material models that properly characterize biological organs are used. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

Analytical solutions for movement of cold water thermal front in a heterogeneous geothermal reservoir

In the present study an analytical model has been presented to describe the transient temperature distribution and advancement of the thermal front generated due to the reinjection of heat depleted water in a heterogeneous geothermal reservoir. One dimensional heat transport equation in porous media with advection and longitudinal heat conduction has been solved analytically using Laplace transform technique in a semi infinite medium. The heterogeneity of the porous medium is expressed by the spatial variation of the flow velocity and the longitudinal effective thermal conductivity of the medium. A simpler solution is also derived afterwards neglecting the longitudinal conduction depending on the situation where the contribution to the transient heat transport phenomenon in the porous media is negligible. Solution for a homogeneous aquifer with constant values of the rock and fluid parameters is also derived with an aim to compare the results with that of the heterogeneous one. The effect of some of the parameters involved, on the transient heat transport phenomenon is assessed by observing the variation of the results with different magnitudes of those parameters. Results prove the heterogeneity of the medium, the flow velocity and the longitudinal conductivity to have great influence and porosity to have negligible effect on the transient temperature distribution. (C) 2013 Elsevier Inc. All rights reserved.