Search NASA⌕ Search

SEARCH · Search NASA

Results for “loading pattern optimization”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Model Study of Blast Furnace Operation with Central Coke Charging

Blast furnace (BF) remains the dominant ironmaking process worldwide. Central coke charging (CCC) operation is a promising technology for stabilizing BF operations, but it needs reliable and quantified process design and control. In this work, a multi-fluid BF model is further developed for quantitatively investigating flow-thermal-chemical phenomena of a BF under CCC operation. This model features the respective chemical reactions in the respective coke and ore layers, and a specific sub-model of layer profile for the burden structure for the CCC operation. The simulation results confirm that the gas flow patterns and cohesive zone’s shape and location under the CCC operation are quite different from the non-CCC operation. Under the CCC operation, the heat is overloaded at the furnace center while the reduction load is much heavier at the periphery regions; the profiles of top gas temperature and gas utilization show bell-shape and inverse-bell-shape patterns, respectively. More importantly, these differences are characterized quantitatively. In this given case, when the CCC opening radius at the throat is 0.35 m, the cohesive zone top opening radius is around 0.50 m, and the isotherms of CCC operation become much steeper (~ 80 deg) than those of non-CCC operation (~ 60 deg) near BF central regions. In addition, it is confirmed that carbon solution-loss reaction rate can be decreased significantly at BF central regions under CCC operation. The model helps to understand CCC operation and provides a cost-effective method for optimizing BF practice.

36 MATERIALS SCIENCE↗
NASA NTRS1980

An assessment of alternative fuel cell designs for residential and commercial cogeneration

A comparative assessment of three fuel cell systems for application in different buildings and geographic locations is presented. The study was performed at the NASA Lewis Center and comprised the fuel cell design, performance in different conditions, and the economic parameters. Applications in multifamily housing, stores and hospitals were considered, with a load of 10kW-1 MW. Designs were traced through system sizing, simulation/evaluation, and reliability analysis, and a computer simulation based on a fourth-order representation of a generalized system was performed. The cells were all phosphoric acid type cells, and were found to be incompatible with gas/electric systems and more favorable economically than the gas/electric systems in hospital uses. The methodology used provided an optimized energy-use pattern and minimized back-up system turn-on.

Wakefield, R. A.↗
DOE OSTI2019

Status Update of a Harmonic Kicker Development for JLEIC

The Jefferson Lab Electron-Ion Collider (JLEIC) is an asymmetric high luminosity ring-ring collider proposed as the next major R&D facility for the nuclear physics community. Both of JLEIC?s electron and ion collider rings have high beam current with gaps serving the pur-poses of beam abort, ion clearing, etc. Those gaps will cause time-varying beam loading in the storage RF cavi-ties, and may result in modulation of cavity RF phase/voltage and cyclic shift of collision point, which could lead to significant luminosity loss. We studied a few approaches to mitigate the RF phase modulation and IP shift, such as correcting the RF phase/voltage modula-tion with traditional LLRF feedback, one-turn feedback (OTFB), or RF feedforward (FF); optimizing the bunch fill pattern to limit the RF phase/voltage modulation to a small fraction of the bunch trains in the collider ring; or matching the RF phase modulation in the two rings. The preliminary results are discussed in this paper.

Guo, Jiquan↗
DOE OSTI2020

Efficient Data Management in Neutron Scattering Data Reduction Workflows at ORNL

Oak Ridge National Laboratory (ORNL) experimental neutron science facilities produce 1.2 TB a day of raw event-based data that is stored using the standard metadata-rich NeXus schema built on top of the HDF5 file format. Performance of several data reduction workflows is largely determined by the amount of time spent on the loading and processing algorithms in Mantid, an open-source data analysis framework used across several neutron sciences facilities around the world. The present work introduces new data management algorithms to address identified input output (I/O) bottlenecks on Mantid. First, we introduce an in-memory binary-tree metadata index that resemble NeXus data access patterns to provide a scalable search and extraction mechanism. Second, data encapsulation in Mantid algorithms is optimally redesigned to reduce the total compute and memory runtime footprint associated with metadata I/O reconstruction tasks. Results from this work show speed ups in wall-clock time on ORNL data reduction workflows, ranging from 11% to 30% depending on the complexity of the targeted instrument-specific data. Nevertheless, we highlight the need for more research to address reduction challenges as experimental data volumes increase.

