柴油、甲醇和水三组元乳化液滴微爆过程的研究

分别采用激光全息摄影技术和高速数字摄影技术观察了柴油、甲醇和水乳化液喷雾在高温高压(773K,3.1MPa)环境中发生微爆现象的瞬间和全过程,证实了微爆现象的存在.由于微爆机理的复杂性,尚难以用数学方法准确描述该过程.实验分析表明:若环境温度处于"最佳温度"范围内,乳化液滴表面首先形成"无水层",液滴内部形成一个水滴的概率很小,可能形成几个相对较大的水滴,只要其中一个较大水滴的蒸汽压力大于液滴的表面张力和环境压力之和,液滴就有可能发生微爆,微爆不仅与液滴直径、组分的质量分数和组分间的沸点差等乳化液的本身特性有关,而且环境温度和压力的影响也不容忽视.该研究可以为乳化液喷雾微爆过程的数学模拟提供参考.

岩石爆炸增渗模型实验及DEM数值模拟研究

岩体中爆炸提高矿石的渗透性可以极大地提高采矿效率，是碎裂岩型矿床预裂浸出法的关键技术。混凝土与岩体都具有脆性材料的特性，因此用浇筑于铁桶中的混凝土试件进行模型实验研究岩体的爆炸增渗效果。铁桶可以提高实验效率，使边界条件更为简单和易于操作，但是，实验和原型的差异需要论证。用自主开发的基于连续介质力学模型的离散元方法，模拟了有铁桶边界的模型实验，验证了数值模拟方法的可行性及有效性，给出了有弹性侧限边界约束的混凝土和较大尺度无反射边界条件的岩石中的爆炸差别。在此基础上，分析了岩石中爆炸造成的岩石破坏规律。计算结果表明：岩石破坏面总面积和破坏区的最大裂缝宽度受药量和岩石的抗拉强度影响，破坏面总面积和裂缝宽度随药量增加而增大，随抗拉强度增大而减小；在药量相同的条件下，实际岩石环境下的岩块破坏程度比铁桶约束的大。模型实验和数值模拟相结合的办法可以对混凝土和岩石的爆炸破坏给出较为合理的结果。

Studies of ion energy and yield from argon clusters heated by intense femtosecond laser pulses

Using time-of-flight spectrometry, the interaction of intense femtosecond laser pulses with argon clusters has been studied by measuring the energy and yield of emitted ions. With two different supersonic nozzles, the dependence of average ion energy (E) over bar on cluster size (n) over bar in a large range of (n) over bar approximate to 3 x 10(3) similar to 3 x 10(6) has been measured. The experimental results indicate that when the cluster size (n) over bar <= 3 x 10(5), the average ion energy (E) over bar proportional to (n) over bar (0.5), Coulomb explosion is the dominant expansion mechanism. Beyond this size, the average ion energy gets saturated gradually, the clusters exhibit a mixed Coulomb-hydrodynamic expansion behavior. We also find that with the increasing gas backing pressure, there is a maximum ion yield, the ion yield decreases as the gas backing pressure is further increased.

Neutron source from thin foil confined by two laser pulses

Neutron production from a thin deuterium-tritium (D-T) foil irradiated by two intense femtosecond laser pulses from opposite sides with zero phase difference is studied analytically and numerically. For the interaction of a laser pulse of amplitude a = 7, focal area 100 mu m(2) and areal density 4.4 x 10(18) cm(-2) with a D-T plasma foil, about 1.17 x 10(21) neutron s(-1) can be obtained, much more than from other methods. The profiles of the ion and electron densities are also calculated.

A plasma channel scheme for the interaction of an intense femtosecond laser pulse with a deuterium cluster jet

We propose a plasma channel scheme to obtain an improved table-top laser driven fusion neutron yield as a result of explosions of large deuterium clusters irradiated by an intense laser pulse. A cylindrical plasma channel is created by two moderate intensity laser prepulses at the edge of a deuterium cluster jet along which an intense main laser pulse propagates several nanoseconds later. With the aid of this plasma channel, the main laser pulse will be allowed to deposit its energy into the central region of the deuterium gas jet where the cluster sizes are larger and the atomic density is higher. The plasma channel formation and its impact on the deuterium ion energy spectrum and the consequent fusion neutron yield have been investigated. The calculated results show that a remarkable increase of the table-top laser driven fusion neutron yield would be expected.

