Wave turbulence buildup in a vibrating plate

We report experimental and numerical results on the buildup of the energy spectrum in wave turbulence of a vibrating thin elastic plate. Three steps are observed: first a short linear stage, then the turbulent spectrum is constructed by the propagation of a front in wave number space and finally a long time saturation due to the action of dissipation. The propagation of a front at the second step is compatible with scaling predictions from the Weak Turbulence Theory.     

On choosing a nonlinear initial iterate for solving the 2-D 3-T heat conduction equations

The 2-D 3-T heat conduction equations can be used to approximately describe the energy broadcast in materials and the energy swapping between electron and photon or ion. To solve the equations, a fully implicit finite volume scheme is often used as the discretization method. Because the energy diffusion and swapping coefficients have a strongly nonlinear dependence on the temperature, and some physical parameters are discontinuous across the interfaces between the materials, it is a challenge to solve the discretized nonlinear algebraic equations. Particularly, as time advances, the temperature varies so greatly in the front of energy that it is difficult to choose an effective initial iterate when the nonlinear algebraic equations are solved by an iterative method. In this paper, a method of choosing a nonlinear initial iterate is proposed for iterative solving this kind of nonlinear algebraic equations. Numerical results show the proposed initial iterate can improve the computational efficiency, and also the convergence behavior of the nonlinear iteration.     

Effects of high activation energies on acoustic timescale detonation initiation

Acoustic timescale Deflagration-to-Detonation Transition (DDT) has been shown to occur through the generation of compression waves emitted by a hot spot or reaction centre where the pressure and temperature increase with little diminution of density. In order to compensate for the multi-scale nature of the physico-chemical processes, previous numerical simulations in this area have been limited to relatively small activation energies. In this work, a computational study investigates the effect of increased activation energy on the time required to form a detonation wave and the change in behaviour of each hot spot as the activation energy is increased. The simulations use a localised spatially distributed thermal power deposition of limited duration into a finite volume of reactive gas to facilitate DDT. The Adaptive Wavelet-Collocation Method is used to solve efficiently the 1-D reactive Euler equations with one-step Arrhenius kinetics. The DDT process as described in previous work is characterised by the formation of hot spots during an initial transient period, explosion of the hot spots and creation of an accelerating reaction front that reaches the lead shock and forms an overdriven detonation wave. Current results indicate that as the activation energy is raised the chemical heat release becomes more temporally distributed. Hot spots that produce an accelerating reaction front with low activation energies change behaviour with increased activation energy so that no accelerating reaction front is created. An acoustic timescale ratio is defined that characterises the change in behaviour of each hot spot.     

On generating initial conditions for turbulence models: the case of Rayleigh–Taylor instability turbulent mixing

A new approach to generate initial conditions for RANS simulations of Rayleigh–Taylor (RT) turbulence is presented. The strategy is to provide profiles of turbulent model variables when it is suitable for the turbulence model to be started, and then use these profiles for the turbulence model initialization. The generation of turbulence model variable profiles is achieved with a two-step process. In the first step, a nonlinear modal model assuming small amplitude initial perturbations, incompressible and inviscid fluids is used to track the growth of modes that exist in a given initial perturbation spectrum, and also modes generated by mode interactions. The amplitude development of each mode represents the penetration distance of the light fluid into the heavy fluid (bubble penetration), for a given mode perturbation. The penetration distance of heavy fluid into the light fluid (spike penetration), for a given mode perturbation, is inferred from the bubble's height by an empirical relation valid for small initial amplitudes, and established by DNS simulations that depend on a nondimensional time, and the density contrast (Atwood number). It is hypothesized that the bubble front position of the RT mixing layer can be approximated by the largest penetration distance among all existing modes. The spike front position is approximated in the same fashion. The nonlinear model is evaluated by comparing the bubble front height evolution predicted by the model against the bubble front height predicted by an incompressible implicit large eddy simulations (ILES) code. Comparisons of results for “top-hat” and two-band initial perturbation spectra at Atwood numbers, A_T =0.3 and A_T =0.5 for the former, and A_T =0.01 and A_T =0.5 for the latter, show reasonable agreement. In the second step, the bubble and spike front positions, their derived velocities, and simplified profiles of the mixture fraction distribution of each fluid between the bubble and spike fronts are used with a two-fluid approximation to derive profiles for the turbulence model variables. When initialized with modal model profiles at start time τ₀, (i.e., the time when the turbulence model variable profiles are inferred from the modal model results), the RANS simulations provide at all times τ>τ₀ profiles that show good agreement with ILES simulations. The procedure for determining the time at which the RANS model should be started is a representative use, other parameters can be used depending on the application. In this paper, for the purpose of demonstration of the full strategy, τ₀ is taken as the time at which the mixing layer growth rate parameter α has reached its asymptotic value in the corresponding ILES simulation.     

