Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porosity”

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 955 records · Page 53

Petrophysical and Geochemical Properties of Columbia River Flood Basalt: Implications for Carbon Sequestration

Abstract This study presents borehole geophysical data and sidewall core chemistry from the Wallula Pilot Sequestration Project in the Columbia River flood basalt. The wireline logging data were reprocessed, core-calibrated and interpreted in the framework of reservoir and seal characterization for carbon dioxide storage. Particular attention is paid to the capabilities and limitations of borehole spectroscopy for chemical characterization of basalt. Neutron capture spectroscopy logging is shown to provide accurate concentrations for up to 8 major and minor elements but has limited sensitivity to natural alteration in fresh-water basaltic reservoirs. The Wallula borehole intersected 26 flows from 7 members of the Grande Ronde formation. The logging data demonstrate a cyclic pattern of sequential basalt flows with alternating porous flow tops (potential reservoirs) and massive flow interiors (potential caprock). The log-derived apparent porosity is extremely high in the flow tops (20%-45%), and considerably overestimates effective porosity obtained from hydraulic testing. The flow interiors are characterized by low apparent porosity (0-8%) but appear pervasively fractured in borehole images. Electrical resistivity images show diverse volcanic textures and provide an excellent tool for fracture analysis, but neither fracture density nor log-derived porosity uniquely correlate with hydraulic properties of the Grande Ronde formation. While porous flow tops in these deep flood basalts may offer reservoirs with high mineralization rates, long leakage migration paths, and thick sections of caprock for CO2 storage, a more extensive multi- well characterization would be necessary to assess lateral variations and establish sequestration capacity in this reservoir.

Petrophysical Geochemical Properties Columbia Rive↗

A Pore‐Scale Investigation of Mineral Precipitation Driven Diffusivity Change at the Column‐Scale

Abstract Mineral precipitation affects the pore structure and thus transport properties of porous media. In this study, we investigated the pore‐scale dynamics of precipitation in diffusion controlled systems and the resulting impacts on the effective diffusivity, using a micro‐continuum reactive transport model. Forty two‐dimensional pore structures representing both idealized and realistic geometries were simulated with consideration of different precipitation scenarios and rates. A homogeneous nucleation scenario reproduced patterns observed in previous experimental study showing mixing‐induced precipitation, and a surface growth scenario captured the pattern for mineral precipitation on a substrate with the same or similar mineral structures. In all cases, local precipitation resulted in the reduction in the average porosity of the domain ( Φ ) until the cessation of diffusive transport and the termination of precipitation. The minimum porosity reached was referred to as the critical porosity ( Φ c ). The effective diffusivity ( D eff ) decreased with Φ and dropped sharply to effectively zero, that is, the critical effective diffusivity ( D eff c ), as Φ c was reached. These pore‐scale dynamics can be captured by a revised relationship that explicitly considers the critical porosity and the corresponding effective diffusivity, and the pre‐exponential coefficient and the exponent of the relationship varied with initial pore structure and the precipitation kinetics. Overall, the homogeneous nucleation scenario results in systematically larger Φ c and coefficients that give rise to a sharper decrease in diffusivity as Φ c is approached, compared to the surface growth scenario. The revised relationship was also implemented at continuum scale and used to examine column scale diffusivity change and reactions.

58 GEOSCIENCES↗

Keyhole fluctuation and pore formation mechanisms during laser powder bed fusion additive manufacturing

Keyhole porosity is a key concern in laser powder-bed fusion (LPBF), potentially impacting component fatigue life. However, some keyhole porosity formation mechanisms, e.g., keyhole fluctuation, collapse and bubble growth and shrinkage, remain unclear. Using synchrotron X-ray imaging we reveal keyhole and bubble behaviour, quantifying their formation dynamics. The findings support the hypotheses that: (i) keyhole porosity can initiate not only in unstable, but also in the transition keyhole regimes created by high laser power-velocity conditions, causing fast radial keyhole fluctuations (2.5-10 kHz); (ii) transition regime collapse tends to occur part way up the rear-wall; and (iii) immediately after keyhole collapse, bubbles undergo rapid growth due to pressure equilibration, then shrink due to metal-vapour condensation. Concurrent with condensation, hydrogen diffusion into the bubble slows the shrinkage and stabilises the bubble size. The keyhole fluctuation and bubble evolution mechanisms revealed here may guide the development of control systems for minimising porosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grain boundary widening controls siderite (FeCO3) replacement of limestone (CaCO3)

