Analysis Of Finite Difference Method Applied In Energy Transfer Of Energetic Materials And Product Formula System Design

Energy transfer is always an important subject for its specificity. Numerical simulation methods draw more and more attentions with the development of computation technology. Simulation methods can be classified as Finite element method (FEM) and finite difference method (FDM). Foreign workers have exploited powerful software based on FEM, such as ANSYS, ABAQUS, and NASTRAN. But it is too large for exploitation of professional formulation system software. Furthermore it is not safe according to investigation of National Security Agency.FDM can overcome these problems. Grid division and software implementation is easier. It is convenient to embed into professional software. For these reasons, analysis of FDM applied in simulation of energy transfer is very necessary.Somebody has already done a lot of work about energy transfer based on Lax-Friedrichs difference scheme. In this thesis, we compare four different finite difference schemes. They are Lax-Friedrichs scheme, upwind scheme, Lax-Wendroff scheme and Beam-Warming scheme. After analysis of truncation error and stability combined with a large number of numerical tests, we finally choose a method based on Beam-Warming scheme and Up-Wind scheme (UW-BW method) as the tool of our professional formulation system. A switch is designed according to proportional distance.In numerical experiments, we firstly simulate three Elemental energetic materials and add these materials into our material library. Secondly, we simulate hybrid high energy materials mixed with two elemental energetic materials. Relational expressions about JWL coefficients and material proportion are summarized. Thirdly, we simulate energetic materials with Al and AP. Relational expressions about JWL coefficients and material proportion are also summarized. These expressions can be used in development of material library. At last, we develop single professional formulation system software, which includes computation of detonation parameters, simulation of different materials based on UW-BW difference method, calculation of overpressure peak by experiential formula.