Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “benchmark testing”

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 919 records · Page 51

Research on Field Emission and Dark Current in ILC Cavities

Field emission and dark current are issues of concern for SRF cavity performance and SRF linac operation. Complete understanding and reliable control of the issue are still needed, especially in full-scale multi-cell cavities. Our work aims at developing a generic procedure for finding an active field emitter in a multi-cell cavity and benchmarking the procedure through cavity vertical test. Our ultimate goal is to provide feedback to cavity preparation and cavity string assembly in order to reduce or eliminate filed emission in SRF cavities. Systematic analysis of behaviors of field emitted electrons is obtained by ACE3P developed by SLAC. Experimental benchmark of the procedure was carried out in a 9-cell cavity vertical test at JLab. The energy spectrum of Bremsstrahlung X-rays is measured using a NaI(Tl) detector. The end-point energy in the X-ray energy spectrum is taken as the highest kinetic electron energy to predict longitudinal position of the active field emitter. Angular location of the field emitter is determined by an array of silicon diodes around irises of the cavity. High-resolution optical inspection was conducted at the predicted field emitter location.

Liu, Kexin↗

Preliminary Two-Phase Terry Turbine Nozzle Models for RCIC Off-Design Operation Conditions

This report presents the effort to extend the single-phase analytical Terry turbine model to cover two-phase off-design conditions. The work includes: (1) adding well-established two-phase choking models – the Isentropic Homogenous Equilibrium Model (IHEM) and Moody’s model, and (2) theoretical development and implementation of a two-phase nozzle expansion model. The two choking models provide bounding cases for the two-phase choking mass flow rate. The new two-phase Terry turbine model uses the choking models to calculate the mass flow rate, the critical pressure at the nozzle throat, and steam quality. In the divergent stage, we only consider the vapor phase with a similar model for the single-phase case by assuming that the liquid phase would slip along the wall with a much slower speed and will not contribute the impulse on the rotor. We also modify the stagnation conditions according to two-phase choking conditions at the throat and the cross-section areas for steam flow at the nozzle throat and at the nozzle exit. The new two-phase Terry turbine model was benchmarked with the same steam nozzle test as for the single-phase model. Better agreement with the experimental data is observed than from the single-phase model. We also repeated the Terry turbine nozzle benchmark work against the Sandia CFD simulation results with the two-phase model for the pure steam inlet nozzle case. The RCIC start-up tests were simulated and compared with the single-phase model. Similar results are obtained. Finally, we designed a new RCIC system test case to simulate the self-regulated Terry turbine behavior observed in Fukushima accidents. In this test, a period inlet condition for the steam quality varying from 1 to 0 is applied. For the high quality inlet period, the RCIC system behaves just like the normal operation condition with a high pump injection flow rate and a nominal steam release rate through the turbine, with the net addition of water to the primary system; for the low quality inlet period, the RCIC turbine shaft work dramatically decreases and results in a much reduced pump injection flow rate, and the mixture flow rate through the turbine increases due to the high liquid phase flow rate. The net effect for this period is net removal of coolant from the primary loop. With the periodic addition and removal of coolant to the primary loop, the self-regulation mode of the RCIC system can be maintained for a quite long time. Both the IHEM and Moody’s models generate similar phenomena; however noticeable differences can be observed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Carbon Back Sputter Modeling for Hall Thruster Testing

Lifetime requirements for electric propulsion devices, including Hall Effect thrusters, are continually increasing, driven in part by NASA's inclusion of this technology in it's exploration architecture. NASA will demonstrate high-power electric propulsion system on the Solar Electric Propulsion Technology Demonstration Mission (SEP TDM). The Asteroid Redirect Robotic mission is one candidate SEP TDM, which is projected to require tens of thousands of thruster life. As thruster life is increased, for example through the use of improved magnetic field designs, the relative influence of facility effects increases. One such effect is the sputtering and redeposition, or back sputter, of facility materials by the high energy thruster plumes. In support of wear testing for the Hall Effect Rocket with Magnetic Shielding (HERMeS) project, the back sputter from a Hall effect thruster plume has been modeled for the NASA Glenn Research Center's Vacuum Facility 5. The predicted wear at a near-worst case condition of 600 V, 12.5 kW was found to be on the order of 1 micron/kh in a fully carbon-lined chamber. A more detailed numerical Monte Carlo code was also modified to estimate back sputter for a detailed facility and pumping configuration. This code demonstrated similar back sputter rate distributions, but is not yet accurately modeling the magnitudes. The modeling has been benchmarked to recent HERMeS wear testing, using multiple microbalance measurements. These recent measurements have yielded values on the order of 1.5 - 2 micron/kh at 600 V and 12.5 kW.