Godoy, William↗
DOE OSTI2021

3D Printed, Low Tortuosity Garnet Framework For Beyond 500 Wh/kg Batteries

In this project, we developed LLZO garnet ink recipes and processes for 3D-printing highly ordered ionically conductive garnet porous structures on dense garnet separators. Using this technique, we are able to fabricate controlled architecture LLZO garnet solid-state electrolyte (SSE) trilayers for application in solid-state lithium batteries. The trilayer comprises a thin dense center layer sandwiched between a 3D-printed patterned porous layer and a random porous layer. The dense layer functions as the ionic separator between the anode and cathode. The random porous layer hosts the lithium-metal anode and provides the structural support. The 3D-printed SSE patterned porous layer hosts the cathode, providing continuous, low tortuosity pathways for fast 3D Li+ transport through the cell while increasing the electrode/electrolyte interface area to decrease the interfacial resistance. Compared to the random porous structure, this ordered patterned structure possesses more vacant space for higher cathode loading without sacrificing ionically conducting capability, thus potentially greatly increasing the cell energy density. For demonstration purposes, we developed two patterns for the 3D-printed porous layer: grids and columns, for hosting sulfur and NMC cathode, respectively. The corresponding two types of cells were fabricated and tested, and have demonstrated achievement of theoretical discharge capacity without cathode calendaring. In addition, we developed a fundamental solid-state ionic and electronic transport model to optimize the 3D-printed structures for maximum energy and power density. The model was validated by experiment and provides the critical design criteria for achieving the >500 Wh/kg energy goal as function of C-rate.

25 ENERGY STORAGE↗
DOE OSTI2026

Time Matters: A Survival Analysis of Public Electric Vehicle Charging Infrastructure Utilization

The rapid adoption of plug-in electric vehicles (PEVs) places significant demands on public charging infrastructure, making it critical to understand and optimize charger utilization. This study provides one of the most comprehensive analyses of charging behavior to date by applying a survival analysis to a dataset of nearly 16 million level 2 (L2) and direct current (DC) fast charger sessions across the United States from 2017 to 2022. Using Kaplan-Meier curves and log rank tests, our analysis reveals statistically significant and distinct duration patterns influenced by charger type, time of day, and day of the week. We find that L2 charging sessions exhibit high variability tied to venue type, whereas DC sessions are more uniform, typically lasting 30-45 min. This study introduces the operational efficiency score (OES), a metric for standardizing the performance evaluation of charging stations. Our findings offer actionable insights for optimizing charger deployment, developing dynamic pricing strategies to reduce vehicle dwell time, and improving load management for grid operators, ultimately enhancing the efficiency and availability of public charging infrastructure.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1999

The Effects of Partial Mechanical Loading and Ibandronate on Skeletal Tissues in the Adult Rat Hindquarter Suspension Model for Microgravity

