Local measurement of nonclassical ion heating during magnetic reconnection

Local ion temperature and flows are measured directly in the well-characterized reconnection layer of a laboratory plasma. The measurements indicate strongly that ions are heated due to reconnection and that more than half of the reconnected field energy is converted to ion thermal energy. Neither classical viscous damping of the observed sub-Alfvenic ion flows nor classical energy exchange with electrons is sufficient to account for the ion heating, suggesting the importance of nonclassical dissipation mechanisms in the reconnection layer.     

On energy transfer and generation of an electric current in the vicinity of an electrode dipped in an electrolyte and strongly heated by a current passing through it

Processes occurring when a metal electrode dipped in an electrolyte is heated by intense evaporation of the electrolyte are considered in terms of a physically rigorous model. Based on the Onsager principle of least energy dissipation rate in nonequilibrium processes, the fractions of thermal energy that are spent on heating and evaporating the electrolyte and on heating the vapor are found. The energy is released within the vapor-gas sheath when an electric current flows between the electrode and electrolyte surface. It is found that the electrolyte vapor temperature exceeds 1300 K. Analytical expressions are derived for the vapor-gas sheath thickness, the electrolyte vapor pressure, and the velocity of the vapor escaping the discharge zone. It is shown that field evaporation of thermally activated negative ions from the electrolyte surface cannot provide an electric current with densities found in experiments but is responsible for the generation of free electrons near the electrolyte surface. These electrons arise when the ions decay via collisions with excited molecules.     

Conditional analysis of temperature and strain rate effects on dissipation structure in turbulent non-premixed jet flames

This work presents results from simultaneous high-resolution temperature and velocity measurements in a series of turbulent non-premixed jet flames. The filtered Rayleigh scattering (FRS)-based temperature measurements demonstrate sufficient signal-to-noise (SNR) and spatial resolution to estimate the smallest scalar length scales and accurately determine dissipation rate fields. A comprehensive set of conditional statistics are used to characterize the small-scale structure, including the dependence of dissipation layer widths on Reynolds number, temperature, and dissipation magnitude. In general, the dissipation layer thickness decrease with increasing Reynolds number and increase with increasing temperature. However, dissipation layer widths show two distinct behaviors with respect to dissipation magnitude. For small dissipation values, increases in magnitude results in broadening of the dissipation layer, while for larger magnitude values of dissipation, the layer widths are thinned, highlighting the complexity of small-scale turbulent mixing. Additionally, measured ratios of the dissipation layer width to the Batchelor length scale are consistent across all Reynolds numbers and agree with previous studies in non-reacting flows. The unique aspect about the current set of measurements is the ability to examine the interaction of dissipation structure with turbulent flow parameters for the first time in turbulent non-premixed flames. Particularly, the strain rate/dissipation relationship is examined and compared to previous studies in non-reacting flows. It is found that the dissipation layers tend to align normal to the principal compressive strain axis and this tendency increases with increasing Reynolds number. For the lowest Reynolds number case, no dependence of the dissipation layer width nor dissipation rate magnitude on strain rate is found. However, for higher Reynolds numbers, a strong dependence of the dissipation layer width and dissipation rate magnitude on the principal compressive strain rate is observed. These results indicate the direct role of the compressive strain rate field on small-scale mixing structure in reacting flows.     

Anomalous Dissipation and Energy Cascade in 3D Inviscid Flows

Adopting the setting for the study of existence and scale locality of the energy cascade in 3D viscous flows in physical space recently introduced by the authors to 3D inviscid flows, it is shown that the anomalous dissipation is – in the case of decaying turbulence – indeed capable of triggering the cascade which then continues ad infinitum, confirming Onsager’s predictions.     

Three-dimensional modelling of dissipative quantum transport in quantum dots and atomistic scale devices using nonHermitian generalized potentials

