CTVD格式数值计算非均质炸药爆轰问题

将高分辨率激波捕捉格式CTVD格式拓展应用到非均质炸药爆轰的数值模拟问题.增加了化学反应率控制方程,引入Lee-Tarver点火成长模型,未反应的固体炸药和化学反应气体产物都使用JWL形式状态方程.数值模拟了非均质固体炸药PBX-9404和TATB的冲击起爆问题.获得了较高的爆轰波分辨率和光滑解区的数值精度,对具有复杂物态方程形式的固体炸药爆轰问题,CTVD格式具有简单实用、高效和高分辨的特点.     

On the State of Ionization,Non-Ideality and Electric Conductivity of the Cathode Plasma of Vacuum Arcs

The cathode spot plasma of a metal vapour are in vacuum has a high density (10²⁵–10²⁸ m³) with a relatively low temperature (1–5 eV), therefore it is strongly non-ideal. The degree of ionization α and the electric conductivity σ of such plasmas are discussed, considering that non-ideality. The burning voltage of the arc U_B is estimated with the help of σ-values, obtained from theories and measurements for non-ideal copper plasmas. It is shown, that a high current density of the cathode spot of j = 10¹² A/m² is not in contradiction to a low burning voltage U_B = 20–30 V, but these values give a consistent picture of the spot processes.     

Macrokinetics of the retention of condensed carbon and detonation diamond in a hermetic explosion chamber

The effect of the thermophysical parameters of a cooling medium on the macrokinetics of secondary physicochemical processes that occur in a hermetic chamber after the explosion of a solid explosive charge is studied. The yields of condensed carbon and the content of the diamond phase in it are mainly determined by the temperature of the medium in the chamber after explosion. The maximum yield of detonation diamond synthesized from a trinitrotoluene-hexogen TG50/50 alloy is equal to ～10% of the initial explosive mass and is achieved when the steady-state temperature of the medium in the chamber does not exceed T_m=550±50 K. As this temperature increases, the yield of detonation diamond decreases approximately in inverse proportion to the temperature, and, at T_m>2800 K, there is virtually no diamond phase in the explosion products. The conversion of condensed carbon due to the presence of the oxygen-containing components of the explosion products (CO₂, H₂O) begins at a temperature of the medium of above 1550±150 K. The decrease in the final energy release that is experimentally detected in calorimetric studies in the case of an explosion of solid explosives with a negative oxygen balance in an inert gaseous medium or in the case where explosives are surrounded by massive shells results from the endothermic conversion of condensed carbon, which absorbs a significant portion of the explosion energy.     

Statistics of inter-particle distances and angles in plasmas

Accurate treatment of the plasma density effects requires a detailed knowledge of the spatial distribution of individual ions around a test ion. In the present work, rigorous expressions are derived for the main 2- and 3-particle spatial distribution functions involving the nearest neighbor (NN) and the next-nearest neighbor (NNN) ions. These expressions, valid for both ideal and nonideal plasmas, present the distributions as functionals of the potentials U _NN and U _NNN at the nearest and next-nearest ion locations. All of the distribution functions except one are derived and discussed in the present work for the first time ever. For utilization of our results in practical calculations, we suggest semi-empirical expressions for U _NN and U _NNN in the ion-ion coupling parameter range 0 ?Γ < 1. In order to test the accuracy of our expressions for U _NN and U _NNN we conduct Molecular Dynamics (MD) simulations. The simulations utilize the pure Coulomb particle-particle interaction potentials, regularized at close range to avoid classical Coulomb collapse, and are free from the assumptions made to find U _NN and U _NNN. Thus, the results of the MD simulations provide an independent test of our theoretical results. Excellent agreement has been found between the results of the theory and of the MD simulations. Finally, we outline the implications of the present findings on the problem of tunneling and charge exchange in dense plasmas. Received 27 October 2000 and Received in final form 30 January 2001     