We report initial data from a suspended rat model that quantitatively relates chronic partial weightbearing to bone loss. Chronic partial weightbearing is our simulation of the effect of limited artificial gravity aboard spacecraft or reduced planetary gravity. Preliminary analysis of bone by PQCT, histomorphometry, mechanical testing and biochemistry suggest that chronic exposure to half of Earth gravity is insufficient to prevent severe bone loss. The effect of episodic full weightbearing activity (Earth Gravity) on rats otherwise at 50% weightbearing was also explored. This has similarity to treatment by an Earth G-rated centrifuge on a spacecraft that normally maintained artificial gravity at half of Earth G. Our preliminary evidence, using the above techniques to analyze bone, indicate that 2 hours daily of full weightbearing was insufficient to prevent the bone loss observed in 50% weightbearing animals. The effectiveness of partial weightbearing and episodic full weightbearing as potential countermeasures to bone loss in spaceflight was compared with treatment by ibandronate. Ibandronate, a long-acting potent bisphosphonate proved more effective in preventing bone loss and associated functionality based upon structure than our first efforts at mechanical countermeasures. The effectiveness of ibandronate was notable by each of the testing methods we used to study bone from gross structure and strength to tissue and biochemistry. These results appear to be independent of generalized systemic stress imposed by the suspension paradigm. Preliminary evidence does not suggest that blood levels of vitamin D were affected by our countermeasures. Despite the modest theraputic benefit of mechanical countermeasures of partial weightbearing and episodic full weightbearing, we know that some appropriate mechanical signal maintains bone mass in Earth gravity. Moreover, the only mechanism that correctly assigns bone mass and strength to oppose regionally specific force applied to bone is mechanical, a process based upon bone strain. Substantial evidence indicates that the specifics of dynamic loading i.e. time-varying forces are critical. Bone strain history is a predictor of the effect that mechanical conditions have on bone structure mass and strength. Using servo-controlled force plates on suspended rats with implanted strain gauges we manipulated impact forces of ambulation in the frequency (Fourier) domain. Our results indicate that high frequency components of impact forces are particularly potent in producing bone strain independent of the magnitude of the peak force or peak energy applied to the leg. Because a servo-system responds to forces produced by the rat's own muscle activity during ambulation, the direction of ground-reaction loads act on bone through the rat's own musculature. This is in distinction to passive vibration of the floor where forces reach bone through the natural filters of soft tissue and joints. Passive vibration may also be effective, but it may or may not increase bone in the appropriate architectural pattern to oppose the forces of normal ambulatory activity. Effectiveness of high frequency mechanical stimulation in producing regional (muscle directed) bone response will be limited by 1. the sensitivity of bone to a particular range of frequencies and 2. the inertia of the muscles, limiting their response to external forces by increasing tension along insertions. We have begun mathematical modeling of normal ambulatory activity. Effectiveness of high frequency mechanical stimulation in producing regional (muscle directed) bone response will be limited by 1. the sensitivity of bone to a particular range of frequencies and 2. the inertia of the muscles, limiting their response to external forces by increasing tension along insertions. We have begun mathematical modeling of the rat forelimb as a transfer function between impact force and bone strain to predict optimal dynamic loading conditions for this system. We plan additional studies of mechanical counter-measures that incorporate improved dynamic loading, features relevant to anticipated evaluation of artificial gravity, exercise regimens and exposure to Martian gravity, The combination of mechanical countermeasures with ibandronate will also be investigated for signs of synergy.

Schultheis, Lester W.↗
DOE OSTI2024

The Gene Fitness Atlas: A Roadmap for Predicting Evolution

We developed a novel, high-throughput microfluidic device design containing “interaction zones” where progeny cell lines compete against each other allowing for accurate analysis of bacterial cell fitness. The goal of the project was to use the device for two applications: 1) gene knockout screening and 2) antibiotic resistance screening. The microfluidic platform was fabricated using photolithography and soft lithography in polydimethylsiloxane (PDMS). E.coli Keio mutants and fluorescent wildtype parent were chosen for the study. Cells were grown overnight and their loading into the devices and seeding in mother machines was optimized. For mutant screening, the least fit mutant and wildtype parent were cultured individually and then added to the microfluidic device. The mother machines which were seeded with mutant and wildtype were imaged through time lapse microscopy and the growth of cells was observed. For antibiotic screening, wildtype E.coli cells which were grown overnight were added to the device and washed with media containing the antibiotic ampicillin. The growth pattern in presence and absence of ampicillin was observed through time lapse microscopy. It was observed that over a period of four hours, both the mutant and the wildtype divided in the mother machine and pushed daughter cells out into the interaction zone. In case of the antibiotic screening experiment, the fluorescent wildtype divided both in the absence and presence of sublethal concentration of ampicillin. This study is a proof of concept demonstration of high- throughput single cell analysis of cells using a novel microfluidics device.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Quantifying agricultural productive use of energy load in Sub-Saharan Africa and its impact on microgrid configurations and costs