In this paper we introduce a phenomenology for inserting dissipation into the single-particle Schrödinger equation for carrier transport by utilizing appropriate nonHermitian additions to the Hamiltonian. The nonHermitian terms are determined by incorporating model particle trapping/de-trapping, momentum gain/loss, energy gain/loss into the quantum continuity equations derived within the Bohm picture and then reconstructing the full Hamiltonian by reversing the Bohm projection. The new phenomenology is designed to obtain quantum velocity flows using the Bohm projection of solutions to the nonHermitian Schrödinger equation for applications in 2D and 3D quantum dots and mesoscopic MOSFETs. For this purpose we introduce a novel fast algorithm to compute the wave function in 2D and 3D based on a two time step iteration and direct integration.     

Statistics of the turbulent/non-turbulent interface in a spatially developing mixing layer

The thin interface separating the inner turbulent region from the outer irrotational fluid is analysed in a direct numerical simulation of a spatially developing turbulent mixing layer. A vorticity threshold is defined to detect the interface separating the turbulent from the non-turbulent regions of the flow, and to calculate statistics conditioned on the distance from this interface. The conditional statistics for velocity are in remarkable agreement with the results for other free shear flows available in the literature, such as turbulent jets and wakes. In addition, an analysis of the passive scalar field in the vicinity of the interface is presented. It is shown that the scalar has a jump at the interface, even stronger than that observed for velocity. The strong jump for the scalar has been observed before in the case of high Schmidt number (Sc). In the present study, such a strong jump is observed for a scalar with Sc ≈ 1. Conditional statistics of kinetic energy and scalar dissipation are presented. While the kinetic energy dissipation has its maximum far from the interface, the scalar dissipation is characterised by a strong peak very close to the interface. Finally, it is shown that the geometric features of the interfaces correlate with relatively large scale structures as visualised by low-pressure isosurfaces.     

Intermittency of energy in rapid granular shear flows

Hard-disk simulations are used for two-dimensional rapid granular shear flows of circular disks between two rotating cylinders. The intermittency effects associated with the rate of the energy dissipation of collisions are studied. The statistics of intermittent signals of energy dissipation reveals that a power law governs the dynamics of rapid shear granular flows. A dynamical system approach based on the Gledzer-Ohkitani-Yamada shell model of turbulence is employed to reproduce signals for energy dissipation that are statistically consistent with those from simulations. The results suggest that rapid granular flows can be analyzed by appropriate turbulent models.     

Three-dimensional finite-difference lattice Boltzmann model and its application to inviscid compressible flows with shock waves

In this paper, a three-dimensional (3D) finite-difference lattice Boltzmann model for simulating compressible flows with shock waves is developed in the framework of the double-distribution-function approach. In the model, a density distribution function is adopted to model the flow field, while a total energy distribution function is adopted to model the temperature field. The discrete equilibrium density and total energy distribution functions are derived from the Hermite expansions of the continuous equilibrium distribution functions. The discrete velocity set is obtained by choosing the abscissae of a suitable Gauss–Hermite quadrature with sufficient accuracy. In order to capture the shock waves in compressible flows and improve the numerical accuracy and stability, an implicit–explicit finite-difference numerical technique based on the total variation diminishing flux limitation is introduced to solve the discrete kinetic equations. The model is tested by numerical simulations of some typical compressible flows with shock waves ranging from 1D to 3D. The numerical results are found to be in good agreement with the analytical solutions and/or other numerical results reported in the literature.     

Quantization of 2D hole gas in conductive hydrogenated diamond surfaces observed by electron field emission

Discrete jumps are observed in the emitted current density (J) versus extraction electric field (E) curves in electron field emission measurements from a conductive, hydrogen-terminated air-exposed diamond surface. These jumps are well reproduced by computations based on the assumption that a 2D nanoscale quantum system with discrete energy levels exists in the diamond near-surface layer. The present results confirm the formation of well-defined quantum states of holes in the 2D surface layer present on hydrogenated air-exposed diamond surfaces.     

基于总能形式的耦合的双分布函数热晶格玻尔兹曼数值方法

