Search

Particle production in cosmology is often efficiently computed in terms of Bogoliubov transforms. Restricting to a particular class of dispersion relationships, we identify a map between the number of particles produced in a special kinematic limit and a Stokes phenomena related topology of analytic continuation of the Bogoliubov coefficient functions. Intuitively, this kinematic limit corresponds to the long wavelength limit although a more precise description depends on the nature of the curved spacetime. To identify the topology, we reformulate the usual Bogoliubov computations as a type of SU(1,1) gauged differential equation and utilize a special gauge together with a discrete symmetry that naturally characterizes the dispersion relationship. Using a dark matter model and a nonzero constant spatial curvature model, we estimate how such topological contributions will arise in physical applications., Comment: 42 pages + appendices and references, 8 figures

We investigate nonlinear energy transfer for channel flows at friction Reynolds numbers of $Re_{\tau}=180$ and $590$. The key feature of the analysis is that we introduce a new variable, which quantifies the energy transferred from a source mode to a recipient mode through explicit examination of nonlinear triadic interactions in streamwise-spanwise wavenumber space. First, we use this variable to quantify the nonlinear energy transfer gain and loss for individual Fourier modes. The nonlinear energy transfer gain and loss cannot be directly obtained from the turbulent kinetic energy (TKE) equation. Second, we quantify the nonlinear energy transfer budgets for three types of structures: streamwise streaks, oblique waves and Tollmien-Schlichting waves. We found that a transverse cascade from streamwise-elongated modes to spanwise-elongated modes exists in all three structures. Third, we quantify the forward and inverse cascades between resolved scales and subgrid scales in the spirit of large-eddy simulation. For the cutoff wavelength range we consider, the forward and inverse cascades between the resolved scales and subgrid scales result in a net forward cascade from the resolved scales to the subgrid scales. The shape of the net forward cascade curve with respect to the cutoff wavelength resembles the net forward cascade predicted by the Smagorinsky eddy viscosity.

One well known method of generating a large blue spectral index for axionic isocurvature perturbations is through a flat direction not having a quartic potential term for the radial partner of the axion field. In this work, we show how one can obtain a large blue spectral index even with a quartic potential term associated with the Peccei-Quinn symmetry breaking radial partner. We use the fact that a large radial direction with a quartic term can naturally induce a conformal limit which generates an isocurvature spectral index of 3. We point out that this conformal representation is intrinsically different from both the ordinary equilibrium axion scenario or massless fields in Minkowski spacetime. Another way to view this limit is as a scenario where the angular momentum of the initial conditions slows down the radial field or as a superfluid limit. Quantization of the non-static system in which derivative of the radial field and the derivative of the angular field do not commute is treated with great care to compute the vacuum state. The parametric region consistent with axion dark matter and isocurvature cosmology is discussed., Comment: 52 pages of main text, 13 figures

To aid in prediction of turbulent boundary layer flows over rough surfaces, a new model is proposed to estimate hydrodynamic roughness based solely on geometric surface information. The model is based on a fluid-mechanics motivated geometric parameter called the wind-shade factor. Sheltering is included using a rapid algorithm adapted from the landscape shadow literature, while local pressure drag is estimated using a piecewise potential flow approximation. Similarly to evaluating traditional surface parameters such as skewness or average slope magnitude, the wind-shade factor is purely geometric and can be evaluated efficiently from knowing the surface elevation map and the mean flow direction. The wind-shade roughness model is applied to over 100 different surfaces available in a public roughness database and some others, and the predicted sandgrain-roughness heights are compared to measured values. Effects of various model ingredients are analyzed, and transitionally rough surfaces are treated by adding a term representing the viscous stress component., Comment: Paper submitted to J. Fluid Mechanics for publication, currently under review. Contains Notebooks but links will work only once it is published