The use of advanced energy technologies for agricultural purposes—such as irrigation, refrigeration, crop processing, and egg incubation—has the potential to increase crop yield, reduce vulnerability to changing precipitation patterns, increase shelf life, strengthen income and employment opportunities in rural areas, and reduce emissions by displacing fossil fuel-based technologies. These productive uses of energy (PUE) in remote areas could potentially be powered by microgrids that additionally serve otherwise unelectrified communities, most of which are located in rural Sub-Saharan Africa. Here, in this paper, we use high-resolution geospatial data to estimate the end-use electricity demand for a range of agricultural PUE across Sub-Saharan Africa, and we share these data in an open-access mapping tool. Next, we use REopt®, a techno-economic optimization model of energy systems, to determine the cost and system sizing implications of incorporating agricultural PUE into microgrid designs in Kenya and Zambia. We estimate the upper bound of agricultural PUE demand for irrigation, milling, shelling, refrigeration, and egg incubation across Sub-Saharan to be 16.8 TWh/yr. We find that incorporating local agricultural PUE into microgrid system designs increases the required system sizing while having minimal impact on the levelized cost of energy of these systems. Our analysis is the first to demonstrate the PUE potential in the agricultural sector at a 10x10-kilometer resolution across Sub-Saharan Africa and to show, at scale, how site-specific PUE can impact the cost and sizing of microgrids that are otherwise deployed to serve local household and community load.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRSDate not supplied

Runway Configuration Management with Offline Reinforcement Learning

Runway configuration management (RCM) is a challenging task, and it affects the efficiency of the National Airspace System (NAS) and airport surface operations significantly. Each airport, depending on the geometry, capacity, local climate patterns, etc. has multiple configurations for the runway usage for arriving and departing flights. Many factors such as the incoming/outgoing traffic load, wind direction and speed, convective weather, cloud ceiling and other environmental factors might affect the choice of a runway configuration at any point in time. However, other factors such as safety measures and regulations, noise abatement, capacity of each configuration, and preference of the air traffic controllers (ATCs) can also play a significant role in selecting the configuration. A sub-optimal selection of the runway configuration, or delay in making configuration changes might result in significant increase in taxi times for aircraft on the surface of the airport, fuel and energy use of the aircraft, and maintenance costs. It can also lead to safety concerns, such as an aircraft performing one or more go-arounds before being able to land. All these factors make RCM an extremely important and challenging decision-making process for the ATCs. The current state of practice sets the runway configuration by the ATCs based on relevant information available at the time including weather, traffic, noise abatement, safety bounds, etc. This makes the decision-making process subjective based on the accuracy of the available information and the bias in human decision making. Unfortunately, this approach yields poor results (e.g., significant delays) if the predicted outcomes are uncertain and their relative impact is not well understood. This is especially evident when the uncertainty increases the size of possible predicted outcomes (combinatorial explosion in possible scenarios) that cannot be handled by human reasoning. On the other hand, an automated approach based on machine intelligence can make use of historical data and search through all (or significant amount of) possible scenarios under uncertainty and make well-informed decisions.

Milad Memarzadeh↗
DOE OSTI2026

Overview of Wendelstein 7-X high-performance operation