Abstract The microstructure of minerals and rocks can significantly alter reaction rates. This study focuses on identifying transport paths in low porosity rocks based on the hypothesis that grain boundary widening accelerates reactions in which one mineral is replaced by another (replacement reaction). We conducted a time series of replacement experiments of three limestones (CaCO 3 ) of different microstructures and solid impurity contents using FeCl 2 . Reacted solids were analyzed using chemical imaging, small angle X-ray and neutron scattering and Raman spectroscopy. In high porosity limestones replacement is reaction controlled and complete replacement was observed within 2 days. In low porosity limestones that contain 1–2% dolomite impurities and are dominated by grain boundaries, a reaction rim was observed whose width did not change with reaction time. Siderite (FeCO 3 ) nucleation was observed in all parts of the rock cores indicating the percolation of the solution throughout the complete core. Dolomite impurities were identified to act as nucleation sites leading to growth of crystals that exert force on the CaCO 3 grains. Widening of grain boundaries beyond what is expected based on dissolution and thermal grain expansion was observed in the low porosity marble containing dolomite impurities. This leads to a self-perpetuating cycle of grain boundary widening and reaction acceleration instead of reaction front propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Resolving experimental biases in the interpretation of diffusion experiments with a user-friendly numerical reactive transport approach

The reactive transport code CrunchClay was used to derive effective diffusion coefficients (D e ), clay porosities (ε), and adsorption distribution coefficients (K D ) from through-diffusion data while considering accurately the influence of unavoidable experimental biases on the estimation of these diffusion parameters. These effects include the presence of filters holding the solid sample in place, the variations in concentration gradients across the diffusion cell due to sampling events, the impact of tubing/dead volumes on the estimation of diffusive fluxes and sample porosity, and the effects of O-ring-filter setups on the delivery of solutions to the clay packing. Doing so, the direct modeling of the measurements of (radio)tracer concentrations in reservoirs is more accurate than that of data converted directly into diffusive fluxes. While the above-mentioned effects have already been described individually in the literature, a consistent modeling approach addressing all these issues at the same time has never been described nor made easily available to the community. A graphical user interface, CrunchEase, was created, which supports the user by automating the creation of input files, the running of simulations, and the extraction and comparison of data and simulation results. While a classical model considering an effective diffusion coefficient, a porosity and a solid/solution distribution coefficient (D e –ε–K D ) may be implemented in any reactive transport code, the development of CrunchEase makes it easy to apply by experimentalists without a background in reactive transport modeling. CrunchEase makes it also possible to transition more easily from a D e –ε–K D modeling approach to a state-of-the-art process-based understanding modeling approach using the full capabilities of CrunchClay, which include surface complexation modeling and a multi-porosity description of the clay packing with charged diffuse layers.

58 GEOSCIENCES↗

Colloidal Three-Dimensional Covalent Organic Frameworks and Their Application as Porous Liquids