Large blue tilted spectral index axionic isocurvature perturbations can be produced when the axion sector is far out of equilibrium during inflation through an initial Peccei-Quinn (PQ) symmetry breaking field displacement along a nearly flat direction in the effective potential. As a companion to a previous work, we present analytic formulae for the blue isocurvature spectrum for the case of the kinetic energy density of the PQ symmetry breaking field being larger than the quartic power of the final spontaneous PQ symmetry breaking scale. It corresponds to a regime in which the nonlinearities of the classical potential become important many times during the formation of the axion isocurvature quantum perturbations leading to interesting resonant behavior. One consequence of this nonlinearity-driven resonance is the chaotic nature of the map that links the underlying Lagrangian parameters to the isocurvature amplitudes. We point out an accidental duality symmetry between the perturbation equations and the background field equations that can be used to understand this. Finally, we present two types of analytic results. The first relies on a computation utilizing an effective potential wherein fast time scale fluctuations have been integrated out. The second is grounded in a functional ansatz, requiring only a limited set of fitting parameters. Both analytic results should be useful for carrying out forecasts and fits to the data., Comment: 64 pages, 19 figures, 1 table

This study presents an experimental dataset documenting the evolution of a turbulent boundary layer downstream of a rough-to-smooth surface transition. To investigate the effect of upstream flow conditions, two groups of experiments are conducted. For the \emph{Group-Re} cases, a nominally constant viscous-scaled equivalent sand grain roughness $k_{s0}^+\approx160$ is maintained on the rough surface, while the friction Reynolds number $Re_{\tau 0}$ ranges from 7100 to 21000. For the \emph{Group-ks} cases, $Re_{\tau 0}\approx14000$ is maintained while $k_{s0}^+$ ranges from 111 to 228. The wall-shear stress on the downstream smooth surface is measured directly using oil-film interferometry to redress previously reported uncertainties in the skin-friction coefficient recovery trends. In the early development following the roughness transition, the flow in the internal layer is not in equilibrium with the wall-shear stress. This conflicts with the common practise of modelling the mean velocity profile as two log laws below and above the internal layer height, as first proposed by Elliott (\textit{Trans. Am. Geophys. Union}, vol. 39, 1958, pp 1048--1054). As a solution to this, the current data are used to model the recovering mean velocity semi-empirically by blending the corresponding rough-wall and smooth-wall profiles. The over-energised large-scale motions leave a strong footprint in the near-wall region of the energy spectrum, the frequency and magnitude of which exhibit dependence on $Re_{\tau 0}$ and $k_{s0}^+$ respectively. The energy distribution in near-wall small scales is mostly unaffected by the presence of the outer flow with rough-wall characteristics, which can be used as a surrogate measure to extract the local friction velocity.

This paper examines recovery of the wall-shear stress of a turbulent boundary layer that has undergone a sudden transition from a rough to a smooth surface. Early works of Antonia and Luxton questioned the reliability of standard smooth-wall methods to measure wall-shear stress in such conditions, and subsequent studies show significant disagreement depending on the approach used to determine the wall-shear stress downstream. Here we address this by utilising a collection of experimental databases at Re_\tau \approx 4100 that have access to both `direct' and `indirect' measures of the wall-shear stress to understand the recovery to equilibrium conditions to the new surface. Our results reveal that the viscous region (z^+\lesssim 4) recovers almost immediately to an equilibrium state with the new wall conditions, however, the buffer region and beyond takes several boundary layer thicknesses before recovering to equilibrium conditions, which is longer than previously thought. A unique direct numerical simulation database of a wall-bounded flow with a rough-to-smooth wall transition is employed to confirm these findings. In doing so, we present evidence that any estimate of the wall-shear stress from the mean velocity profile in the buffer region or further away from the wall tends to underestimate its magnitude in the near vicinity of the rough-to-smooth transition, and this is likely to be partly responsible for the large scatter of recovery lengths to equilibrium conditions reported in the literature. Our results also reveal that the smaller energetic scales in the near-wall region recover to an equilibrium state associated with the new wall conditions within one boundary layer thickness downstream of the transition, while the larger energetic scales exhibit an over-energised state for several boundary layer thicknesses downstream of the transition.