The Wendelstein 7-X (W7-X) stellarator has completed two consecutive experimental campaigns OP 2.2 (Sep.-Dec. 2024) and OP 2.3 (Feb.-May 2025) under a new operational strategy enabling more than one year of uninterrupted device availability. This approach, supported by exceptionally high subsystem reliability, allowed sustained high-efficiency plasma operations with up to 80–100 discharges per day across a broad range of magnetic configurations. Several key technical upgrades-most notably the first operation of a 1.5 MW class steady-state gyrotron, a new steady-state pellet injector, and advanced real-time feedback control systems significantly enhanced heating, fueling, and plasma control capabilities. Together, these improvements enabled major advances in long-pulse performance, high-β operation, and confinement optimization. Long-pulse discharges achieved 1.8 GJ of injected energy under fully detached divertor conditions, while reduced-field scenarios facilitated record volume-averaged β values approaching 3%. High-performance plasmas with centrally peaked density profiles, created via neutral beam injection (NBI) or sustained pellet fueling, demonstrated strongly reduced turbulent transport and stellarator-record fusion triple products. Complementary studies of power exhaust and divertor heat loads revealed the role of scrape-off-layer drift physics in shaping strike-line patterns under attached conditions. Together, the results from OP 2.2 and OP 2.3 significantly expand the operational space of W7-X and strengthen its role as a leading platform for steady-state stellarator research and reactor-relevant plasma scenarios.

Wendelstein 7-X↗
NASA NTRS2017

Optimal Control Surface Layout for an Aeroservoelastic Wingbox

This paper demonstrates a technique for locating the optimal control surface layout of an aeroservoelastic Common Research Model wingbox, in the context of maneuver load alleviation and active utter suppression. The combinatorial actuator layout design is solved using ideas borrowed from topology optimization, where the effectiveness of a given control surface is tied to a layout design variable, which varies from zero (the actuator is removed) to one (the actuator is retained). These layout design variables are optimized concurrently with a large number of structural wingbox sizing variables and control surface actuation variables, in order to minimize the sum of structural weight and actuator weight. Results are presented that demonstrate interdependencies between structural sizing patterns and optimal control surface layouts, for both static and dynamic aeroelastic physics.

Stanford, Bret K.↗
NASA NTRS2003

Reliability Sensitivity Analysis and Design Optimization of Composite Structures Based on Response Surface Methodology

This report discusses the development and application of two alternative strategies in the form of global and sequential local response surface (RS) techniques for the solution of reliability-based optimization (RBO) problems. The problem of a thin-walled composite circular cylinder under axial buckling instability is used as a demonstrative example. In this case, the global technique uses a single second-order RS model to estimate the axial buckling load over the entire feasible design space (FDS) whereas the local technique uses multiple first-order RS models with each applied to a small subregion of FDS. Alternative methods for the calculation of unknown coefficients in each RS model are explored prior to the solution of the optimization problem. The example RBO problem is formulated as a function of 23 uncorrelated random variables that include material properties, thickness and orientation angle of each ply, cylinder diameter and length, as well as the applied load. The mean values of the 8 ply thicknesses are treated as independent design variables. While the coefficients of variation of all random variables are held fixed, the standard deviations of ply thicknesses can vary during the optimization process as a result of changes in the design variables. The structural reliability analysis is based on the first-order reliability method with reliability index treated as the design constraint. In addition to the probabilistic sensitivity analysis of reliability index, the results of the RBO problem are presented for different combinations of cylinder length and diameter and laminate ply patterns. The two strategies are found to produce similar results in terms of accuracy with the sequential local RS technique having a considerably better computational efficiency.

Rais-Rohani, Masoud↗
DOE OSTI2025

Understanding the structural mechanics of ligated DNA crystals via molecular dynamics simulation