Developing solid porous materials into free-flowing liquids with permanent porosity is a promising strategy for overcoming certain limitations of conventional sorbent materials employed in gas storage and separation applications. The ability to control the pore size and chemical functionalities of organic frameworks gives these particular nanoporous materials distinct advantages over other small cage-like molecules or hollow particles that have been developed into porous liquids. Here, we describe the synthesis of a 3D imine-linked colloidal covalent organic framework (COF)-based porous liquid, designed for efficient size-exclusion of solvent, as well as for long-term stability. By tethering ionic liquids to the colloid surface, the colloids can be dried, purified, and resuspended in a variety of solvents without irreversible aggregation typically observed of COF colloids. Colloid size could be controlled between 50 and 400 nm, with surface areas as high as 800 m2/g. The 3D intertwining morphology of the colloids had pore sizes ranging from 5 to 14 Å, allowing them to efficiently sizeexclude bulky ionic liquids. The COF colloids were stable towards flocculation in an ionic liquid for > 1 year. Permanent porosity was confirmed with a combination of 19F NMR measurements and gas sorption techniques. CO2 and CH4 uptake in these porous liquids increased more than 10 and 20-fold, respectively, over non-porous, neat ionic liquid control samples. The work not only advances the state of COF-based colloid science but also represents a practical advance towards developing more robust, tunable framework-based porous liquid materials for a host of gas storage and separation applications.

Mow, Rachel E.↗

Microstructure characterization of electric field assisted sintering (EFAS) sintered metallic and ceramic materials using local thermal diffusivity measurement

Electric Field Assisted Sintering (EFAS, also referred to as spark plasma sintering) is a powerful technology for the consolidation of powder materials. The high heating rate during the sintering process is critical for minimizing energy consumption, but it can also cause microstructure heterogeneities in sintered parts, such as spatially varied porosity. The examination of localized porosity usually requires the use of a scanning electron microscope with a carefully prepared surface. In this paper, photothermal radiometry is used to measure local thermal diffusivity and extract localized porosity of EFAS-sintered parts by using a percolation-threshold model. Applying this approach, we identified the radial position-dependent porosity variation in EFAS parts, which is likely formed due to the large temperature gradient during the sintering process. This approach has a unique advantage because it can measure samples with minimal or no surface preparation, enabling the possibility of in situ characterization in EFAS with proper system modification. Necessary modifications on the measurement approach for EFAS deployment and in situ characterization are also discussed.

36 MATERIALS SCIENCE↗

Macroscopic modeling of gas permeability in hierarchical micro/nanoporous media: A unified characterization of rarefaction using Klinkenberg theory and equivalent diameter

Estimating gas transport through a hierarchical micro/nanoporous system is challenging due to non-equilibrium gas dynamics. The primary difficulty lies in determining the rarefaction level, because identifying a representative flow dimension in a complex porous system with multiple pore scales is not straightforward. Our study performed a pore-level analysis for gas permeability in dual-scale porous media with varying porosity, throat size, and secondary pore size under different rarefaction conditions. We found that secondary porosity negatively affects permeability due to increased friction forces, with this influence growing as the secondary pore size and porosity increase until the secondary pore becomes comparable to the throat. However, rarefaction reduces the effects of secondary pores due to boundary slip. Traditional Knudsen number (Kn) calculations based on Darcy-defined height failed to accurately describe the rarefaction effects on gas permeability. Instead, we introduced an equivalent diameter to calculate the Kn, which provided an accurate normalization of apparent gas permeability independent of pore geometry. Furthermore, the extended Kozeny–Carman–Klinkenberg model developed in our previous study successfully yielded a macroscopic model for apparent gas permeability in hierarchical micro/nanoporous systems as a function of the traditional Darcy height and porosity.

04 OIL SHALES AND TAR SANDS↗

A Flux Detection Probe to Quantify Dynamic Groundwater-Surface Water Exchange in the Hyporheic Zone