双分布函数热晶格玻尔兹曼数值方法在微尺度热流动系统中得到广泛的应用. 本文基于晶格玻尔兹曼平衡分布函数低阶Hermite展开式, 创新性地提出了包含黏性热耗散和压缩功的耦合的双分布函数热晶格玻尔兹曼数值方法, 将能量场内温度的变化以动量源的形式引入晶格波尔兹曼动量演化方程, 实现了能量场与动量场之间的耦合. 研究了考虑黏性热耗散和压缩功的和不考虑的两种热自然对流模型, 重点分析了不同瑞利数和普朗特数下流场内的流动情况以及温度、速度和平均努赛尔数的变化趋势. 本文实验结果与文献结果一致, 验证了本文数值方法的可行性和准确性. 研究结果表明: 随着瑞利数和普朗特数的增大, 方腔内对流传热作用逐渐增强, 边界处形成明显的边界层; 考虑黏性热耗散和压缩功的模型对流作用相对增强, 黏性热耗散和压缩功对自然对流的影响在微尺度流动过程中不能忽略.     

Dissipative interaction and anomalous surface absorption of bulk phonons at a two-dimensional defect in a solid

We predict an extreme sensitivity to the dissipative losses of the resonant interaction of bulk phonons with a 2D defect in a solid. We show that the total resonant reflection of the transverse phonon at the 2D defect, described earlier without an account for dissipation, occurs only in the limit of extremely weak dissipation and is changed into almost total transmission by relatively weak bulk absorption. Anomalous surface absorption of the transverse phonon, when one half of the incident acoustic energy is absorbed at the 2D defect, is predicted for the case of “intermediate” bulk dissipation.     

Giant enhancement of noncontact friction between closely spaced bodies by dielectric films and two-dimensional systems

The effect of an external bias voltage and spatial variations of the surface potential on the damping of cantilever vibrations in an atomic force microscope (AFM) is considered. The damping is due to an electrostatic friction that arises due to dissipation of the energy of an electromagnetic field generated in the sample by oscillating static charges induced on the surface of the AFM probe tip by the bias voltage or spatial variations of the surface potential. A similar effect appears when the tip is oscillating in an electrostatic field created by charged defects present in the dielectric sample. The electrostatic friction is compared to the van der Waals (vdW) friction between closely spaced bodies, which is caused by a fluctuating electromagnetic field related to the quantum and thermal fluctuations of current density inside the bodies. It is shown that the electrostatic friction and the vdW friction can be strongly enhanced in the presence of dielectric films or two-dimensional (2D) structures—such as a 2D electron system or an incommensurate layer of adsorbed ions exhibiting acoustic oscillations—on the probe tip and sample surfaces. It is also shown that the damping of cantilever oscillations caused by the electrostatic friction in the presence of such 2D structures can have the same order of magnitude and the same dependence on the distance as observed in experiment by Stipe et al. [Phys. Rev. Lett. 87, 096801 (2001)]. At small distances, the vdW friction can be large enough to be measured in experiment. In interpreting the experimental data that obey a quadratic dependence on the bias voltage, one can reject a phonon mechanism according to which the friction depends on the fourth power of the voltage.     

Inverse energy cascade in three-dimensional isotropic turbulence

We study the statistical properties of homogeneous and isotropic three-dimensional (3D) turbulent flows. By introducing a novel way to make numerical investigations of Navier-Stokes equations, we show that all 3D flows in nature possess a subset of nonlinear evolution leading to a reverse energy transfer: from small to large scales. Up to now, such an inverse cascade was only observed in flows under strong rotation and in quasi-two-dimensional geometries under strong confinement. We show here that energy flux is always reversed when mirror symmetry is broken, leading to a distribution of helicity in the system with a well-defined sign at all wave numbers. Our findings broaden the range of flows where the inverse energy cascade may be detected and rationalize the role played by helicity in the energy transfer process, showing that both 2D and 3D properties naturally coexist in all flows in nature. The unconventional numerical methodology here proposed, based on a Galerkin decimation of helical Fourier modes, paves the road for future studies on the influence of helicity on small-scale intermittency and the nature of the nonlinear interaction in magnetohydrodynamics.     

Multipolar Structure of Equilibrium Shear Flow Field in Toroidal Plasmas

The multipolar velocity field structures are investigated by 2D momentum conservation equation with 3D equilibrium sheared flows in the full toroidal system. Numerical results show that the non-existence of radial velocity field in equilibrium surfaces is suitable only for the zero-order term of our 2D simulation. The non-zero-order radial velocity field is still preserved, even when converted to conventional magnetic surface coordinates. The distribution of velocity field vectors of the order of 1, 2, and 3 are presented respectively in 2, 4, and 6 polar fields with the local vortex structure. The excitation mechanisms of these velocity vortices are the coupling effects of the magneto-fluid structure patterns and the toroidal effects. These results can help us understand the complexity of core physics, the transverse transport across magnetic field by the radial plasma flow and the formation of velocity vortices.     

A robust numerical scheme for highly compressible magnetohydrodynamics: Nonlinear stability,implementation and tests

The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate that the theoretical stability properties are valid in practice with negligible compromises to accuracy. The result is a highly robust scheme with state-of-the-art efficiency. The scheme’s robustness is due to entropy stability, positivity and properly discretised Powell terms. The implementation takes the form of a modification of the MHD module in the FLASH code, an adaptive mesh refinement code. We compare the new scheme with the standard FLASH implementation for MHD. Results show comparable accuracy to standard FLASH with the Roe solver, but highly improved efficiency and stability, particularly for high Mach number flows and low plasma β. The tests include 1D shock tubes, 2D instabilities and highly supersonic, 3D turbulence. We consider turbulent flows with RMS sonic Mach numbers up to 10, typical of gas flows in the interstellar medium. We investigate both strong initial magnetic fields and magnetic field amplification by the turbulent dynamo from extremely high plasma β. The energy spectra show a reasonable decrease in dissipation with grid refinement, and at a resolution of 512³ grid cells we identify a narrow inertial range with the expected power law scaling. The turbulent dynamo exhibits exponential growth of magnetic pressure, with the growth rate higher from solenoidal forcing than from compressive forcing. Two versions of the new scheme are presented, using relaxation-based 3-wave and 5-wave approximate Riemann solvers, respectively. The 5-wave solver is more accurate in some cases, and its computational cost is close to the 3-wave solver.     

On nonlinear periodic capillary-fluctuation wave flow in a thin liquid film on a solid substrate

Analytical calculation of a nonlinear periodic wave flow on the free surface of a charged layer of an ideal incompressible conducting liquid resting on a solid substrate is carried out for the case when fluctuation-induced forces (the dispersion component of the wedging pressure) have a decisive effect on the system. It is shown that wave flows emerge in the liquid in calculations of the second order of smallness in the wave amplitude, which is assumed to be small compared with the thickness of the liquid layer. These flows result from nonlinear interaction as nonlinear corrections to the waves set at the zero time. The field of fluctuation-induced forces displaces these flows toward the periphery of the area of influence of these forces. This effect takes place both in the presence of an external electric field near the free surface and in its absence. The sign and value of the nonlinear corrections depend on whether an electric field is present near the free surface of the liquid. In the presence of an electric field, the curvature of the crest of the nonlinear waves increases; in its absence, the curvature decreases.     

Electronic structure of Ag adsorbed on Si(111); experiment and theory

Experimental results (low energy electron loss spectroscopy) and band structure calculations relating to the early stages of Ag growth on a Si(111) surface are presented. Crystallography and thermal desorption kinetics studies of this interface, previously published, gave rise to the following conclusions. At room temperature and below 200°C, two-dimensional (2D) (111) epitaxial layers develop on top of a first ordered layer (√3 × √3), while at higher temperatures three-dimensional (3D) clusters develop over this first layer. Low energy electron loss experiments were performed at various surface coverages θ. They display different evolutions according to the growth modes. For the 2D epitaxial growth, one observes the disappearance of the peaks characteristic of a Si surface and the onset of Ag induced peaks located at 7.1 and 4.6 eV at completion of the √3 layer. These peaks narrow and shift to the bulk Ag excitation energies at 7.5 and 4 eV when a second Ag layer is deposited. In order to explain these results, we present a theoretical calculation of the electronic density of states of the interface using a tight binding approximation. This calculation accounts for the development of the Ag d band from the √3 coverage range to the (111) epitaxial Ag planes. The evolution of the spectra when θ is increased is discussed in view of these results.     