Moduli Space of BPS Walls in Supersymmetric Gauge Theories

Existence and uniqueness of the solution are proved for the ‘master equation’ derived from the BPS equation for the vector multiplet scalar in the U(1) gauge theory with N _F charged matter hypermultiplets with eight supercharges. This proof establishes that the solutions of the BPS equations are completely characterized by the moduli matrices divided by the V-equivalence relation for the gauge theory at finite gauge couplings. Therefore the moduli space at finite gauge couplings is topologically the same manifold as that at infinite gauge coupling, where the gauged linear sigma model reduces to a nonlinear sigma model. The proof is extended to the U(N _C) gauge theory with N _F hypermultiplets in the fundamental representation, provided the moduli matrix of the domain wall solution is U(1)-factorizable. Thus the dimension of the moduli space of U(N _C) gauge theory is bounded from below by the dimension of the U(1)-factorizable part of the moduli space. We also obtain sharp estimates of the asymptotic exponential decay which depend on both the gauge coupling and the hypermultiplet mass differences.     

Monte-Carlo algorithm for a two-dimensional Hubbard model of high-T c superconductivity

A Monte-Carlo procedure is given for the two-dimensional (2-D) Hubbard model using the Suzuki-Trotter transformation. The resulting three-dimensional (3-D) classical model does not have the usual problems with negative transition probabilities in the large-U limit (U-repulsive interactions). Numerical simulations based on the algorithm described are expected to be of importance for the theory of high-T _c superconductivity.     

Some extensions on the soliton solutions for the Novikov equation with cubic nonlinearity

In this paper, by using the bifurcation method of dynamical systems, we derive the traveling wave solutions of the nonlinear equation UU_τyy ? U_yU_τy + U²U_τ + 3U_y = 0. Based on the relationship of the solutions between the Novikov equation and the nonlinear equation, we present the parametric representations of the smooth and nonsmooth soliton solutions for the Novikov equation with cubic nonlinearity. These solutions contain peaked soliton, smooth soliton, W-shaped soliton and periodic solutions. Our work extends some previous results.     

Methods to determine surface energies of natural fibres: a review

To tailor the interaction across composite interfaces especially for the development of green composites, i.e. composites made completely from renewable materials, information about the fibre surfaces is required. We review the current state of the art of methods to determine the surface tension of natural fibres and discuss the advantages and disadvantages of techniques used. Although numerous techniques have been employed to characterise surface tension of natural fibres, it seems that commonly used wetting techniques are very much more affected by the non-ideal character of natural fibres. Inverse Gas Chromatography (IGC) is a much better suited technique to determine the surface energetic properties of natural fibres than wetting techniques. The surface tension of natural reinforcements, determined using IGC, was reported for nanosized bacterial cellulose as well as bamboo, cornhusk, flax, hemp and sisal, covering a wide range of cellulose content. The effect of methods to separate/extract fibres from the plants as well as of a few surface modification procedures on the fibre surface properties is also reviewed. The dispersive part of the natural fibre surface tension γ ^d _S varies from 32 to 61 mJ/m². The fibre surface tension increases with increasing cellulose content of natural fibres. We also found that a higher basicity (Donor Number, K _B to Acceptor Number, K _Aratio) was observed for fibres containing more cellulose. This may be reflective of higher crystalline cellulose content in the surfaces of the fibres, as only the ether linkage of the cellulose is labile for hydrogen bonding.     

Ternary Z2 × Z3 Graded Algebras and Ternary Dirac Equation

The wave equation generalizing the Dirac operator to the Z₃-graded case is introduced, whose diagonalization leads to a sixth-order equation. It intertwines not only quark and anti-quark state as well as the u and d quarks, but also the three colors, and is therefore invariant under the product group Z₂ × Z₂ × Z₃. The solutions of this equation cannot propagate because their exponents always contain non-oscillating real damping factor. We show how certain cubic products can propagate nevertheless. The model suggests the origin of the color SU(3) symmetry and of the SU(2) × U(1) that arise automatically in this model, leading to the full bosonic gauge sector of the Standard Model.

     