Propulsion↗

Carbon Back Sputter Modeling for Hall Thruster Testing

In support of wear testing for the Hall Effect Rocket with Magnetic Shielding (HERMeS) program, the back sputter from a Hall effect thruster plume has been modeled for the NASA Glenn Research Centers Vacuum Facility 5. The predicted wear at a near-worst case condition of 600 V, 12.5 kW was found to be on the order of 3 4 mkhour in a fully carbon-lined chamber. A more detailed numerical monte carlo code was also modified to estimate back sputter for a detailed facility and pumping configuration. This code demonstrated similar back sputter rate distributions, but is not yet accurately modeling the magnitudes. The modeling has been benchmarked to recent HERMeS wear testing, using multiple microbalance measurements. These recent measurements have yielded values, on the order of 1.5- 2 microns/khour.

Facility↗

CNRS Molten Salt Reactor Benchmark Analysis Using Griffin-Pronghorn Coupled Multi-physics Code System

Molten salt reactors (MSRs) with flowing fuel have the unique feature of utilizing the fuel salt for heat generation and extraction at the same time, since the fuel salt is circulating through the whole primary loop. This movement of fuel salt results in a partial decay of the delayed neutron precursors (DNPs) outside the core and the corresponding redistribution of the DNPs in the active region of the core. To capture this effect, neutronics and thermal-hydraulics (T/H) codes need to be able to handle the movement of DNPs and their decay [1, 2]. Idaho National Laboratory (INL) is actively working on developing the neutronics and T/H codes to model the MSRs with flowing fuel. As one of the efforts, the neutronics code Griffin [3] and the T/H code Pronghorn [4] that are built upon the MOOSE framework [5] have been extended to handle the flowing fuel involving the drift of DNPs, and this capability is being tested and verified with MSR benchmark problems including the CNRS benchmark problem [6]. Additionally, a multi-physics analysis of the molten salt reactor experiment (MSRE) is being conducted to verify and validate the codes against available experimental data [7]. This work presents the analysis of the CNRS benchmark problem using the Griffin-Pronghorn coupled multi-physics code system of INL. Both steady-state and transient problems are analyzed, and the results are compared with those of other participants presented in [6]. The following section provides a description of the benchmark, along with input parameters and observables for each step. Then, the results for each step are presented in comparison with other participants’ results, followed by summary and conclusion.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bringing HPE Slingshot 11 support to Open MPI

The Cray HPE Slingshot 11 network is used on the new exascale systems arriving at the U.S. Department of Energy (DoE) laboratories (e.g., Frontier, Aurora, Perlmutter). As such, the support of this network is an important capability to meet the needs of exascale applications. Here, this article highlights recent work to develop supporting infrastructure to enable Open MPI to efficiently support these new platforms. A key component of this effort involves development of a new Open Fabrics Interface (OFI) provider, LinkX. We discuss the design and development of enhancements that take advantage of the new Slingshot 11 network and AMD GPUs. We include performance data from tests on the Frontier supercomputer using synthetic communication benchmarks, and the vendor provided MPI as a baseline for comparison. The tests demonstrate full functionality of Open MPI on the system and initial results show favorable performance when compared to the highly tuned vendor implementation.

97 MATHEMATICS AND COMPUTING↗

Benchmark specifications and data requirements for EBR-II shutdown heat removal tests SHRT-17 and SHRT-45R