Acceleration of ions with a nonadiabatic linearly polarized laser pulse

Acceleration of ions during interaction of a nonadiabatic laser pulse (i.e., a pulse with a sharp front) with a nanofilm is considered. If the amplitude of such a pulse is large enough, all electrons are removed from the target at the beginning of interaction and an energy of the most energetic ions follows approximately a parabolic law with time. Two main physical mechanisms limiting the maximal ions energy are identified and investigated in detail with the help of two-dimensional (2D) particle-in-cell (PIC) simulations. The first effect is a compensation of the ions charge due to the longitudinal return of the electrons to their initial position. The second effect is the compensation of the ions charge in the laser spot due to the transverse motion of the electrons from the periphery of the target. The theory for both effects is developed and a good agreement with the 2D PIC results is established. This theory allows predicting the optimal parameters for ions acceleration.     

Growth rates of the buoyancy-driven instability of an autocatalytic reaction front in a narrow cell

Experimental studies were performed on the buoyancy-driven instability of an autocatalytic reaction front in a quasi-2D cell. The unstable density stratification at an ascending front leads to convection that results in a fingerlike front deformation. The growth rates of the spatial modes of the instability are determined at the initial stage. A stabilization is found at higher wave numbers, while the system is unstable against low wave number perturbations. Whereas comparison with a reported model governed by Hele-Shaw flow fails, a two-dimensional Navier-Stokes model yields more satisfactory results. Still, present deviations suggest the presence of an additional mechanism that suppresses the growth.     

Dynamic behaviors of supersonic granular media under vertical vibration

We present experimental study of vibrofluidized granular materials by high speed photography. Statistical results present the averaged dynamic behaviors of granular materials in one cycle, including the variations of height, velocity and mechanical energy of the center of mass. Furthermore, time-space distribution of granular temperature which corresponds to the random kinetic energy shows that a temperature peak forms in the compression period and propagates upward with a steepened front. The Mach number in the steepened front is found to be greater than unity, indicating a shock propagating in the supersonic granular media.     

液电脉冲激波特性分析

基于相应的数学模型来表征液电脉冲激波的产生和传播过程,搭建了液电式碎岩综合试验平台,分析了典型的激波特性的仿真和实测结果。给出了不同充电电压下液电脉冲激波特性的仿真结果,并分析了充电电压对激波特性的影响。结果表明:充电电压为11 kV时,激波的压力峰值为2.67 MPa,激波能量为27.30 J,波前时间为2.16μs,激波加载速率为1.24 MPa/μs,电能转化为激波能量的效率为13.35%;提高电容充电电压,激波压力峰值和激波能量增大,波前时间减少,激波加载速率增加,但电能转化为激波能量的效率降低。利用建模分析的方法,可以根据放电回路参数预测液电脉冲激波特性,从而为进一步研究激波破碎岩石的形态和效果提供理论依据。     

Reflexionsversuche mit elektromagnetisch erzeugten Verdichtungsstößen hoher Mach-Zahlen

In a 80 cm long quartz-tube of 33 mm diameter plasma is generated by theta-pinch mechanism with conical coil and accelerated in the direction along the axis of the tube. A luminous front is observed and its velocity has been studied in, in dependence of coil distance and initial pressure. In Helium in a pressure range from 0·3 to 5 Torr it has been found by reflection experiments that before the luminous front a discontinuity in the pressure is moving ahead of the luminous front. The distanceΔx between this discontinuity and luminous front increases with increasing distance from the coil and with increasing initial pressure for instance up to 10 cm at a distance of about 50 cm and an initial pressure of 5 Torr. The results are compared to theoretical considerations.     

A new definition of boundary point between early reflections and late reverberation in room impulse responses

The early reflections in the room impulse response are usually defined as those observed within the initial 80 ms after the arrival of the direct sound, after which time the sound field is called reverberant. This number was chosen from measurements of other functions in a limited number of halls. In order to give an objective foundation to this time separation and to establish a physical indicator for it, a new method is proposed that defines a "transition time t(L)," which is the time at which the energy correlation between the direct plus initial sound and the subsequent decaying sound first achieves a specified low value. For various halls this number is shown and its relevance as a new parameter is discussed.     

超短脉冲超强激光与固体靶相互作用中高能离子的产生

 用2D3V PIC粒子模拟方法得到了超短脉冲超强激光与固体靶相互作用中高能离子产生的图像，并对其机理进行了研究。在靶前后表面都观察到了高能离子的产生，并诊断了离子能谱。模拟结果表明，在靶前表面所产生的高能离子，角分布较大，在向靶内输运过程中会损失能量；在靶后表面产生的高能离子，定向性很好，能获得很高的能量。模拟得到的离子能量和实验观测结果在量级上相符。     

Vacuum discharge instability at laser initiation of a cathode spot

The dynamics of fast (a current rise time of ≤10¹¹ A/s) laser-induced vacuum discharges with moderate amplitudes of the current (≤10 kA) and voltage (≤20 kV), as well as medium storage energy (20 J), is studied. It is shown experimentally that the initial conditions specified by the energy and duration of laser radiation are a decisive factor governing the discharge dynamics. Two types of beam-plasma instabilities separated in space and time are discovered, and their occurrence conditions are analyzed. The first type of instability, observed early in the discharge, is associated with pinch structures at the front of the cathode jet expanding into a vacuum. The second type arises either at the peak or at the trailing edge of the current pulse and is accompanied by generation of hard (with an energy of ≥100 keV) bremsstrahlung from the anode region. The increase of the hard component energy over the current source potential is attributed to breaking due to plasma erosion.     

中能重离子碰撞中前平衡核子发射的同位旋效应

利用同位旋相关的输运理论,研究了中能重离子碰撞中丰中子和缺中子碰撞系统在较宽能区范围内前平衡发射核子的同位旋效应.结果表明除低能区外,在100MeV/u以上核子–核子碰撞动力学效应起主要作用的能区,前平衡发射的中子–质子比仍然对介质中核子–核子碰撞截面的同位旋效应不灵敏,然而对于对称势的改变非常灵敏.故从低能量到较高能量的较大能区内,前平衡发射的中子–质子比的理论值与实验值的比较是提取对称势形式和强度的灵敏探针.     

X射线散射对面密度测量结果影响的数值模拟研究

定性分析了X射线散射对面密度测量结果的影响,给出了面密度测量结果偏差的评估表达式。利用蒙特卡罗方法数值模拟了不同照相布局、不同光源能量和不同材料类型情形下X射线散射对测量结果的影响。模拟结果表明:在无防护窗情形下,被测物自身散射将导致测量结果偏小;在有防护窗情形下,前窗散射对测量结果的影响不大,后窗散射对测量结果的影响与接收介质离被测物的距离有关,具体表现为距离越小影响越大;在一定光子能量范围内,光源能量越低散射导致的测量偏差越大;被测物的原子序数和密度越大,散射导致的测量偏差越大。     

On the propagation of a light pulse in the “transparency window” of a medium with population inversion