A new probe was designed to quantify groundwater-surface water exchange in the hyporheic zone under dynamic stage condition. Current methods focus on either vertical pore water velocity or Darcy flux measurements. Both parameters must be understood to evaluate residence time and mass flux of constituents. Furthermore, most instruments are not well suited for monitoring instantaneous velocity or flux under dynamic exchange conditions. For this reason, the flux detection probe (FDP) was designed that employs electrogeophysical measurements to estimate in situ sediment porosity, which can be used to convert pore water velocity to Darcy flux. Dynamic pore water velocity is obtained by monitoring fluid conductivity and temperature along the FDP probe. Pressure sensors deployed at the top and bottom of the probe provide the additional information necessary to estimate vertical permeability. This study focuses on the use of a geophysical method to estimate pore water velocity, porosity, and permeability within a controlled soil column where simulated river water displaces simulated groundwater. The difference between probe derived and theoretical pore water velocity using natural tracers such as electrical conductivity and temperature was -4.9 and 3.9% for downward flow and 1.1 and 12.8% for upward flow, respectively. The difference in porosity calculated from mass and volume packed in the soil column and probe measure porosity ranged between -3.2% and 1.5%. Also, the calculated hydraulic conductivity differed from probe derived values by -8.9%.

54 ENVIRONMENTAL SCIENCES↗

Foaming during Nuclear Waste Melter Feeds Conversion to Glass: Application of Evolved Gas Analysis

During the final stages of batch-to-glass conversion in a waste-glass melter, gases evolving in the cold cap produce primary foam, the formation and collapse of which control the glass production rate via its effect on heat transfer to the reacting batch. We performed quantitative evolved gas analysis (EGA) for several HLW melter feeds with temperatures ranging from 100 °C to 1150 °C, the whole temperature span in a cold cap. EGA results were supplemented with visual observation of batch-to-glass transition using the feed expansion tests. Upon heating, most of the gases—mainly H2O, CO2, NO, NO2, N2 and O2—evolve at temperatures below 700 °C and escape directly to the atmosphere through open porosity. However, as open porosity closes when enough glass-forming melt appears at 720 °C, the residual gas evolution leads to the formation of primary foam. We found that primary foaming is mostly caused by the decomposition of residual carbonates, though oxygen evolution from iron-redox reaction can also play a role. We also show that the gas evolution shifts to a higher temperature when the heating rate increases. The implications for the mathematical modeling of foam layer in the cold cap are presented.

Hujova, Miroslava↗

Computational Evaluation of Thermal Response of Open-Cell Foam With Circular Pore (Computional Evaluation of Thermal Response of Open-Cell Foam with Circular Pore)

The evaluation of effective material properties in heterogeneous materials (e.g., composites or multicomponent structures) critically relevant to a wide spectrum of applications, including nuclear power, electronic packaging, flame retardants, hypersonics, and gas turbine power. The work described in this paper is centered around the numerical assessment of the thermal behavior of porous materials obtained from finite element thermal modeling and simulation. Here, two-dimensional, steady state analyses were performed on unit cells with centered, circular pores using a second order accurate Galerkin finite element method (FEM). The effective thermal conductivities of the porous systems were examined, encompassing a range of porosities from 4.9% to 60.1%. The geometries of the models were generated based on ordered circular pores for each modeled porosity level. The system response quantity (SRQ) under investigation was the dimensionless effective thermal conductivity across the unit cell. The dimensionless effective thermal conductivity was compared across all simulated cases, producing a trend between porosity and effective thermal conductivity. In the presented investigation, the method of manufactured solutions (MMS) was used to perform code verification, and the grid convergence index (GCI) was employed to estimate discretization uncertainty as solution verification. Code verification concluded an approximately second order accurate Galerkin FEM solver. It was found that the introduction of porosity to the unit cell material structure reduces effective thermal conductivity, as anticipated. Numerical results obtained in this study are compared to an analytical solution and to a sample of empirical data.

36 MATERIALS SCIENCE↗

Additive Manufacturing of Porous Ceramics With Foaming Agent