Effect of laser spot size on fusion neutron yield in laser-deuterium cluster interactions

The effect of the laser spot size on the neutron yield of table-top nuclear fusion from explosions of a femtosecond intense laser pulse heated deuterium clusters is investigated by using a simplified model, in which the cluster size distribution and the energy attenuation of the laser as it propagates through the cluster jet are taken into account. It has been found that there exists a proper laser spot size for the maximum fusion neutron yield for a given laser pulse and a specific deuterium gas cluster jet. The proper spot size, which is dependent on the laser parameters and the cluster jet parameters, has been calculated and compared with the available experimental data. A reasonable agreement between the calculated results and the published experimental results is found.

甲烷团簇纯库仑爆炸的离子平均动能

利用一简单的经典静电模型研究了甲烷团簇纯库仑爆炸情况下产生的离子平均动能。研究表明，甲烷团簇爆炸后离子获得的平均动能和离子的初始平均静电势能的比值，与团簇的尺寸大小无关。这意味着在甲烷团簇纯库仑爆炸近似下，不必使用分子动力学模拟，离子获得的动能可以通过这一比值以及离子的初始静电势能进行估算。给出了不同碳离子价态下的离子平均动能和其初始平均静电势能的比值。

Part I - The radiation of elastic waves from a spherical cavity in a half space. Part II - Precision determination of focal depths and epicenters of earthquakes.

Part I: The dynamic response of an elastic half space to an explosion in a buried spherical cavity is investigated by two methods. The first is implicit, and the final expressions for the displacements at the free surface are given as a series of spherical wave functions whose coefficients are solutions of an infinite set of linear equations. The second method is based on Schwarz's technique to solve boundary value problems, and leads to an iterative solution, starting with the known expression for the point source in a half space as first term. The iterative series is transformed into a system of two integral equations, and into an equivalent set of linear equations. In this way, a dual interpretation of the physical phenomena is achieved. The systems are treated numerically and the Rayleigh wave part of the displacements is given in the frequency domain. Several comparisons with simpler cases are analyzed to show the effect of the cavity radius-depth ratio on the spectra of the displacements.

Part II: A high speed, large capacity, hypocenter location program has been written for an IBM 7094 computer. Important modifications to the standard method of least squares have been incorporated in it. Among them are a new way to obtain the depth of shocks from the normal equations, and the computation of variable travel times for the local shocks in order to account automatically for crustal variations. The multiregional travel times, largely based upon the investigations of the United States Geological Survey, are confronted with actual traverses to test their validity.

It is shown that several crustal phases provide control enough to obtain good solutions in depth for nuclear explosions, though not all the recording stations are in the region where crustal corrections are considered. The use of the European travel times, to locate the French nuclear explosion of May 1962 in the Sahara, proved to be more adequate than previous work.

A simpler program, with manual crustal corrections, is used to process the Kern County series of aftershocks, and a clearer picture of tectonic mechanism of the White Wolf fault is obtained.

Shocks in the California region are processed automatically and statistical frequency-depth and energy depth curves are discussed in relation to the tectonics of the area.

Regional variations in upper mantle structure beneath North America

Several types of seismological data, including surface wave group and phase velocities, travel times from large explosions, and teleseismic travel time anomalies, have indicated that there are significant regional variations in the upper few hundred kilometers of the mantle beneath continental areas. Body wave travel times and amplitudes from large chemical and nuclear explosions are used in this study to delineate the details of these variations beneath North America.

As a preliminary step in this study, theoretical P wave travel times, apparent velocities, and amplitudes have been calculated for a number of proposed upper mantle models, those of Gutenberg, Jeffreys, Lehman, and Lukk and Nersesov. These quantities have been calculated for both P and S waves for model CIT11GB, which is derived from surface wave dispersion data. First arrival times for all the models except that of Lukk and Nersesov are in close agreement, but the travel time curves for later arrivals are both qualitatively and quantitatively very different. For model CIT11GB, there are two large, overlapping regions of triplication of the travel time curve, produced by regions of rapid velocity increase near depths of 400 and 600 km. Throughout the distance range from 10 to 40 degrees, the later arrivals produced by these discontinuities have larger amplitudes than the first arrivals. The amplitudes of body waves, in fact, are extremely sensitive to small variations in the velocity structure, and provide a powerful tool for studying structural details.

