Adaptive Executions of Multi-Physics Coupled Applications on Batch Grids

Long running multi-physics coupled parallel applications have gained prominence in recent years. The high computational requirements and long durations of simulations of these applications necessitate the use of multiple systems of a Grid for execution. In this paper, we have built an adaptive middleware framework for execution of long running multi-physics coupled applications across multiple batch systems of a Grid. Our framework, apart from coordinating the executions of the component jobs of an application on different batch systems, also automatically resubmits the jobs multiple times to the batch queues to continue and sustain long running executions. As the set of active batch systems available for execution changes, our framework performs migration and rescheduling of components using a robust rescheduling decision algorithm. We have used our framework for improving the application throughput of a foremost long running multi-component application for climate modeling, the Community Climate System Model (CCSM). Our real multi-site experiments with CCSM indicate that Grid executions can lead to improved application throughput for climate models.

Development of novel multiple quantummethodologies for the analyses of complex protonNMR spectra of scalar coupled spins

One of the significant advancements in Nuclear Magnetic Resonance spectroscopy (NMR) in combating the problem of spectral complexity for deriving the structure and conformational information is the incorporation of additional dimension and to spread the information content in a two dimensional space. This approach together with the manipulation of the dynamics of nuclear spins permitted the designing of appropriate pulse sequences leading to the evolution of diverse multidimensional NMR experiments. The desired spectral information can now be extracted in a simplified and an orchestrated manner. The indirect detection of multiple quantum (MQ) NMR frequencies is a step in this direction. The MQ technique has been extensively used in the study of molecules aligned in liquid crystalline media to reduce spectral complexity and to determine molecular geometries. Unlike in dipolar coupled systems, the size of the network of scalar coupled spins is not big in isotropic solutions and the MQ 1H detection is not routinely employed,although there are specific examples of spin topology filtering. In this brief review, we discuss our recent studies on the development and application of multiple quantum correlation and resolved techniques for the analyses of proton NMR spectra of scalar coupled spins.

Repulsive fermions in optical lattices: Phase separation versus coexistence of antiferromagnetism and d-wave superfluidity

We investigate a system of fermions on a two-dimensional optical square lattice in the strongly repulsive coupling regime. In this case, the interactions can be controlled by laser intensity as well as by Feshbach resonance. We compare the energetics of states with resonating valence bond d-wave superfluidity, antiferromagnetic long-range order, and a homogeneous state with coexistence of superfluidity and antiferromagnetism. Using a variational formalism, we show that the energy density of a hole e(hole)(x) has a minimum at doping x = x(c) that signals phase separation between the antiferromagnetic and d-wave paired superfluid phases. The energy of the phase-separated ground state is, however, found to be very close to that of a homogeneous state with coexisting antiferromagnetic and superfluid orders. We explore the dependence of the energy on the interaction strength and on the three-site hopping terms and compare with the nearest-neighbor hopping t-J model.

Seismic hazard map of Coimbatore using subsurface fault rupture

This study presents the future seismic hazard map of Coimbatore city, India, by considering rupture phenomenon. Seismotectonic map for Coimbatore has been generated using past earthquakes and seismic sources within 300 km radius around the city. The region experienced a largest earthquake of moment magnitude 6.3 in 1900. Available earthquakes are divided into two categories: one includes events having moment magnitude of 5.0 and above, i.e., damaging earthquakes in the region and the other includes the remaining, i.e., minor earthquakes. Subsurface rupture character of the region has been established by considering the damaging earthquakes and total length of seismic source. Magnitudes of each source are estimated by assuming the subsurface rupture length in terms of percentage of total length of sources and matched with reported earthquake. Estimated magnitudes match well with the reported earthquakes for a RLD of 5.2% of the total length of source. Zone of influence circles is also marked in the seismotectonic map by considering subsurface rupture length of fault associated with these earthquakes. As earthquakes relive strain energy that builds up on faults, it is assumed that all the earthquakes close to damaging earthquake have released the entire strain energy and it would take some time for the rebuilding of strain energy to cause a similar earthquake in the same location/fault. Area free from influence circles has potential for future earthquake, if there is seismogenic source and minor earthquake in the last 20 years. Based on this rupture phenomenon, eight probable locations have been identified and these locations might have the potential for the future earthquakes. Characteristic earthquake moment magnitude (M-w) of 6.4 is estimated for the seismic study area considering seismic sources close to probable zones and 15% increased regional rupture character. The city is divided into several grid points at spacing of 0.01 degrees and the peak ground acceleration (PGA) due to each probable earthquake is calculated at every grid point in city by using the regional attenuation model. The maximum of all these eight PGAs is taken for each grid point and the final PGA map is arrived. This map is compared to the PGA map developed based on the conventional deterministic seismic hazard analysis (DSHA) approach. The probable future rupture earthquakes gave less PGA than that of DSHA approach. The occurrence of any earthquake may be expected in near future in these eight zones, as these eight places have been experiencing minor earthquakes and are located in well-defined seismogenic sources.

