SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2022

On the operational characteristics and economic value of pumped thermal energy storage

Pumped thermal energy storage (PTES) systems use an electrically-driven heat pump to store electricity in the form of thermal energy, and subsequently dispatch the stored thermal energy to generate electricity using a thermodynamic heat engine. Optimal day-ahead operational scheduling and annual value of a PTES system based on Joule-Brayton thermodynamic cycles and two-tank molten salt hot thermal storage is evaluated in this work. Production cost models, which simultaneously optimize commitment and dispatch schedules for an entire set of generators to minimize the cost of satisfying electricity demand, are employed to determine system-optimal operation and day-ahead energy value of the PTES system within each of six hypothetical near-future grid scenarios intended to approximately represent the U.S. Western Interconnection or the Texas Interconnection. Sensitivity to grid scenario (including the contribution of variable renewable energy sources), thermal storage capacity, relative heat pump and heat engine capacities, and startup/shutdown cycling costs are evaluated. PTES energy value and heat engine annual capacity factor increase strongly as the contribution of variable renewable resources increases, heat pump capacity increases relative to heat engine capacity, or PTES cycling costs decrease. Grid scenarios in which the contribution of variable renewable energy is dominated by solar photovoltaics (PV) vs. wind produce inherently different PTES operational patterns. Annual PTES energy value within PV-dominated scenarios increased with storage capacity only up to approximately seven hours of full-load discharge capacity, whereas that within wind-dominated scenarios exhibited a continual increase with storage duration up to at least 16 hours.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Coordinated PEV Charging for Distribution System Management

Plug-in electric vehicles (PEVs) can help to reduce worldwide dependence on petroleum and carbon emissions. In the past a few years, PEVs have received considerable attention as an eco-friendly and cost-effective alternative to conventional gasoline vehicles. PEVs consume higher power from the grid during charging compared to conventional residential loads. Therefore, the emerging fleet of PEVs will introduce a considerable amount of additional load on power system. On the other hand, majority of PEVs are parked for more than 90% of the time, making them ideal for providing various services through Electric Vehicle Supply Equipment (EVSE), a.k.a. charging stations which are PEVs’ connection points to the power system. In this project, coordinated PEV charging are developed to provide distribution system services, considering PEV characteristics, realistic travel pattern, charging behaviors, and EVSE power rating and availability. The main efforts and contributions are summarized as follows. • Mobility models and EVSE charging capability/availability are indispensable for evaluating PEV charging coordination strategies and developing PEV charging control. Many existing studies utilize simplified mobility models assuming that the entire vehicle fleet returns home in the evening and is parked at home until the next morning. Some other PEV models better represent diversified home arrival and departure time, but cannot capture varying charging flexibility and capability at different locations, and therefore are not appropriate to study the impacts of public EVSE on PEV load and utilization. In this work, a mobility and charging flexibility model is proposed to better represent the temporal availability and varying charging capability from PEV onboard charger and EVSE. • Based on the mobility model, a generalized optimization method is proposed to evaluate different PEV charging coordination strategies. In existing literature, different algorithms and methods need to be designed to evaluate each charging strategy from both vehicle owners’ and system’s perspective. With the optimization method proposed in this work, different charging control strategies can be studied and compared by only updating objective functions. Moreover, optimization tricks are provided to convert the optimization problems to equivalent linear programming problems, which can be efficiently solved with existing solvers. • An innovative scheduling and control framework is proposed to enable smart PEV charging for grid services while meeting PEV owners’ travel needs. In the proposed framework, each set of EVSE is equipped with a controller that estimates charging power and energy flexibility based on vehicle characteristics, EVSE power rating, battery energy state, and upcoming trip information. With the simplified charging flexibility model received from each EVSE controller, the central coordinator determines the optimal power allocation for a look-ahead time window for the given grid services. The proposed charging control can help reduce the computational complexity and communication requirement compared with existing methods. It is also scalable to the expanding PEV fleet and robust to uncertainties in upcoming vehicle trips and future system condition. • The proposed charging strategy evaluation and smart charging control methods are applied to one of the prototypical feeders developed by the Pacific Northwest National Laboratory (PNNL) for case studies. Detailed trip information extracted from the National Household Travel Survey (NHTS) is used to represent travel and parking patterns. Different scenarios are simulated to evaluate the performance of the proposed methods and understand how PEV mobility and public EVSE availability affect PEV load and potential grid services.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2024

