Stabilized SPH-based simulations of impact dynamics using acceleration-corrected artificial viscosity

Artificial viscosity in SPH-based computations of impact dynamics is a numerical artifice that helps stabilize spurious oscillations near the shock fronts and requires certain user-defined parameters. Improper choice of these parameters may lead to spurious entropy generation within the discretized system and make it over-dissipative. This is of particular concern in impact mechanics problems wherein the transient structural response may depend sensitively on the transfer of momentum and kinetic energy due to impact. In order to address this difficulty, an acceleration correction algorithm was proposed in Shaw and Reid (''Heuristic acceleration correction algorithm for use in SPH computations in impact mechanics'', Comput. Methods Appl. Mech. Engrg., 198, 3962-3974) and further rationalized in Shaw et al. (An Optimally Corrected Form of Acceleration Correction Algorithm within SPH-based Simulations of Solid Mechanics, submitted to Comput. Methods Appl. Mech. Engrg). It was shown that the acceleration correction algorithm removes spurious high frequency oscillations in the computed response whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. In this paper, we aim at gathering further insights into the acceleration correction algorithm by further exploring its application to problems related to impact dynamics. The numerical evidence in this work thus establishes that, together with the acceleration correction algorithm, SPH can be used as an accurate and efficient tool in dynamic, inelastic structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Microstructure and dry sliding wear resistance evaluation of plasma nitrided austenitic stainless steel type AISI 316LN against different sliders

In this work, Plasma Nitriding was carried out at a temperature of 570 degrees C on nuclear grade austenitic stainless steel type AISI 316 LN (316LN SS) in a gas mixture of 20% N-2-80% H-2 to improve the surface hardness and thereby sliding wear resistance. The Plasma Nitride (PN) treated surface has been characterized by Vickers microhardness measurements, Scanning Electron Microscopic (SEM) examination, X-ray Diffraction (XRD) and sliding wear assessment. The average thickness of the PN layer was found to be 70 mu m. Microhardness measurements showed a significant increase in the hardness from 210 HV25g (unnitrided sample) to 1040 HV25g (Plasma Nitrided sample). The XRD reveals that PN layer consists of CrN, Fe4N and Fe3N phases along with austenite phase. The tribological parameters such as the friction coefficient and wear mechanism have been evaluated at ambient conditions for PN treated ring (PN ring) vs. ASTM A453 grade 660 pin (ASTM pin), PN ring vs. Nickel based alloy hard faced pin (Colmonoy pin), PN ring vs. 316LN SS pin and 316LN SS ring vs. 316LN SS pin. The wear tracks have been analyzed by SEM, Energy Dispersive X-ray Analysis (EDX) and Optical Profilometry. The untreated 316LN SS ring vs. 316LN SS pin produced severe wear and was characterized by a combination of delamination and adhesion wear mechanism, whereas wear mechanism of the PN rings reveals mild abrasion and a transfer layer from pin materials. (C) 2012 Elsevier B.V. All rights reserved.

Quasiuniversal Transient Behavior of a Nonequilibrium Mott Insulator Driven by an Electric Field

We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field and that for a broad range of field strengths a long-lived transient prethermalized state emerges. This long-lived transient regime implies that thermal equilibrium may be out of reach of the time scales accessible in present cold atom experiments but shows that an interesting new quasiuniversal transient state exists in nonequilibrium governed by a thermalized kinetic energy but not a thermalized potential energy. In addition, when the field strength is equal in magnitude to the interaction between atoms, the system undergoes a rapid thermalization, characterized by a different quasiuniversal behavior of the current and spectral function for different values of the hopping. DOI: 10.1103/PhysRevLett.109.260402

Scaling of pressure spectrum in turbulent boundary layers

Scaling of pressure spectrum in zero-pressure-gradient turbulent boundary layers is discussed. Spatial DNS data of boundary layer at one time instant (Re-theta = 4500) are used for the analysis. It is observed that in the outer regions the pressure spectra tends towards the -7/3 law predicted by Kolmogorov's theory of small-scale turbulence. The slope in the pressure spectra varies from -1 close to the wall to a value close to -7/3 in the outer region. The streamwise velocity spectra also show a -5/3 trend in the outer region of the flow. The exercise carried out to study the amplitude modulation effect of the large scales on the smaller ones in the near-wall region reveals a strong modulation effect for the streamwise velocity, but not for the pressure fluctuations. The skewness of the pressure follows the same trend as the amplitude modulation coefficient, as is the case for the velocity. In the inner region, pressure spectra were seen to collapse better when normalized with the local Reynolds stress (-(u'v') over bar) than when scaled with the local turbulent kinetic energy (q(2) = (u'(2)) over bar + (v'(2)) over bar + (w'(2)) over bar)