Argonne National Laboratory’s (ANL) Experimental Breeder Reactor II (EBR-II) was a liquid metal reactor (LMR) with a sodium-bonded metallic fuel core that contributed very favorably to the reactor’s negative reactivity feedback. On June 20, 1984, a severe loss-of-flow test in the Shutdown Heat Removal Test (SHRT) series demonstrated the effectiveness of natural circulation in the EBR-II reactor. This test was SHRT-17 and at the beginning of the test the primary pumps were tripped at the same time as a full control rod insertion.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Matrix-Free High-Performance Saddle-Point Solvers for High-Order Problems in \(\boldsymbol{H}(\operatorname{\textbf{div}})\)

Here, this work describes the development of matrix-free GPU-accelerated solvers for high-order finite element problems in H(div). The solvers are applicable to grad-div and Darcy problems in saddle-point formulation, and have applications in radiation diffusion and porous media flow problems, among others. Using the interpolation–histopolation basis, efficient matrix-free preconditioners can be constructed for the (1, 1)-block and Schur complement of the block system. With these approximations, block-preconditioned MINRES converges in a number of iterations that is independent of the mesh size and polynomial degree. The approximate Schur complement takes the form of an M-matrix graph Laplacian and therefore can be well-preconditioned by highly scalable algebraic multigrid methods. High-performance GPU-accelerated algorithms for all components of the solution algorithm are developed, discussed, and benchmarked. Numerical results are presented on a number of challenging test cases, including the “crooked pipe” grad-div problem, the SPE10 reservoir modeling benchmark problem, and a nonlinear radiation diffusion test case.

97 MATHEMATICS AND COMPUTING↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Benchmark Accelerated Aging of Harvested Hypalon/Epr and Cspe/Xlpe Power and I&C Cable in Nuclear Power Plants

As part of the Light Water Reactor and Sustainability (LWRS) program in the U.S. Department of Energy (DOE) Office of Nuclear Energy, material aging and degradation research is currently geared to support the long-term operation of existing nuclear power plants (NPPs) as they move beyond their initial 40 year licenses. The goal of this research is to provide information so that NPPs can develop aging management programs (AMPs) to address replacement and monitoring needs as they look to operate for 20 years, and in some cases 40 years, beyond their initial operating lifetimes. For cable insulation and jacket materials that support instrument, control, and safety systems, accelerated aging data are needed to determine priorities in cable aging management programs. Before accelerated thermal and radiation aging of harvested, representative cable insulation and jacket materials, the benchmark performance of a new test capability at Oak Ridge National Laboratory (ORNL) was evaluated for temperatures between 70 and 135°C, dose rates between 100 and 500 Gy/h, and accumulated doses up to 20 kGy, Samples that were characterized and are representative of current materials in use were harvested from the Callaway NPP near Fulton, Missouri, and the San Onofre NPP north of San Diego, California. From the Callaway NPP, a multiconductor control rod cable manufactured by Boston Insulated Wire (BIW), with a Hypalon/ chorolosulfonated polyethylene (CSPE) jacket and ethylene-propylene rubber (EPR) insulation, was harvested from the auxiliary space during a planned outage in 2013. This cable was placed into service when the plant was started in 1984. From the San Onofre NPP, a Rockbestos Firewall III (FRIII) cable with a Hypalon/ CSPE jacket with cross-linked polyethylene (XLPE) insulation was harvested from an on-site, climate-controlled storage area. This conductor, which was never placed into service, was procured around 2007 in anticipation of future operation that did not occur. Benchmark aging for both jacket and insulation material was carried out in air at a temperature of 125°C or in uniform 140 Gy/h gamma irradiation over a period of 60 days. Their mechanical properties over the course of their exposures were compared with reference data from comparable cable jacket/insulation compositions and aging conditions. For both accelerated thermal and radiation aging, it was observed that the mechanical properties for the Callaway BIW control rod cable were consistent with those previously measured. However, for the San Onofre Rockbestos FRIII, there was an observable functional difference for accelerated thermal aging at 125°C. Details on possible sources for this difference and plans for resolving each source are given in this paper.

aging↗

BENCHMARK ACCELERATED AGING OF HARVESTED HYPALON/EPR AND CSPE/XLPE POWER AND I&C CABLE IN NUCLEAR POWER PLANTS