The macro-porous ceramics has promising durability and thermal insulation performance. As porous ceramics find more and more applications across many industries, a cost-effective and scalable additive manufacturing technique for fabricating macro-porous ceramics is highly desirable. Herein, we reported a facile additive manufacturing approach to fabricate porous ceramics and control the printed porosity. Several printable ceramic inks were prepared, and the foaming agent was added to generate gaseous bubbles in the ink, followed by the direct ink writing and the ambient-pressure and room-temperature drying to create the three-dimensional geometries. In this work, a set of experimental studies were performed to optimize the printing quality. The results revealed the optimal process parameters for printing the foamed ceramic ink with a high spatial resolution and fine surface quality. Varying the concentration of the foaming agent enables the controllability of the structural porosity. The maximum porosity can reach 85%, with a crack-free internal porous structure. The tensile tests showed that the printed macro-porous ceramics possessed enhanced durability with the addition of fiber. With a high-fidelity three-dimensional (3D) printing process and the precise controllability of the porosity, we showed that the printed samples exhibited a remarkably low thermal conductivity and durable mechanical strength.

36 MATERIALS SCIENCE↗

Manufacturing Cathodes via Dry-Processing for Lithium-Ion Batteries

Conventional lithium-ion battery (LIB) electrodes are prepared through a wet slurry process with n-methyl pyrrolidone solvent, especially for cathodes. The wet slurry process encounters several disadvantages such as binder migration, electrode cracking in thick electrodes, energy intense heat-dry NMP solvent removal, and costly NMP recovery. The cost and energy consumption of coating and drying of electrode are about 11.5 % and > 46 % in LIB manufacturing, respectively. Thereby, it is essential to develop a facile roll-to-roll solvent-free LIB electrode processing for reducing the cost and energy consumption. Recently, the Maxwell-type dry processing (DP) shines new lights on LIB manufacturing, which mainly bases on dry mixing (DM) of electrode component powder followed by calendering into electrode films and laminating onto current collectors, realizing the rapid manufacturing of LIB electrodes in a powder-to-film manner for industries. This report shares some recent progress on the DP from our group. We aim to further advance the manufacturing science of DP by correlating the processing conditions with electrode properties and performance. Particularly, we investigate the effect of DM, and compression on the polytetrafluoroethylene (PTFE) binder fiberization, porosity, mechanical properties, electrical conductivity and electrochemical behaviors of electrodes. The DM study suggests that PTFE fiberization heavily relies on the degree of DM. Insufficient DM results in poor PTFE fiberization while outrageous DM damages the formed PTFE fibers. Both negatively affect the mechanical behaviors of the electrodes and their rate capability. However, moderate DM is highly beneficial. In addition, our study of the porosity impact reveals that LiNi0.8Mn0.1Co0.1O2 (NMC) secondary particles can be broken into primary particles due to compression, especially at low porosity. Those fractured NMC secondary particles exhibit lower modulus. We propose that a moderate porosity of around 32% favors the electronic conductivity, charge transfer impedance and rate capability. The study of the cathodic electrolyte interphase layer of PTFE-based DPed electrode confirms that side reactions of PTFE binder due to the formation of LiF in LiClO4-based electrolyte.

Tao, Runming↗

Modeling Battery Performance Due to Intercalation Driven Volume Change in Porous Electrodes

Simulations are presented that result from incorporating dimensional and porosity changes in porous electrodes caused by volume changes in the active material during intercalation into a detailed lithium-ion battery model. Porosity and dimensional changes in an electrode can significantly affect the resistance of the battery during cycling, which in turn alters the reaction distributions in the porous electrodes. In addition, volume changes generate stresses in the electrode which can lead to premature failure of the battery. Material conservation equations are coupled with the mechanical properties of porous electrodes to link dimensional and porosity changes to stresses and the resulting resistances that occur during the intercalation processes. Through the use of porous rock mechanics, porosity and strain gradients can be predicted based on state of discharge and discharge rate. Several different battery casings and discharge rates are examined and operating curves are predicted.

25 ENERGY STORAGE↗

Additive Manufacturing of Controlled Porous Elastomeric Nanocomposites for Enhanced Sensing Function

