Radiatively inefficient accretion flow simulations with cooling: implications for black hole transients

We study the effects of optically thin radiative cooling on the structure of radiatively inefficient accretion flows (RIAFs). The flow structure is geometrically thick, and independent of the gas density and cooling, if the cooling time is longer than the viscous time-scale (i.e. t(cool) greater than or similar to t(visc)). For higher densities, the gas can cool before it can accrete and forms the standard geometrically thin, optically thick Shakura-Sunyaev disc. For usual cooling processes (such as bremsstrahlung), we expect an inner hot flow and an outer thin disc. For a short cooling time the accretion flow separates into two phases: a radiatively inefficient hot coronal phase and a cold thin disc. We argue that there is an upper limit on the density of the hot corona corresponding to a critical value of t(cool)/t(ff)( similar to 10-100), the ratio of the cooling time and the free-fall time. Based on our simulations, we have developed a model for transients observed in black hole X-ray binaries (XRBs). An XRB in a quiescent hot RIAF state can transition to a cold blackbody-dominated state because of an increase in the mass accretion rate. The transition from a thin disc to a RIAF happens because of mass exhaustion due to accretion; the transition happens when the cooling time becomes longer than the viscous time at inner radii. Since the viscous time-scale for a geometrically thin disc is quite long, the high-soft state is expected to be long-lived. The different time-scales in black hole transients correspond to different physical processes such as viscous evolution, cooling and free fall. Our model captures the overall features of observed state transitions in XRBs.

Numerical simulations of bubble formation from submerged needles under non-uniform direct current electric field

In several chemical and space industries, small bubbles are desired for efficient interaction between the liquid and gas phases. In the present study, we show that non-uniform electric field with appropriate electrode configurations can reduce the volume of the bubbles forming at submerged needles by up to three orders of magnitude. We show that localized high electric stresses at the base of the bubbles result in slipping of the contact line on the inner surface of the needle and subsequent bubble formation occurs with contact line inside the needle. We also show that for bubble formation in the presence of highly non-uniform electric field, due to high detachment frequency, the bubbles go through multiple coalescences and thus increase the apparent volume of the detached bubbles. (C) 2013 AIP Publishing LLC.

Effect of SST Variation on ITCZ in APE Simulations

The effect of meridional variation of sea surface temperature (SST) on tropical atmospheric circulation is analyzed using Aqua-planet Experiment (APE) simulations. The meridional SST gradient around the narrow SST peak in CONTROL simulation favours a strong and single equatorial Intertropical Convergence Zone (ITCZ, defined by the maximum of zonally averaged total precipitation) in all APE models. In contrast, flat equatorial SST peak (FLAT simulation) favours split/double ITCZs flanking the SST maximum, in the majority of the APE models. Although there is reasonable agreement for SST sensitivity of ITCZ among the APE models in CONTROL, there exists disparity among them in FLAT case. Similarly, while the total and convective precipitation responses are consistent among the models, the large-scale precipitation response shows considerable inter-model variations in FLAT case. The APE intercomparison indicates that the occurrence and positioning of the ITCZ are primarily related to boundary layer moisture convergence as a response to the meridional variation of SST. Furthermore, the meridional gradient of tropospheric temperature is found to be an important factor that can influence the positioning of ITCZ. FLAT SST distribution is found to be similar to the observed distribution over the Indian region during summer season. Models that yield double ITCZs in this case simulate an easterly jet over the equatorial region (similar to 15 degrees equatorward of the ITCZ). This is analogous to the Tropical Easterly Jet (TEJ), which is a unique feature observed over the Indian region during summer monsoon season, with its core at 12 degrees N, equatorward of the seasonal convergence zone centered along 25 degrees N. In these models, positive meridional temperature gradient and the associated easterly shear in the atmosphere strengthened by moisture convergence penetrate up to the upper troposphere, with which TEJ is in thermal wind balance.

Molecular dynamic simulations of elastomer structure and its influence on anisotropic stress under time-varying strain

We investigate the evolution of polymer structure and its influence on uniaxial anisotropic stress under time-varying uniaxial strain, and the role of external control variables such as temperature, strain rate, chain length, and density, using molecular dynamics simulation. At temperatures higher than glass transition, stress anisotropy in the system is reduced even though the bond stretch is greater at higher temperatures. There is a significant increase in the stress level with increasing density. At higher densities, the uncoiling of the chains is suppressed and the major contribution to the deformation is by internal deformation of the chains. At faster rates of loading stress anisotropy increases. The deformation mechanism is mostly due to bond stretch and bond bending rather than overall shape and size. Stress levels increase with longer chain length. There is a critical value of the functionality of the cross-linkers beyond which the uniaxial stress developed increases caused primarily by bond stretching due to increased constraint on the motion of the monomers. Stacking of the chains in the system also plays a dominant role in the behaviour in terms of excluded volume interactions. Low density, high temperature, low values of functionality of cross-linkers, and short chain length facilitate chain uncoiling and chain slipping in cross-linked polymers.