Hysteresis of nanocrystalline R–Fe–B

The energetic model of ferromagnetic hysteresis calculates the magnetic state of materials by minimizing the total energy function for statistical domain behavior. The approach shows good agreement with the magnetization curves of mechanically alloyed Pr₉Fe₈₅B₆ powder, heat treated at different temperatures.     

惰性物质等离子体物态方程研究

对高温高压下惰性等离子体的电离度和物态方程，给出了一种基于Thomas-Feimi(TF)统计模型的简化计算新方法，即首先对TF模型电离势的数值结果进行函数逼近，得出近似计算电离势的简单解析函数；在局部热动平衡情况下，假定离子数密度n(Z^*)为Z^*的连续函数，再由Debye-Hückel修正的 Saha 方程，得出了一个便于数值求解的电离度的近似计算公式，从而建立了一种惰性等离子体物态方程的简化模型，并对氦、氖、氩三种惰性物质等离子体进行了计算.计算结果与其他文献计算结果和实验值均符合很好.所提出的简单模型也适用于计算混合物物态方程，在高温高密度强电离等离子体领域将有更为广阔的应用前景.     

Cosmological constraints on unparticle dark matter

In unparticle dark matter (unmatter) models, the equation of state of the unmatter is given by p=ρ/(2d _U+1), where d _U is the scaling factor. Unmatter with such equations of state would have a significant impact on the history of the expansion of the universe. Using type Ia supernovae (SNIa), the baryon acoustic oscillation (BAO) measurements and the shift parameter of the cosmic microwave background (CMB) to place constraints on such unmatter models, we find that if only the SNIa data are used, the constraints are weak. However, with the BAO and CMB shift parameter data added, strong constraints can be obtained. For the ΛUDM model, in which unmatter is the sole dark matter, we find that d _U>60 at 95% C.L. For comparison, in most unparticle physics models it is assumed that d _U<2. For the ΛCUDM model, in which unmatter co-exists with cold dark matter, we found that the unmatter can at most make up a few percent of the total cosmic density if d _U<10; thus it cannot be the major component of dark matter.     

Equation of State for Neutron Star Matter with the Lowest Order Constrained Variational Method and Kaon Condensation

We investigate the equation of state of kaon condensed neutron star matter using the lowest order constrained variational (LOCV) method. In this paper, we employ the AV₁₈ and AV₁₄ potentials for nucleon-nucleon interaction and use the chiral model for kaon-nucleon interaction. According to our results, the kaon condensation occurs at some density values and softens the equation of state as other techniques have shown.     

Ferromagnetism in the Hubbard model on the square lattice: Improved instability criterion for the Nagaoka state

The instability of the fully polarized ferromagnetic state (Nagaoka state) with respect to single spin flips is re-examined for the Hubbard model on the square lattice with a large family of variational wave functions which include correlation effects of the majority spins in the vicinity of the flipped spin. We find a critical hole density of δ_cr = 0.251 for U = ∞ and a critical coupling of U_cr = 77.7t. Both values improve previous variational results considerably.     

Single-particle potential in a chiral approach to nuclear matter including short-range NN-terms

We extend a recent chiral approach to nuclear matter of Lutz et al.Phys. Lett. B 474, 7 (2000)) by calculating the underlying (complex-valued) single-particle potential U(p, k _f) + iW(p, k _f). The potential for a nucleon at the bottom of the Fermi sea, U(0, k _f0) = - 20.0 MeV, comes out as much too weakly attractive in this approach. Even more seriously, the total single-particle energy does not rise monotonically with the nucleon momentum p, implying a negative effective nucleon mass at the Fermi surface. Also, the imaginary single-particle potential, W(0, k _f0) = 51.1 MeV, is too large. More realistic single-particle properties together with a good nuclear-matter equation of state can be obtained if the short-range contributions of non-pionic origin are treated in mean-field approximation (i.e. if they are not further iterated with 1π-exchange). We also consider the equation of state of pure neutron matter ˉE_n(k _n) and the asymmetry energy A(k _f) in that approach. The downward bending of these quantities above nuclear-matter saturation density seems to be a generic feature of perturbative chiral pion-nucleon dynamics. Received: 19 December 2002 / Accepted: 11 February 2003 / Published online: 15 April 2003     