As part of the Light Water Reactor and Sustainability (LWRS) program in the U.S. Department of Energy (DOE) Office of Nuclear Energy, material aging and degradation research is currently geared to support the long-term operation of existing nuclear power plants (NPPs) as they move beyond their initial 40 year licenses. The goal of this research is to provide information so that NPPs can develop aging management programs (AMPs) to address replacement and monitoring needs as they look to operate for 20 years, and in some cases 40 years, beyond their initial operating lifetimes. For cable insulation and jacket materials that support instrument, control, and safety systems, accelerated aging data are needed to determine priorities in cable aging management programs. Before accelerated thermal and radiation aging of harvested, representative cable insulation and jacket materials, the benchmark performance of a new test capability at Oak Ridge National Laboratory (ORNL) was evaluated for temperatures between 70 and 135 C, dose rates between 100 and 500 Gy/h, and accumulated doses up to 20 kGy, Samples that were characterized and are representative of current materials in use were harvested from the Callaway NPP near Fulton, Missouri, and the San Onofre NPP north of San Diego, California. From the Callaway NPP, a multiconductor control rod cable manufactured by Boston Insulated Wire (BIW), with a Hypalon/ chorolosulfonated polyethylene (CSPE) jacket and ethylene-propylene rubber (EPR) insulation, was harvested from the auxiliary space during a planned outage in 2013. This cable was placed into service when the plant was started in 1984. From the San Onofre NPP, a Rockbestos Firewall III (FRIII) cable with a Hypalon/ CSPE jacket with cross-linked polyethylene (XLPE) insulation was harvested from an on-site, climate-controlled storage area. This conductor, which was never placed into service, was procured around 2007 in anticipation of future operation that did not occur. Benchmark aging for both jacket and insulation material was carried out in air at a temperature of 125 C or in uniform 140 Gy/h gamma irradiation over a period of 60 days. Their mechanical properties over the course of their exposures were compared with reference data from comparable cable jacket/insulation compositions and aging conditions. For both accelerated thermal and radiation aging, it was observed that the mechanical properties for the Callaway BIW control rod cable were consistent with those previously measured. However, for the San Onofre Rockbestos FRIII, there was an observable functional difference for accelerated thermal aging at 125 C. Details on possible sources for this difference and plans for resolving each source are given in this paper.

Duckworth, Robert C↗

Facile and Scalable Synthesis of Metal- and Nitrogen-Doped Carbon Nanotubes for Efficient Electrochemical CO 2 Reduction

Metal- and nitrogen-doped carbon (M–N–C) is a promising material to catalyze electrochemical CO 2 reduction reaction (CO 2 RR). However, most M–N–C catalysts in the literature require complicated synthesis procedures and produce small quantities per batch, limiting the commercialization potential. In this work, we developed a simple and scalable synthesis method to convert metal-impurity-containing commercial carbon nanotubes (CNTs) and nitrogen-containing organic precursors into M–N–C via one-step moderate-temperature (650 °C) pyrolysis without any other treatment nor the need to add metal precursors. Batches of catalysts in varied mass up to 10 g (150 mL in volume) per batch were synthesized, and repeatable catalytic performances were demonstrated. To the best of our knowledge, the 10 g batch is one of the largest batches of CO2RR catalysts synthesized in the literature while requiring minimal synthesis steps. The catalyst possessed single-atomic iron–nitrogen (Fe–N) sites, enabling a high performance of >95% CO product selectivity at a high current density of 400 mA/cm 2 and high stability for 45 h at 100 mA/cm 2 in a flow cell testing. The catalyst outperformed a benchmark noble-metal nanoparticle catalyst and achieved longer stability than many other reported M–N–C catalysts in the literature. The scalable and cost-effective synthesis developed in this work paves a pathway toward practical CO 2 RR applications. The direct utilization of metal impurities from raw CNTs for efficient catalyst synthesis with minimal treatment is a green and sustainable engineering approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Using FIPD and OPTD to Benchmark Metallic Fuel Performance

This report serves as an introduction, tutorial, and benchmark specification for out-of-pile tests on metallic fuel. It introduces a new user to the EBR-II legacy fuel performance test program and the fast reactor fuel performance databases built to preserve the records. It then details the information stored in each database and how to find it. A benchmark specification is included for a small set of out-of-pile tests on U-10Zr fuel to function as a tutorial demonstrating how the legacy fuel performance data sets stored in the FIPD and OPTD databases can be used together to benchmark fuel performance models for steady-state and transient performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical studies of electron cyclotron wave current drive on HL-2A tokamak