Calculation of Chemical Potentials and Occupancies in Clathrate Hydrates through Monte Carlo Molecular Simulations

The flexibility of the water lattice in clathrate hydrates and guest-guest interactions has been shown in previous studies to significantly affect the values of the thermodynamic properties, such as chemical potentials and free energies. Here we describe methods for computing occupancies, chemical potentials, and free energies that account for the flexibility of water lattice and guest-guest interactions in the hydrate phase. The methods are validated for a wide variety of guest molecules, such as methane, ethane, carbon dioxide, and tetrahydrodfuran by comparing the predicted occupancy values of guest molecules with those obtained from isothermal isobaric semigrand Monte Carlo simulations. The proposed methods extend the van der Waals and Platteuw theory for clathrate hydrates, and the Langmuir constant is calculated based on the structure of the empty hydrate lattice. These methods in combination with development of advanced molecular models for water and guest molecules should lead to a more thermodynamically consistent theory for clathrate hydrates.

Simulations and experiments in punching spring-steel devices with sub-millimeter features

This work demonstrates the feasibility of mesoscale (100 μm to mm) punching of multiple holes of intricate shapes in metals. Analytical modeling, finite element (FE)simulation, and experimentations are used in this work. Two dimensional FE simulations in ABAQUS were done with an assumed material modeling and plane-strain condition. A known analytical model was used and compared with the ABAQUS simulation results to understand the effects of clearance between the punch and the die. FE simulation in ABAQUS was done for different clearances and corner radii at punch, die, and holder. A set of punches and dies were used to punch out a miniature spring-steel gripper. Comparison of compliant grippers manufactured by wire-cut electro discharge machining(EDM) and punching shows that realizing sharp interior and re-entrant corners by punching is not easy to achieve. Punching of circular holes with 5 mm and 2.5 mm diameter is achieved. The possibility of realizing meso-scale parts with complicated shapes through punching is demonstrated in this work; and some strategies are suggested for improvement.

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved.

Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of {10 (1) over bar2}-type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load-displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin-twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.

Synthesis, DNA binding, docking and photocleavage studies of quinolinyl chalcones

A series of simple quinoline-chalcone conjugates have been synthesized by Claisen-Schmidt condensation reactions of substituted acetophenones with 2-chloro-3-formyl-quinoline and evaluated for their nucleolytic activity. The structures of the synthesized quinoline-chalcone conjugates were confirmed by IR, H-1 NMR, C-13 NMR and mass spectral analyses. Most of the prepared compounds showed significant DNA binding and photocleavage activities. The incorporation of an electron-donating group into ring A caused a moderate increase in the DNA binding and photocleavage activities. Compounds 3c and 3d exhibited promising DNA photocleavage against pUC 19 DNA with 85% inhibition at 100 mu M concentration. A structure-activity relationship analysis of these compounds was performed; compounds 3c and 3d are potential candidates for future drug discovery and development.

Riemann shock tube: 1D normal shocks in air, simulations and experiments

This paper presents numerical simulation of the evolution of one-dimensional normal shocks, their propagation, reflection and interaction in air using a single diaphragm Riemann shock tube and validate them using experimental results. Mathematical model is derived for one-dimensional compressible flow of viscous and conducting medium. Dimensionless form of the mathematical model is used to construct space-time finite element processes based on minimization of the space-time residual functional. The space-time local approximation functions for space-time p-version hierarchical finite elements are considered in higher order GRAPHICS] spaces that permit desired order of global differentiability of local approximations in space and time. The resulting algebraic systems from this approach yield unconditionally positive-definite coefficient matrices, hence ensure unique numerical solution. The evolution is computed for a space-time strip corresponding to a time increment Delta t and then time march to obtain the evolution up to any desired value of time. Numerical studies are designed using recently invented hand-driven shock tube (Reddy tube) parameters, high/low side density and pressure values, high- and low-pressure side shock tube lengths, so that numerically computed results can be compared with actual experimental measurements.

Synthesis, molecular docking and anti-mycobacterial evaluation of new imidazo1,2-a]pyridine-2-carboxamide derivatives