Isocurvature perturbations with a blue power spectrum are one of the natural targets for the future large scale structure observations which are probing shorter length scales with greater accuracy. We present a Fisher forecast for the Euclid and MegaMapper (MM) experiments in their ability to detect blue isocurvature perturbations. We construct the theoretical predictions in the EFTofLSS and bias expansion formalisms at quartic order in overdensities which allows us to compute the power spectrum at one loop order and bispectrum at tree level and further include theoretical error at the next to leading order for the covariance determination. We find that Euclid is expected to provide at least a factor of few improvement on the isocurvature spectral amplitude compared to the existing Planck constraints for large spectral indices while MM is expected to provide about 1 to 1.5 order of magnitude im provement for a broad range of spectral indices. We find features that are specific to the blue isocurvature scenario including the leading parametric degeneracy being with the Laplacian bias and a UV sensitive bare sound speed parameter., Comment: v2: 45 pages (32+13), 9 figures, 2 tables, minor corrections (results same as in v1)

Previous numerical investigations of gravitational particle production during the coherent oscillation period of inflation displayed unexplained fluctuations in the spectral density of the produced particles. We argue that these features are due to the quantum interference of the coherent scattering reactions that produce the particles. We provide accurate analytic formulae to compute the particle production amplitude for a conformally-coupled scalar field, including the interference effect in the kinematic region where the production can be interpreted as inflaton scattering into scalar final states via graviton exchange., Comment: 19 pages + appendices and references, 3 figures, v2 - references added, matches JHEP

Previous computations of strongly blue tilted axionic isocurvature spectra were computed in the parametric region in which the lightest time-dependent mass is smaller than the Hubble expansion rate during inflation, leading to an overdamped time evolution. Here we present the strongly blue tilted axionic isocurvature spectrum in an underdamped time evolution parametric regime. Somewhat surprisingly, there exist parametric regions with a strong resonant spectral behavior that leads to a rich isocurvature spectral shape. We focus on computing this resonant spectrum analytically in a large parametric region amenable to such computations. Because the spectrum is sensitive to nonperturbative classical field dynamics, a wide variety of analytic techniques are used including a time-space effective potential obtained by integrating out high frequency fluctuations., Comment: 120 pages, 25 figures, added references and better typesetting

Analytic and numerical techniques are presented for computing gravitational production of scalar particles in the limit that the inflaton mass is much larger than the Hubble expansion rate at the end of inflation. These techniques rely upon adiabatic invariants and time modeling of a typical inflaton field which has slow and fast time variation components. A faster computation time for numerical integration is achieved via subtraction of slowly varying components that are ultimately exponentially suppressed. The fast oscillatory remnant results in production of scalar particles with a mass larger than the inflationary Hubble expansion rate through a mechanism analogous to perturbative particle scattering. An improved effective Boltzmann collision equation description of this particle production mechanism is developed. This model allows computation of the spectrum using only adiabatic invariants, avoiding the need to explicitly solve the inflaton equations of motion., Comment: 37 pages, 5 figures

Direct numerical simulation is used to study turbulent flow over irregular rough surfaces in the periodic minimal channel configuration. The generation of irregular rough surface is based on a random algorithm, in which the power spectrum of the roughness height function along with its probability density function can be directly prescribed. The hydrodynamic properties of the roughness are investigated and compared to those obtained from full-size DNS for 12 roughness topographies with systematically varied PDF and PS at four roughness height. The comparison confirms the viability of the minimal channel approach for characterization of rough surfaces providing excellent agreement in roughness function and zero-plane displacement across various types of roughness and different regimes. Results also indicates that different realizations of roughness, with a fixed PS and PDF, translate to similar values of roughness function with a small scatter. In addition to the global flow properties, the distribution of time-averaged surface force exerted by the roughness onto the fluid is examined and compared to the roughness height distribution for different cases. It is shown that the surface force distribution has an anisotropic structure with spanwise-elongated coherent regions. The anisotropy translates into a very small streamwise integral length scale, which weakly depends on the considered roughness topography, while the larger spanwise integral length scale shows a stronger dependence on roughness characteristics. It is also shown that the sheltering model describes well the spatial distribution of the surface force. Finally, existing roughness correlations are assessed using the present dataset. It was shown that the most correlations can reproduce the values of equivalent sand-grain roughness from DNS within+-30% error while none of the correlations shows a superior predictive accuracy., Comment: 34 pages, 18 figures, 3 tables, submitted to Journal of Fluid Mechanics