DNA self-assembly is a highly programmable method to construct arbitrary architectures based on sequence complementarity. Among various constructs, DNA crystals are macroscopic crystalline materials formed by assembling motifs via sticky end association. Due to their high structural integrity and size ranging from tens to hundreds of micrometers, DNA crystals offer unique opportunities to study the structural properties and deformation behaviors of DNA assemblies. For example, enzymatic ligation of sticky ends can selectively seal nicks resulting in more robust structures with enhanced mechanical properties. However, the research efforts have been mostly on experiments involving different motif designs, structural optimization, or new synthesis methods, while their mechanics are not yet fully understood. The complex properties of DNA crystals are difficult to study via experiments alone, and numerical simulation can complement and aid the experiments. The coarse-grained molecular dynamics (MD) simulation is a powerful tool that can probe the mechanics of DNA assemblies. Here, we investigate DNA crystals made of four different motif lengths with various ligation patterns (full ligation, major directions, connectors, and in-plane) using oxDNA, an open-source, coarse-grained MD platform. We found that several distinct deformation stages emerge in response to mechanical loading and that the number and the location of ligated nucleotides can significantly modulate structural behaviors. These findings should be useful for predicting crystal properties and thus improving the design.

DNA crystal↗
DOE OSTI2021

Addressing Load Imbalance in Bioinformatics and Biomedical Applications: Efficient Scheduling across Multiple GPUs

Computational bioinformatics and biomedical applications frequently contain heterogeneously sized units of work or tasks, for instance due to variability in the sizes of biological sequences and molecules. Variable-sized workloads lead to load imbalances in parallel implementations which detract from efficiency and performance. Many modern computing resources now have multiple graphics processing units(GPUs) per computer for acceleration. These multiple GPU resources need to be used efficiently through balancing of workloads across the GPUs. OpenMP is a portable directive-based parallel programming API used ubiquitously in bioscience applications to program CPUs; recently, the use of OpenMP directives for GPU acceleration has become possible. Here, motivated by experiences with imbalanced loads in GPU-accelerated bioinformatics applications, we address the load balancing problem using OpenMP task-to-GPU scheduling combined with OpenMP GPU offloading for multiply heterogeneous workloads – loads with both variable input sizes, and simultaneously, variable convergence rates for algorithms with a stochastic component – scheduled across multiple GPUs. We aim to develop strategies which are both easy to use and have lower overheads, and may be incorporated incrementally in existing programs which already make use of OpenMP for CPU-based threading in order to make use of multi-GPU computers. We test different combinations of input size variability and convergence rate variability, and characterize the effects of these different scenarios on the performance of scheduling strategies across multiple GPUs with OpenMP. We present several dynamic scheduling solutions for different parallel patterns, explore optimizations, and provide publicly available example computational kernels to make these strategies easy to use in programs. This work will enable application developers to efficiently and easily use multiple GPUs for imbalanced workloads found in bioinformatics and biomedical applications.

Thavappiragasam, Mathialakan↗
DOE OSTI2022

Feasibility of lead-bismuth-cooled fast reactor with fuel shuffling strategy

This paper investigates the feasibility of fuel shuffling of the Advanced Nitride-fueled Simplified 100 MWe (ANTS-100e) Lead-bismuth-cooled Fast Reactor (LFR). The preliminary design has the active core of 235 cm in diameter and 100 cm high. The core is composed of 138 hexagonal fuel assemblies enclosed with a 15-15Ti steel duct. The core adopts the onion zoning with the two zones of low-enriched uranium for a better breeding capability to maximize the first cycle length to 10 effective full-power years (EFPYs). The fuel shuffling study is performed by using the UNIST in-house Monte Carlo code MCS. The neutronic parameters such as the multiplication factor and burnup are being investigated. The fuel shuffling strategy is evaluated through several loading patterns for flattening radial power distribution. It is also found that the core with fuel shuffling could lengthen the core lifetime by more than 5 EFPYs without introducing fresh fuel. Thus, it appears that the ANTS-100e can be principally designed to have the second cycle with another 10 EFPYs if introducing fresh fuel into the core. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Analyzing the Validity of Brazilian Testing Using Digital Image Correlation and Numerical Simulation Techniques