The propagation of a light pulse through a layer of a substance with population inversion and, correspondingly, with an anomalous dispersion in the transparency range is considered. It is shown that, when the wave packet approaches the layer of the substance with a strong anomalous dispersion, two additional wave packets arise at the rear boundary of the layer. One of these packets has a negative energy and moves in the layer from its rear to the front boundary, while the other packet moves away from the rear boundary of the layer into a vacuum. As the initial wave packet comes into contact with the front boundary of the layer of the substance, it meets the packet of the negative energy. As a result, the two packets annihilate one another, after which only one packet remains, which moves away from the layer.     

神龙一号直线感应加速器X光转换靶破坏诊断

 利用能量约450 keV、焦斑直径1~4 mm的低能X光对神龙一号直线感应加速器束靶作用后钽靶的破坏进行诊断,利用增强型电荷耦合器件（ICCD）对诊断过程记录,得到束靶作用后数μs时间内钽靶材料密度的变化。结果表明：在束靶作用后约1 μs内靶材料密度基本没有变化,且该时间段内ICCD相机没有观察到有靶前钽靶材料的微粒喷射。     

Detailed measurements of local scalar front structures in stagnation-type turbulent premixed flames

Local scalar front structures of OH mole fraction, reaction progress variable, and its three-dimensional gradient have been measured in stagnation-type turbulent premixed flames. The reaction progress variable front is observed to change with increasing turbulence from parallel iso-scalar contours but reduced progress variable gradients, called the lamella-like front, to disrupted non-parallel iso-contours that deviate substantially from those of wrinkled laminar flamelets, called the non-flamelet front. This transition is attributed to the different scales of interaction between the flame internal structure and a spectrum of turbulence extending from the integral scale to the Kolmogorov scale. The lamella-like front pattern occurs when the length scales of interaction are smaller than the laminar flame thickness but the time scales are greater than the flame residence time. The non-flamelet front pattern occurs when the length scales of interaction are greater than the laminar flame thickness but the time scales are smaller than the flame residence time. This difference corresponds to the change of combustion regime from complex-strain flame front to turbulent flame front on a revised regime diagram. A correlation is also proposed for the turbulent flame brush thickness as a function of turbulent Reynolds number and heat release parameter. The heat release parameter is considered to arise from the non-passive effects of flame-surface wrinkling.     

Numerical experiments of the planar equal-mass three-body problem: the effects of rotation

The planar three-body problem with angular momentum is numerically and systematically studied as a generalization of the free-fall problem (i.e., the three-body problem with zero initial velocities). The initial conditions in the configuration space exhaust all possible forms of a triangle, whereas the initial conditions in the momentum space are chosen so that position vectors and momentum vectors are orthogonal. Numerical results are organized according to the value of virial ratio k defined as the ratio of the total kinetic energy to the total potential energy. Final motions are mapped in the initial value space. Several interesting features are found. Among others, binary collision curves seem to spiral into the Lagrange point, and for large k, binary collision curves connect the Lagrange point and the Euler point. The existence of a lunar periodic orbit and a periodic orbit of petal-type is suggested. The number of escape orbits as a function of the escape time is analyzed for different k. The behavior of this number for different time and k shows most remarkably the effects of rotation of triple systems. The number of escape orbits increases exponentially for k相似文献   

飞秒激光-金属薄膜靶相互作用中靶前后超热电子能谱的比较

采用电子谱仪测量了飞秒激光-金属薄膜靶相互作用中靶前和靶后产生的超热电子能谱.结果显示：靶前超热电子能谱的峰出现在约430 keV处，靶后超热电子能谱的峰出现在约175 keV处；靶前超热电子的有效温度分别为218 keV和425 keV，靶后超热电子能谱出现“软化”现象，其有效温度分别为96 keV和347 keV.靶前和靶后超热电子能谱明显不同是由于超热电子输运穿越过密等离子体和冷材料的靶，并在靶后建立Debye鞘，鞘电场使靶后超热电子能谱峰向低能端移动，鞘电场和自生磁场导致靶后超热电子能谱产生“软化”，估算出的鞘电场小于激光电场.