Wall-roughness induces extra drag in wall-bounded turbulent flows. Mapping any given roughness geometry to its fluid dynamic behaviour has been hampered by the lack of accurate and direct measurements of skin-friction drag. Here the Taylor-Couette (TC) system provides an opportunity as it is a closed system and allows to directly and reliably measure the skin-friction. However, the wall-curvature potentially complicates the connection between the wall friction and the wall roughness characteristics. Here we investigate the effects of a hydrodynamically fully rough surface on highly turbulent, inner cylinder rotating, TC flow. We find that the effects of a hydrodynamically fully rough surface on TC turbulence, where the roughness height k is three orders of magnitude smaller than the Obukhov curvature length Lc (which characterizes the effects of curvature on the turbulent flow, see Berghout et al. arXiv: 2003.03294, 2020), are similar to those effects of a fully rough surface on a flat plate turbulent boundary layer (BL). Hence, the value of the equivalent sand grain height ks, that characterizes the drag properties of a rough surface, is similar to those found for comparable sandpaper surfaces in a flat plate BL. Next, we obtain the dependence of the torque (skin-friction drag) on the Reynolds number for given wall roughness, characterized by ks, and find agreement with the experimental results within 5 percent. Our findings demonstrate that global torque measurements in the TC facility are well suited to reliably deduce wall drag properties for any rough surface., Comment: 18 pages, 13 figures

We conduct direct simulations of turbulent channels imposing different virtual origins for all three velocities using Robin, slip-like boundary conditions to study the effect of displacing the origins perceived by different flow components, a mechanism common to small-textured surfaces. We also explore this effect using opposition control. For riblets, Luchini et al. (1991) proposed that their effect could be reduced to the offset between the streamwise- and spanwise-velocity origins, the latter being the origin perceived by turbulence. Later results on superhydrophobic surfaces suggest that the apparent wall-normal-velocity origin could also play a role. Our results support that the relevant parameter is the offset between the origins perceived by the mean flow and by turbulence. The former is determined by the streamwise slip length, but the latter depends on both the wall-normal and spanwise slip lengths. The slip on the streamwise velocity fluctuations, in turn, has a negligible effect on the virtual origin for turbulence, at least in the drag-reduction regime. This suggests that the origin perceived by the quasi-streamwise vortices, which induce the cross velocities at the surface, is key in determining the virtual origin for turbulence, while that perceived by the near-wall streaks, which are associated with the streamwise velocity fluctuations, plays a secondary role. In this framework, the changes in turbulent quantities typically reported in the flow-control literature are shown to be only a result of the choice of origin and are absent when referenced to the origin experienced by turbulence. Other than this shift in origin, we show that turbulence remains essentially smooth-wall-like, and a simple expression can predict the virtual origin for turbulence. The effect on drag can be reproduced a priori by introducing the virtual origins into a smooth-wall eddy-viscosity framework., Comment: 38 pages, 23 figures. Changes made to address reviewer comments

Taylor--Couette (TC) flow is the shear-driven flow between two coaxial independently rotating cylinders. In recent years, high-fidelity simulations and experiments revealed the shape of the streamwise and angular velocity profiles up to very high Reynolds numbers. However, due to curvature effects, so far no theory has been able to correctly describe the turbulent streamwise velocity profile for given radius ratio, as the classical Prandtl--von K\'arm\'an logarithmic law for turbulent boundary layers over a flat surface at most fits in a limited spatial region. Here we address this deficiency by applying the idea of a Monin--Obukhov curvature length to turbulent TC flow. This length separates the flow regions where the production of turbulent kinetic energy is governed by pure shear from that where it acts in combination with the curvature of the streamlines. We demonstrate that for all Reynolds numbers and radius ratios, the mean streamwise and angular velocity profiles collapse according to this separation. We then derive the functional form of the velocity profile. Finally, we match the newly derived angular velocity profile with the constant angular momentum profile at the height of the boundary layer, to obtain the dependence of the torque on the Reynolds number, or, in other words, of the generalized Nusselt number (i.e., the dimensionless angular velocity transport) on the Taylor number.

