An economical method for direct numerical simulation studies of transitional round jets

A pseudo-spectral method based on Fourier expansions in a Cartesian coordinate system is shown to be an economical method for direct numerical simulation studies of transitional round jets, Several characteristics of the solutions are presented to establish the validity of the solutions in spite of the unnatural choices. We show that neither periodicity, nor the use of a Cartesian system have adversely affected the simulations, Instead, there are benefits in terms of ease of computing and lack of the usual restrictions due to grid structure near the jet axis. By computing the simultaneous evolution of passive scalers, the process of reaction in round jet burners, between a fuel-laden jet and an ambient oxidizer, was also simulated. Some typical solutions are shown and then the results of analysis of these data are summarized. (C) 2001 Elsevier Science Ltd, All rights reserved.

Direct numerical simulation of turbulent spots

A group of high-order finite-difference schemes for incompressible flow was implemented to simulate the evolution of turbulent spots in channel flows. The long-time accuracy of these schemes was tested by comparing the evolution of small disturbances to a plane channel flow against the growth rate predicted by linear theory. When the perturbation is the unstable eigenfunction at a Reynolds number of 7500, the solution grows only if there are a comparatively large number of (equispaced) grid points across the channel. Fifth-order upwind biasing of convection terms is found to be worse than second-order central differencing. But, for a decaying mode at a Reynolds number of 1000, about a fourth of the points suffice to obtain the correct decay rate. We show that this is due to the comparatively high gradients in the unstable eigenfunction near the walls. So, high-wave-number dissipation of the high-order upwind biasing degrades the solution especially. But for a well-resolved calculation, the weak dissipation does not degrade solutions even over the very long times (O(100)) computed in these tests. Some new solutions of spot evolution in Couette flows with pressure gradients are presented. The approach to self-similarity at long times can be seen readily in contour plots.

Reconstructing Solid Model from 2D Scanned Images of Biological Organs for Finite Element Simulation

This work presents a methodology to reconstruct 3D biological organs from image sequences or other scan data using readily available free softwares with the final goal of using the organs (3D solids) for finite element analysis. The methodology deals with issues such as segmentation, conversion to polygonal surface meshes, and finally conversion of these meshes to 3D solids. The user is able to control the detail or the level of complexity of the solid constructed. The methodology is illustrated using 3D reconstruction of a porcine liver as an example. Finally, the reconstructed liver is imported into the commercial software ANSYS, and together with a cyst inside the liver, a nonlinear analysis performed. The results confirm that the methodology can be used for obtaining 3D geometry of biological organs. The results also demonstrate that the geometry obtained by following this methodology can be used for the nonlinear finite element analysis of organs. The methodology (or the procedure) would be of use in surgery planning and surgery simulation since both of these extensively use finite elements for numerical simulations and it is better if these simulations are carried out on patient specific organ geometries. Instead of following the present methodology, it would cost a lot to buy a commercial software which can reconstruct 3D biological organs from scanned image sequences.

An improved Monte Carlo scheme for simulation of synthesis of nanoparticles in reverse micelles

An improved Monte Carlo technique is presented in this work to simulate nanoparticle formation through a micellar route. The technique builds on the simulation technique proposed by Bandyopadhyaya et al. (Langmuir 2000, 16, 7139) which is general and rigorous but at the same time very computation intensive, so much so that nanoparticle formation in low occupancy systems cannot be simulated in reasonable time. In view of this, several strategies, rationalized by simple mathematical analyses, are proposed to accelerate Monte Carlo simulations. These are elimination of infructuous events, removal of excess reactant postreaction, and use of smaller micelle population a large number of times. Infructuous events include collision of an empty micelle with another empty one or with another one containing only one molecule or only a solid particle. These strategies are incorporated in a new simulation technique which divides the entire micelle population in four classes and shifts micelles from one class to other as the simulation proceeds. The simulation results, throughly tested using chi-square and other tests, show that the predictions of the improved technique remain unchanged, but with more than an order of magnitude decrease in computational effort for some of the simulations reported in the literature. A post priori validation scheme for the correctness of the simulation results has been utilized to propose a new simulation strategy to arrive at converged simulation results with near minimum computational effort.

