Source author record

Baidurya Bhattacharya

Baidurya Bhattacharya appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2016arXiv

Heat pump without particle transport or external work on the medium achieved by differential thermostatting of the phase space

We propose a new mechanism that enables heat flow from a colder region to a hotter region without necessitating either particle transport or external work on the conductor, thereby bypassing the compressor part of a classical heat pump cycle. Our mechanism relies on thermostatting the kinetic and configurational temperatures of the same particle differently. We keep the two ends of a conductor, which in the present study is a single dimensional $ϕ^4$ chain, at the same kinetic temperature $T_0$, but at different configurational temperatures - one end hotter and the other end colder than $T_0$. While external energy is needed within the thermostatted regions to achieve this differential thermostatting, no external work is performed on the system itself. We show that the mechanism satisfies the statistical form of the second law of thermodynamics (the fluctuation theorem). The proposed mechanism reveals two interesting findings - (i) contrary to traditional thermodynamics where only the kinetic temperature is thought to govern heat conduction, configurational temperature can also play an important role, and (ii) the relative temperature difference between the kinetic and configurational variables governs the direction of heat flow. The challenge, however, is in developing experimental techniques to thermostat the kinetic and configurational variables of the same particle at different values.

preprint2015arXiv

Effect of randomly occurring Stone-Wales defects on mechanical properties of carbon nanotubes using atomistic simulation

While CNTs are found to have ultra high stiffness and strength, an enormous scatter is also observed in available laboratory results. This paper studies the effects of randomly distributed Stone Wales (SW or 5 7 7 5) defects on the mechanical properties of single walled nanotubes (SWNTs) using the technique of atomistic simulation (AS). A Matern hard core random field applied on a finite cylindrical surface is used to describe the spatial distribution of the Stone Wales defects. We simulate a set of displacement controlled tensile loading up to fracture of SWNTs with (6,6) armchair and (10, 0) zigzag configurations and aspect ratio around 6. A modified Morse potential is adopted to model the interatomic forces. We found that fracture invariably initiates from a defect if one is present; for a defect-free tube the crack initiates at quite random locations. The force-displacement curve typically behaves almost linearly up to about half way, although there is no obvious yield point. Three mechanical properties - stiffness, ultimate strength and ultimate strain - are calculated from the simulated force and displacement time histories. The randomness in mechanical behavior resulting only from initial velocity distribution was found to be insignificant at room temperature. The mean values of stiffness, ultimate strength and ultimate strain of the tube decrease as the average number of defects increases although the coefficients of variation do not show such monotonic trend. The introduction of an additional defect has the most pronounced effect on the randomness in mechanical properties when the tube is originally defect free. We also found that, for a given mean number of defects in the tube, the zigzag configuration has less strength and less ultimate strain on the average, but more uncertainty in its stiffness and ultimate strain, compared with the armchair tube.

preprint2015arXiv

Fracture resistance of zigzag single walled carbon nanotubes

Brittle fracture is one of the important failure modes of Single-Walled Carbon Nanotube (SWNT) due to mechanical loading. In this paper, the fracture resistance of zigzag SWNTs with preexisting defects is calculated using fracture mechanics concepts based on atomistic simulations. The problem of unstable crack growth at finite temperature, presumably caused by lattice trapping effect, is circumvented by computing the strain energy release rate through a series of displacement-controlled tensile loading of SWNTs (applied through moving the outermost layer of atoms at one end at constant strain rate of 9.4x10-4/ps) with pre-existing crack-like defects of various lengths. The strain energy release rate, G, is computed for (17,0), (28,0) and (35,0) SWNTs (each with aspect ratio 4) with pre-existing cracks up to 29.5Å long. The fracture resistance, Gc, is determined as a function of crack length for each tube at three different temperatures (1K, 300K and 500K). A significant dependence of Gc on crack length is observed reminiscent of the rising R curve behavior of metals at the macroscale: for the zigzag nanotubes Gc increases with crack length at small length, and tends to reach a constant value if the tube diameter is large enough. We suspect that the lattice trapping effect plays the role of crack tip plasticity at the atomic scale. For example, at 300 Kelvin, Gc for the (35,0) tube with aspect ratio 4 converges to 6 Joule/m2 as the crack length exceeds 20 Angstrom. This value is comparable with the fracture toughness of graphite and Silicon. The fracture resistance of the tubes is found to decrease significantly as the temperature increases. To study the length effects, the computations are repeated for zigzag nanotubes with the same three chiralities but with aspect ratio 8 at 1K.

preprint2015arXiv

Performance metrics in a hybrid MPI-OpenMP based molecular dynamics simulation with short-range interactions

We discuss the computational bottlenecks in molecular dynamics (MD) and describe the challenges in parallelizing the computation intensive tasks. We present a hybrid algorithm using MPI (Message Passing Interface) with OpenMP threads for parallelizing a generalized MD computation scheme for systems with short range interatomic interactions. The algorithm is discussed in the context of nanoindentation of Chromium films with carbon indenters using the Embedded Atom Method potential for Cr Cr interaction and the Morse potential for Cr C interactions. We study the performance of our algorithm for a range of MPIthread combinations and find the performance to depend strongly on the computational task and load sharing in the multicore processor. The algorithm scaled poorly with MPI and our hybrid schemes were observed to outperform the pure message passing scheme, despite utilizing the same number of processors or cores in the cluster. Speed-up achieved by our algorithm compared favourably with that achieved by standard MD packages.

preprint2015arXiv

The Extremal Index and the Maximum of a Dependent Stationary Pulse Load Process Observed above a High Threshold

Observing a load process above high thresholds, modeling it as a pulse process with random occurrence times and magnitudes, and extrapolating life-time maximum or design loads from the data is a common task in structural reliability analyses. In this paper, we consider a stationary live load sequence that arrive according to a dependent point process and allow for a weakened mixing-type dependence in the load pulse magnitudes that asymptotically decreases to zero with increasing separation in the sequence. Inclusion of dependence in the model eliminates the unnecessary conservatism introduced by the i.i.d. (independent and identically distributed) assumption often made in determining maximum live load distribution. The scale of fluctuation of the loading process is used to identify clusters of exceedances above high thresholds which in turn is used to estimate the extremal index of the process. A Bayesian updating of the empirical distribution function, derived from the distribution of order statistics in a dependent stationary series, is performed. The pulse arrival instants are modeled as a Cox process goverened by a stationary lognormal intensity. An illustrative example utilizes in-service peak strain data from ambient traffic collected on a high volume highway bridge, and analyzes the asymptotic behavior of the maximum load.

preprint2014arXiv

Approximating the entire spectrum of nonequilibrium steady state distributions using relative entropy: An application to thermal conduction

We show that distribution functions of nonequilibrium steady states (NESS) evolving under a slowly varying protocol can be accurately obtained from limited data and the closest known detailed state of the system. In this manner, one needs to perform only a few detailed experiments to obtain the nonequilibrium distribution function for the entire gamut of nonlinearity. We achieve this by maximizing the relative entropy functional (MaxRent), which is proportional to the Kullback-Leibler distance from a known density function, subject to constraints supplied by the problem definition and new measurements. MaxRent is thus superior to the principle of maximum entropy (MaxEnt), which maximizes Shannon's informational entropy for estimating distributions but lacks the ability of incorporating additional prior information. The MaxRent principle is illustrated using a toy model of $ϕ^4$ thermal conduction consisting of a single lattice point. An external protocol controlled position-dependent temperature field drives the system out of equilibrium. Two different thermostatting schemes are employed: the Hoover-Holian deterministic thermostat (which produces multifractal dynamics under strong nonlinearity) and the Langevin stochastic thermostat (which produces phase space-filling dynamics). Out of the 80 possible states produced by the protocol, we assume that 4 states are known to us in detail, one of which is used as input into MaxRent at a time. We find that MaxRent accurately approximates the phase space density functions at all values of the protocol even when the known distribution is far away. MaxEnt, however, is unable to capture the fine details of the phase space distribution functions. We expect this method to be useful in other external protocol driven nonequilibrium cases as well.

preprint2014arXiv

Improving the ergodic characteristics of thermostats using higher order temperatures

Most deterministic schemes for controlling temperature by kinetic variables (NH thermsotat), configurational variables (BT thermostat) and all phase space variables (PB thermostat) are non-ergodic for systems with a few degrees of freedom. While for the NH thermostat ergodicity has been achieved by controlling the higher order moments of kinetic energy, the issues of nonergodicity of BT and PB thermostats still persist. In this paper, we propose a family of modifications for improving their ergodic characteristics. To do so, we introduce two new measures of configurational temperature (second and third order) based on the generalized temperature - curvature relationship. The equations of motion for the existing thermostats are modified by controlling the relevant higher order temperature variables. The ergodic characteristics of the proposed modifications are tested using a single harmonic oscillator. For the PB thermostat, the fastest route to ergodicity is by controlling the usual configurational and kinetic temperatures along with the $2^{\text{nd}}$ order configurational temperature. For the BT thermostat, controlling the usual configurational temperature along with the $2^{\text{nd}}$ order configurational temperature is sufficient to make it ergodic. Our method also provides a new ergodic generalization of NH thermostat that reduces to kinetic moments method for a single harmonic oscillator.

preprint2014arXiv

Non-Ergodicity of Nose-Hoover chain thermostat in computationally achievable time

The widely used Nose-Hoover chain (NHC) thermostat in molecular dynamics simulations is generally believed to impart the canonical distribution as well as quasi- (i.e., space filling) ergodicity on the thermostatted physical system (PS). Working with the standard single harmonic oscillator, we prove analytically that the two chain Nose-Hoover thermostat with unequal thermostat masses approach the standard Nose-Hoover dynamics and hence the PS loses its canonical and quasi-ergodic nature. We also show through numerical simulations over substantially long times that for certain Poincare sections, for both the equal and unequal thermostat mass cases, the bivariate distribution function of position and momentum (${x,p}$) and of reservoir degrees of freedom ($ξ,η$) lose their Gaussian nature. Further, the 4-dimensional $x-p-ξ-η$ extended phase space exhibits two holes of non-zero measure. The NHC thermostat therefore does not generate the canonical distribution or preserve quasi-ergodicity for the PS.