A turbulence characteristic length scale for compressible flows

The current RANS models are generally established and calibrated under incompressible condition and these kinds of models could succeed in predicting many features of incompressible flows. However, these models extended to the high-speed, compressible flows are always less accurate. In the paper, a compressible von Kármán length scale is proposed for compressible flows considering the variable densities. It contains no empirical coefficients and is based on phenomenological theory. In the turbulent kinetic equation, the extra unclosed terms induced by non-constant densities are treated as dissipation terms and the equation is closed algebraically via the introduction of the von Kármán length scale. The original and the proposed von Kármán length scale lead to two different kinds of SAS (scale adaption simulation) models, KDO (turbulence kinetic energy dependent only) and CKDO (compressible KDO), respectively. Compressible mixing layer with significant compressibility is studied within standard k–?, k–ω, KDO turbulence models and their compressible versions. The compressibility effects such as the reduced mixing layer thickness, growth rate and turbulence intensity can be reproduced by CKDO. The new length scale can improve the performances of the model in predicting the mixing layer thickness, stream-wise velocity and Reynolds shear stresses when the convective Mach number is 0.8. Besides, the new length scale also leads to accurate computed growth rate when the convective Mach number ranges from 0.1 to 1.0.     

Temperature gradient spectra and temperature dissipation rate in a turbulent convective flow

Determining mixing coefficients in oceanographic flows relies on the form of temperature gradient spectra in turbulent water flows at large wavenumbers. Several recent investigations concluded that these spectra are best described by the functional form proposed by Kraichnan rather than by the Batchelor form, more commonly used in oceanography. In this study, we provide additional support for this conclusion using laboratory measurements of the temperature gradient spectra in a Rayleigh–Bénard convective cell, in order to avoid difficulties inherent in oceanographic field measurements. The range of Rayleigh numbers in experiments is between Ra = 3×10⁷ and Ra = 5×10⁹. In addition to a traditional method of traversing thermistors, a novel optical technique recently introduced for oceanic measurements was used to obtain the spectra; comparison between these two methods serves as a validation test for the new optical technique. The temperature measurements were also augmented by 2D particle image velocimetry (PIV) observations. The measured dependence of the Nusselt number on the Rayleigh number followed Nu ∝ Ra^0.29 at Pr = 6 and was consistent with the literature data. We observed the temperature dissipation rate to vary by an order of magnitude over a horizontal transect at Ra > 10⁹. The temperature dissipation spectra obtained by both methods were in agreement over the Ra interval considered. The location of the temperature dissipation peaks was also consistent with PIV measured energy dissipation rates. Our data suggest increasing importance of top/bottom boundaries for the momentum and the temperature dissipation with increasing Ra number. Applied to oceanic upper ocean convection, our results imply that most of the dissipation occurs close to the air–sea boundary. Thus, attempts to parameterise or measure air–sea turbulent convective fluxes have to reflect the dominant role of near boundary dissipation at large Ra.     

Formation of relief on a metallic target surface under the action of compression-plasma flows

The results of experimental and theoretical investigations of relief formation on the surface of a steel target (grade St 3 steel, GOST (State Standard) 380) during treatment by compression-plasma flows are represented. The density of energy absorbed by the target varied in the range from 15 to 25 J/cm², the pulse duration was 100 μs, and the pulse number was N = 1, 3, 5, 7. The experiment revealed the expansion of boundaries of the central area (the area on which the plasma flow is incident normally to the surface) with increasing pulse number. This is explained by the more uniform surface treatment at a greater pulse number. It is shown that to describe relief formation in the central area there is a need to take into account the pressure of the plasma flow on the target surface, apart from surface tension forces and energy dissipation due to viscosity.