Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow analysis”

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 1,081 records · Page 60

Spectral Analysis of Regular Material Point Method and its Application to Study High Pressure Reverse Osmosis Membrane Compaction and Embossing

Material Point Method (MPM) is gaining widespread interest in applied continuum mechanics. The fact that all the continuum properties are stored on the particles (or material points) and the governing equations are solved on these material points makes MPM extremely suited to problems involving severe material deformations, such as crack propagation, soil movement, and fluid flows. Despite its popularity, only a few studies have focused on the numerical properties of MPM. This presentation introduces a global spectral analysis of the regular material point method. Contrary to previous studies, the analysis focuses on the numerical properties of the method in the spectral space. The amplification factor is derived as a function of the non- dimensional wave numbers. It provides insights into the stability and dissipative properties of the method for various CFL and Fourier numbers. The effect of the grid shape functions, number of particles per cell and their locations inside the grid cell are also analyzed. The EXAGOOP MPM solver (https://github.com/NREL/Exagoop.git) is developed at the National Renewable Energy Laboratory as a part of the NAWI UHPRO project and is based on the AMReX framework. A single-level, uniform cartesian grid is used as the background mesh, while the particle class in AMReX is used to manage the material point operations. Linear hat and B-splines are used as grid shape functions, while the time integration is performed using explicit Euler time integration. EXAGOOP is both CPU and GPU compatible and has been demonstrated to work well on multiple compute architectures. The performance of EXAGOOP on various computing architectures is presented along with its application to study compaction and embossing of high-pressure reverse osmosis membranes. The MPM solution accurately reproduces the membrane deformation. The deformed pore size and structure simulated using MPM also agree well with experimental SEM images.

material point method↗

Power-Domain Impedance Theory for the Analysis and Mitigation of Interarea Oscillations: Preprint

This paper presents an impedance-based theory for the analysis and mitigation of interarea oscillations in power systems. Because low-frequency power system oscillations primarily manifest in phasor quantities, including active and reactive power flows, and the magnitude and frequency of bus voltages, a new-type of impedance, called power-domain impedance, is defined in terms of these phasor quantities. The power-domain impedance provides an intuitive framework for the analysis of interarea oscillations without requiring internal details of generators, and it is ideally suited for developing oscillation damping controls in inverter-based resources. The power-domain impedance theory is demonstrated in this paper using PSCAD simulations of a two-area system. The theory is also used for designing damping control in synchronous generators and a wind power plant in the two-area system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady geothermal field, Nevada, USA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a throughgoing network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density.

58 GEOSCIENCES↗

Code Validation of SAM Using Forced and Natural Circulation Data from NACIE-UP Benchmark

Heavy liquid metals (HLMs) are promising candidates as coolants of Generation IV fast reactors due to their thermophysical properties. In the last decade, experimental work has been proposed as part of research and development efforts to develop such systems. In this context, researchers from the Brasimone Research Center have conducted many experiments using the Natural Circulation Experiment Upgrade (NACIE-UP) facility to study the thermofluid dynamic behavior of HLMs in rod bundle configurations with or without wire wrappers. This facility consists of a rectangular loop operated with lead-bismuth eutectic. Sensors across the loop monitor relevant parameters, i.e. temperatures, heat transfer, and flow conditions. Here, in the present work, we carefully select published data from NACIE-UP to validate the System Analysis Module (SAM), a modern system analysis code developed at Argonne National Laboratory. We developed one SAM model using specifications of the facility geometry and materials existing in relevant papers and reports. On top of that, these references provided the boundary conditions for simulating natural circulation and forced convection experiments in either steady or transient conditions. The SAM model simulates five test cases with diverse operating conditions. Ultimately, the code is proven to predict temperatures and mass flow rates that closely match the experiments. The discrepancies between numerical predictions and diverse transients are limited to a few degrees Celsius, showcasing that SAM is well suited for analyzing nuclear systems relying on HLM coolants.

advanced reactors↗

Intermittent Frontogenesis in the Alboran Sea

Abstract We present a phenomenological description and dynamical analysis of the Alboran fronts using a realistic simulation at submesoscale resolution. The study is focused on east Alboran fronts emerging within relatively strong flows that separate from the Spanish coast into the basin interior. Despite modest lateral shifting associated with the position of the Alboran anticyclonic gyres and variations in intensity, these fronts present a similar structure and dynamical configuration as the climatological Almeria–Oran front. The statistical analysis of our solution shows that strained-induced frontogenesis is a recurrent submesoscale mechanism associated with these fronts, and the process is assessed in terms of the advective Lagrangian frontogenetic tendencies associated with buoyancy and velocity horizontal gradients. Intermittency in their strength and patterns is indicative of high variability in the occurrence of active frontogenesis in association with the secondary (overturning) circulation across the frontal gradient. As a result, we find many episodes with strong surface fronts that do not have much associated downwelling. Frontogenesis and the associated secondary circulation are further explored during two particular frontal events, both showing strong downwelling of (1) cm s −1 extending down into the pycnocline. A frontogenetic contribution of turbulent vertical momentum mixing to the secondary circulation is identified in the easternmost region during the cold season, when the dynamics are strongly influenced by the intrusion of the salty Northern Current. The background vertical velocity fields observed during the analyzed events indicate other currents in the submesoscale range, including tidal and topographic internal waves.

Capó, Esther↗

Normalizing flows for likelihood-free inference with fusion simulations

Fluid-based scrape-off layer transport codes, such as UEDGE, are heavily utilized in tokamak analysis and design, but typically require user-specified anomalous transport coefficients to match experiments. Determining the uniqueness of these parameters and the uncertainties in them to match experiments can provide valuable insights to fusion scientists. Here, we leverage recent work in the area of likelihood-free inference (‘simulation-based inference’) to train a neural network, which enables accurate statistical inference of the anomalous transport coefficients given experimental plasma profile input. UEDGE is treated as a black-box simulator and runs multiple times with anomalous transport coefficients sampled from priors, and the neural network is trained on these simulations to emulate the posterior. The neural network is trained as a normalizing flow model for density estimation, allowing it to accurately represent complicated, high-dimensional distribution functions. With a fixed simulation budget, we compare a single-round procedure to a multi-round approach that guides the training simulations toward a specific target observation. Finally, we discuss the future possibilities for use of amortized models, which train on a wide range of simulations and enable fast statistical inference for results during experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Steady-state off-design modeling of the supercritical carbon dioxide recompression cycle for concentrating solar power applications with two-tank sensible-heat storage

Concentrating solar power researchers are evaluating the potential of the supercritical carbon dioxide recompression cycle to improve the thermal efficiency and decrease the capital costs of next-generation systems. This analysis investigates the steady-state off-design performance of a recompression cycle integrated with two-tank sensible-heat thermal energy storage as the ambient temperature and heat-transfer fluid (HTF) inlet conditions to the cycle change. This paper presents off-design component models and then cycle convergence and control models to maximize net cycle power output while constraining the high-side pressure and air-cooler fan power to their respective design values and fixing the cycle HTF outlet temperature to its design value. Results show that inventory control and air-cooler fan power are important control parameters that can be optimized to maximize off-design power output as the ambient temperature and HTF mass flow rate diverge from design. The high-side pressure and fan power constraints cause the cycle net power to degrade when the ambient temperature is warmer than the design value, and this study calculates a maximum mass flow rate for hot days above which the cycle cannot achieve the design HTF outlet temperature. Finally, the analysis shows that optimizing compressor shaft speeds can improve cold-day performance by around 1.5 percentage points and part-load performance by up to 2.5 percentage points versus a baseline case with no active compressor control.

14 SOLAR ENERGY↗

SOFE2023 Presentation: "Tritium Safety and Sensitivity Analysis of Tritium Extraction eXperiment (TEX) PbLi Loop Using MELCOR-TMAP"

The efficient extraction of tritium from the breeder blanket (BB) in a fusion reactor is vital for maintaining the tritium economy for a sustainable and self-sufficient fusion power plant. The vacuum permeator (VP) design is one of several promising technology platforms for extracting tritium from the BB with the advantages of being capable of operating continuously at high temperatures with the aim of achieving high tritium extraction efficiencies. The Tritium Extraction eXperiment (TEX), under construction at the Safety and Tritium Applied Research (STAR) Facility at Idaho National Laboratory (INL), is designed to test the viability of VP technology by introducing atomic hydrogen (H, D, and T) from a source permeator (SP) into a forced convection lead lithium eutectic (PbLi) loop and determining the hydrogen extraction rates from the vanadium VP. Presented here is a computational analysis of the hydrogen inventory and permeation behavior in TEX using the hydrogen transport and accident analysis code MELCOR-TMAP developed at INL. The results of a sensitivity analysis based on MELCOR-TMAP simulations of hydrogen extraction efficiencies from the TEX test section using varying PbLi flow rates, hydrogen source rates from the SP, and material property uncertainties are also reported. Operation scenarios for TEX in consideration of the analyses based on results from MELCOR-TMAP are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Improving air leakage prediction of buildings using the fan pressurization method with the Weighted Line of Organic Correlation

In many countries, the fan pressurization method is the most frequently chosen approach for measuring the air leakage of houses. The measurements are usually performed at pressures that far exceed pressures to which buildings are exposed to under normal conditions. A fit of these tests to the power-law formulation allows an extrapolation to data points outside the measured pressure range. With the Ordinary Least Square (OLS) fitting method, the pressure exponent and flow coefficient can be determined. However, the measurement results are highly sensitive to uncertainties induced by external factors like changing wind conditions during the tests, which is neglected by OLS. This may lead to errors in the prediction of flows at lower pressures. The Weighted Line of Organic Correlation (WLOC) is an alternative approach and takes measurement uncertainty into account. In this paper, a statistical analysis of an extensive data set of pressurization measurements has been performed. Both regression techniques have been compared for almost 7500 fan pressurization measurements of six houses in 109 different house leak configurations. The variability in predicting pressure exponent and flow coefficient for both WLOC and OLS regression was analyzed using probability density functions. It was found that the Weighted Line of Organic Correlation significantly decreases the uncertainty in predicting pressure exponent, flow coefficient, and other low-pressure air leakage metrics compared to the Ordinary Least Square fitting. Here, the authors highly recommend the implementation of WLOC in current measurement standards and test equipment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Evaluating large scale aqueous organic redox flow battery performance with a hybrid numerical and machine learning framework

Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low-capacity degradation in 10 cm$^2$ cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier to commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm$^2$ DHP-based AORFB by combining a physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. Such combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.

25 ENERGY STORAGE↗

Exploring Classification of Topological Priors With Machine Learning for Feature Extraction

In many scientific endeavors, increasingly abstract representations of data allow for new interpretive methodologies and conceptualization of phenomena. For example, moving from raw imaged pixels to segmented and reconstructed objects allows researchers new insights and means to direct their studies toward relevant areas. Thus, the development of new and improved methods for segmentation remains an active area of research. With advances in machine learning and neural networks, scientists have been focused on employing deep neural networks such as U-Net to obtain pixel-level segmentations, namely, defining associations between pixels and corresponding/referent objects and gathering those objects afterward. Topological analysis, such as the use of the Morse-Smale complex to encode regions of uniform gradient flow behavior, offers an alternative approach: first, create geometric priors, and then apply machine learning to classify. This approach is empirically motivated since phenomena of interest often appear as subsets of topological priors in many applications. Using topological elements not only reduces the learning space but also introduces the ability to use learnable geometries and connectivity to aid the classification of the segmentation target. Here, in this article, we describe an approach to creating learnable topological elements, explore the application of ML techniques to classification tasks in a number of areas, and demonstrate this approach as a viable alternative to pixel-level classification, with similar accuracy, improved execution time, and requiring marginal training data.

97 MATHEMATICS AND COMPUTING↗

Effects of surface runoff and infiltration partition methods on hydrological modeling: A comparison of four schemes in two watersheds in the Northeastern US

The ratio according to which the precipitation/snowmelt is divided into surface and subsurface water is of primary importance to watershed scale hydrological modeling. The present study compared four different surface runoff and infiltration partition methods in the Richards-equation-based SWAT model (RSWAT), i.e., daily curve number method (DC-RSWAT), hourly curve number method (HC-RSWAT), Green-Ampt method (GA-RSWAT), and effective infiltration capacity method (EIC-RSWAT). These four versions of RSWAT and a daily-curve-number-based SWAT (DC-SWAT) were used to simulate daily flow rate (from 2001 to 2015) at the outlets of the Tuckahoe Creek Watershed (TCW) and Greensboro Watershed (GW) of the Chesapeake Bay in Maryland. Parameter sensitivity analysis, calibration, and prediction uncertainty analysis were conducted. The results indicate that surface runoff and infiltration partition methods have pronounced impacts on model performance. Comparative analysis showed that DC-RSWAT outperformed other versions of RSWAT for both calibration (from 2001 to2007) and validation (from 2008 to 2015) periods in the two test watersheds. In general, DC-RSWAT achieved more accurate deterministic flow simulation and reliable uncertainty estimation than other versions. As such, daily-curve-number method is recommended for surface runoff and infiltration partition in RSWAT for simulating stream flow.

hydrology, modeling, streamflow↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗

Involute Working Group – FSI Analysis of Fuel Plates Using Finite Volume and Finite Element Methods

The three involute plate research reactors RHF, HFIR, and FRM II have expressed an interest in using computational software to carry their steady-state safety analysis. Since these tools represent a significant departure from the methods used currently (one-dimensional), the acceptability of the new approach by regulators requires thorough verification and validation of these tools. Therefore, Argonne National Laboratory and the three involute-plate reactors formed an informal group called the Involute Working Group aiming at qualifying computational tools to perform steady-state safety analysis. The present report focuses on a comparison of finite volume and finite element methods to model solids in fluid-structure interaction problems with the goal to estimate the coolant flow-induced fuel plate deflections obtained with the two methods. The finite volume method will be obsoleted in STARCCM+ by the end of 2021, nevertheless, this evaluation is important because the method was used by ANL researchers to model the response of the fuel plates, despite its drawbacks, which are discussed in the report. It was essential to check how those estimates compare to the results obtained with the finite element method that is considered superior for structural analysis. Various geometries, i.e., flat, cylindrical and circle-involute fuel plates, as well as coolant flow speed, were considered. The comparison shows that, independently of the plate geometry, the finite volume method significantly underestimates the deflection as compared to finite element method for coarser meshes. When the discretization is developed as a result of a mesh sensitivity study using finite element method, the result obtained using finite volume method can be a few times smaller than the corresponding finite element method solution. A code-to-code comparison , between STAR-CCM+ and LS-DYNA was included in the analysis. Within the LS-DYNA models, two types of finite element formulations were used: solid and shell finite elements. Mesh sensitivity study showed that both approaches converge to a similar value that was obtained with STAR-CCM+ finite element solver. The evaluation of the computational solvers was extended by adding two benchmark cases from the STAR-CCM+ Verification Suite and presented in the Appendix A. The selected cases are: (1) bending of a cantilever beam under external load, and (2) cylindrical shell deformation analysis, known in the literature as ‘Scordelis-Lo roof’. The problems were solved with finite volume, and finite element methods, and the results confirmed the previously discussed findings. The analysis shows that the finite element solver is superior to the finite volume solver in terms of representation of model geometry and estimating the structural behavior of fuel plates. Depending on the ratio of the load to the flexibility of the plate, the finite volume solver can greatly under- or overestimate the structural response if a very carefully selected mesh is not used.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Energy Analysis of Combi Heat Pump System Configurations for Space Conditioning and Domestic Hot Water Heating in Residential Buildings

Combi heat pump systems, also referred to multifunctional variable refrigerant flow heat recovery (MF-VRFHR) systems, are specifically designed for residential applications to manage both space conditioning and domestic hot water (DHW). They have attracted attention due to their potential for energy conservation through heat recovery. The incorporation of a hot water tank introduces various system configurations, each characterized by distinct pros and cons related to energy efficiency, system stability, and maintenance. Despite this, a critical gap exists as the specific energy performance remains unquantified under diverse operational modes (e.g., heating mode and heat recovery mode). This paper aims to bridge this gap by conducting a comprehensive comparative analysis of two prevalent system configurations while considering feasible proposed control logics. Configuration 1 integrates a separate hot water tank and a refrigerant-to-water heat exchanger (HEX), also known as a Hydro Kit while Configuration 2 incorporates a refrigerant-wrapped hot water tank. To facilitate this analysis, we developed high-fidelity system models for both configurations in Modelica, capturing system dynamics and detailed control sequences effectively. These system models were built upon the TIL library for HVAC equipment components and the Buildings library for residential building thermal load calculations. The validation of the simulation testbed utilized data from experiments conducted in the PNNL lab home for Configuration 1. To establish the simulation testbed for Configuration 2, we extended the modeling setup derived from Configuration 1. This extension specifically involved substituting the separate hot water tank and Hydro Kit with a refrigerant-wrapped hot water tank of similar sizing sourced from an actual product. The simulation analysis of heating-only and heat recovery modes reveals that Configuration 2 not only saves energy and maintains warmer tank temperatures but also demonstrates faster water heating capabilities. This is attributed to decreased energy loss and improved heat transfer. The study encompasses a wide range of scenarios, considering diverse thermal loads and water usage patterns across heating and heat recovery modes. Overall, the comprehensive results indicate that Configuration 2 achieves energy savings ranging from 3.5% to 12.2% compared to Configuration 1, depending on factors such as water usage patterns, thermal loads, and operational modes.

Configuration, Comparison, Multi-functional, Resid↗

Towards a heat- and mass-balanced kinetic model of TATB decomposition

We report TATB (1,3,5-triamino-2,4,6-trinitrobenzene) was thermally degraded by two small-scale analytical methods – simultaneous differential scanning calorimetry and thermogravimetric analysis (SDT) and a hot-stage microscope with Fourier Transform Infrared (FTIR) analysis capabilities. SDT used ramped heating, isothermal soaking, and thermal pretreatment at various conditions. The heat flow and mass loss were monitored during various treatment conditions to derive chemical decomposition kinetics and Arrhenius parameters. FTIR experiments used isothermal heating, and changes were monitored spectroscopically. Solid samples generated at specific conditions were collected from both methods and were analyzed by DMSO extraction followed by chemical speciation by optical and mass spectrometric methods. Characterization provided the following reaction insights: TATB decreases in a sigmoidal pattern in isothermally heated samples. Other soluble products gradually increase in concentration and then abruptly decline in concentration during the second exotherm, such as diamino-dinitro-benzofurazan and amino-nitro-benzodifurazan. FTIR showed gradual changes in the amino and nitro functionality, shifting positions and decreasing intensity for the first 40 min. Then the solid gradually appeared more like an amorphous C with N incorporated, similar to previous studies on thermally degraded TATB-type materials. Extracted residues (DMSO-soluble components removed) examined by FTIR showed an abrupt change in chemical composition between 40- and 45-min isothermal treatment, indicating early forming solids are different than later forming residues. A reliable mass- and energy-balanced global reaction network must include at least two autocatalytic reactions, either in parallel or series, and at least one must have an explicit initiation reaction having a low activation energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prompt D 0 , D + , and D* + production in Pb–Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5.02 TeV

The production of prompt D 0 , D + , and D* + mesons was measured at midrapidity (|y| < 0.5) in Pb–Pb collisions at the centre-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV with the ALICE detector at the LHC. The D mesons were reconstructed via their hadronic decay channels and their production yields were measured in central (0–10%) and semicentral (30–50%) collisions. The measurement was performed up to a transverse momentum (p T ) of 36 or 50 GeV/c depending on the D meson species and the centrality interval. For the first time in Pb–Pb collisions at the LHC, the yield of D 0 mesons was measured down to p T = 0, which allowed a model-independent determination of the p T -integrated yield per unit of rapidity (dN/dy). A maximum suppression by a factor 5 and 2.5 was observed with the nuclear modification factor (R AA ) of prompt D mesons at p T = 6–8 GeV/c for the 0–10% and 30–50% centrality classes, respectively. The D-meson RAA is compared with that of charged pions, charged hadrons, and J/ψ mesons as well as with theoretical predictions. The analysis of the agreement between the measured R AA , elliptic (v 2 ) and triangular (v 3 ) flow, and the model predictions allowed us to constrain the charm spatial diffusion coefficient D s . Furthermore the comparison of R AA and v 2 with different implementations of the same models provides an important insight into the role of radiative energy loss as well as charm quark recombination in the hadronisation mechanisms.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