The electron cyclotron wave (ECW) current drive (CD) for the HL-2A tokamak is investigated numerically with a new ray-tracing and Fokker-Planck code. The code is benchmarked with other well-tested linear and quasilinear codes and is then used to study the electron cyclotron current drive on the HL-2A tokamak. The wave propagation, power deposition, and driven-current profiles are presented. The effect of electron trapping is also assessed. It is found that quasilinear effects are negligible at the present ECW power levels and that when both waves are injected at an angle of 20° on the plasma equatorial plane, the CD efficiency for the HL-2A saturates at ∼0.029 × 10{sup 20 }A/W/m{sup 2} and ∼0.020 × 10{sup 20 }A/W/m{sup 2} for the 0.5 MW/68 GHz first harmonic ordinary (O1) and 1 MW/140 GHz second harmonic extraordinary (X2) modes, respectively. The effects of the plasma density, temperature, and wave-launching position on the driven current are also investigated analytically and numerically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics modelling in PyLith: poroelasticity

SUMMARY PyLith, a community, open-source code for modelling quasi-static and dynamic crustal deformation with an emphasis on earthquake faulting, has recently been updated with a flexible multiphysics implementation. We demonstrate the versatility of the multiphysics implementation by extending the code to model fully coupled continuum poromechanics. We verify the newly incorporated physics using standard benchmarks for a porous medium saturated with a slightly compressible fluid. The benchmarks include the 1-D consolidation problem as outlined by Terzaghi, Mandel’s problem for the 2-D case, and Cryer’s problem for the 3-D case. All three benchmarks have been added to the PyLith continuous integration test suite. We compare the closed form analytical solution for each benchmark against solutions generated by our updated code, and lastly, demonstrate that the poroelastic material formulation may be used alongside the existing fault implementation in PyLith.

Geochemistry & Geophysics↗

PMT-004 Validation Testing for Thermal Scattering Laws and Some Interesting Temperature-Dependent Implications [Slides]

This presentation finds that PMT-004 has four new benchmark cases highly sensitive to PE TSL (2 cases) and PMMA TSL (2 Cases). Additionally discovered was that the PE cases were well predicted using MCNP6.2 and ENDF/B-VIII.0. PMMA cases overpredicted by approximately 0.6-0.7% in k eff likely due to the lack of integral benchmarks sensitive to PMMA. The NCSU ENDF/B-VIII.1 PMMA TSL is an improvement over VIII.0. This was Accepted at last ICSBEP meeting pending comment resolution, should be published in 2024 version of the Handbook. It was found that temperature had a large impact on reactivity of the critical configuration. This has Implications for validation work for thermal cases- need to adjust TSL data to correct temperature as it can have hundreds of pcm effects for a few °C. Future thermal experiments should try and measure reactivity at multiple temperatures to add in data testing and benchmark adjustment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A robust interface finite element formulation for modeling brittle material failure problems

Abstract Failure of many brittle materials and structures can be modeled using interface‐oriented finite elements combined with intrinsic cohesive zone models. The discontinuous Galerkin (DG) finite element method provides an innovative framework for modeling brittle crack propagation with zero‐thickness interface elements, which can accommodate extrinsic cohesive laws to avoid the artificial compliance required in intrinsic cohesive models. However, robust formulations and implementations of DG methods are critical in alleviating the well‐known convergence issues for both crack nucleation and propagation with reduced instability. This paper presents a robust interface element formulation by modifying the incomplete interior penalty Galerkin (IIPG) method, which successfully avoids the initial element interface penetration across elements that occurs prior to crack nucleation, and thereby greatly reduces the instability issue as cracks open. We further verified and validated our implementation by using a bar tension test and a beam fracturing benchmark. The robustness of our proposed interface element method was demonstrated by a micromechanics fiber/matrix debonding problem with 64 fibers embedded in a bulk matrix.

Liu, Ruijie↗