Source author record

Krishna Kumar

Krishna Kumar 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

32works
17topics
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

32 published item(s)

preprint2023arXiv

Quantification of geogrid lateral restraint using transparent sand and deep learning-based image segmentation

An experimental technique is presented to quantify the lateral restraint provided by a geogrid embedded in granular soil at the particle level. Repeated load triaxial tests were done on transparent sand specimens with geosynthetic inclusions simulating geogrids. Particle outlines on laser illuminated planes through the specimens were segmented using a deep learning-based segmentation algorithm. The particle outlines were characterized in terms of Fourier shape descriptors and tracked across sequentially captured images. The accuracy of the particle displacement measurements was validated against Digital Image Correlation (DIC) measurements. In addition, the method's resolution and repeatability is presented. Based on the measured particle displacements and rotations, a state boundary line between probable and improbable particle motions was identified for each test. The size of the zone of probable motions could be used to quantify the lateral restraint provided by the inclusions. Overall, the tests results revealed that the geosynthetic inclusions restricted both particle displacements and rotations. However, the particle displacements were found to be restrained more significantly than the rotations. Finally, a unique relationship was found between the magnitude of the permanent strains of the specimens and the size of the zone of probable motions.

preprint2022arXiv

A Frequency-Velocity CNN for Developing Near-Surface 2D Vs Images from Linear-Array, Active-Source Wavefield Measurements

This paper presents a frequency-velocity convolutional neural network (CNN) for rapid, non-invasive 2D shear wave velocity (Vs) imaging of near-surface geo-materials. Operating in the frequency-velocity domain allows for significant flexibility in the linear-array, active-source experimental testing configurations used for generating the CNN input, which are normalized dispersion images. Unlike wavefield images, normalized dispersion images are relatively insensitive to the experimental testing configuration, accommodating various source types, source offsets, numbers of receivers, and receiver spacings. We demonstrate the effectiveness of the frequency-velocity CNN by applying it to a classic near-surface geophysics problem, namely, imaging a two-layer, undulating, soil-over-bedrock interface. This problem was recently investigated in our group by developing a time-distance CNN, which showed great promise but lacked flexibility in utilizing different field-testing configurations. Herein, the new frequency-velocity CNN is shown to have comparable accuracy to the time-distance CNN while providing greater flexibility to handle varied field applications. The frequency-velocity CNN was trained, validated, and tested using 100,000 synthetic near-surface models. The ability of the proposed frequency-velocity CNN to generalize across various acquisition configurations is first tested using synthetic near-surface models with different acquisition configurations from that of the training set, and then applied to experimental field data collected at the Hornsby Bend site in Austin, Texas, USA. When fully developed for a wider range of geological conditions, the proposed CNN may ultimately be used as a rapid, end-to-end alternative for current pseudo-2D surface wave imaging techniques or to develop starting models for full waveform inversion.

preprint2022arXiv

A machine learning approach to predicting pore pressure response in liquefiable sands under cyclic loading

Shear stress history controls the pore pressure response in liquefiable soils. The excess pore pressure does not increase under cyclic loading when shear stress amplitude is lower than the peak prior amplitude -- the shielding effect. Many sophisticated constitutive models fail to capture the shielding effect observed in the cyclic liquefaction experiments. We develop a data-driven machine learning model based on the LSTM neural network to capture the liquefaction response of soils under cyclic loading. The LSTM model is trained on 12 laboratory cyclic simple shear tests on Nevada sand in loose and dense conditions subjected to different cyclic simple shear loading conditions. The LSTM model features include the relative density of soil and the previous stress history to predict the pore water pressure response. The LSTM model successfully replicates the pore pressure response for three cyclic simple test results considering the shielding and density effects.

preprint2022arXiv

Effect of color reconnection and hadronic re-scattering on underlying events in pp collisions at LHC energies

Underlying events dominate most of the hadronic activity in p$-$p collisions and are spanned from perturbative to non-perturbative QCD, having a sensitivity ranging from the multi-scale to very low-x scale physics. A detailed understanding of such events plays a crucial role in the accurate understanding of Standard Model ()SM and Beyond Standard Model physics. The underlying event activities has been studied within the framework of Pythia 8 Monte Carlo model, considering the underlying events observables mean charged particle multiplicity density , $\langle d^{2}N /dηdϕ\rangle$ and mean scalar $p_T$ sum, $\langle d^{2} \sum p_{T} /dηdϕ\rangle$ as a function of leading charged particle in towards, away, and transverse region of p$-$p collisions at $\sqrt{s}$ = 2.76, 7 and 13 TeV. The towards, away, and transverse regions have been defined on an azimuthal plane relative to leading particle in p$-$p collisions. The energy dependence of underlying events and their activities in the central and forward region has also been studied. The effect of hadronic re-scattering, color reconnection, and rope hadronization mechanism implemented in Pythia 8 has been studied in details to gain insight into the different processes contributing to underlying events in soft sector.