Molecular Dynamics Study of Phonon Screening in Graphene

Phonon interaction with electrons or phonons or with structural defects result in a phonon mode conversion. The mode conversion is governed by the frequency wave-vector dispersion relation. The control over phonon mode or the screening of phonon in graphene is studied using the propagation of amplitude modulated phonon wave-packet. Control over phonon properties like frequency and velocity opens up several wave guiding, energy transport and thermo-electric applications of graphene. One way to achieve this control is with the introduction of nano-structured scattering in the phonon path. Atomistic model of thermal energy transport is developed which is applicable to devices consisting of source, channel and drain parts. Longitudinal acoustic phonon mode is excited from one end of the device. Molecular dynamics based time integration is adopted for the propagation of excited phonon to the other end of the device. The amount of energy transfer is estimated from the relative change of kinetic energy. Increase in the phonon frequency decreases the kinetic energy transmission linearly in the frequency band of interest. Further reduction in transmission is observed with the tuning of channel height of the device by increasing the boundary scattering. Phonon mode selective transmission control have potential application in thermal insulation or thermo-electric application or photo-thermal amplification.

Turbulence and cooling in galaxy cluster cores

We study the interplay between turbulent heating, mixing, and radiative cooling in an idealized model of cool cluster cores. Active galactic nuclei (AGN) jets are expected to drive turbulence and heat cluster cores. Cooling of the intracluster medium (ICM) and stirring by AGN jets are tightly coupled in a feedback loop. We impose the feedback loop by balancing radiative cooling with turbulent heating. In addition to heating the plasma, turbulence also mixes it, suppressing the formation of cold gas at small scales. In this regard, the effect of turbulence is analogous to thermal conduction. For uniform plasma in thermal balance (turbulent heating balancing radiative cooling), cold gas condenses only if the cooling time is shorter than the mixing time. This condition requires the turbulent kinetic energy to be a parts per thousand(3) the plasma internal energy; such high velocities in cool cores are ruled out by observations. The results with realistic magnetic fields and thermal conduction are qualitatively similar to the hydrodynamic simulations. Simulations where the runaway cooling of the cool core is prevented due to mixing with the hot ICM show cold gas even with subsonic turbulence, consistent with observations. Thus, turbulent mixing is the likely mechanism via which AGN jets heat cluster cores. The thermal instability growth rates observed in simulations with turbulence are consistent with the local thermal instability interpretation of cold gas in cluster cores.

Isotropic finite volume discretization

Finite volume methods traditionally employ dimension by dimension extension of the one-dimensional reconstruction and averaging procedures to achieve spatial discretization of the governing partial differential equations on a structured Cartesian mesh in multiple dimensions. This simple approach based on tensor product stencils introduces an undesirable grid orientation dependence in the computed solution. The resulting anisotropic errors lead to a disparity in the calculations that is most prominent between directions parallel and diagonal to the grid lines. In this work we develop isotropic finite volume discretization schemes which minimize such grid orientation effects in multidimensional calculations by eliminating the directional bias in the lowest order term in the truncation error. Explicit isotropic expressions that relate the cell face averaged line and surface integrals of a function and its derivatives to the given cell area and volume averages are derived in two and three dimensions, respectively. It is found that a family of isotropic approximations with a free parameter can be derived by combining isotropic schemes based on next-nearest and next-next-nearest neighbors in three dimensions. Use of these isotropic expressions alone in a standard finite volume framework, however, is found to be insufficient in enforcing rotational invariance when the flux vector is nonlinear and/or spatially non-uniform. The rotationally invariant terms which lead to a loss of isotropy in such cases are explicitly identified and recast in a differential form. Various forms of flux correction terms which allow for a full recovery of rotational invariance in the lowest order truncation error terms, while preserving the formal order of accuracy and discrete conservation of the original finite volume method, are developed. Numerical tests in two and three dimensions attest the superior directional attributes of the proposed isotropic finite volume method. Prominent anisotropic errors, such as spurious asymmetric distortions on a circular reaction-diffusion wave that feature in the conventional finite volume implementation are effectively suppressed through isotropic finite volume discretization. Furthermore, for a given spatial resolution, a striking improvement in the prediction of kinetic energy decay rate corresponding to a general two-dimensional incompressible flow field is observed with the use of an isotropic finite volume method instead of the conventional discretization. (C) 2014 Elsevier Inc. All rights reserved.

