The paper presents an approximate algorithm for modeling a stationary discrete random process with marginal and bivariate distributions of its consecutive components in the form of a mixture of two Gaussian distributions. The algorithm is based on a combination of the conditional distribution method and the rejection method. An example of application of the proposed algorithm for simulating time series of daily maximum air temperatures is given.
The density of shock-compressed liquid krypton ρ ≈ 9 g/cm3 and the temperature T ≈ 55 000 K on the main Hugoniot adiabat in the pressure range of Р ≈ 140–255 GPa have been measured with hemispherical shock wave generators. Using spherical devices, gaseous krypton has been compressed to the density ρ ≈ 20 g/cm3 by the pressure Р ≈ 2700 GPa at the temperature T ≈ 120 000 K. The experiment has been carried out at the X-ray diffraction complex RGK-B-L (Russian Federal Nuclear Center All-Russian Scientific Research Institute of Experimental Physics) consisting of BIM 234.3000 betatrons with a boundary energy of 60 MeV used in the multipulse bremsstrahlung generation regime with a multichannel optoelectronic detection system for X-ray images. The designs of the experimental devices have been described, and the thermodynamic parameters of the krypton plasma reached in these devices have been estimated. The obtained data have been analyzed and compared to previously reported data. The compression of krypton to the density ρ ≈ 20 g/cm3 by the pressure Р ≈ 2700 GPa is currently a record achievement.
Two devices intended for copper cylindrical liner gasdynamic acceleration to velocities of 5–7 km/s using the chemicals explosion energy have been investigated. It has been demonstrated that the acceleration of quasi-isentropically and isentropically loaded liners under the conditions of high-level dynamics, symmetry of deposition, and suppression of shock-induced dusting is feasible.
Comparative experimental studies of the shock-induced particle ejection (“dusting”) from the free rough (Rz20) surface of copper and lead liners exposed to one or two successive shock waves separated in time by 0.2 μs have been carried out for the first time. This situation usually occurs in cumulative systems for the compression of the plasma by cylindrical or spherical liners shock or quasi-isentropically accelerated by explosion products. Using pulsed X-ray diffraction, laser optical recording, piezoelectric pressure sensors, and heterodyne interferometry, a qualitative picture has been studied and the quantitative characteristics of particle ejection from the free surface such as the velocities of the free surface and the particle flux front and the density (mass) distribution of particle flux in the direction of its motion, which are necessary for more accurate determination of features and the development of more appropriate models of the effect, have been evaluated.
The quasi-isentropic compressibility of a strongly nonideal helium plasma in the pressure range 250–600 GPa is experimentally studied in devices with cylindrical geometry. The temperature at the front of a cylindrical shock wave in helium ( T ≈ 10 000 K) and the flight speed of the inner cascade ( W ≈ 3.5 km/s), in the cavity of which the maximum compressed plasma density is achieved, are measured. Data on the compression of a nonideal helium plasma to a density ρ ≈ 3 g/cm 3 at an approximately constant final temperature of 21000 K are obtained. The trajectories of the metallic shells compressing the plasma are detected using high-power pulsed X-ray sources with a boundary electron energy of up to 60 MeV. The helium plasma density is determined using the radii of the shells measured at the time of their “stop.” The compressed plasma pressure is obtained using gasdynamic calculations. Comparative theoretical calculations of the quasi-isentropic compression parameters have been carried out using the following two theoretical models: the traditional chemical plasma model (SAHA code) and an ab initio quantum molecular dynamics (QMD) approach. No anomaly of the experimental data in the pressure range of the plasma phase transition theoretically assumed in helium is detected.
The quasi-isentropic compressibility of a strongly nonideal helium plasma in the pressure range 250–600 GPa is experimentally studied in devices with cylindrical geometry. The temperature at the front of a cylindrical shock wave in helium (T ≈ 10 000 K) and the flight speed of the inner cascade (W ≈ 3.5 km/s), in the cavity of which the maximum compressed plasma density is achieved, are measured. Data on the compression of a nonideal helium plasma to a density ρ ≈ 3 g/cm3 at an approximately constant final temperature of 21000 K are obtained. The trajectories of the metallic shells compressing the plasma are detected using high-power pulsed X-ray sources with a boundary electron energy of up to 60 MeV. The helium plasma density is determined using the radii of the shells measured at the time of their “stop.” The compressed plasma pressure is obtained using gasdynamic calculations. Comparative theoretical calculations of the quasi-isentropic compression parameters have been carried out using the following two theoretical models: the traditional chemical plasma model (SAHA code) and an ab initio quantum molecular dynamics (QMD) approach. No anomaly of the experimental data in the pressure range of the plasma phase transition theoretically assumed in helium is detected.
In the development of special explosion-proof chambers that must meet strict requirements for strength reliability, an important issue is the choice of the material of the load-bearing shell subjected to pulsed (dynamic and shock-wave) loads. As a rule, these structures are made from industrial low-alloy steel pipes of various standard sizes. This always raises the question of choosing the steel grade, especially at the stage of design-basis justification of their explosion resistance, since the dynamic strength characteristics of the pipe material are generally unknown. This paper is the first to present the results of analysis of the static, dynamic, and shock-wave compressive and tensile strengths of 17G1S, 09G2S, 10G2FBYu, and K60 strength class pipe steels. In addition, comparative data are given on the explosion resistance of pipes of 09G2S and 10G2FBYu steels at a strain rate of (2-5)· 10^2 s ^-1 .
Two devices intended for copper cylindrical liner gasdynamic acceleration to velocities of 5–7 km/s using the chemicals explosion energy have been investigated. It has been demonstrated that the acceleration of quasi-isentropically and isentropically loaded liners under the conditions of high-level dynamics, symmetry of deposition, and suppression of shock-induced dusting is feasible.
A new iterative method for modeling of non-Gaussian random vectors with given marginal distributions and a covariance matrix is proposed in this paper. The algorithm is compared with another iterative algorithm for modeling of non-Gaussian vectors, based on reordering of a sample of independent random variables with given marginal distributions. Our numerical studies show that both algorithms are equivalent in terms of the accuracy of reproduction of a given covariance matrix, but the offered algorithm turns out to be more efficient in terms of memory usage and, in many cases, is faster than the other one.
A nonideal helium plasma has been compressed to a density of ρ ≈ 14 g/cm 3 at a pressure of $$P \approx 20{\kern 1pt} $$ TPa (200 Mbar) in a spherical two-cascade device. The design of the device has been presented and the thermodynamic parameters of the helium plasma reached in it have been estimated. Data on the dynamics of the external and internal cascades in the device used to choose the calculation method have been obtained in a preliminary experiment with a hemispherical prototype. The experiment has been performed on the X-ray complex at the Russian Federal Nuclear Center All-Russian Research Institute of Experimental Physics, which includes BIM 234.3000 betatrons with a boundary energy of 60 MeV used in the multipulse bremsstrahlung generation regime with a multichannel optoelectronic system of recording X-ray images. A high-current linear accelerator LIU-R-T and an image detector assembly consisting of photochromic ADC screens have been used in addition to betatrons to detect the compressed shell cavity.
An Erratum to this paper has been published: https://doi.org/10.1134/S106377612230001X
We present data on the quasi-isentropic compressibility of a strongly coupled (nonideal) plasma of a deuterium–helium mixture in the range of pressures 150–250 GPa in devices with a cylindrical geometry and a new result obtained in an experiment with “pure” deuterium at a pressure ~200 GPa. The trajectory of the plasma-compressing metallic shells was recorded using powerful pulsed X-ray sources with a boundary electron energy up to 60 MeV. The densities of the plasma of deuterium and its mixture with helium were determined from the measured radii of the shells at the instant of their “stopping.” We derived the pressure of the compressed plasma based on gasdynamic computations including the real characteristics of the experimental devices. The data obtained confirm the conclusion previously reached at VNIIEF about a phase transition in the plasma of compressed deuterium in the range of pressures 150–160 GPa and suggest that this phase transition is retained in the mixture of deuterium with helium.
The electric response and processes occurring in piezoelectric elements subjected to pulsed pressure are usually described for two extreme cases: when the time of pressure variation is much longer than the time of propagation of an acoustic wave through a piezoelectric element, when the pressure in the piezoelectric element can be treated as constant (thin sensor regime), and under loading by a rectangular shock wave, when the shock front divides the piezoelectric element into the compressed and uncompressed zones (thick sensor regime). In the former case, the voltage across the piezoelectric element and electric fields emerging in it directly depend on the electric load, and the induced charge is proportional to the applied pressure (in the linear region). In the latter case, the current generated by the piezoelectric element is proportional to the shock load pressure, and the fields in the bulk of the piezoelectric element appear even in the case of short-circuiting of its electrodes. In this study, we consider the electric response of piezoelectrics to the action of pressure, which noticeably varies during the time commensurate with the time of propagation of the shock wave through the piezoelectric element. Such a situation takes place, for example, when a high-velocity particle flow formed when the shock wave emerging on the free surface of the metal plate (dusting) [1–4] is incident on the piezoelectric sensor. The calculations based on the constructed mathematical model show that under an increasing pressure on the piezoelectric element, electric fields with magnitudes depending on the rate of pressure increase and nonuniform over the piezoelectric element thickness appear in its bulk, and the induced electric charge is proportional to the average pressure in the piezoelectric element. Under certain conditions, these fields can attain values leading to the emergence of breakdown in the piezoelectric element and to distortion of generated signals. We consider experimentally observed cases of manifestation of breakdown effects under the action of rapidly increasing pressures in piezoelectric elements.