preprint2022arXiv

Hybrid Finite Element and Material Point Method to Simulate Granular Column Collapse from Failure Initiation to Runout

The performance evaluation of a potentially unstable slope involves two key components: the initiation of the slope failure and the post-failure runout. The Finite Element Method (FEM) excels at modeling the initiation of instability but quickly loses accuracy in modeling large-deformation problems due to mesh distortion. Hence, the FEM is unable to accurately model post-failure slope runout. Hybrid Eulerian-Lagrangian methods, such as the Material Point Method (MPM), offer a promising alternative for solving large-deformation problems, because particles can move freely across a background mesh, allowing for large deformation without computational issues. However, the use of moving material points in MPM for integration rather than the fixed Gauss points of the FEM reduces the accuracy of MPM in predicting stress distribution and thus failure initiation. We have created a hybrid method by initiating a failure simulation in FEM and subsequently transferring the coordinates, velocities, and stresses to MPM particles to model the runout behavior, combining the strength of both methods. We demonstrate the capability of the hybrid approach by simulating the collapse of a frictional granular column, comparing it to an empirical solution, and evaluating a suitable time to transfer from FEM to MPM by trialing multiple iterations with transfers at different stages of the collapse.

preprint2022arXiv

In situ visualization of regional-scale natural hazards with Galaxy and Material Point Method

Visualizing regional-scale landslides is essential to conveying the threat of natural hazards to stakeholders and policymakers. Traditional visualization techniques are restricted to post-processing a limited subset of simulation data and are not scalable to rendering regional-scale models. In situ visualization is a technique of rendering simulation data in real-time, i.e., rendering visuals in tandem while the simulation is running. This study develops a scalable N:M interface architecture to visualize regional-scale landslides. We demonstrate the scalability of the architecture by simulating the long runout of the 2014 Oso landslide using the Material Point Method coupled with the Galaxy ray tracing engine rendering 4.2 million material points as spheres. In situ visualization has an amortized runtime increase of 2% compared to non-visualized simulations. The developed approach can achieve in situ visualization of regional-scale landslides with minimal impact on the simulation process.

preprint2022arXiv

Minority Report: A Graph Network Oracle for In Situ Visualization

In situ visualization techniques are hampered by a lack of foresight: crucial simulation phenomena can be missed due to a poor sampling rate or insufficient detail at critical timesteps. Keeping a human in the loop is impractical, and defining statistical triggers can be difficult. This paper demonstrates the potential for using a machine-learning-based simulation surrogate as an oracle to identify expected critical regions of a large-scale simulation. These critical regions are used to drive the in situ analysis, providing greater data fidelity and analysis resolution with an equivalent I/O budget to a traditional in situ framework. We develop a distributed asynchronous in situ visualization by integrating TACC Galaxy with CB-Geo MPM for material point simulation of granular flows. We employ a PyTorch-based 3D Graph Network Simulator (GNS) trained on granular flow problems as an oracle to predict the dynamics of granular flows. Critical regions of interests are manually tagged in GNS for in situ rendering in MPM.

preprint2021arXiv

An HPC-Based Hydrothermal Finite Element Simulator for Modeling Underground Response to Community-Scale Geothermal Energy Production

Geothermal heat, as renewable energy, shows great advantage with respect to its environmental impact due to its significantly lower CO2 emissions than conventional fossil fuel. Open and closed-loop geothermal heat pumps, which utilize shallow geothermal systems, are an efficient technology for cooling and heating buildings, especially in urban areas. Integrated use of geothermal energy technologies for district heating, cooling, and thermal energy storage can be applied to optimize the subsurface for communities to provide them with multiple sustainable energy and community resilience benefits. The utilization of the subsurface resources may lead to a variation in the underground environment, which might further impact local environmental conditions. However, very few simulators can handle such a highly complex set of coupled computations on a regional or city scale. We have developed high-performance computing (HPC) based hydrothermal finite element (FE) simulator that can simulate the subsurface and its hydrothermal conditions at a scale of tens of km. The HPC simulator enables us to investigate the subsurface thermal and hydrologic response to the built underground environment (such as basements and subways) at the community scale. In this study, a coupled hydrothermal simulator is developed based on the open-source finite element library deal.II. The HPC simulator was validated by comparing the results of a benchmark case study against COMSOL Multiphysics, in which Aquifer Thermal Energy Storage (ATES) is modeled and a process of heat injection into ATES is simulated. The use of an energy pile system at the Treasure Island redevelopment site (San Francisco, CA, USA) was selected as a case study to demonstrate the HPC capability of the developed simulator. The simulator is capable of modeling multiple city-scale geothermal scenarios in a reasonable amount of time.

preprint2021arXiv

Frequency power spectra of global quantities in magnetoconvection

We present the results of direct numerical simulations of power spectral densities for kinetic energy, convective entropy and heat flux for unsteady Rayleigh-Bénard magnetoconvection in the frequency space. For larger values of frequency, the power spectral densities for all the global quantities vary with frequency $f$ as $f^{-2}$. The scaling exponent is independent of Rayleigh number, Chandrasekhar's number and thermal Prandtl number.

preprint2021arXiv

Using Convolutional Neural Networks to Develop Starting Models for 2D Full Waveform Inversion

Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamentally change engineering site characterization by enabling the recovery of high resolution 2D/3D maps of subsurface stiffness. Yet, the accuracy of FWI remains quite sensitive to the choice of the initial starting model due to the complexity and non-uniqueness of the inverse problem. In response, we present the novel application of convolutional neural networks (CNNs) to transform an experimental seismic wavefield acquired using a linear array of surface sensors directly into a robust starting model for 2D FWI. We begin by describing three key steps used for developing the CNN, which include: selection of a network architecture, development of a suitable training set, and performance of network training. The ability of the trained CNN to predict a suitable starting model for 2D FWI was compared against other commonly used starting models for a classic near-surface imaging problem; the identification of an undulating, two-layer, soil-bedrock interface. The CNN developed during this study was able to predict complex 2D subsurface images of the testing set directly from their seismic wavefields with an average mean absolute percent error of 6%. When compared to other common approaches, the CNN approach was able to produce starting models with smaller seismic image and waveform misfits, both before and after FWI. The ability of the CNN to generalize to subsurface models which were dissimilar to the ones upon which it was trained was assessed using a more complex, three-layered model. While the predictive ability of the CNN was slightly reduced, it was still able to achieve seismic image and waveform misfits comparable to the other commonly used starting models. This study demonstrates that CNNs have great potential as a tool for developing good starting models for FWI ...

preprint2020arXiv

Rhombic Patterns Near a Bicritical Point in Periodically Forced Surface Waves

We present here a study of selection of rhombic patterns close to a bicritical point at the onset of primary surface instability in viscous fluids under two-frequency vertical vibration. Rhombic patterns appear to be natural at the primary instability in the form of a bicritical point if the ratio of driving frequencies is selected properly. We present two different patterns which may be accessible in a Faraday experiment.

preprint2020arXiv

Thermal flux in unsteady Rayleigh-Bénard magnetoconvection

We present results of numerical investigation on thermal flux in Rayleigh-Bénard magnetoconvection in the presence of a uniform vertical magnetic field. We have studied thermal flux in different viscous fluids with a range of Prandtl number ($0.1 \le \mathrm{Pr} < 6.5$) and a range of Chandrasekhar number ($50 \le \mathrm{Q} \le 2.5 \times 10^4$). The power spectral density of the Nusselt number varies with frequency $f$ approximately as $f^{-2}$. The probability distribution function of the fluctuating part of the Nusselt number is nearly normal distribution with slight asymmetric tails. For a fixed value the Rayleigh number $\mathrm{Ra}$, the time averaged Nusselt number $\langle \mathrm{Nu} (\mathrm{Q})\rangle$ decreases logarithmically with Chandrasekhar number for $\mathrm{Q} > \mathrm{Q}_c$, which depends on $\mathrm{Ra}$ and $\mathrm{Pr}$. The reduced Nusselt number $\mathrm{Nu_r}$ $=$ $\langle \mathrm{Nu}(\mathrm{Q})\rangle/{\langle \mathrm{Nu}(0)\rangle}$ rises sharply, reaches a maximum slightly above unity and then start decreasing very slowly to unity as the value of a dimensionless parameter $\sqrt{\mathrm{Ra/(Q~Pr)}}$ is raised. The probability distribution function of the local thermal flux in the vertical direction is found to be asymmetric and non-Gaussian with a cusp at its maximum.

preprint2019arXiv

Instability of a horizontal water half-cylinder under vertical vibration

We present the results of an experimental investigation on parametrically driven waves in a water half-cylinder on a rigid horizontal plate, which is sinusoidally vibrated in the vertical direction. As the forcing amplitude is raised above a critical value, stationary waves are excited in the water half-cylinder. Parametrically excited subharmonic waves are non-axisymmetricand qualitatively different from the axisymmetric Savart-Plateau-Rayleigh waves in a vertical liquid cylinder or jet. Depending on the driving frequency, stationary waves of different azimuthal wave numbers are excited. A linear theory is also supplemented, which captures the observed dispersion relations quantitatively.

preprint2016arXiv

Effects of a small magnetic field on homoclinic bifurcations in a low-Prandtl-number fluid

Effects of a uniform magnetic field on homoclinic bifurcations in Rayleigh-Bénard convection in a fluid of Prandtl number $Pr = 0.01$ are investigated using direct numerical simulations (DNS). A uniform magnetic field is applied either in the vertical or in the horizontal direction. For a weak vertical magnetic field, the possibilities of both forward and backward homoclinic bifurcations are observed leading to a spontaneous merging of two limit cycles into one as well as a spontaneous breaking of a limit cycle into two for lower values of the Chandrasekhar's number ($Q\leq 5$). A slightly stronger magnetic field makes the convective flow time independent giving the possibility of stationary patterns at the secondary instability. For horizontal magnetic field, the $x\leftrightharpoons y$ symmetry is destroyed and neither a homoclinic gluing nor a homoclinic breaking is observed. Two low-dimensional models are also constructed: one for a weak vertical magnetic field and another for a weak horizontal magnetic field. The models qualitatively capture the features observed in DNS and help understanding the unfolding of bifurcations close to the onset of magnetoconvection.

preprint2016arXiv

Numerical evidence for a chiral spin liquid in the XXZ antiferromagnetic Heisenberg model on the kagome lattice at $m=\frac{2}{3}$ magnetization

We perform an exact diagonalization study of the spin-$1/2$ XXZ Heisenberg antiferromagnet on the kagome lattice at finite magnetization $m = \frac{2}{3}$ with an emphasis on the XY point ($J_z = 0$), and in the presence of a small chiral term. Recent analytic work by Kumar, Sun and Fradkin [Phys. Rev. B 90, 174409 (2014)] on the same model, using a newly developed flux attachment transformation, predicts a plateau at this value of the magnetization described by a chiral spin liquid (CSL) with a spin Hall conductance of $σ_{xy} = \frac{1}{2}$. Such a state is topological in nature, has a ground state degeneracy and exhibits fractional excitations. We analyze the degeneracy structure in the low energy manifold, identify the candidate topological states and use them to compute the modular matrices and Chern numbers all of which strongly agree with expected theoretical behavior for the $σ_{xy} = \frac{1}{2}$ CSL. We argue that the evidence suggests the CSL is robust even in the limit of zero external chirality.

preprint2016arXiv

Spectra and probability distributions of thermal flux in turbulent Rayleigh-Bénard convection

The spectra of turbulent heat flux $\mathrm{H}(k)$ in Rayleigh-Bénard convection with and without uniform rotation are presented. The spectrum $\mathrm{H}(k)$ scales with wave number $k$ as $\sim k^{-2}$. The scaling exponent is almost independent of the Taylor number $\mathrm{Ta}$ and Prandtl number $\mathrm{Pr}$ for higher values of the reduced Rayleigh number $r$ ($ > 10^3$). The exponent, however, depends on $\mathrm{Ta}$ and $\mathrm{Pr}$ for smaller values of $r$ ($<10^3$). The probability distribution functions of the local heat fluxes are non-Gaussian and have exponential tails.

preprint2015arXiv

A discretized Chern-Simons gauge theory on arbitrary graphs

In this paper, we show how to discretize the abelian Chern-Simons gauge theory on generic planar lattices/graphs (with or without translational symmetries) embedded in arbitrary 2D closed orientable manifolds. We find that, as long as a one-to-one correspondence between vertices and faces can be defined on the graph such that each face is paired up with a neighboring vertex (and vice versa), a discretized Chern-Simons theory can be constructed consistently. We further verify that all the essential properties of the Chern-Simons gauge theory are preserved in the discretized setup. In addition, we find that the existence of such a one-to-one correspondence is not only a sufficient condition for discretizing a Chern-Simons gauge theory but, for the discretized theory to be nonsingular and to preserve some key properties of the topological field theory, this correspondence is also a necessary one. A specific example will then be provided, in which we discretize the Chern-Simons gauge theory on a tetrahedron.

preprint2015arXiv

A model for Rayleigh-Bénard magnetoconvection

A model for three-dimensional Rayleigh-Bénard convection in low-Prandtl-number fluids near onset with rigid horizontal boundaries in the presence of a uniform vertical magnetic field is constructed and analyzed in detail. The kinetic energy $K$, the convective entropy $Φ$ and the convective heat flux ($Nu-1$) show scaling behaviour with $ε= r-1$ near onset of convection, where $r$ is the reduced Rayleigh number. The model is also used to investigate various magneto-convective structures close to the onset. Straight rolls, which appear at the primary instability, become unstable with increase in $r$ and bifurcate to three-dimensional structures. The straight rolls become periodically varying wavy rolls or quasiperiodically varying structures in time with increase in $r$ depending on the values of Prandtl number $Pr$. They become irregular in time, with increase in $r$. These standing wave solutions bifurcate first to periodic and then quasiperiodic traveling wave solutions, as $r$ is raised further. The variations of the critical Rayleigh number $Ra_{os}$ and the frequency $ω_{os}$ at the onset of the secondary instability with $Pr$ are also studied for different values of Chandrasekhar's number $Q$.

preprint2015arXiv

Chiral spin liquids on the kagome Lattice

We study the nearest neighbor $XXZ$ Heisenberg quantum antiferromagnet on the kagome lattice. Here we consider the effects of several perturbations: a) a chirality term, b) a Dzyaloshinski-Moriya term, and c) a ring-exchange type term on the bowties of the kagome lattice, and inquire if they can suppport chiral spin liquids as ground states. The method used to study these Hamiltonians is a flux attachment transformation that maps the spins on the lattice to fermions coupled to a Chern-Simons gauge field on the kagome lattice. This transformation requires us to consistently define a Chern-Simons term on the kagome lattice. We find that the chirality term leads to a chiral spin liquid even in the absence of an uniform magnetic field, with an effective spin Hall conductance of $\sxy = \frac{1}{2}$ in the regime of $XY$ anisotropy. The Dzyaloshinkii-Moriya term also leads a similar chiral spin liquid but only when this term is not too strong. An external magnetic field also has the possibility of giving rise to additional plateaus which also behave like chiral spin liquids in the $XY$ regime. Finally, we consider the effects of a ring-exchange term and find that, provided its coupling constant is large enough, it may trigger a phase transition into a chiral spin liquid by the spontaneous breaking of time-reversal invariance.

preprint2015arXiv

Measuring the weak mixing angle with the P2 experiment at MESA

The P2 experiment in Mainz aims to measure the weak mixing angle in electron- proton scattering to a precision of 0.13 %. In order to suppress uncertainties due to proton structure and contributions from box graphs, both a low average momentum transfer $Q^2$ of $4.5\cdot 10^{-3}$ GeV$^2/c^2$ and a low beam energy of 155 MeV are chosen. In order to collect the enormous statistics required for this measurement, the new Mainz Energy Recovery Superconducting Accelerator (MESA) is being constructed. These proceedings describe the motivation for the measurement, the experimental and accelerator challenges and how we plan to tackle them.

preprint2014arXiv

Chern-Simons theory of the magnetization plateaus of the spin-1/2 quantum XXZ Heisenberg model on Kagome Lattice

Frustrated spin systems on Kagome lattices have long been considered to be a promising candidate for realizing exotic spin liquid phases. Recently, there has been a lot of renewed interest in these systems with the discovery of materials such as Volborthite and Herbertsmithite that have Kagome like structures. In the presence of an external magnetic field, these frustrated systems can give rise to magnetization plateaus of which the plateau at $m=\frac{1}{3}$ is considered to be the most prominent. Here we study the problem of the antiferromagnetic spin-1/2 quantum XXZ Heisenberg model on a Kagome lattice by using a Jordan-Wigner transformation that maps the spins onto a problem of fermions coupled to a Chern-Simons gauge field. This mapping relies on being able to define a consistent Chern-Simons term on the lattice. Using a recently developed method to rigorously extend the Chern-Simons term to the frustrated Kagome lattice we can now formalize the Jordan-Wigner transformation on the Kagome lattice. We then discuss the possible phases that can arise at the mean-field level from this mapping and focus specifically on the case of $\frac{1}{3}$-filling ($m=\frac{1}{3}$ plateau) and analyze the effects of fluctuations in our theory. We show that in the regime of $XY$ anisotropy the ground state at the $1/3$ plateau is equivalent to a bosonic fractional quantum Hall Laughlin state with filling fraction $1/2$ and that at the $5/9$ plateau it is equivalent to the first bosonic Jain daughter state at filling fraction $2/3$.

preprint2014arXiv

Early Development of UVM based Verification Environment of Image Signal Processing Designs using TLM Reference Model of RTL

With semiconductor industry trend of smaller the better, from an idea to a final product, more innovation on product portfolio and yet remaining competitive and profitable are few criteria which are culminating into pressure and need for more and more innovation for CAD flow, process management and project execution cycle. Project schedules are very tight and to achieve first silicon success is key for projects. This necessitates quicker verification with better coverage matrix. Quicker Verification requires early development of the verification environment with wider test vectors without waiting for RTL to be available. In this paper, we are presenting a novel approach of early development of reusable multi-language verification flow, by addressing four major activities of verification like Early creation of Executable Specification, Early creation of Verification Environment, Early development of test vectors and Better and increased Re-use of blocks. Although this paper focuses on early development of UVM based Verification Environment of Image Signal Processing designs using TLM Reference Model of RTL, same concept can be extended for non-image signal processing designs. Main Keywords are SystemVerilog, SystemC, Transaction Level Modeling, Universal Verification Methodology (UVM), Processor model, Universal Verification Component (UVC), Reference Model.

preprint2014arXiv

Entropy and energy spectra in low-Prandtl-number convection with rotation

We present results for entropy and kinetic energy spectra computed from direct numerical simulations for low-Prandtl-number ($Pr < 1$) turbulent flow in Rayleigh-Bénard convection with uniform rotation about a vertical axis. The simulations are performed in a three-dimensional periodic box for a range of Taylor number ($ 0 \leq Ta \leq 10^8$) and reduced Rayleigh number $r = Ra/Ra_{\circ} (Ta, Pr)$ ($1.0 \times 10^2 \le r \le 5.0 \times 10^3$). The Rossby number $Ro$ varies in the range $1.34 \le Ro \le 73$. The entropy spectrum $E_θ(k)$ shows bi-splitting into two branches for lower values of wave number $k$. The entropy in the lower branch scales with $k$ as $k^{-1.4\pm 0.1}$ for $r > 10^3$ for the rotation rates considered here. The entropy in the upper branch also shows scaling behavior with $k$, but the scaling exponent decreases with increasing $Ta$ for all $r$. The energy spectrum $E_v(k)$ is also found to scale with the wave number $k$ as $k^{-1.4\pm 0.1}$ for $r > 10^3$. The scaling exponent for the energy spectrum and the lower branch of the entropy spectrum vary between $-1.7$ to $-2.4$ for lower values of $r$ ($< 10^3$). We also provide some simple arguments based on the variation of the Kolmogorov picture to support the results of simulations.

preprint2014arXiv

Frequency spectra of turbulent thermal convection with uniform rotation

The frequency spectra of the entropy and kinetic energy along with the power spectrum of the thermal flux are computed from direct numerical simulations for turbulent Rayleigh-Bénard convection with uniform rotation about a vertical axis in low-Prandtl-number fluids ($\mathrm{Pr} < 0.6$). Simulations are done for convective Rossby numbers $\mathrm{Ro} \ge 0.2$. The temporal fluctuations of these global quantities show two scaling regimes: (i) $ω^{-2}$ at higher frequencies for all values of $\mathrm{Ro}$ and (ii) $ω^{-γ_1}$ at intermediate frequencies with $γ_1 \approx 4$ for $\mathrm{Ro} > 1$, while $4 < γ_1 < 6.6$ for $0.2 \le \mathrm{Ro} < 1$.

preprint2014arXiv

Oscillatory instability and fluid patterns in low-Prandtl-number Rayleigh-Bénard convection with uniform rotation

We present the results of direct numerical simulations of flow patterns in a low-Prandtl-number ($Pr = 0.1$) fluid above the onset of oscillatory convection in a Rayleigh-Bénard system rotating uniformly about a vertical axis. Simulations were carried out in a periodic box with thermally conducting and stress-free top and bottom surfaces. We considered a rectangular box ($L_x \times L_y \times 1$) and a wide range of Taylor numbers ($750 \le Ta \le 5000$) for the purpose. The horizontal aspect ratio $η= L_y/L_x$ of the box was varied from $0.5$ to $10$. The primary instability appeared in the form of two-dimensional standing waves for shorter boxes ($0.5 \le η< 1$ and $1 < η< 2$). The flow patterns observed in boxes with $η= 1$ and $η= 2$ were different from those with $η< 1$ and $1 < η< 2$. We observed a competition between two sets of mutually perpendicular rolls at the primary instability in a square cell ($η= 1$) for $Ta < 2700$, but observed a set of parallel rolls in the form of standing waves for $Ta \geq 2700$. The three-dimensional convection was quasiperiodic or chaotic for $750 \le Ta < 2700$, and then bifurcated into a two-dimensional periodic flow for $Ta \ge 2700$. The convective structures consisted of the appearance and disappearance of straight rolls, rhombic patterns, and wavy rolls inclined at an angle $ϕ= \fracπ{2} - \arctan{(η^{-1})}$ with the straight rolls.

preprint2014arXiv

Rayleigh-Bénard convection with uniform vertical magnetic field

We present the results of direct numerical simulations of Rayleigh-Bénard convection in the presence of a uniform vertical magnetic field near instability onset. We have done simulations in boxes with square as well as rectangular cross-sections in the horizontal plane. We have considered horizontal aspect ratio $η= L_y/L_x =1$ and $2$. The onset of the primary and secondary instabilities are strongly suppressed in the presence of the vertical magnetic field for $η=1$. The Nusselt number $\mathrm{Nu}$ scales with Rayleigh number $\mathrm{Ra}$ close to the primary instability as $[\mathrm{\{Ra - Ra_c (Q)\}/Ra_c (Q)} ]^{0.91}$, where $\mathrm{Ra_c (Q)}$ is the threshold for onset of stationary convection at a given value of the Chandrasekhar number $\mathrm{Q}$. $\mathrm{Nu}$ also scales with $\mathrm{Ra/Q}$ as $(\mathrm{Ra/Q})^μ$. The exponent $μ$ varies in the range $0.39 \le μ\le 0.57$ for $\mathrm{Ra/Q} \ge 25$. The primary instability is stationary as predicted by Chandrasekhar. The secondary instability is temporally periodic for $\mathrm{Pr}=0.1$ but quasiperiodic for $\mathrm{Pr} = 0.025$ for moderate values of $\mathrm{Q}$. Convective patterns for higher values of $\mathrm{Ra}$ consist of periodic, quasiperiodic and chaotic wavy rolls above onset of the secondary instability for $η=1$. In addition, stationary as well as time dependent cross-rolls are observed, as $\mathrm{Ra}$ is further raised. The ratio $r_{\circ}/\mathrm{Pr}$ is independent of $\mathrm{Q}$ for smaller values of $\mathrm{Q}$. The delay in onset of the oscillatory instability is significantly reduced in a simulation box with $η=2$. We also observe inclined stationary rolls for smaller values of $\mathrm{Q}$ for $η=2$.

preprint2014arXiv

Zero-Prandtl-number convection with slow rotation

We present the results of our investigations of the primary instability and the flow patterns near onset in zero-Prandtl-number Rayleigh-Bénard convection with uniform rotation about a vertical axis. The investigations are carried out using direct numerical simulations of the hydrodynamic equations with stress-free horizontal boundaries in rectangular boxes of size $(2π/k_x) \times (2π/k_y) \times 1$ for different values of the ratio $η= k_x/k_y$. The primary instability is found to depend on $η$ and $Ta$. Wavy rolls are observed at the primary instability for smaller values of $η$ ($1/\sqrt{3} \le η\le 2$ except at $η= 1$) and for smaller values of $Ta$. We observed Küppers-Lortz (KL) type patterns at the primary instability for $η= 1/\sqrt{3}$ and $ Ta \ge 40$. The fluid patterns are found to exhibit the phenomenon of bursting, as observed in experiments [Bajaj et al. Phys. Rev. E {\bf 65}, 056309 (2002)]. Periodic wavy rolls are observed at onset for smaller values of $Ta$, while KL-type patterns are observed for $ Ta \ge 100$ for $η=\sqrt{3}$. In case of $η= 2$, wavy rolls are observed for smaller values of $Ta$ and KL-type patterns are observed for $25 \le Ta \le 575$. Quasi-periodically varying patterns are observed in the oscillatory regime ($Ta > 575$). The behavior is quite different at $η= 1$. A time dependent competition between two sets of mutually perpendicular rolls is observed at onset for all values of $Ta$ in this case. Fluid patterns are found to burst periodically as well as chaotically in time. It involved a homoclinic bifurcation. We have also made a couple of low-dimensional models to investigate bifurcations for $η= 1$, which is used to investigate the sequence of bifurcations.

preprint2013arXiv

Fundamental Symmetry Tests with Nucleons, Nuclei, and Atoms: A Snowmass Report

Present and prospective fundamental symmetry tests with nucleons, nuclei and atoms are probing for possible new physics at the TeV scale and beyond. These ongoing and proposed table-top as well as accelerator-based experiments are thus a vital component of the Intensity Frontier. At the same time, these tests provide increasingly sophisticated probes of long-distance strong interactions that are responsible for the structure of nucleons and nuclei. In this community report, some of the most compelling opportunities with nucleons, nuclei and atoms are summarized, drawing largely on input received from the nuclear and atomic physics communities. In particular, this report includes many contributions submitted to two recent Intensity Frontier Workshops.

preprint2013arXiv

Pattern dynamics near homoclinic bifurcation in Rayleigh-Bénard convection

We report for the first time the pattern dynamics in the vicinity of an inverse homoclinic bifurcation in an extended dissipative system. We observe, in direct numerical simulations of three dimensional Rayleigh-Bénard convection, a spontaneous breaking of a competition of two mutually perpendicular sets of oscillating cross rolls to one of two possible sets of oscillating cross rolls as the Rayleigh number is raised above a critical value. The time period of the cross-roll patterns diverges, and shows scaling behavior near the bifurcation point. This is an example of a transition from nonlocal to local pattern dynamics near an inverse homoclinic bifurcation. We also present a simple four-mode model that captures the pattern dynamics quite well.

preprint2013arXiv

Role of uniform horizontal magnetic field on convective flow

The effect of uniform magnetic field applied along a fixed horizontal direction in Rayleigh-Bénard convection in low-Prandtl-number fluids has been studied using a low dimensional model. The model shows the onset of convection (primary instability) in the form of two dimensional stationary rolls in the absence of magnetic field, when the Rayleigh number $R$ is raised above a critical value $R_c$. The flow becomes three dimensional at slightly higher values of Rayleigh number via wavy instability. These wavy rolls become chaotic for slightly higher values of $R$ in low-Prandtl-number ($P_r$) fluids. A uniform magnetic field along horizontal plane strongly affects all kinds of convective flows observed at higher values of $R$ in its absence. As the magnetic field is raised above certain value, it orients the convective rolls in its own direction. Although the horizontal magnetic field does not change the threshold for the primary instability, it affects the threshold for secondary (wavy) instability. It inhibits the onset of wavy instability. The critical Rayleigh number $R_o (Q,P_r)$ at the onset of wavy instability, which depends on Chandrasekhar's number $Q$ and $P_r$, increases monotonically with $Q$ for a fixed value of $P_r$. The dimensionless number $R_o (Q, P_r)/(R_c Q P_r)$ scales with $Q$ as $Q^{-1}$. A stronger magnetic field suppresses chaos and makes the flow two dimensional with roll pattern aligned along its direction.

preprint2013arXiv

Theory Uncertainty in Extracting the Proton's Weak Charge: White Paper

We review the state-of-the-art and address open questions relative to the theory uncertainty of the $γ-Z$ box contribution to the $A_{PV}$ measurement within the QWEAK experiment at Jefferson Lab. This white paper summarizes the contributions by participants and discussion sessions on this topic within the MITP Workshop on Precision Electroweak Physics held in Mainz, Germany, September 23 - October 11, 2013 \url{http://indico.cern.ch/conferenceDisplay.py?confId=248072

preprint2012arXiv

Physics Opportunities with the 12 GeV Upgrade at Jefferson Lab

This white paper summarizes the scientific opportunities for utilization of the upgraded 12 GeV Continuous Electron Beam Accelerator Facility (CEBAF) and associated experimental equipment at Jefferson Lab. It is based on the 52 proposals recommended for approval by the Jefferson Lab Program Advisory Committee.The upgraded facility will enable a new experimental program with substantial discovery potential to address important topics in nuclear, hadronic, and electroweak physics.