New anti-tubercular agents, imidazo1,2-a]pyridine-2-carboxamide derivatives (5a-q) have been designed and synthesized. The structural considerations of the designed molecules were further supported by the docking study with a long-chain enoyl-acyl carrier protein reductase (InhA). The chemical structures of the new compounds were characterized by IR, H-1 NMR, C-13 NMR, HRMS and elemental analysis. In addition, single crystal X-ray diffraction has also been recorded for compound 5f. Compounds were evaluated in vitro against Mycobacterium tuberculosis H37Rv, and cytotoxicity against HEK-293T cell line. Amongst the tested compounds 5j, 5l and 5q were emerged as good anti-tubercular agents with low cytotoxicity. The structure-anti TB activity relationship of these derivatives was explained by molecular docking. (C) 2014 Elsevier Masson SAS. All rights reserved.

NUMERICAL SIMULATIONS OF LIQUID JET BREAK UP IN A CROSSFLOW

Atomization is the process of disintegration of a liquid jet into ligaments and subsequently into smaller droplets. A liquid jet injected from a circular orifice into cross flow of air undergoes atomization primarily due to the interaction of the two phases rather than an intrinsic break up. Direct numerical simulation of this process resolving the finest droplets is computationally very expensive and impractical. In the present study, we resort to multiscale modelling to reduce the computational cost. The primary break up of the liquid jet is simulated using Gerris, an open source code, which employs Volume-of-Fluid (VOF) algorithm. The smallest droplets formed during primary atomization are modeled as Lagrangian particles. This one-way coupling approach is validated with the help of the simple test case of tracking a particle in a Taylor-Green vortex. The temporal evolution of the liquid jet forming the spray is captured and the flattening of the cylindrical liquid column prior to breakup is observed. The size distribution of the resultant droplets is presented at different distances downstream from the location of injection and their spatial evolution is analyzed.

A stabilization approach for mesh-free simulations of systems developing shocks or extreme strain localizations

A new stabilization scheme, based on a stochastic representation of the discretized field variables, is proposed with a view to reduce or even eliminate unphysical oscillations in the mesh-free numerical simulations of systems developing shocks or exhibiting localized bands of extreme deformation in the response. The origin of the stabilization scheme may be traced to nonlinear stochastic filtering and, consistent with a class of such filters, gain-based additive correction terms are applied to the simulated solution of the system, herein achieved through the element-free Galerkin method, in order to impose a set of constraints that help arresting the spurious oscillations. The method is numerically illustrated through its Applications to inviscid Burgers' equations, wherein shocks may develop as a result of intersections of the characteristics, and to a gradient plasticity model whose response is often characterized by a developing shear band as the external load is gradually increased. The potential of the method in stabilized yet accurate numerical simulations of such systems involving extreme gradient variations in the response is thus brought forth. (C) 2014 Elsevier Ltd. All rights reserved.

Modeling of solar radiation management: a comparison of simulations using reduced solar constant and stratospheric sulphate aerosols

The climatic effects of Solar Radiation Management (SRM) geoengineering have been often modeled by simply reducing the solar constant. This is most likely valid only for space sunshades and not for atmosphere and surface based SRM methods. In this study, a global climate model is used to evaluate the differences in the climate response to SRM by uniform solar constant reduction and stratospheric aerosols. Our analysis shows that when global mean warming from a doubling of CO2 is nearly cancelled by both these methods, they are similar when important surface and tropospheric climate variables are considered. However, a difference of 1 K in the global mean stratospheric (61-9.8 hPa) temperature is simulated between the two SRM methods. Further, while the global mean surface diffuse radiation increases by similar to 23 % and direct radiation decreases by about 9 % in the case of sulphate aerosol SRM method, both direct and diffuse radiation decrease by similar fractional amounts (similar to 1.0 %) when solar constant is reduced. When CO2 fertilization effects from elevated CO2 concentration levels are removed, the contribution from shaded leaves to gross primary productivity (GPP) increases by 1.8 % in aerosol SRM because of increased diffuse light. However, this increase is almost offset by a 15.2 % decline in sunlit contribution due to reduced direct light. Overall both the SRM simulations show similar decrease in GPP (similar to 8 %) and net primary productivity (similar to 3 %). Based on our results we conclude that the climate states produced by a reduction in solar constant and addition of aerosols into the stratosphere can be considered almost similar except for two important aspects: stratospheric temperature change and the consequent implications for the dynamics and the chemistry of the stratosphere and the partitioning of direct versus diffuse radiation reaching the surface. Further, the likely dependence of global hydrological cycle response on aerosol particle size and the latitudinal and height distribution of aerosols is discussed.

Simulations using meshfree methods

In this paper, attempt is made to solve a few problems using the Polynomial Point Collocation Method (PPCM), the Radial Point Collocation Method (RPCM), Smoothed Particle Hydrodynamics (SPH), and the Finite Point Method (FPM). A few observations on the accuracy of these methods are recorded. All the simulations in this paper are three dimensional linear elastostatic simulations, without accounting for body forces.