Sensors manufactured using piezoresistive elastomeric nanocomposites have a wide range of applications in many fields such as structural health monitoring, robotics, and biomedical industries. Additive manufacturing of these sensors, achieved using Direct Ink Writing (DIW), has provided numerous advantages including increased sensitivity and the fabrication of complex geometry. However, solid sensors manufactured either additively or using conventional methods, display a greater level of hysteresis, especially at strains higher than 10%, which limits their sensitivity. The solid sensors also display a greater level of piezoresistive relaxation where the resistance change with time either in a static or dynamic loading environment. Porosity can be introduced in the form of varying infill patterns and densities to reduce bulk material effects. In this paper, a lattice structure with four infill densities are investigated, and the samples are printed using a grid infill pattern. The fabricated samples are characterized using a scanning electron microscope (SEM) to validate the microstructural features and layer bonding. The pressure-sensing capability of each sensor is investigated using cyclic compression loading at various strain levels. Sensing results show that an increase in sensitivity, as well as a decrease in hysteresis with an increase in porosity, compared to solid samples (100% in-fill). The highest porosity was measured to be 62.3% with a corresponding 112% increase in the gauge factor compared to the solid bulk. These results indicate that introducing porosity using DIW is a sensible strategy to improve the piezoresistive performance of nanocomposites and to allow for the tunability of sensing capacity in pressure sensors

Sukati, Musa↗

Investigation and modeling of the degradation of zirconium-based fuel cladding during corrosion in steam and air-steam mixtures at high temperatures - 2016-0033

This paper presents the main results of a study conducted to quantify and to model the degradation state of zirconium-based fuel claddings submitted to severe-accident conditions in a nuclear reactor core: high temperatures and either pure steam or an air-steam mixture. Due to the progressive thickening of a dense and protective ZrO{sub 2} layer, the oxidation kinetics of zirconium-based claddings in steam at high temperatures typical of severe nuclear accidents are generally cubic or parabolic. However, for some temperature domains, this oxide layer may crack, becoming porous and no longer protective. In these 'breakaway' conditions, the oxidation kinetics change from (sub)parabolic to linear or even accelerated. In addition, the temperature increase can lead core materials to melt and to relocate down to the vessel lower head, threatening its integrity. If it fails, and for specific conditions, air ingress may take place into the reactor. Hence, oxygen and nitrogen both react with zirconium-based claddings successively through the oxidation of zirconium (forming a ZrO{sub 2} layer), nitriding of zirconium (forming zirconium nitride particles), and the oxidation of zirconium nitride (forming ZrO{sub 2} and releasing nitrogen). These self-sustained chemical reactions enhance the deterioration of zirconium-based claddings and their ZrO{sub 2} layers, inducing an increase in their open porosity. To quantify this porosity, a series of two-step experiments was conducted. First, ZIRLO{sup TM} cladding samples were isothermally oxidized in pure steam or in a 50:50 mol% air-steam mixture at several different temperatures and durations. The main thermal effects on reaction kinetics and the high impact of air on the cladding degradation were all confirmed by experimental results. Second, pioneering porosimetry measurements by mercury intrusion were realized for the first time on such corroded cladding samples. In both atmospheres, it was pointed out that 1,200 and 1,250 K led to particularly porous oxide layers, especially due to strong breakaway effects. Moreover, we confirmed that the presence of air strongly enhances the oxide cracking: cladding samples were more porous when oxidized in the air-steam mixture than under pure steam. Finally, we observed that in all conditions, the open porous volume fraction of ZIRLO claddings continuously rose during their corrosion process. Hence, for each experimental condition, porosity correlations were determined through linear regressions, and porosity increase rates were deduced by derivation versus time and validated against porosimetry results of cladding samples corroded in transient nonisothermal) conditions. (authors)

36 MATERIALS SCIENCE↗

Sandy Soil Microaggregates: Rethinking Our Understanding of Hydraulic Function

