Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA complete equation of state for non-ideal condensed phase explosives
67
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The objective of this work is to improve the robustness and accuracy of numerical simulations of both ideal and non-ideal explosives by introducing temperature dependence in mechanical equations of state for reactants and products. To this end, we modify existing mechanical equations of state to appropriately approximate the temperature in the reaction zone. Mechanical equations of state of Mie-Gr{u}neisen form are developed with extensions, which allow the temperature to be evaluated appropriately, and the temperature equilibrium condition to be applied robustly. Furthermore the snow plow model is used to capture the effect of porosity on the reactants equation of state. We apply the methodology to predict the velocity of compliantly confined detonation waves. Once reaction rates are calibrated for unconfined detonation velocities, simulations of confined rate sticks and slabs are performed, and the experimental detonation velocities are matched without further parameter alteration, demonstrating the predictive capability of our simulations. We apply the same methodology to both ideal (PBX9502, a high explosive with principal ingredient TATB) and non-ideal (EM120D, an ANE or ammonium nitrate based emulsion) explosives.
rate research
Read More
The effects of mid-range repulsion in Lattice Boltzmann models on the coalescence/breakup behaviour of single-component, non-ideal fluids are investigated. It is found that mid-range repulsive interactions allow the formation of spray-like, multi-droplet configurations, with droplet size directly related to the strength of the repulsive interaction. The simulations show that just a tiny ten-percent of mid-range repulsive pseudo-energy can boost the surface/volume ratio of the phase- separated fluid by nearly two orders of magnitude. Drawing upon a formal analogy with magnetic Ising systems, a pseudo-potential energy is defined, which is found to behave like a quasi-conserved quantity for most of the time-evolution. This offers a useful quantitative indicator of the stability of the various configurations, thus helping the task of their interpretation and classification. The present approach appears to be a promising tool for the computational modelling of complex flow phenomena, such as atomization, spray formation and micro-emulsions, break-up phenomena and possibly glassy-like systems as well.
Truncated Taylor expansions of smooth flow maps are used in Hamiltons principle to derive a multiscale Lagrangian particle representation of ideal fluid dynamics. Numerical simulations for scattering of solutions at one level of truncation are found to produce solutions at higher levels. These scattering events to higher levels in the Taylor expansion are interpreted as modeling a cascade to smaller scales.
We propose the rhoLENT method, an extension of the unstructured Level Set / Front Tracking (LENT) method, based on the collocated Finite Volume equation discretization, that recovers exact numerical stability for the two-phase momentum convection with a range of density ratios, namely $rho^-/rho^+in [1, 10000]$. We provide the theoretical basis for the numerical inconsistency in the collocated finite volume equation discretization of the single-field two-phase momentum convection. The cause of the numerical inconsistency lies in the way the cell-centered density is computed in the new time step ($rho_c^{n+1}$). Specifically, if $rho_c^{n+1}$ is computed from the approximation of the fluid interface at $t^{n+1}$, and it is not computed by solving a mass conservation equation (or its equivalent), the two-phase momentum convection term will automatically be inconsistently discretized. We provide the theoretical justification behind using the auxiliary mass conservation equation to stabilize flows with strong density ratios. The evaluation of the face-centered (mass flux) density we base on the fundamental principle of mass conservation, used to model the single-field density, contrary to the use of different weighted averages of cell-centered single-field densities and alternative reconstructions of the mass flux density by other contemporary methods. Implicit discretization of the two-phase momentum convection term is achieved, removing the CFL stability criterion. Numerical stability is demonstrated in terms of the relative $L_infty$ velocity error norm with realistic viscosity and strong surface tension forces. The stabilization technique in the rhoLENT method is also applicable to other two-phase flow simulation methods that utilize the collocated unstructured Finite Volume Method to discretize single-field two-phase Navier-Stokes Equations.
The smoothed particle hydrodynamics (SPH) method is a useful numerical tool for the study of a variety of astrophysical and planetlogical problems. However, it turned out that the standard SPH algorithm has problems in dealing with hydrodynamical instabilities. This problem is due to the assumption that the local density distribution is differentiable. In order to solve this problem, a new SPH formulation, which does not require the differentiability of the density, have been proposed. This new SPH method improved the treatment of hydrodynamical instabilities. This method, however, is applicable only to the equation of state (EOS) of the ideal gas. In this paper, we describe how to extend the new SPH method to non-ideal EOS. We present the results of various standard numerical tests for non-ideal EOS. Our new method works well for non-ideal EOS. We conclude that our new SPH can handle hydrodynamical instabilities for an arbitrary EOS and that it is an attractive alternative to the standard SPH.
We solve the Navier-Stokes equations with two simultaneous forcings. One forcing is applied at a given large-scale and it injects energy. The other forcing is applied at all scales belonging to the inertial range and it injects helicity. In this way we can vary the degree of turbulence helicity from non helical to maximally helical. We find that increasing the rate of helicity injection does not change the energy flux. On the other hand the level of total energy is strongly increased and the energy spectrum gets steeper. The energy spectrum spans from a Kolmogorov scaling law $k^{-5/3}$ for a non-helical turbulence, to a non-Kolmogorov scaling law $k^{-7/3}$ for a maximally helical turbulence. In the later case we find that the characteristic time of the turbulence is not the turnover time but a time based on the helicity injection rate. We also analyse the results in terms of helical modes decomposition. For a maximally helical turbulence one type of helical mode is found to be much more energetic than the other one, by several orders of magnitude. The energy cascade of the most energetic type of helical mode results from the sum of two fluxes. One flux is negative and can be understood in terms of a decimated model. This negative flux is however not sufficient to lead an inverse energy cascade. Indeed the other flux involving the least energetic type of helical mode is positive and the largest. The least energetic type of helical mode is then essential and cannot be neglected.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