Eco-driving Profile Optimization by Dynamic Programming for Battery Electric Vehicles

Although full automation has not yet been achieved, automated vehicles are a valid research area. Not only would automated vehicles provide ultimate driver convenience, but they would maximize energy efficiency by eliminating undesired human driving behaviors and optimally controlling the powertrain. From the perspective of control related to energy saving, speed profile optimization is important for improving system efficiency and satisfying passenger demands. This study employs Dynamic Programming (DP) to solve the constrained optimal problem for travel time, distance, and speed limit by exploring all possible control options. The solutions obtained by DP demonstrate consistent control patterns combining four control modes-acceleration, cruising, coasting, and braking, with cruising or coasting being selective depending on the boundary conditions. Further, this study introduces DP-based simulation results and attempts to provide comprehensive interpretations of the optimal policy by analyzing the essential factors that affect the control problem, including boundary conditions, road load, and powertrain characteristics. Based on these interpretations, the control concepts can be explained as the optimal policy selecting the best control option based on system efficiency and boundary conditions. The results of DP are compared with a human-like driver model to show that the optimal speed profiles can effectively reduce energy consumption.

Autonomous vehicles↗

DOE OSTI2023

Designing Low Tortuosity Electrodes through Pattern Optimization for Fast-Charging

The development of fast-charging technologies is crucial for expediting the progress and promotion of electric vehicles. In addition to innovative material exploration, reduction in the tortuosity of electrodes is a favored strategy to enhance the fast-charging capability of lithium-ion batteries by optimizing the ion-transfer kinetics. To realize the industrialization of low-tortuosity electrodes, a facile, cost-effective, highly controlled, and high-output continuous additive manufacturing roll-to-roll screen printing technology is proposed to render customized vertical channels within electrodes. Extremely precise vertical channels are fabricated by applying the as-developed inks, using LiNi 0.6 Mn 0.2 Co 0.2 O 2 as the cathode material. Additionally, the relationship between the electrochemical properties and architecture of the channels, including the pattern, channel diameter, and edge distance between channels, is revealed. The optimized screen-printed electrode exhibited a seven-fold higher charge capacity (72 mAh g -1 ) at a current rate of 6 C and superior stability compared with that of the conventional bar-coated electrode (10 mAh g -1 , 6 C) at a mass loading of 10 mg cm -2 . Finally, this roll-to-roll additive manufacturing can potentially be applied to various active materials printing to reduce electrode tortuosity and enable fast charging in battery manufacturing.

42 ENGINEERING↗

DOE OSTI2018

Optimal feed flow sequence for multi-effect distillation system integrated with supercritical carbon dioxide Brayton cycle for seawater desalination

Multi-effect distillation is often integrated with a Rankine power cycle for cogeneration (simultaneous production of power and desalinated water). Such integration increases the condenser operating pressure for the Rankine cycle (and increases the heat-rejection temperature) to produce desalinated water, resulting in decreased power-plant efficiency. The supercritical carbon dioxide Brayton cycle has a higher efficiency compared to the Rankine cycle. The waste heat rejected from supercritical carbon dioxide Brayton cycle is at sufficiently hot temperature to have feasible energy integration with multi-effect distillation system. The paper introduces the novel concept of cogeneration without being a parasitic load to the power cycle. For the illustrative example considered, integrating 4-effect distillation system with a 115 MWe power plant can produce 3041 m3 of distillate per day at 1.06 $/m3, at a constant power plant efficiency of 49.2%. The pattern in which feed enters the desalination system often dictates its energy consumption. The two most commonly used feed configurations for multi-effect distillation system are parallel/cross feed and forward feed. With steam as a heat source (i.e., a latent heat source), parallel/cross feed is the most energy-efficient feed configuration. The other objective of this paper is to identify the optimal feed flow configuration for multi-effect distillation system integrated with a supercritical carbon dioxide power cycle (i.e., a sensible heat source). For the simplified network, forward feed is the best feed configuration, which yields a 7.5% increase in distillate production at 2.6% reduced distillate cost. Additionally, different methods for reducing the brine discharge are studied, which can help to achieve zero liquid discharge. Result show increasing the maximum brine concentration gives superior results compared to brine recycling. The system modelling is done using the principle of process integration, and an analytical methodology for cogeneration is derived.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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↗

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↗

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↗

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↗

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↗