Basis tensor gauge theory (BTGT) is a vierbein analog reformulation of ordinary gauge theories in which the vierbein field describes the Wilson line. After a brief review of the BTGT, we clarify the Lorentz group representation properties associated with the variables used for its quantization. In particular, we show that starting from an SO(1,3) representation satisfying the Lorentz-invariant U(1,3) matrix constraints, BTGT introduces a Lorentz frame choice to pick the Abelian group manifold generated by the Cartan subalgebra of u(1,3) for the convenience of quantization even though the theory is frame independent. This freedom to choose a frame can be viewed as an additional symmetry of BTGT that was not emphasized before. We then show how an $S_4$ permutation symmetry and a parity symmetry of frame fields natural in BTGT can be used to construct renormalizable gauge theories that introduce frame dependent fields but remain frame independent perturbatively without any explicit reference to the usual gauge field.

Assessing tumor tissue heterogeneity via ultrasound has recently been suggested for predicting early response to treatment. The ultrasound backscattering characteristics can assist in better understanding the tumor texture by highlighting local concentration and spatial arrangement of tissue scatterers. However, it is challenging to quantify the various tissue heterogeneities ranging from fine-to-coarse of the echo envelope peaks in tumor texture. Local parametric fractal features extracted via maximum likelihood estimation from five well-known statistical model families are evaluated for the purpose of ultrasound tissue characterization. The fractal dimension (self-similarity measure) was used to characterize the spatial distribution of scatterers, while the Lacunarity (sparsity measure) was applied to determine scatterer number density. Performance was assessed based on 608 cross-sectional clinical ultrasound RF images of liver tumors (230 and 378 demonstrating respondent and non-respondent cases, respectively). Crossvalidation via leave-one-tumor-out and with different k-folds methodologies using a Bayesian classifier were employed for validation. The fractal properties of the backscattered echoes based on the Nakagami model (Nkg) and its extend four-parameter Nakagami-generalized inverse Gaussian (NIG) distribution achieved best results - with nearly similar performance - for characterizing liver tumor tissue. Accuracy, sensitivity and specificity for the Nkg/NIG were: 85.6%/86.3%, 94.0%/96.0%, and 73.0%/71.0%, respectively. Other statistical models, such as the Rician, Rayleigh, and K-distribution were found to not be as effective in characterizing the subtle changes in tissue texture as an indication of response to treatment. Employing the most relevant and practical statistical model could have potential consequences for the design of an early and effective clinical therapy., Comment: 31 pages, 7 figures, 3 tables, journal article

Heat and momentum transfer in wall-bounded turbulent flow, coupled with the effects of wall-roughness, is one of the outstanding questions in turbulence research. In the standard Rayleigh-B\'enard problem for natural thermal convection, it is notoriously difficult to reach the so-called ultimate regime in which the near-wall boundary layers are turbulent. Following the analyses proposed by Kraichnan [Phys. Fluids vol 5., pp. 1374-1389 (1962)] and Grossmann & Lohse [Phys. Fluids vol. 23, pp. 045108 (2011)], we instead utilize recent direct numerical simulations of forced convection over a rough wall in a minimal channel [MacDonald, Hutchins & Chung, J. Fluid Mech. vol. 861, pp. 138--162 (2019)] to directly study these turbulent boundary layers. We focus on the heat transport (in dimensionless form, the Nusselt number $Nu$) or equivalently the heat transfer coefficient (the Stanton number $C_h$). Extending the analyses of Kraichnan and Grossmann & Lohse, we assume logarithmic temperature profiles with a roughness-induced shift to predict an effective scaling of $Nu \sim Ra^{0.42}$, where $Ra$ is the dimensionless temperature difference, corresponding to $C_h \sim Re^{-0.16}$, where $Re$ is the centerline Reynolds number. This is pronouncedly different from the skin-friction coefficient $C_f$, which in the fully rough turbulent regime is independent of $Re$, due to the dominant pressure drag. In rough-wall turbulence the absence of the analog to pressure drag in the temperature advection equation is the origin for the very different scaling properties of the heat transfer as compared to the momentum transfer. This analysis suggests that, unlike momentum transfer, the asymptotic ultimate regime, where $Nu\sim Ra^{1/2}$, will never be reached for heat transfer at finite $Ra$., Comment: Accepted (In press) as a Rapid Communication in Physical Review Fluids

Basis tensor gauge theory is a vierbein analog reformulation of ordinary gauge theories in which the difference of local field degrees of freedom has the interpretation of an object similar to a Wilson line. Here we present a non-Abelian basis tensor gauge theory formalism. Unlike in the Abelian case, the map between the ordinary gauge field and the basis tensor gauge field is nonlinear. To test the formalism, we compute the beta function and the two-point function at the one-loop level in non-Abelian basis tensor gauge theory and show that it reproduces the well-known results from the usual formulation of non-Abelian gauge theory.

Highly turbulent Taylor-Couette flow with spanwise-varying roughness is investigated experimentally and numerically (direct numerical simulations (DNS) with an immersed boundary method (IBM)) to determine the effects of the spacing and axial width $s$ of the spanwise varying roughness on the total drag and {on} the flow structures. We apply sandgrain roughness, in the form of alternating {rough and smooth} bands to the inner cylinder. Numerically, the Taylor number is $\mathcal{O}(10^9)$ and the roughness width is varied between $0.47\leq \tilde{s}=s/d \leq 1.23$, where $d$ is the gap width. Experimentally, we explore $\text{Ta}=\mathcal{O}(10^{12})$ and $0.61\leq \tilde s \leq 3.74$. For both approaches the radius ratio is fixed at $\eta=r_i/r_o = 0.716$, with $r_i$ and $r_o$ the radius of the inner and outer cylinder respectively. We present how the global transport properties and the local flow structures depend on the boundary conditions set by the roughness spacing $\tilde{s}$. Both numerically and experimentally, we find a maximum in the angular momentum transport as function of $\tilde s$. This can be atributed to the re-arrangement of the large-scale structures triggered by the presence of the rough stripes, leading to correspondingly large-scale turbulent vortices., Comment: 20 pages, 7 figures, draft for JFM