Effect of Rayleigh numbers on the evolution of double-diffusive salt fingers

This is a transient two-dimensional numerical study of double-diffusive salt fingers in a two-layer heat-salt system for a wide range of initial density stability ratio (R-rho 0) and thermal Rayleigh numbers (Ra-T similar to 10(3) - 10(11)). Salt fingers have been studied for several decades now, but several perplexing features of this rich and complex system remain unexplained. The work in question studies this problem and shows the morphological variation in fingers from low to high thermal Rayleigh numbers, which have been missed by the previous investigators. Considerable variations in convective structures and evolution pattern were observed in the range of Ra-T used in the simulation. Evolution of salt fingers was studied by monitoring the finger structures, kinetic energy, vertical profiles, velocity fields, and transient variation of R-rho(t). The results show that large scale convection that limits the finger length was observed only at high Rayleigh numbers. The transition from nonlinear to linear convection occurs at about Ra-T similar to 10(8). Contrary to the popular notion, R-rho(t) first decrease during diffusion before the onset time and then increase when convection begins at the interface. Decrease in R-rho(t) is substantial at low Ra-T and it decreases even below unity resulting in overturning of the system. Interestingly, all the finger system passes through the same state before the onset of convection irrespective of Rayleigh number and density stability ratio of the system. (C) 2014 AIP Publishing LLC.

Low frequency modulation of jets in quasigeostrophic turbulence

Quasigeostrophic turbulence on a beta-plane with a finite deformation radius is studied numerically, with particular emphasis on frequency and combined wavenumber-frequency domain analyses. Under suitable conditions, simulations with small-scale random forcing and large-scale drag exhibit a spontaneous formation of multiple zonal jets. The first hint of wave-like features is seen in the distribution of kinetic energy as a function of frequency; specifically, for progressively larger deformation scales, there are systematic departures in the form of isolated peaks (at progressively higher frequencies) from a power-law scaling. Concomitantly, there is an inverse flux of kinetic energy in frequency space which extends to lower frequencies for smaller deformation scales. The identification of these peaks as Rossby waves is made possible by examining the energy spectrum in frequency-zonal wavenumber and frequency-meridional wavenumber diagrams. In fact, the modified Rhines scale turns out to be a useful measure of the dominant meridional wavenumber of the modulating Rossby waves; once this is fixed, apart from a spectral peak at the origin (the steady jet), almost all the energy is contained in westward propagating disturbances that follow the theoretical Rossby dispersion relation. Quite consistently, noting that the zonal scale of the modulating waves is restricted to the first few wavenumbers, the energy spectrum is almost entirely contained within the corresponding Rossby dispersion curves on a frequency-meridional wavenumber diagram. Cases when jets do not form are also considered; once again, there is a hint of Rossby wave activity, though the spectral peaks are quite muted. Further, the kinetic energy scaling in frequency domain follows a -5/3 power-law and is distributed much more broadly in frequency-wavenumber diagrams. (C) 2015 AIP Publishing LLC.

Identification of dominant modes in random dynamical and aeroelastic systems

Identification of dominant modes is an important step in studying linearly vibrating systems, including flow-induced vibrations. In the presence of uncertainty, when some of the system parameters and the external excitation are modeled as random quantities, this step becomes more difficult. This work is aimed at giving a systematic treatment to this end. The ability to capture the time averaged kinetic energy is chosen as the primary criterion for selection of modes. Accordingly, a methodology is proposed based on the overlap of probability density functions (pdf) of the natural and excitation frequencies, proximity of the natural frequencies of the mean or baseline system, modal participation factor, and stochastic variation of mode shapes in terms of the modes of the baseline system - termed here as statistical modal overlapping. The probabilistic descriptors of the natural frequencies and mode shapes are found by solving a random eigenvalue problem. Three distinct vibration scenarios are considered: (i) undamped arid damped free vibrations of a bladed disk assembly, (ii) forced vibration of a building, and (iii) flutter of a bridge model. Through numerical studies, it is observed that the proposed methodology gives an accurate selection of modes. (C) 2015 Elsevier Ltd. All rights reserved.