Learning polyhedral classifiers using logistic function

In this paper we propose a new algorithm for learning polyhedral classifiers. In contrast to existing methods for learning polyhedral classifier which solve a constrained optimization problem, our method solves an unconstrained optimization problem. Our method is based on a logistic function based model for the posterior probability function. We propose an alternating optimization algorithm, namely, SPLA1 (Single Polyhedral Learning Algorithm1) which maximizes the loglikelihood of the training data to learn the parameters. We also extend our method to make it independent of any user specified parameter (e.g., number of hyperplanes required to form a polyhedral set) in SPLA2. We show the effectiveness of our approach with experiments on various synthetic and real world datasets and compare our approach with a standard decision tree method (OC1) and a constrained optimization based method for learning polyhedral sets.

Coupled Wavenumbers in an Infinite Flexible Fluid-Filled Circular Cylindrical Shell: Comparison between Different Shell Theories

Analytical expressions are found for the wavenumbers in an infinite flexible in vacuo I fluid-filled circular cylindrical shell based on different shell-theories using asymptotic methods. Donnell-Mushtari theory (the simplest shell theory) and four higher order theories, namely Love-Timoshenko, Goldenveizer-Novozhilov, Flugge and Kennard-simplified are considered. Initially, in vacuo and fluid-coupled wavenumber expressions are presented using the Donnell-Mushtari theory. Subsequently, the wavenumbers using the higher order theories are presented as perturbations on the Donnell-Mushtari wavenumbers. Similarly, expressions for the resonance frequencies in a finite shell are also presented, using each shell theory. The basic differences between the theories being what they are, the analytical expressions obtained from the five theories allow one to see how these differences propagate into the asymptotic expansions. Also, they help to quantify the difference between the theories for a wide range of parameter values such as the frequency range, circumferential order, thickness ratio of the shell, etc.

Density matrix renormalization group algorithm for Bethe lattices of spin-1/2 or spin-1 sites with Heisenberg antiferromagnetic exchange

A density matrix renormalization group (DMRG) algorithm is presented for the Bethe lattice with connectivity Z = 3 and antiferromagnetic exchange between nearest-neighbor spins s = 1/2 or 1 sites in successive generations g. The algorithm is accurate for s = 1 sites. The ground states are magnetic with spin S(g) = 2(g)s, staggered magnetization that persists for large g > 20, and short-range spin correlation functions that decrease exponentially. A finite energy gap to S > S(g) leads to a magnetization plateau in the extended lattice. Closely similar DMRG results for s = 1/2 and 1 are interpreted in terms of an analytical three-site model.

A coupled VOF-IBM-enthalpy approach for modeling motion and growth of equiaxed dendrites in a solidifying melt

A coupled methodology for simulating the simultaneous growth and motion of equiaxed dendrites in solidifying melts is presented. The model uses the volume-averaging principles and combines the features of the enthalpy method for modeling growth, immersed boundary method for handling the rigid solid-liquid interfaces, and the volume of fluid method for tracking the advection of the dendrite. The algorithm also performs explicit-implicit coupling between the techniques used. A two-dimensional framework with incompressible and Newtonian fluid is considered. Validation with available literature is performed and dendrite growth in the presence of rotational and buoyancy driven flow fields is studied. It is seen that the flow fields significantly alter the position and morphology of the dendrites. (C) 2012 Elsevier Inc. All rights reserved.

Coupled phonons, magnetic excitations, and ferroelectricity in AlFeO3: Raman and first-principles studies

We determine the nature of coupled phonons and magnetic excitations in AlFeO3 using inelastic light scattering from 5 to 315 K covering a spectral range from 100 to 2200 cm(-1) and complementary first-principles density functional theory-based calculations. A strong spin-phonon coupling and magnetic ordering-induced phonon renormalization are evident in (1) anomalous temperature dependence of many modes with frequencies below 850 cm(-1), particularly near the magnetic transition temperature T-c approximate to 250 K, and (2) distinct changes in band positions of high-frequency Raman bands between 1100 and 1800 cm(-1); in particular, a broad mode near 1250 cm(-1) appears only below T-c, attributed to the two-magnon Raman scattering. We also observe weak anomalies in the mode frequencies similar to 100 K due to a magnetically driven ferroelectric phase transition. Understanding of these experimental observations has been possible on the basis of first-principles calculations of the phonons' spectrum and their coupling with spins.

Nonlocal continuum mechanics formulation for axial, flexural, shear and contraction coupled wave propagation in single walled carbon nanotubes

This paper presents the effect of nonlocal scaling parameter on the coupled i.e., axial, flexural, shear and contraction, wave propagation in single-walled carbon nanotubes (SWCNTs). The axial and transverse motion of SWCNT is modeled based on first order shear deformation theory (FSDT) and thickness contraction. The governing equations are derived based on nonlocal constitutive relations and the wave dispersion analysis is also carried out. The studies shows that the nonlocal scale parameter introduces certain band gap region in all wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. Explicit expressions are derived for cut-off and escape frequencies of all waves in SWCNT. It is also shown that the cut-off frequencies of shear and contraction mode are independent of the nonlocal scale parameter. The results provided in this article are new and are useful guidance for the study and design of the next generation of nanodevices that make use of the coupled wave propagation properties of single-walled carbon nanotubes.

Synchronization of globally coupled two-state stochastic oscillators with a state-dependent refractory period

We present a model of identical coupled two-state stochastic units, each of which in isolation is governed by a fixed refractory period. The nonlinear coupling between units directly affects the refractory period, which now depends on the global state of the system and can therefore itself become time dependent. At weak coupling the array settles into a quiescent stationary state. Increasing coupling strength leads to a saddle node bifurcation, beyond which the quiescent state coexists with a stable limit cycle of nonlinear coherent oscillations. We explicitly determine the critical coupling constant for this transition.

How good are the simulations of tropical SST-rainfall relationship by IPCC AR4 atmospheric and coupled models?

The failure of atmospheric general circulation models (AGCMs) forced by prescribed SST to simulate and predict the interannual variability of Indian/Asian monsoon has been widely attributed to their inability to reproduce the actual sea surface temperature (SST)-rainfall relationship in the warm Indo-Pacific oceans. This assessment is based on a comparison of the observed and simulated correlation between the rainfall and local SST. However, the observed SSTconvection/rainfall relationship is nonlinear and for this a linear measure such as the correlation is not an appropriate measure. We show that the SST-rainfall relationship simulated by atmospheric and coupled general circulation models in IPCC AR4 is nonlinear, as observed, and realistic over the tropical West Pacific (WPO) and the Indian Ocean (IO). The SST-rainfall pattern simulated by the coupled versions of these models is rather similar to that from the corresponding atmospheric one, except for a shift of the entire pattern to colder/warmer SSTs when there is a cold/warm bias in the coupled version.

Asymptotic expansions for the coupled wavenumbers in an infinite orthotropic flexible fluid-filled cylindrical shell

Analytical expressions are found for the coupled wavenumbers in flexible, fluid-filled, circular cylindrical orthotropic shells using the asymptotic methods. These expressions are valid for arbitrary circumferential orders. The Donnell-Mushtari shell theory is used to model the shell and the effect of the fluid is introduced through the fluid-loading parameter mu. The orthotropic problem is posed as a perturbation on the corresponding isotropic problem by defining a suitable orthotropy parameter epsilon, which is a measure of the degree of orthotropy. For the first study, an isotropic shell is considered (by setting epsilon = 0) and expansions are found for the coupled wavenumbers using a regular perturbation approach. In the second study, asymptotic expansions are found for the coupled wavenumbers in the limit of small orthotropy (epsilon << 1). For each study, isotropy and orthotropy, expansions are found for small and large values of the fluid-loading parameter mu. All the asymptotic solutions are compared with numerical solutions to the coupled dispersion relation and the match is seen to be good. The differences between the isotropic and orthotropic solutions are discussed. The main contribution of this work lies in extending the existing literature beyond in vacuo studies to the case of fluid-filled shells (isotropic and orthotropic).

Coupled polynomial field approach for elimination of flexure and torsion locking phenomena in the Timoshenko and Euler-Bernoulli curved beam elements

The curvature related locking phenomena in the out-of-plane deformation of Timoshenko and Euler-Bernoulli curved beam elements are demonstrated and a novel approach is proposed to circumvent them. Both flexure and Torsion locking phenomena are noticed in Timoshenko beam and torsion locking phenomenon alone in Euler-Bernoulli beam. Two locking-free curved beam finite element models are developed using coupled polynomial displacement field interpolations to eliminate these locking effects. The coupled polynomial interpolation fields are derived independently for Timoshenko and Euler-Bernoulli beam elements using the governing equations. The presented of penalty terms in the couple displacement fields incorporates the flexure-torsion coupling and flexure-shear coupling effects in an approximate manner and produce no spurious constraints in the extreme geometric limits of flexure, torsion and shear stiffness. the proposed couple polynomial finite element models, as special cases, reduce to the conventional Timoshenko beam element and Euler-Bernoulli beam element, respectively. These models are shown to perform consistently over a wide range of flexure-to-shear (EI/GA) and flexure-to-torsion (EI/GJ) stiffness ratios and are inherently devoid of flexure, torsion and shear locking phenomena. The efficacy, accuracy and reliability of the proposed models to straight and curved beam applications are demonstrated through numerical examples. (C) 2012 Elsevier B.V. All rights reserved.