Most of eastern North America, including the Canadian Shield has a Pn velocity of about 8.1 km/sec, with a nearly abrupt increase in compressional velocity by ~ 0.3 km/sec near at a depth varying regionally between 60 and 90 km. Variations in the structure of this part of the mantle are significant even within the Canadian Shield. The low-velocity zone is a minor feature in eastern North America and is subject to pronounced regional variations. It is 30 to 50 km thick, and occurs somewhere in the depth range from 80 to 160 km. The velocity decrease is less than 0.2 km/sec.

Consideration of the absolute amplitudes indicates that the attenuation due to anelasticity is negligible for 2 hz waves in the upper 200 km along the southeastern and southwestern margins of the Canadian Shield. For compressional waves the average Q for this region is > 3000. The amplitudes also indicate that the velocity gradient is at least 2 x 10^-3 both above and below the low-velocity zone, implying that the temperature gradient is < 4.8°C/km if the regions are chemically homogeneous.

In western North America, the low-velocity zone is a pronounced feature, extending to the base of the crust and having minimum velocities of 7.7 to 7.8 km/sec. Beneath the Colorado Plateau and Southern Rocky Mountains provinces, there is a rapid velocity increase of about 0.3 km/sec, similar to that observed in eastern North America, but near a depth of 100 km.

Complicated travel time curves observed on profiles with stations in both eastern and western North America can be explained in detail by a model taking into account the lateral variations in the structure of the low-velocity zone. These variations involve primarily the velocity within the zone and the depth to the top of the zone; the depth to the bottom is, for both regions, between 140 and 160 km.

The depth to the transition zone near 400 km also varies regionally, by about 30-40 km. These differences imply variations of 250 °C in the temperature or 6 % in the iron content of the mantle, if the phase transformation of olivine to the spinel structure is assumed responsible. The structural variations at this depth are not correlated with those at shallower depths, and follow no obvious simple pattern.

The computer programs used in this study are described in the Appendices. The program TTINV (Appendix IV) fits spherically symmetric earth models to observed travel time data. The method, described in Appendix III, resembles conventional least-square fitting, using partial derivatives of the travel time with respect to the model parameters to perturb an initial model. The usual ill-conditioned nature of least-squares techniques is avoided by a technique which minimizes both the travel time residuals and the model perturbations.

Spherically symmetric earth models, however, have been found inadequate to explain most of the observed travel times in this study. TVT4, a computer program that performs ray theory calculations for a laterally inhomogeneous earth model, is described in Appendix II. Appendix I gives a derivation of seismic ray theory for an arbitrarily inhomogeneous earth model.

Sputtering of insulators in the electronic stopping region

Using track detectors we have measured sputtering yields induced by MeV light ions incident on a uranium containing glass, UO₂ and UF₄. No deviation from the behavior predicted by the Sigmund theory was detected in the glass or the UO₂. The same was true for UF₄ bombarded with ⁴He at 1 MeV and with ¹⁶O and ²⁰Ne at 100 keV. In contrast to this, 4.75 MeV ¹⁹F(+2) sputters uranium from UF₄ with a yield of 5.6 ± 1.0, which is about 3 orders of magnitude larger than expected from the Sigmund theory. The energy dependence of the yield indicates that it is generated by electronic rather than nuclear stopping processes. The yield depends on the charge state of the incident fluorine but not on the target temperature. We have also measured the energy spectrum of the uranium sputtered from the UF₄. Ion explosions, thermal spikes, chemical rearrangement and induced desorption are considered as possible explanations for the anomalous yields.

High frequency wave propagation in the earth: theory and observation

The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:

1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.

2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.

3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.

These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.

New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = Δ/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.

Resolução numérica de escoamentos compressíveis empregando um método de partículas livre de malhas e o processamento em paralelo (CUDA)

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem).