Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explosive flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Gas Transport in the Insensitive High Explosive PBX 9502

Small-scale laboratory experiments were performed to analyze gas transport mechanics in PBX 9502, an insensitive high explosive (IHE) composition. Two independent investigations are reported here. First, gas permeametry and dynamic pycnometry techniques were used to measure the molecular flow (Knudsen) coefficient and the internally interconnected void fraction of pristine and thermally damaged samples in two uniaxial pressing orientations. The permeability of PBX 9502 was found to be unmeasurably low (of the same order of magnitude as PTFE) with gas transport being diffusion-dominated. Secondly, a pressure vessel experiment was developed to measure quasi-static and dynamic gas generation as the explosive was heated to self-ignition (cookoff). The gas generation results and the permeability/diffusivity findings provide evidence that PBX 9502 remains impermeable until seconds prior to self-ignition. At ignition, internal void-pressure drives macro-scale cracking and the sample becomes uniformly incorporated in the following deconsolidating deflagration. These results are discussed within the context of previous observations of pressure-dependent cookoff behavior of PBX 9502 and provide a more complete description of thermal damage evolution in this explosive composition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Material Model Parameters Optimization in Liquid Mercury Target Dynamics Simulation With Machine Learning Surrogates

A pulsed spallation target is subjected to very short (∼0.7μs) but intense loads (23.3 kJ) from repeated proton pulses, which knock away neutrons from the mercury atoms’ nuclei for a wide range application in physics, engineering, medicine, petroleum exploration, biology, chemistry, etc. The effect of this pulsed loading on the stainless-steel target module which contains the flowing mercury target material is difficult to predict not only due to its short but intense explosive-like physical reaction, but also the nonlinear material behavior of the liquid mercury in the structure. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at high power level stage. However, prediction of the resultant loading on the target is more difficult when helium gas is intentionally injected into the mercury. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. A parameter sensitivity study was firstly employed to understand their impact on the simulation strains. The investigated parameters included E, μ, γ, σ, n, VFgas, and gas cumulative volume curve control parameters a and b. Verification and validation results from sparse polynomial expansions (SPE) method and directional Gaussian smoothing (DGS) optimization show that the surrogate model had training error of ∼7% and validation error of ∼15%, indicating that machine learning methods and surrogate models can help optimize the uncertain parameters in the complex 2-phase material model. This approach is expected to fill the knowledge gap between unknown liquid-gas mixture material model and measured vessel strain responses.

Lin, Lianshan↗

Preliminary results of geologic and remote sensing studies of Rima Mozart

In order to better understand the processes responsible for the formation of lunar sinuous rilles, a study of Rima Mozart was conducted using a variety of geologic, photographic, and remote sensing data. The apparent source of this rille is located in a highlands unit of known composition and it is hypothesized that thermal and mechanical erosion played an important role in the formation of Rima Mozart. Excellent photographic, topographic, multispectral, and radar data exist for this rille. The preliminary results of an analysis of this data are presented. Photographic data indicates the presence of two volcanic source vents for Rima Mozart: Kathleen and Ann. It is suggested that Rima Mozart, like many other lunar sinuous rilles, was most likely formed by a combination of events. Rima Mozart does follow a pre-existing, dominant NW/SE structural trend suggesting the influence of structural features on the rille, however, the tectonic influence is not the sole source for the formation of the rille, as suggested by the presence of the two source vents and the spatter around Ann. It is suggested that the rille formation began with an explosive eruption at Kathleen which later calmed down to a pulsating, high volume, low-viscosity lava flow. The rapid effusion rate of the magma as well as its high temperature and turbid nature helped carve the sinuous rille into the fractured and structurally weak Apennine Bench Formation underneath. Similar eruptions and subsequent flows were also created at Ann and joined to the main channel by a NE-trending secondary rille.

Coombs, C. R.↗

Effects of highly intermittent magnetic fields and plasma flows on solar wind dynamics

Filamentary magnetic fields and intermittent mass flows with highly variant physical parameters as observed in coronal holes provide, from a theoretical point of view, natural conditions for strongly nonlinear dynamics. The presence of sheared mass flows along fine scale magnetic structures results in strong nonlinear instability, most important of which is the explosive instability. We specify the physical conditions for several different manifestations of the onset of explosive instability and its further evolution: (1) fully developed explosive instability - explosive release of the energy; (2) shock formation - stabilization of instability by small scale spatial inhomogeneities leads to formation of subsequent shocks having a number of peculiarities that is determined by the interplay of thermal and viscous losses (for example, in predominance of thermal losses the isothermal jump occurs); and (3) solitary waves - stabilization of explosive instability by nonlinear dispersion effect leads to formation of a 'gas' of solitons which are later either damped away with characteristic time and energy input or evolve to solitons with explosively growing amplitudes. Each scenario is completely determined by the physical parameters of the medium, thus producing a quite uneven distribution of energy in a coronal hole and, respectively, an uneven outward propagation of the energy flux.

Ryutova, M.↗

Spontaneous Hot Flow Anomalies at Mars and Venus

We report the first observations of Spontaneous Hot Flow Anomalies (SHFAs) at Venus and Mars, demonstrating their existence in the foreshocks of other planets beyond Earth. Using data from the ESA Venus Express and the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft, we present magnetic and plasma observations from events at both planets, exhibiting properties similar to “classical” Hot Flow Anomalies, with bounding shock-like compressive regions and a hot and diffuse core. However, these explosive foreshock transients were observed without any attendant interplanetary magnetic field discontinuity, consistent with SHFAs observed at Earth and our hybrid simulations.

Collinson, Glyn A.↗

Symposium /International/ on Combustion, 17th, Leeds University, Leeds, England, August 20-25, 1978, Proceedings

The Symposium focused on deflagration to detonation transition, coal combustion, turbulent-combustion interactions, kinetics, furnace combustion, inhibition and ignition, flame structure and chemistry, combustion studies, measurement techniques, fire and explosion, engine combustion, soot, and propellants and explosives. Papers were presented on numerical modeling of the deflagration-to-detonation transition, the interaction between turbulence and combustion, turbulent flame propagation in premixed gases, spray evaporation in recirculating flow, dissociation of nitric oxide in shock waves, pollutant emissions from partially mixed turbulent flames, energy transfer and quenching rates of laser-pumped electronically excited alkalis in flames, a study of flammability limits using counterflow flames, the unified theory of explosions with fuel consumption, and the dynamics and radiant intensity of large hydrogen flames.

Source record↗

Effect of Reynolds Number in Turbulent-Flow Range on Flame Speeds of Bunsen Burner Flames

The effect of flow conditions on the geometry of the turbulent Bunsen flame was investigated. Turbulent flame speed is defined in terms of flame geometry and data are presented showing the effect of Reynolds number of flow in the range of 3000 to 35,000 on flame speed for burner diameters from 1/4 to 1 1/8 inches and three fuels -- acetylene, ethylene, and propane. The normal flame speed of an explosive mixture was shown to be an important factor in determining its turbulent flame speed, and it was deduced from the data that turbulent flame speed is a function of both the Reynolds number of the turbulent flow in the burner tube and of the tube diameter.

Bollinger, Lowell M↗

Quantifying uncertainty for deep learning based forecasting and flow-reconstruction using neural architecture search ensembles

Classical problems in computational physics such as data-driven forecasting and signal reconstruction from sparse sensors have recently seen an explosion in deep neural network (DNN) based algorithmic approaches. However, most DNN models do not provide uncertainty estimates, which are crucial for establishing the trustworthiness of these techniques in downstream decision making tasks and scenarios. In recent years, ensemble-based methods have achieved significant success for the uncertainty quantification in DNNs on a number of benchmark problems. However, their performance on real-world applications remains under-explored. In this work, we present an automated approach to DNN discovery and demonstrate how this may also be utilized for ensemble-based uncertainty quantification. Specifically, we propose the use of a scalable neural and hyperparameter architecture search for discovering an ensemble of DNN models for complex dynamical systems. We highlight how the proposed method not only discovers high-performing neural network ensembles for our tasks, but also quantifies uncertainty seamlessly. This is achieved by using genetic algorithms and Bayesian optimization for sampling the search space of neural network architectures and hyperparameters. Subsequently, a model selection approach is used to identify candidate models for an ensemble set construction. Afterwards, a variance decomposition approach is used to estimate the uncertainty of the predictions from the ensemble. We demonstrate the feasibility of this framework for two tasks — forecasting from historical data and flow reconstruction from sparse sensors for the sea-surface temperature. In conclusion, we demonstrate superior performance from the ensemble in contrast with individual high-performing models and other benchmarks.

Deep ensembles↗

VIPIR: A High-Throughput Drop-Weight Impact Instrument for Imaging the Initiation and Propagation of Reactions in Energetic Materials

Characterizing the handling safety and sensitivity of explosives has been a challenging area of study for over 60 years. Historically one of the most accessible and widely utilized experiments has been the drop-weight impact test, which involves dropping a weight on a small sample sandwiched between two anvils. Because this experiment generally only utilizes sound thresholds to determine whether or not a sample reacted, the physical and chemical properties governing sensitivity remain convolved. Better understanding of chemical and material characteristics is needed to give the chemistry and engineering communities a predictive tool to determine the handling sensitivity of explosives prior to pursuing expensive and potentially hazardous synthesis and formulation operations. Here, we are developing a high throughput drop tower instrument capable of imaging the deformation and flow of energetic materials during impact and the resulting thermal ignition and propagation events. This instrument is based on previous designs but has been modified for higher throughput and tailorable modifications in the future. Herein, we present key design features that enable high-speed visible and thermal imaging of explosive initiation by sub-shock impacts, as well as preliminary results in which ignition sites were observed in an erythritol tetranitrate sample.

47 OTHER INSTRUMENTATION↗

Volcano evolution on Mars

The diversity of volcanic activity on Mars throughout geologic time was one of the major factors that has controlled the spatial distribution of surface mineralogies. The traditional view of Martian volcanism is one in which effusive activity has dominated the entire preserved geologic history of the planet, with the minor exception of phreatomagnetic activity and associated volcano ground-ice interactions. However, two lines of evidence have caused reconsidering of this view, and have led to the possible role of explosive volcanism on Mars. First, detailed analysis of high resolution Viking Orbiter images has provided good evidence for explosive activity on Hecates Tholus and Alba Patera. Secondly, the problems believed to exist in associating explosive volcanism with silicic magmas on Mars, and the consequent unusual magmatic evolutionary trend for Martian volcanoes from silica-rich to silica-poor, may now be circumvented by the consideration of basatic plinian activity similar in kind to terrestrial eruptions such as the 1886 Tarawera eruption. The morphologic evidence for an early phase of explosive activity on Mars is briefly reviewed, and a model is presented for the emplacement of ash-flow deposits on Martian volcanoes. The volcanoes Alba Patera and Olympus Mons are considered in this context, along with some of the older Martian tholi and paterae

Mouginis-Mark, Pete↗

Effusive silicic volcanism in the Central Andes: The Chao dacite and other young lavas of the Altiplano-Puna Volcanic Complex

The largest known Quaternary silicic lava body in the world is Cerro Chao in north Chile, a 14-km-long coulee with a volume of at least 26 cu km. It is the largest of a group of several closely similar dacitic lavas erupted during a recent (less than 100,000 year old) magmatic episode in the Altiplano-Puna Volcanic Complex (APVC; 21-24 deg S) of the Centra; Andean Volcanic Zone. The eruption of Chao proceeded in three phases. Phase 1 was explosive and produced approximately 1 cu km of coarse, nonwelded dacitic pumice deposits and later block and ash flows that form an apron in front of the main lava body. Phase 2 was dominantly effusive and erupted approximately 22.5 cu km of magma in the form of a composite coulee covering approximately 53 sq km with a 400-m-high flow front and a small cone of poorly expanded pumice around the vent. The lava is homogeneous with rare flow banding and vesicular tops and selvages. Ogives (flow ridges) reaching heights of 30 m form prominent features on its surface. Phase 3 produced a 6-km-long, 3-km-wide flow that emanated from a collapsed dome. Ogives are subdued, and the lava is glassier than that produced in previous phases. All the Chao products are crystal-rich high-K dacites and rhyodacites with phenocrysts of plagioclase, quartz, hornblende, biotite, sphene, rare snidine, and oxides. Phenocryst contents reach 40-60 vol % (vesicle free) in the main phase 2 lavas but are lower in the phase 1 (20-25%) and phase 3 (approximately 40%) lavas. Ovoid andesitic inclusions with vesicular interiors and chilled margins up to 10 cm are found in the later stages of phase 2 and compose up to 5% of the phase 3 lava. There is little evidence for preeruptive zonation of the magma body in composition, temperature (approximately 840 C), fO2 (19(exp -11), or water content, so we propose that eruption of the Chao complex was driven by intrusion of fresh, hot andesitic magma into a crystallizing and largely homogeneous body of dacitic magma. Morphological measurements suggest that the Chao lavas had internal plastic viscosities of 10(exp 10) to 10(exp 12) Pa s, apparent viscosities of 10(exp 9) Pa s, surface viscosities of 10(exp 15) to 10(exp 24) Pa s, and a yield strength of 8 x 10(exp 5) Pa. These estimates indicate that Chao would have exhibited largely similar rheological properties to other silicic lava extrusions, notwithstanding its high phenocryst content. We suggest that Chao's anomalous size is a function of both the relatively steep local slope (20 deg to 3 deg) and the available volume of magma. The eruption duration for Chao's emplacement is thought to have been about 100 to 150 years, with maximum effusion rates of about 25 cu m/s for short periods. Four other lavas in the vicinity with volumes of approximately 5 cu km closely resemble Chao and are probably comagnetic. The suite as a whole shares a petrologic and chemical similarity with the voluminous regional Tertiary to Pleistocene ignimbrites of the APVC and may be derived from a zone of silicic magmatism that is thought to have been active since the late Tertiary. Chao and the other young lavas may represent either the waning of this system or a new episode fueled by intrusions of mafic magma.

De Silva, S. L.↗

Space-Time Conservation Element and Solution Element Method Being Developed

The engineering research and design requirements of today pose great computer-simulation challenges to engineers and scientists who are called on to analyze phenomena in continuum mechanics. The future will bring even more daunting challenges, when increasingly complex phenomena must be analyzed with increased accuracy. Traditionally used numerical simulation methods have evolved to their present state by repeated incremental extensions to broaden their scope. They are reaching the limits of their applicability and will need to be radically revised, at the very least, to meet future simulation challenges. At the NASA Lewis Research Center, researchers have been developing a new numerical framework for solving conservation laws in continuum mechanics, namely, the Space-Time Conservation Element and Solution Element Method, or the CE/SE method. This method has been built from fundamentals and is not a modification of any previously existing method. It has been designed with generality, simplicity, robustness, and accuracy as cornerstones. The CE/SE method has thus far been applied in the fields of computational fluid dynamics, computational aeroacoustics, and computational electromagnetics. Computer programs based on the CE/SE method have been developed for calculating flows in one, two, and three spatial dimensions. Results have been obtained for numerous problems and phenomena, including various shock-tube problems, ZND detonation waves, an implosion and explosion problem, shocks over a forward-facing step, a blast wave discharging from a nozzle, various acoustic waves, and shock/acoustic-wave interactions. The method can clearly resolve shock/acoustic-wave interactions, wherein the difference of the magnitude between the acoustic wave and shock could be up to six orders. In two-dimensional flows, the reflected shock is as crisp as the leading shock. CE/SE schemes are currently being used for advanced applications to jet and fan noise prediction and to chemically reacting flows.

Chang, Sin-Chung↗

Non-Contact Vapor Detection of Illicit Drugs via Atmospheric Flow Tube-Mass Spectrometry

Real-time, non-contact detection of illicit drugs is a desirable goal for the interdiction of these controlled substances, but the relatively low vapor pressures of such species present a challenge for trace vapor detection technologies. The introduction of atmospheric flow tube-mass spectrometry (AFT-MS), which has previously been demonstrated to detect gas-phase analytes at low parts-per-quadrillion levels for explosives and organophosphorus compounds, also enables the potential for non-contact drug detection. With AFT-MS, direct vapor detection of cocaine and methamphetamine from ~5 µg residues at room temperature is demonstrated herein. Furthermore, thermal desorption of low- to sub-picogram levels of cocaine, methamphetamine, fentanyl, and heroin is observed via AFT-MS using a carrier flow rate of several L/min of air. These low levels can permit non-contact sampling through collection of vapor, effectively preconcentrating the analyte before desorption and analysis. Quantitative evaluation of the thermal desorption approach has yielded limits of detection (LODs) on order of 10 fg for cocaine and fentanyl, 100 fg for methamphetamine, and 1.6 pg for heroin. The LOD for heroin was lowered to 200 fg by using tributyl phosphate as a dopant to form a proton-bound heterodimer with heroin. When used with AFT-MS, the intentional formation of specific drug-dopant adducts has the potential to enhance detection limits and selectivity of additional drug species. Species that are prone to form adducts present a challenge to analysis, but that difficulty can be overcome by the intentional addition of a dopant. Molecules unlikely to form adducts will remain essentially unimpacted, but the adduct-forming species will interact with the dopant to compress the analyte signal into a single peak. This approach would be valuable in the application of non-contact screening for illicit substances via vapor collection followed by thermal desorption for analysis.

methamphetamine, cocaine, heroin, fentanyl, atmosp↗

The Effect of Clearance Distribution on the Performance of a Compression-ignition Engine with a Precombustion Chamber

The clearance distribution in a precombustion chamber cylinder head was varied so that for a constant compression ratio of 13.5 the spherical auxiliary chambers contained 20, 35, 50, and 70 per cent of the total clearance volume. Each chamber was connected to the cylinder by a single circular passage, flared at both ends, and of a cross-sectional area proportional to the chamber volume, thereby giving the same calculated air-flow velocity through each passage. Results of engine-performance tests are presented with variations of power, fuel consumption, explosion pressure, rate of pressure rise, ignition lag, heat loss to the cooling water, and motoring characteristics. For good performance the minimum auxiliary chamber volume, with the cylinder head design used, was 35 per cent of the total clearance volume; for larger volumes the performance improves but slightly. With the auxiliary chamber that contained 35 percent of the clearance volume there were obtained the lowest explosion pressures, medium rates of pressure rise, and slightly less than the maximum power. For all clearance distributions an increase in engine speed decreased the ignition lag in seconds and increased the rate of pressure rise.

Moore, C. S.↗

Modeling the AC Power Flow Equations with Optimally Compact Neural Networks: Application to Unit Commitment

Nonlinear power flow constraints render a variety of power system optimization problems computationally intractable. Emerging research shows, however, that the nonlinear AC power flow equations can be successfully modeled using neural networks. These neural networks can be exactly transformed into mixed integer linear programs and embedded inside challenging optimization problems, thus replacing nonlinearities that are intractable for many applications with tractable piecewise linear approximations. Such approaches, though, suffer from an explosion of the number of binary variables needed to represent the neural network. Accordingly, this paper develops a technique for training an "optimally compact'' neural network, i.e., one that can represent the power flow equations with a sufficiently high degree of accuracy while still maintaining a tractable number of binary variables. We demonstrate the use of this neural network as an approximator of the nonlinear power flow equations by embedding it in the AC unit commitment problem, transforming the problem from a mixed integer nonlinear program into a more manageable mixed integer linear program. We use the 14-, 57-, and 89-bus networks as test cases and compare the AC-feasibility of commitment decisions resulting from the neural network, DC, and linearized power flow approximations. Our results show that the neural network model outperforms both the DC and linearized power flow approximations when embedded in the unit commitment problem. The neural network formulation most often selects a feasible unit commitment schedule, and furthermore, it only s

AC power flow↗

Effects of Model Uncertainties in Underground Chemical Explosions on Far-field Results

We used the CTH shock physics code to simulate the explosion of an 18-t chemical explosive at a depth of 250 m. We used the CTH in the two-dimensional axisymmetric (cylindrical) geometry (2DC) and most simulations included fully tamped explosions in wet tuff. Our study focused on parametric studies of three of the traditional strength models available in CTH, namely, geologic-yield, elastic perfectly-plastic von Mises, and Johnson-Cook strength (flow stress) models. We processed CTH results through a code that generates Reduced Displacement Potential (RDP) histories for each simulation. Since RDP is the solution of the linear wave equation in spherical coordinates, it is mainly valid at far-enough distance from the explosion the elastic radius. Among various parameters examined, we found the yield strength to have the greatest effect on the resulting RDP, where the peak RDP reduces almost linearly in log-log space as the yield strength increases. Moreover, an underground chemical explosion results in a cavity whose final diameter is inversely proportional to the material yield strength, i.e., as the material's yield strength increases the resulting final cavity radius decreases. Additionally, we found the choice of explosive material (COMP-C4 versus COMP-B) has minor effects on the peak RDP, where denser COMP-C4 shows higher peak RDP than the less dense COMP-B by a factor of ~1.1. In addition to wet tuff, we studied explosions in dry tuff, salt, and basalt, for a single strength model and yield strength value. We found wet tuff has the highest peak RDP value, followed by dry tuff, salt, and basalt. 2DC simulations of explosions in 11 m radius spherical, hemispherical, and cylindrical cavities showed the RDP signals have much lower magnitude than tamped explosions, where the cavity explosions mimicked nearly decoupled explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Magnetic reconnection in the era of exascale computing and multiscale experiments

Astrophysical plasmas have the remarkable ability to preserve magnetic topology, which inevitably gives rise to the accumulation of magnetic energy within stressed regions including current sheets. This stored energy is often released explosively through the process of magnetic reconnection, which produces a reconfiguration of the magnetic field, along with high-speed flows, thermal heating and nonthermal particle acceleration. Either collisional or kinetic dissipation mechanisms are required to overcome the topological constraints, both of which have been predicted by theory and validated with in situ spacecraft observations or laboratory experiments. However, major challenges remain in understanding magnetic reconnection in large systems, such as the solar corona, where the collisionality is weak and the kinetic scales are vanishingly small in comparison with macroscopic scales. The plasmoid instability or formation of multiple plasmoids in long, reconnecting current sheets is one possible multiscale solution for bridging this vast range of scales, and new laboratory experiments are poised to study these regimes. In conjunction with these efforts, we anticipate that the coming era of exascale computing, together with the next generation of observational capabilities, will enable new progress on a range of challenging problems, including the energy build-up and onset of reconnection, partially ionized regimes, the influence of magnetic turbulence and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct observation of diamond formation in a shock-compressed high explosive

Understanding the formation timescale and structure of carbonaceous reaction products is critical for modeling the high-pressure equation-of-state of organic materials. We use the National Ignition Facility to shock-compress polycrystalline TATB (C 6 H 6 N 6 O 6 ) samples to ~70–130 GPa and ~4000–5500 K, employing in situ nanosecond X-ray diffraction to probe reaction products and velocimetry to measure transmitted compression wave profiles. Our diffraction data is consistent with the formation of diamond over timescales less than ~60 ns. This represents carbon condensation from a molecular explosive on timescales three times faster than previously reported and the earliest observation of diamond produced from reacting TATB. Reactive flow simulations with explicit chemistry reproduce the observed temporal structure within wave profiles to inform the distribution of P-T states. These findings provide direct evidence of ultrafast diamond formation in a reactive system at extreme conditions and provide new constraints for models of shock and detonation chemistry.

Clarke, Samantha M. [Lawrence Livermore National L↗