Response Dynamics of Recirculation Structures in Coaxial Nonpremixed Swirl-Stabilized Flames Subjected to Acoustic Forcing

This paper reports the time-mean and phase-locked response of nonreacting as well as reacting flow field in a coaxial swirling jet/flame (nonpremixed). Two distinct swirl intensities plus two different central pipe flow rates at each swirl setting are investigated. The maximum response is observed at the 105 Hz mode in the range of excitation frequencies (0-315 Hz). The flow/flame exhibited minimal response beyond 300 Hz. It is seen that the aspect ratio change of inner recirculation zone (IRZ) under nonreacting conditions (at responsive modes) manifests as a corresponding increase in the time-mean flame aspect ratio. This is corroborated by similar to 25% decrease in the IRZ transverse width in both flame and cold flow states. In addition, 105 Hz excited states are found to shed high energy regions (eddies) asymmetrically when compared to dormant 315 Hz pulsing frequency. The kinetic energy (KE) of the flow field is subsequently reduced due to acoustic excitation and a corresponding increase (similar to O (1)) in fluctuation intensity is witnessed. The lower swirl intensity case is found to be more responsive than the high swirl case as in the former flow state the resistance offered by IRZ to incoming acoustic perturbations is lower due to inherently low inertia. Next, the phase-locked analysis of flow and flame structure is employed to further investigate the phase dependence of flow/flame response. It is found that the asymmetric shifting of IRZ mainly results at 270 deg acoustic forcing. The 90 deg phase angle forcing is observed to convect the IRZ farther downstream in both swirl cases as compared to other phase angles. The present work aims primarily at providing a fluid dynamic view point to the observed nonpremixed flame response without considering the confinement effects.

Deformation pathways and breakup modes in acoustically levitated bicomponent droplets under external heating

Controlled breakup of droplets using heat or acoustics is pivotal in applications such as pharmaceutics, nanoparticle production, and combustion. In the current work we have identified distinct thermal acoustics-induced deformation regimes (ligaments and bubbles) and breakup dynamics in externally heated acoustically levitated bicomponent (benzene-dodecane) droplets with a wide variation in volatility of the two components (benzene is significantly more volatile than dodecane). We showcase the physical mechanism and universal behavior of droplet surface caving in leading to the inception and growth of ligaments. The caving of the top surface is governed by a balance between the acoustic pressure field and the restrictive surface tension of the droplet. The universal collapse of caving profiles for different benzene concentration (<70% by volume) is shown by using an appropriate time scale obtained from force balance. Continuous caving leads to the formation of a liquid membrane-type structure which undergoes radial extension due to inertia gained during the precursor phase. The membrane subsequently closes at the rim and the kinetic energy leads to ligament formation and growth. Subsequent ligament breakup is primarily Rayleigh-Plateau type. The breakup mode shifts to diffusional entrapment-induced boiling with an increase in concentration of the volatile component (benzene >70% by volume). The findings are portable to any similar bicomponent systems with differential volatility.

Quantitation and characterization of glutathionyl haemoglobin as an oxidative stress marker in chronic renal failure by mass spectrometry

