Fragment spatial distribution of prismatic casing under internal explosive loading

Non-cylindrical casings filled with explosives have undergone rapid development in warhead design and explosion control.The fragment spatial distribution of prismatic casings is more complex than that of traditional cylindrical casings.In this study,numerical and experimental investigations into the fragment spatial distribution of a prismatic casing were conducted.A new numerical method,which adds the Lagrangian marker points to the Eulerian grid,was proposed to track the multi-material interfaces and material dynamic fractures.Physical quantity mappings between the Lagrangian marker points and Eulerian grid were achieved by their topological relationship.Thereafter,the fragment spatial distributions of the prismatic casing with different fragment sizes,fragment shapes,and casing geometries were obtained using the numerical method.Moreover,fragment spatial distribution experiments were conducted on the prismatic casing with different fragment sizes and shapes,and the experimental data were compared with the numerical results.The effects of the fragment and casing geometry on the fragment spatial distributions were determined by analyzing the numerical results and experimental data.Finally,a formula including the casing geometry parameters was fitted to predict the fragment spatial distribution of the prismatic casing under internal explosive loading.

High-Precision Direct Method for the Radiative Transfer Problems

It is the main aim of this paper to investigate the numerical methods of the radiative transfer equation. Using the five-point formula to approximate the differential part and the Simpson formula to substitute for integral part respectively, a new high-precision numerical scheme, which has 4-order local truncation error, is obtained. Subsequently, a numerical example for radiative transfer equation is carried out, and the calculation results show that the new numerical scheme is more accurate.

A Continuous Minimization Method for Solving Optimal Control Problem

In this paper,a kind of continuous minimization method for solving optimal control problem which transforms the problem into solving the initial value problem of ordinary differential equations is presented so that some new algorithms for solving optimal control problem are obtained by constrcting suitable differential equations and choosing some numerical integration formulas, and prove the every accumulation point developed by the methods is a stationary control.We also point out that existing optimization algorithms for computing optimal control,in general,do not gauratee to converge to stationary control.At last,some numerical examples show the method is elective.

ON CALCULATIONS OF THE DRAG COEFFICIENT Cd AND THE FALL VELOCITY ω OF SPHERICAL BODIES

A host of authors have proposed some theoretical and experimental formulas in hydromechanics concerning the calculation of the drag coefficient Cd of spherical bodies. But all of the existing Cd formulas hold true only at small Reynolds numbers and are restricted within certain flowing range.As regards the fall velocity ω of spherical bodies, there is yet no formula applicable to each flowing range and to a direct expression and calculation of the fall velocity ω.In view of these, from N-S equations, and meanwhile based on measured data and complicated calculations, the author has developed and proposed the following results:(1) The drag coefficient (2) The dimensionless fall velocity where Es, Ω* and constants etc. are indicated in detail in this paper.Through laborious calculation in lgRe<5 larger range, the verification proves that our results well agree with the measured data. And the leading features of formulas of this paper are: (1) simple in form, (2) convenient for general use, (3) preferable on the part of the precision and applicability.Finally, to introduce this process and to illustrate the temperature effects on the fall velocity ω, some examples are discussed in this paper.