Everywhere invariant spaces of metrics and isometries

We define an everywhere invariant space of metricsU to be one for which the set of diffeomorphisms, which leaveU invariant, contains all the isometries of the individual metrics inU. We also generalize Killing's equation to a new equation, the invariance equation, which has as solutions those vector fields which leaveU invariant. By combining these ideas we give a new method for finding the isometries of a given metric.     

Bound State Effects in Quantum Transport Theory

A quantum kinetical equation for the one-particle density operator for inhomogeneous multi component non-ideal gases is derived taking into account retardation effects and external fields. It is shown that the interactions give rise to additional contributions to the quantities in the transport equation. In some terms (e.g., internal energy, pressure, effective external force) the interaction effects are strong in the presence of bound states and yield contributions of the order exp (—E₁₀/kT) where E₁₀ is the ground state energy of bound pairs of particles. Other terms (e.g. friction forces) are rather insensitive with respect to the presence of bound states. An interpretation of these effects in terms of the mass action law is given using the principle of equivalence between bound states and particles.     

First principles study of the magnetism driven by cation defects in CeO2: the important role of O2p states

The magnetism driven by cation defects in undoped CeO2 bulk and thin films is studied by the density functional theory corrected for on-site Coulomb interactions （DFT＋U） with U = 5 eV for the Ce4f states and U = 7 eV for the O2p states. It is found that the Ce vacancies can induce a magnetic moment of the -4 gB/supercell, which arises mainly from the 2p hole state of the nearest neighbouring O atom （-1μB on per oxygen） to the Ce vacancy. The effect of the methodology is investigated, indicating that U = 7 eV for the O2p state is necessary to obtain the localized O2p hole state in defective ceria with cation vacancies.     

Analytic Advantages of Spherically Symmetric Step-Function Potentials in the Dirac Equation with Scalar and Fourth Component of Vector Potential

The root mean square radii of the particle orbits are calculated (semi)analyticallyfor every bound state, using the Dirac equation with a scalar potentialU _s andfourth component of a vector potentialU _v in the case of a spherically symmetricstep-function shape with the same radiusR for these potentials. In addition, a(semi)analytic expression of the expectation value of the corresponding potentialenergy operator is derived. For the above quantities, expressions of the energyeigenvalues in terms of the potential parameters are needed and approximateformulas may be used in certain cases. This study emphasizes the analyticadvantages of the relativistic, spherically symmetric step-function potential model.Its applicability is discussed in connection with a problem of physical interest,namely that of the motion of a particle in hypernuclei.     

KCrF3中的轨道有序及其成因

强关联体系中的轨道有序及其成因一直是凝聚态物理研究的热点问题.轨道有序对于巨磁阻和 超导材料的研究有非常重要的地位.利用第一性原理计算研究了KCrF₃的四方相和立方相中的轨道有序 及其成因.在四方相中, GGA和GGA+U两种方法计算结果都表明其基态是A型反铁磁和G型轨道有序. 对于立方结构, GGA方法得出铁磁半金属态是基态,而GGA+U(U_eff = 3.0 eV)得到的基态是A型 反铁磁绝缘体. 光电导测量是少数能从实验上观察到轨道有序的方法之一,因此计算了其光电导,并结合投影态密度讨论 了KCrF₃中的轨道有序.最后找到了其轨道有序的成因:电子强关联效应,而非电-声子相互作用是其 轨道 有序的物理根源.