Progress in roughness research, mapping any given roughness geometry to its fluid dynamic behaviour, has been hampered by the lack of accurate and direct measurements of skin-friction drag, especially in open systems. The Taylor--Couette (TC) system has the benefit of being a closed system, but its potential for characterizing irregular, realistic, 3-D roughness has not been previously considered in depth. Here, we present direct numerical simulations (DNSs) of TC turbulence with sand grain roughness mounted on the inner cylinder. The model proposed by Scotti (\textit{Phys. Fluids}, vol. 18, 031701, 2006) has been improved to simulate a random rough surface of monodisperse sand grains, which is characterized by the equivalent sand grain height $k_s$. Taylor numbers range from $Ta = 1.0\times 10^7$(corresponding to $Re_\tau = 82$) to $Ta = 1.0\times 10^9$($Re_\tau = 635$). We focus on the influence of the roughness height $k_s^+$ in the transitionally rough regime, through simulations of TC with rough surfaces, ranging from $k_s^+=5$ up to $k_s^+ = 92$, where the superscript `$+$' indicates non-dimensionalization in viscous units. We find that the downwards shift of the logarithmic layer, due to transitionally rough sand grains exhibits remarkably similar behavior to that of the Nikuradse (\textit{VDI-Forschungsheft} 361, 1933) data of sand grain roughness in pipe flow, regardless of the Taylor number dependent constants of the logarithmic layer.

Through a mechanism similar to perturbative particle scattering, particles of mass $m_\chi$ larger than the Hubble expansion rate $H_\mathrm{inf}$ during inflation can be gravitationally produced at the end of inflation without the exponential suppression powers of $\exp(-m_\chi/H_\mathrm{inf})$. Here we develop an analytic formalism for computing particle production for such massive particles. We apply our formalism to specific models that have been previously been studied only numerically, and we find that our analytical approximations reproduce those numerical estimates well., Comment: v2: 24 pages, 1 figure. Refs added. Clarified discussion of time scales at Eq. (6.11)

Previous numerical studies on homogeneous Rayleigh-B\'enard convection, which is Rayleigh-B\'enard convection (RBC) without walls, and therefore without boundary layers, have revealed a scaling regime that is consistent with theoretical predictions of bulk-dominated thermal convection. In this so-called asymptotic regime, previous studies have predicted that the Nusselt number ($Nu$) and the Reynolds number ($Re$) vary with the Rayleigh number ($Ra$) according to $Nu\sim Ra^{1/2}$ and $Re\sim Ra^{1/2}$ at small Prandtl number ($Pr$). In this study, we consider a flow that is similar to RBC but with the direction of temperature gradient perpendicular to gravity instead of parallel; we refer to this configuration as vertical natural convection (VC). Since the direction of the temperature gradient is different in VC, there is no exact relation for the average kinetic dissipation rate, which makes it necessary to explore alternative definitions for $Nu$, $Re$ and $Ra$ and to find physical arguments for closure, rather than making use of the exact relation between $Nu$ and the dissipation rates as in RBC. Once we remove the walls from VC to obtain the homogeneous setup, we find that the aforementioned $1/2$-power-law scaling is present, similar to the case of homogeneous RBC. When focussing on the bulk, we find that the Nusselt and Reynolds numbers in the bulk of VC too exhibit the $1/2$-power-law scaling. These results suggest that the $1/2$-power-law scaling may even be found at lower Rayleigh numbers if the appropriate quantities in the turbulent bulk flow are employed for the definitions of $Ra$, $Re$ and $Nu$. From a stability perspective, at low- to moderate-$Ra$, we find that the time-evolution of the Nusselt number for homogenous vertical natural convection is unsteady, which is consistent with the nature of the elevator modes reported in previous studies on homogeneous RBC.

In thermal convection for very large Rayleigh numbers ($Ra$), the thermal and viscous boundary layers (BL) undergo a transition from a classical state to an ultimate state. In the former state, the BL thicknesses follow a laminar-like Prandtl-Blasius-Polhausen scaling, whereas in the latter, the BLs are turbulent with log-corrections in the sense of Prandtl and von K\'arm\'an. Here, we report evidence of this transition via changes in the BL structure of vertical natural convection (VC), which is a buoyancy driven flow between differentially heated vertical walls. The dataset spans $Ra$-values from $10^5$ to $10^9$ and Prandtl number value of 0.709. For this $Ra$ range, the VC flow exhibits classical state behaviour in a global sense. Yet, with increasing $Ra$, we observe that near-wall higher-shear patches occupy increasingly larger fractions of the wall-areas, which suggest that the BLs are undergoing a transition from the classical state to the ultimate shear-dominated state. The presence of streaky structures-reminiscent of the near-wall streaks in canonical wall-bounded turbulence-further supports the notion of this transition. Within the higher-shear patches, conditionally averaged statistics yield a log-variation in the local mean temperature profiles, in agreement with the log-law of the wall for mean temperature, and a $Ra^{0.37}$ effective power-law scaling of the local Nusselt number, consistent with the logarithmically corrected 1/2-power law scaling predicted for ultimate thermal convection for very large $Ra$. Collectively, the results from this study indicate that turbulent and laminar-like BL coexist in VC at moderate to high $Ra$ and this transition from the classical state to the ultimate state manifests as increasingly larger shear-dominated patches, consistent with the findings reported for Rayleigh-B\'enard convection and Taylor-Couette flows.

Results from direct numerical simulations of vertical natural convection at Rayleigh numbers $1.0\times 10^5$-$1.0\times 10^9$ and Prandtl number $0.709$ support a generalised applicability of the Grossmann-Lohse (GL) theory, which was originally developed for horizontal natural (Rayleigh-B{\'e}nard) convection. In accordance with the GL theory, it is shown that the boundary-layer thicknesses of the velocity and temperature fields in vertical natural convection obey laminar-like Prandtl-Blasius-Pohlhausen scaling. Specifically, the normalised mean boundary-layer thicknesses scale with the $-1/2$-power of a wind-based Reynolds number, where the "wind" of the GL theory is interpreted as the maximum mean velocity. Away from the walls, the dissipation of the turbulent fluctuations, which can be interpreted as the "bulk" or "background" dissipation of the GL theory, is found to obey the Kolmogorov-Obukhov-Corrsin scaling for fully developed turbulence. In contrast to Rayleigh-B{\'e}nard convection, the direction of gravity in vertical natural convection is parallel to the mean flow. The orientation of this flow presents an added challenge because there no longer exists an exact relation that links the normalised global dissipations to the Nusselt, Rayleigh and Prandtl numbers. Nevertheless, we show that the unclosed term, namely the global-averaged buoyancy flux that produces the kinetic energy, also exhibits both laminar and turbulent scaling behaviours, consistent with the GL theory. The present results suggest that, similar to Rayleigh-B{\'e}nard convection, a pure power-law relationship between the Nusselt, Rayleigh and Prandtl numbers is not the best description for vertical natural convection and existing empirical relationships should be recalibrated to better reflect the underlying physics.

In turbulent Rayleigh-B\'{e}nard (RB) convection, a transition to the so-called ultimate regime, in which the boundary layers (BL) are of turbulent type, has been postulated. Indeed, at very large Rayleigh number $Ra \approx 10^{13}-10^{14}$ a transition in the scaling of the global Nusselt number $Nu$ (the dimensionless heat transfer) and the Reynolds number with $Ra$ has been observed in experiments and very recently in direct numerical simulations (DNS) of two-dimensional (2D) RB. In this paper we analyse the local scaling properties of the lateral temperature structure functions in the BLs of this simulation of 2D RB, employing extended self-similarity (ESS) (i.e., plotting the structure functions against each other, rather than only against the scale) in the spirit of the attached eddy hypothesis, as we have recently introduced for velocity structure functions in wall turbulence (Krug \textit{et al.}, \textit{J. Fluid Mech.}, vol. 830, 2017, pp. 797-819). We find no ESS scaling below the transition and in the near wall region. However, beyond the transition and for large enough wall distance $z^+ > 100$, we find clear ESS behaviour, as expected for a scalar in a turbulent boundary layer. In striking correspondence to the $Nu$ scaling, the ESS scaling region is negligible at $Ra = 10^{11}$ and well developed at $Ra = 10^{14}$, thus providing strong evidence that the observed transition in the global Nusselt number at $Ra \approx 10^{13}$ indeed is the transition from a laminar type BL to a turbulent type BL. Our results further show that the relative slopes for scalar structure functions in the ESS scaling regime are the same as for their velocity counterparts extending their previously established universality. The findings are confirmed by comparing to scalar structure functions in 3D turbulent channel flow., Comment: 10pages, 5 figures

Basis tensor gauge theory (BTGT) is a reformulation of ordinary gauge theory that is an analog of the vierbein formulation of gravity and is related to the Wilson line formulation. To match ordinary gauge theories coupled to matter, the BTGT formalism requires a continuous symmetry that we call the BTGT symmetry in addition to the ordinary gauge symmetry. After classically interpreting the BTGT symmetry, we construct using the BTGT formalism the Ward identities associated with the BTGT symmetry and the ordinary gauge symmetry. As a way of testing the quantum stability and the consistency of the Ward identities with a known regularization method, we explicitly renormalize the scalar QED at one-loop using dimensional regularization using the BTGT formalism., Comment: 23 pages

It is known that if the Peccei-Quinn symmetry breaking field is displaced from its minimum during inflation, the axion isocurvature spectrum is generically strongly blue tilted with a break transition to a flat spectrum. We fit this spectrum (incorporated into the "vanilla" $\Lambda$-CDM cosmological model) to the Planck and BOSS DR11 data and find a mild hint for the presence of axionic blue-tilted isocurvature perturbations. We find the best fit parameter region is consistent with all of the dark matter being composed of QCD axions in the context of inflationary cosmology with an expansion rate of order $10^{8}$ GeV, the axion decay constant of order $10^{13}$ GeV, and the initial misalignment angle of order unity. Intriguingly, isocurvature with a spectral break may at least partially explain the low-$\ell$ vs. high-$\ell$ anomalies seen in the CMB data., Comment: 25 pages