Characterizing the mechanical behavior of rocks plays a crucial role to optimize the fracturing process in unconventional reservoirs. However, due to the intrinsic anisotropy and heterogeneity in unconventional resources, fracture process prediction remains the most significant challenge for sustainable and economic hydrocarbon production. During the deformation tracking under compression, deploying conventional methods (strain gauge, extensometer, etc.) is insufficient to measure the deformation since the physical attachment of the device is restricted to the size of the sample, monitoring limited point-wise deformation, producing difficulties in data retrieval, and a tendency to lose track in failure points, etc. Where conventional methods are limited, the application of digital image correlation (DIC) provides detailed and additional information of strain evolution and fracture patterns under loading. DIC is an image-based optical method that records an object with a camera and monitors the random contrast speckle pattern painted on the facing surface of the specimen. To overcome the existing limitations, this paper presents numerical modeling of Brazilian disc tests under quasi-static conditions to understand the full-field deformation behaviors and finally, it is validated by DIC. As the direct tensile test has limitations in sample preparation and test execution, the Brazilian testing principle is commonly used to evaluate indirectly the tensile strength of rocks. The two-dimensional numerical model was built to predict the stress distribution and full-field deformation on Brazilian disc under compression based on the assumptions of a homogenous, isotropic and linear elastic material. The uniaxial compression test was conducted using the DIC technique to determine the elastic properties of Spider Berea sandstone, which were used as inputs for the simulation model. The model was verified by the analytical solution and compared with the digital image correlation. The numerical simulation results showed that the solutions matched reasonably with the analytical solutions where the maximum deviation of stress distribution was obtained as 14.59%. The strain evolution (normal and shear strains) and displacements along the central horizontal and vertical planes were investigated in three distinguishable percentages of peak loads (20%, 40%, and 90%) to understand the deformation behaviors in rock. The simulation results demonstrated that the strain evolution contours consistently matched with DIC generated contours with a reasonable agreement. The changes in displacement along the central horizontal and vertical planes showed that numerical simulation and DIC generated experimental results were repeatable and matched closely. In terms of validation, Brazilian testing to measure the indirect tensile strength of rocks is still an issue of debate. The numerical model of fracture propagation supported by digital image correlation from this study can be used to explain the fracturing process in the homogeneous material and can be extended to non-homogeneous cases by incorporating heterogeneity, which is essential for rock mechanics field applications.

58 GEOSCIENCES↗
DOE OSTI2018

Cohesive finite element modeling of the delamination of HTPB binder and HMX crystals under tensile loading

Accurately modeling the mechanical behavior of the polymer binders and the degradation of interfaces between binder and crystal is important to science-based understanding of the macro-scale response of polymer bonded explosives. The paper presents a description of relatively a simple bi-crystal HMX-HTPB specimen and associated tensile loading experiment including computed tomography imaging, the pertinent constitutive theory, and details of numerical simulations used to infer the behavior of the material during the delamination process. Within this work, mechanical testing and direct numerical simulation of this relatively simple bi-crystal system enabled reasonable isolation of binder-crystal interface delamination, in which the effects of the complicated thermomechanical response of explosive crystals were minimized. Cohesive finite element modeling of the degradation and delamination of the interface between a modified HTPB binder and HMX crystals was used to reproduce observed results from tensile loading experiments on bi-crystal specimens. Several comparisons are made with experimental measurements in order to identify appropriate constitutive behavior of the binder and appropriate parameters for the cohesive traction-separation behavior of the crystal-binder interface. This research demonstrates the utility of directly modeling the delamination between binder and crystal within crystal-binder-crystal tensile specimen towards characterizing the behavior of these interfaces in a manner amenable to larger scale simulation of polycrystalline PBX materials. One critical aspect of this approach is micro computed tomography imaging conducted during the experiments, which enabled comparison of delamination patterns between the direct numerical simulation and actual specimen. In addition to optimizing the cohesive interface parameters, one important finding from this investigation is that understanding and representing the strain-hardening plasticity of HTPB binder is important within the context of using a cohesive traction-separation model for the delamination of a crystal-binder system.

42 ENGINEERING↗