Objectives: Glutathionyl haemoglobin (GS-Hb) belonging to the class of glutathionylated proteins has been investigated as a possible marker of oxidative stress in different chronic diseases. The purpose of this study was to examine whether glutathionyl haemoglobin can serve as an oxidative stress marker in non-diabetic chronic renal failure patients on different renal replacement therapies (RRT) through its quantitation, and characterization of the specific binding site of glutathione in haemoglobin molecule by mass spectrometric analysis. Design and methods: The study group consisted of non-diabetic chronic renal failure patients on renal replacement therapy (RRT): hemodialysis (HD), continuous ambulatory peritoneal dialysis (CAPD) and renal allograft transplant (Txp) patients. Haemoglobin samples of these subjects were analyzed by liquid chromatography electrospray ionization mass spectrometry for GS-Hb quantitation. Characterization of GS-Hb was done by tandem mass spectrometry. Levels of erythrocyte glutathione (GSH) and lipid peroxidation (as thiobarbituric acid reacting substances) were measured spectrophotometrically, while glycated baernoglobin (HbA1c) was measured by HPLC. Results: GS-Hb levels were markedly elevated in the dialysis group and marginally in the transplant group as compared to the controls. GS-Hb levels correlated positively with lipid peroxidation and negatively with the erythrocyte glutathione levels in RRT groups indicating enhanced oxidative stress. De novo sequencing of the chymotryptic fragment of GS-Hb established that glutathione is attached to Cys-93 of the beta globin chain. Mass spectrometric quantitation of total glycated haemoglobin showed good agreement with HbA1c estimation by conventional HPLC method. Conclusions: Glutathionyl haemoglobin can serve as a clinical marker of oxidative stress in chronic debilitating therapies like RRT. Mass spectrometry provides a reliable analytical tool for quantitation and residue level characterization of different post-translational modifications of haemoglobin. (c) 2007 The Canadian Society of Clinical Chemists. Published by Elsevier Inc. All rights reserved.

Detection of alpha(2u)-globulin and its bound putative pheromones in the preputial gland of the Indian commensal rat (Rattus rattus) using mass spectrometry

The role of pheromones and pheromone-binding proteins in the laboratory rat has been extensively investigated. However, we have previously reported that the preputial gland of the Indian commensal rat produces a variety of pheromonal molecules and preputial glands would seem to be the predominant source for pheromonal communication. The presence of pheromone-binding proteins has not yet been identified in the preputial gland of the Indian commensal rat; therefore, the experiments were designed to unravel the alpha(2u)-globulin (alpha 2u) and its bound volatiles in the commensal rat. Total preputial glandular proteins were first fractionated by sodium dodecyl sulfate/polyacrylamide gel electrophoresis (SDS-PAGE) and subsequently analyzed by mass spectrometry. Further, we purified alpha 2u and screened for the presence of bound pheromonal molecules with the aid of gas chromatography/mass spectrometry (GC/MS). A novel alpha 2u was identified with a high score and this protein has not been previously described as present in the preputial gland of Indian commensal rats.This novel alpha 2u was then characterized by tandem mass spectrometry (MS/MS). Peptides with m/z values of 969, 1192, 1303 and 1876 were further fragmented with the aid of MS/MS and generated de novo sequences which provided additional evidence for the presence of alpha 2u in the preputial gland. Finally, we identified the presence of farnesol 1 and 2 bound to alpha 2u. The present investigation confirms the presence of alpha 2u (18.54 kDa) in the preputial gland of the Indian commensal rat and identifies farnesol 1 and 2 as probably involved in chemo-communication by the Indian commensal rat.Copyright (C) 2010 John Wiley & Sons, Ltd.

Sequencing of T-superfamily conotoxins from Conus virgo: Pyroglutamic acid identification and disulfide arrangement by MALDI mass spectrometry

De novo mass spectrometric sequencing of two Conus peptides, Vi1359 and Vi1361, from the vermivorous cone snail Conus virgo, found off the southern Indian coast, is presented. The peptides, whose masses differ only by 2 Da, possess two disulfide bonds and an amidated C-terminus. Simple chemical modifications and enzymatic cleavage coupled with matrix assisted laser desorption ionization (MALDI) mass spectrometric analysis aided in establishing the sequences of Vi1359, ZCCITIPECCRI-NH2, and Vi1361, ZCCPTMPECCRI-NH2, Which differ only at residues 4 and 6 (Z = pyroglutamic acid). The presence of the pyroglutamyl residue at the N-terminus was unambiguously identified by chemical hydrolysis of the cyclic amide, followed by esterification. The presence of Ile residues in both the peptides was confirmed from high-energy collision induced dissociation (CID) studies, using the observation Of W-n- and d(n)-ions as a diagnostic. Differential cysteine labeling, in conjunction with MALDI-MS/MS, permitted establishment of disulfide connectivity in both peptides as Cys2-Cys9 and Cys3-Cys10. The cysteine pattern clearly reveals that the peptides belong to the class of T-superfamily conotoxins, in particular the T-1 superfamily.