Finite element simulation of BAW propagation in inhomogeneous plate due to piezoelectric actuation

A set of finite elements (FEs) is formulated to analyze wave propagation through inhomogeneous material when subjected to mechanical, thermal loading or piezo-electric actuation. Elastic, thermal and electrical properties of the materials axe allowed to vary in length and thickness direction. The elements can act both as sensors and actuators. These elements are used to model wave propagation in functionally graded materials (FGM) and the effect of inhomogeneity in the wave is demonstrated. Further, a surface acoustic wave (SAW) device is modeled and wave propagation due to piezo-electric actuation from interdigital transducers (IDTs) is studied.

Simulation of the infrared spectra of N-methylformamide and N-methylacetamide by the extended molecular mechanics method

The infrared spectra of the matrix isolated species of N-methylformamide (NMF) and N-methylacetamide (NMA) and their N-deuterated molecules have been simulated by the extended molecular mechanics method using an empirical force field which includes charges and charge fluxes as coulombic potential parameters. The structural parameters and dipole. moments of NMF and NMA have. also been computed in satisfactory agreement with the experiment. Good agreement between experimental and calculated vibrational frequencies and infrared absorption band intensities for NMF and NMA and their deuterated molecules has been obtained. The vibrational assignments of NMF and NMA are-discussed taking also into account the infrared absorption intensities.

Impact of satellite observed microwave SST on the simulation of tropical cyclones

The impact of realistic representation of sea surface temperature (SST) on the numerical simulation of track and intensity of tropical cyclones formed over the north Indian Ocean is studied using the Weather Research and Forecast (WRF) model. We have selected two intense tropical cyclones formed over the Bay of Bengal for studying the SST impact. Two different sets of SSTs were used in this study: one from TRMM Microwave Imager (TMI) satellite and other is the weekly averaged Reynold's SST analysis from National Center for Environmental Prediction (NCEP). WRF simulations were conducted using the Reynold's and TMI SST as model boundary condition for the two cyclone cases selected. The TMI SST which has a better temporal and spatial resolution showed sharper gradient when compared to the Reynold's SST. The use of TMI SST improved the WRF cyclone intensity prediction when compared to that using Reynold's SST for both the cases studied. The improvements in intensity were mainly due to the improved prediction of surface latent and sensible heat fluxes. The use of TMI SST in place of Reynold's SST improved cyclone track prediction for Orissa super cyclone but slightly degraded track prediction for cyclone Mala. The present modeling study supports the well established notion that the horizontal SST gradient is one of the major driving forces for the intensification and movement of tropical cyclones over the Indian Ocean.

Simulation Blocks for TOSSIM-T2"

We develop several hardware and software simulation blocks for the TinyOS-2 (TOSSIM-T2) simulator. The choice of simulated hardware platform is the popular MICA2 mote. While the hardware simulation elements comprise of radio and external flash memory, the software blocks include an environment noise model, packet delivery model and an energy estimator block for the complete system. The hardware radio block uses the software environment noise model to sample the noise floor.The packet delivery model is built by establishing the SNR-PRR curve for the MICA2 system. The energy estimator block models energy consumption by Micro Controller Unit(MCU), Radio,LEDs, and external flash memory. Using the manufacturer’s data sheets we provide an estimate of the energy consumed by the hardware during transmission, reception and also track several of the MCUs states with the associated energy consumption. To study the effectiveness of this work, we take a case study of a paper presented in [1]. We obtain three sets of results for energy consumption through mathematical analysis, simulation using the blocks built into PowerTossim-T2 and finally laboratory measurements. Since there is a significant match between these result sets, we propose our blocks for T2 community to effectively test their application energy requirements and node life times.