Core Ideas Sandy soil microaggregates are comprised of a solid sand core encrusted with fines. Thick films control a significant proportion of water retention in coarse sandy soil. Three water pools drain sequentially: capillary water, thick films, and adsorptive films. Film seepage rates should be considered when irrigating sandy soil. This study investigated the peculiar structure of microaggregates in coarse sandy soils that exhibit only external porosity and investigated their control on soil hydrology. The microstructure underpins a hydrologic existence that differs from finer textured soils where aggregates have internal porosity. Understanding the impact of these microaggregates on soil hydrology will permit improved agricultural irrigation management and estimates associated with ecosystem capacity and resiliency. Microstructure was investigated using a digital microscope, and aspects of the structure were quantified by sedimentation and computed microtomography. Sandy soil microaggregates were observed to be comprised of a solid sand‐grain core that is coated with fines, presumably cemented by organic media. This microstructure leads to three distinct water pools during drainage: capillary water, followed by thick films (1–20 μm) enveloping the outer surfaces of the crusted microaggregates, followed by adsorbed thin films (<1 μm). The characteristics of the thick films were investigated using an analytical model. These films may provide as much as 10 to 40% saturation in the range of plant‐available water. Using lubrication theory, it was predicted that thick film drainage follows a power law function with an exponent of 2. Thick films may also have a role in the geochemical evolution of soils and in ecosystem function because they provide contiguous water and gas phases at relatively high moisture contents. And, because the rough outer crust of these microaggregates can provide good niches for microbial activity, biofilm physics will dominate thick film processes, and consequently hydrologic, biologic, and geochemical functions for coarse sandy soils.

Paradiś, Ashley↗

Modeling the thermal deformation of TATB-based explosives. Part 1: Thermal expansion of “neat-pressed” polycrystalline TATB

We detail a modeling approach to simulate the anisotropic thermal expansion of polycrystalline (1,3,5-triamino-2,4,6-trinitrobenzene) TATB-based explosives that utilizes microstructural information including porosity, crystal aspect ratio, and processing-induced texture. This report, the first in a series, focuses on nonlinear thermal expansion of “neat-pressed” polycrystalline TATB specimens which do not contain any binder; additional complexities related to polymeric binder and irreversible ratcheting behavior are briefly discussed, however detailed investigation of these aspects are deferred to subsequent reports. In this work we have, for the first time, developed a mesoscale continuum model relating the thermal expansion of polycrystal TATB specimens to their microstructural characteristics. A self-consistent homogenization procedure is used to relate macroscopic thermoelastic response to the constitutive behavior of single-crystal TATB. The model includes a representation of grain aspect ratio, porosity, and crystallographic texture attributed to the consolidation process. A quantitative model is proposed to describe the evolution of preferred orientation of graphitic planes in TATB during consolidation and an algorithm constructed to develop a discrete representation of the associated orientation distribution function. Analytical and numerical solutions using this model are shown to produce textures consistent with previous measurements and characterization for isostatic and uniaxial “die-pressed” specimens. Predicted thermal strain versus temperature for textured specimens are shown to be in agreement with corresponding experimental measurements. Using the developed modeling approach, several simulations have been run to investigate the influence of microstructure on macroscopic thermal expansion behavior. Results from these simulations are used to identify qualitative trends. Implications of the identified trends are discussed in the context of thermal deformation of engineered components whose consolidation process is generally more complex than isostatic or die-pressed specimens. Finally, an envisioned application of the modeling approach to simulating thermal expansion of weapon systems and components is outlined along with necessary future work to introduce the effects of binder and ratcheting behavior. Key conclusions from this work include the following. Both porosity and grain aspect ratio have an influence on the thermal expansion of polycrystal TATB considering realistic material variability. Thepreferred orientation of the single crystal TATB [001] poles within a polycrystal gives rise to pronounced anisotropy of the macroscopic thermal expansion. The extent of this preferred orientation depends on the magnitude of deformation, and consequently, is expected to vary spatially throughout manufactured components much like porosity. The modeling approach presented here has utility toward bringing spatially variable microstructural features into macroscale system engineering modelsAbstract Not Provided

36 MATERIALS SCIENCE↗