Search NASA⌕ Search

SEARCH · Search NASA

Results for “displacement functions”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Finite-temperature many-body perturbation theory for electrons: Algebraic recursive definitions, second-quantized derivation, linked-diagram theorem, general-order algorithms, and grand canonical and canonical ensembles

A comprehensive and detailed account is presented for the finite-temperature many-body perturbation theory for electrons that expands in power series all thermodynamic functions on an equal footing. Algebraic recursions in the style of the Rayleigh–Schrödinger perturbation theory are derived for the grand potential, chemical potential, internal energy, and entropy in the grand canonical ensemble and for the Helmholtz energy, internal energy, and entropy in the canonical ensemble, leading to their sum-over-states analytical formulas at any arbitrary order. For the grand canonical ensemble, these sum-over-states formulas are systematically transformed to sum-over-orbitals reduced analytical formulas by the quantum-field-theoretical techniques of normal-ordered second quantization and Feynman diagrams extended to finite temperature. It is found that the perturbation corrections to energies entering the recursions have to be treated as a nondiagonal matrix, whose off-diagonal elements are generally nonzero within a subspace spanned by degenerate Slater determinants. They give rise to a unique set of linked diagrams—renormalization diagrams—whose resolvent lines are displaced upward, which are distinct from the well-known anomalous diagrams of which one or more resolvent lines are erased. A linked-diagram theorem is introduced that proves the size-consistency of the finite-temperature many-body perturbation theory at any order. General-order algorithms implementing the recursions establish the convergence of the perturbation series toward the finite-temperature full-configuration-interaction limit unless the series diverges. As a result, the normal-ordered Hamiltonian at finite temperature sheds light on the relationship between the finite-temperature Hartree–Fock and first-order many-body perturbation theories.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Microbial production of advanced biofuels

Concerns over climate change have necessitated a rethinking of our transportation infrastructure. One possible alternative to carbon-polluting fossil fuels is biofuels produced by engineered microorganisms that use a renewable carbon source. Two biofuels, ethanol and biodiesel, have made inroads in displacing petroleum-based fuels, but their uptake has been limited by the amounts that can be used in conventional engines and by their cost. Further, advanced biofuels that mimic petroleum-based fuels are not limited by the amounts that can be used in existing transportation infrastructure but have had limited uptake due to costs. In this Review, we discuss engineering metabolic pathways to produce advanced biofuels, challenges with substrate and product toxicity with regard to host microorganisms and methods to engineer tolerance, and the use of functional genomics and machine learning approaches to produce advanced biofuels and prospects for reducing their costs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Kramers nodal line in the charge density wave state of YTe 3 and the influence of twin domains

Recent studies have focused on the relationship between charge density wave (CDW) collective electronic ground states and nontrivial topological states. YTe 3 , a nonmagnetic quasi-two-dimensional chalcogenide, has been reported to exhibit a CDW state below 334 K. Using angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT), we establish that YTe 3 is a CDW-induced Kramers nodal line (KNL) metal, a recently proposed topological state of matter. Scanning tunneling microscopy and low energy electron diffraction reveal two orthogonal domains, each with a unidirectional CDW and a similar wave vector (𝐪 CDW ). When the influence of twin domains is considered, the effective band structure (EBS) computations that utilize DFT-calculated bands using a noncentrosymmetric structure determined by x-ray crystallography, show excellent agreement with ARPES. The noncentrosymmetry of YTe 3 is established by Raman spectroscopy. The Fermi surface and ARPES intensity plots show weak shadow bands displaced by 𝐪 CDW from the main bands. Furthermore, these are linked to CDW modulation, as the EBS calculation confirms. Bilayer split main and shadow bands suggest the existence of crossings, according to theory and experiment. DFT bands, including spin-orbit coupling, indicate existence of a KNL along the Σ direction from multiple crossings of bands dispersing perpendicular to it. Additionally, doubly degenerate bands are only found along the KNL at all energies, with some bands dispersing through the Fermi level.

Angle-resolved photoemission spectroscopy↗
DOE OSTI2022

Revisiting W–ZrC interfaces: A first principles study

We investigate W–ZrC interfaces using first-principles calculations based on the density functional theory. There have been theoretical and experimental studies exploring W–ZrC interfaces, however, the debate regarding the most stable interface continues to persist. In this study, we systematically simulated various W–ZrC interfaces merging W and ZrC surfaces with different orientations. Subsequently, we evaluated their stabilities and explained the corresponding stabilities in terms of the nature of bonding and charge-transfer processes at the interface. We find ZrC(111)–W(110) is the most stable interface with higher adhesive energy than the other interfaces. The additional stability associated with the ZrC(111)–W(110) results from significant interface reconstruction. Three layers of W and ZrC adjacent to the interface are involved in the charge-transfer process leading to stronger ionic bonds in ZrC(111)–W(110) as compared to the other potential candidate: ZrC(100)–W(100). The C and W atoms are found to be displaced from their symmetric position during the reconstruction process at the interface to facilitate stronger bonds with shorter W–C and W–Zr bonds in ZrC(111)–W(110) as compared to ZrC(100)–W(100). This leads to stronger covalent bonds in ZrC(111)–W(110) than that in ZrC(100)–W(100). Therefore, we conclude that the stronger covalent and ionic forces in ZrC(111)–W(110) than those in ZrC(100)–W(100) are responsible for making ZrC(111)–W(110) to be the most stable interface. In conclusion, this study addresses the long-standing question of the most stable W–ZrC interface and derives a number of implications for other W-transition metal carbide interfaces which are potential candidates for improving the mechanical properties of plasma facing materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

An optical-input Maximum Likelihood Estimation feedback system demonstrated on tokamak horizontal equilibrium control

A readily parallelized Maximum Likelihood Estimation (MLE) algorithm with linear computational complexity is demonstrated in real time using only measurements from an extreme ultraviolet (EUV) diagnostic to control the horizontal position of a tokamak plasma. A set of trial emissivity profiles are parameterized by the control quantity of interest (R m ), and the MLE is identified from the profile which minimizes the signal reconstruction residual. The algorithm depends on an empirically determined likelihood function with exponential form. EUV emission (λ ≈ 15eV-1keV) is captured in a poloidal plane by four 16-channel AXUV diodes mounted at different poloidal angles with radial and angular resolution sufficient to discern plasma equilibrium evolution in HBT-EP. Calculations of the plasma major radius by the system are consistent within diagnostic uncertainty for the majority of the discharge with those of: a weighted average of vertical soft X-ray or EUV chords, magnetic sensors, and an equilibrium reconstruction. The feedback system corrects for a horizontal displacement of the major radius equal to 20% of the plasma minor radius by adjusting the vertical field produced from 40 in-vessel control coils in real time. The MLE calculation is performed on a GPU in a 15 μs cycle, with similar performance in this application to a simple weighted average of vertical chords. Finally, results demonstrate horizontal position control using magnetic actuators and an optical observer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Data-Efficient Dimensionality Reduction and Surrogate Modeling of High-Dimensional Stress Fields

Tensor datatypes representing field variables like stress, displacement, velocity, etc., have increasingly become a common occurrence in data-driven modeling and analysis of simulations. Numerous methods [such as convolutional neural networks (CNNs)] exist to address the meta-modeling of field data from simulations. As the complexity of the simulation increases, so does the cost of acquisition, leading to limited data scenarios. Modeling of tensor datatypes under limited data scenarios remains a hindrance for engineering applications. Here, in this article, we introduce a direct image-to-image modeling framework of convolutional autoencoders enhanced by information bottleneck loss function to tackle the tensor data types with limited data. The information bottleneck method penalizes the nuisance information in the latent space while maximizing relevant information making it robust for limited data scenarios. The entire neural network framework is further combined with robust hyperparameter optimization. We perform numerical studies to compare the predictive performance of the proposed method with a dimensionality reduction-based surrogate modeling framework on a representative linear elastic ellipsoidal void problem with uniaxial loading. The data structure focuses on the low-data regime (fewer than 100 data points) and includes the parameterized geometry of the ellipsoidal void as the input and the predicted stress field as the output. The results of the numerical studies show that the information bottleneck approach yields improved overall accuracy and more precise prediction of the extremes of the stress field. Additionally, an in-depth analysis is carried out to elucidate the information compression behavior of the proposed framework.

artificial intelligence↗
DOE OSTI2024

Search for light long-lived neutral particles from Higgs boson decays via vector-boson-fusion production from pp collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

A search is reported for long-lived dark photons with masses between 0.1 GeV and 15 GeV, from exotic decays of Higgs bosons produced via vector-boson-fusion. Events that contain displaced collimated Standard Model fermions reconstructed in the calorimeter or muon spectrometer are probed. This search uses the full LHC Run 2 (2015–2018) data sample collected in proton–proton collisions at $\sqrt{s}=13$ TeV, corresponding to an integrated luminosity of 139 fb –1 . Dominant backgrounds from Standard Model processes and non-collision sources are estimated using data-driven techniques. The observed event yields in the signal regions are consistent with the expected background. Upper limits on the Higgs boson to dark photon branching fraction are reported as a function of the dark photon mean proper decay length or of the dark photon mass and the coupling between the Standard Model and the potential dark sector. This search is combined with previous ATLAS searches obtained in the gluon–gluon fusion and WH production modes. A branching fraction above 10% is excluded at 95% CL for a 125 GeV Higgs boson decaying into two dark photons for dark photon mean proper decay lengths between 173 and 1296 mm and mass of 10 GeV.

Aad, G. (ORCID:0000000266654934)↗
DOE OSTI2025

University of Miami G-band Vapor Radiometer Calibration (UMGVR_CAL) Field Campaign Report

Cold-air outbreak (CAO) clouds in the Arctic are commonly mixed-phase (MP); however, the partitioning of the amount of ice and water in CAO clouds and precipitation is not always well observed. Understanding how cloud phases partition as a function of cloud life cycle is important for predicting snowfall rates, convective life cycle, and intensity at weather timescales. The partitioning into liquid versus ice also has radiative impacts that are consequential for climate. These concerns motivated the incorporation of an airborne G-band vapor radiometer (GVR) into an National Science Foundation-supported aircraft campaign named the Cold Air outbreak Experiment in the Sub-Arctic Region (CAESAR). The GVR is an upward-pointing passive microwave radiometer using four frequencies centered around the 183.31 GHz water vapor absorption line, displaced by +- 1, 3, 7, and 14 GHz. For context, The U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility operates a surface-based GVR at its North Slope of Alaska (NSA) site. The same GVR has been used previously for a field campaign in the southeast Pacific, where an offset was noticed between brightness temperatures (Tbs) measured under clear skies compared to those calculated from a radiative transfer model. To account for any calibration offsets to the Tbs, a request was made to DOE to allow the GVR to operate at ARM’s Southern Great Plains (SGP) observatory and enable comparisons between its measurements and those available at SGP. This request was granted, titled ‘UMGVR_CAL’, short for ‘UMGVR_Calibration’. From October 30 to November 10, 2023, the GVR was deployed to the ARM SGP site to take advantage of their regular, nearby radiosonde launches under clear-sky conditions. The latter were determined using the SGP total sky imagery data. The GVR brightness temperatures in these clear-sky conditions were compared to those calculated by a radiative transfer code (PAMTRA) based on the SGP radiosondes. During the campaign, four suitable clear-sky episodes could be used for the GVR calibration. While few in number, these proved to be enough to satisfy our goal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Strain‐Driven Stabilization of a Room‐Temperature Chiral Multiferroic with Coupled Ferroaxial and Ferroelectric Order

Noncollinear ferroic materials are sought after as testbeds to explore the intimate connections between topology and symmetry, which result in electronic, optical, and magnetic functionalities not observed in collinear ferroic materials. For example, ferroaxial materials have rotational structural distortions that break mirror symmetry and induce chirality. When ferroaxial order is coupled with ferroelectricity arising from a broken inversion symmetry, it offers the prospect of electric-field-control of the ferroaxial distortions and opens up new tunable functionalities. However, chiral multiferroics, especially ones stable at room temperature, are rare. A strain-stabilized, room-temperature chiral multiferroic phase in single crystals of BaTiS 3 is reported here. Using first-principles calculations, the stabilization of this multiferroic phase having P6 3 space group for biaxial tensile strains exceeding 1.5% applied on the basal ab-plane of the room temperature P6 3 cm phase of BaTiS 3 is predicted. The chiral multiferroic phase is characterized by rotational distortions of TiS 6 octahedra around the long c-axis and polar displacement of Ti atoms along the c-axis. Furthermore, the ferroaxial and ferroelectric distortions and their domains in P6 3 -BaTiS 3 are directly resolved using atomic resolution scanning transmission electron microscopy. Landau-based phenomenological modeling predicts a strong coupling between the ferroelectric and the ferroaxial order making P6 3 -BaTiS 3 an attractive test bed for achieving electric-field-control of chirality.

Chalcogenide↗
DOE OSTI2021

Phase transitions and potential ferroelectricity in noncentrosymmetric KNaNbOF 5

Increasing the number of known ferroelectrics requires expanding the chemical design space and mechanisms producing ferroelectricity. To that end, we examine the displacive, noncentrosymmetric-to-centrosymmetric phase transition in the oxyfluoride KNaNbOF 5 using ab initio calculations and Landau theory. Here we predict an intermediate Pnma phase occurs in the transition from the known low- and high-temperature phases, Pna2 1 and Cmcm, respectively, which we subsequently confirm using high-resolution and in situ 19 F magic-angle spinning nuclear magnetic resonance. Using the intermediate structure, we then calculate the monodomain polarization switching barrier in KNaNbOF 5 to be ~ 93 meV per formula unit, comparable to existing ferroelectrics. The reversal of the polarization is due to uncompensated antipolar displacements of oxygen and fluorine and does not require oxygen and fluorine site exchange as in other nonswitchable heteroanionic materials, which makes KNaNbOF 5 a viable oxyfluoride for a switchable electric polarization.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evidence for In-Plane Electrical Polarization in 3R-β’-In 2 Se 3 Thin Films Grown by Molecular Beam Epitaxy

β-In 2 Se 3 has been identified as a potential ferroelectric. This work describes the growth of β-In 2 Se 3 via molecular beam epitaxy along with a description of its properties, including the search for switchable polarization in 3R-β’-In2Se3. The thin-film morphology, crystal structure, and phase maps of β-In 2 Se 3 on Si(111) and Al 2 O 3 (0001) were evaluated by atomic force microscopy, X-ray diffraction, and Raman spectroscopy as a function of the atomic Se/In flux ratio and growth temperature. Smooth β-In 2 Se 3 thin films were successfully realized on Si(111) at a substrate temperature of 150 °C using a Se/In flux ratio of 5.7, as well as on Al 2 O 3 (0001) at 450 °C using a Se/In flux ratio of 5.5. Scanning transmission electron microscopy (STEM) confirms the 3R polytype of β-In 2 Se 3 in films on both substrates, with a minor disorder associated with the 2H polytype at the interface. Indications of in-plane Se atom displacements characteristic of the 3R-β’-In 2 Se 3 polytype were found by STEM and second harmonic generation analysis in films on Al 2 O 3 (0001), but not in films on Si(111). Furthermore, attempts at electrical polarization switching did not produce compelling evidence for ferroelectricity. Instead, electrical transport measurements demonstrated locally varying anisotropic responses with the applied electric fields along different in-plane directions, with some hysteresis associated with the trapping of charges.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Size dependent lattice pseudosymmetry for frustrated decahedral nanoparticles

Geometric frustration—where geometry prevents simultaneous satisfaction of local interactions—generates pseudosymmetry and emergent behaviors across physical and biological systems. At the nanoscale, pseudosymmetric features in crystalline materials manifest as local strain and distortion, but how they depend on particle size and control structural stability remains unclear. Here, we report the first study of a size-dependent crossover in pseudosymmetry in multi-twinned gold nanoparticles (NPs), combining four-dimensional scanning transmission electron microscopy with nanoscale strain mapping grounded in continuum solid mechanics. Analysis of more than 20 decahedral NPs (20–55 nm) reveals pronounced heterogeneity in multiple modes of in-plane strain and displacement field in small NPs as five tetrahedral grains close the geometric gap, without extended defects. With increasing particle size, strain fields homogenize across grains and local phases shift from predominantly low-symmetry body-centered tetragonal motifs at small sizes to face-centered cubic character approaching the bulk limit. We identify a crossover particle size of ~35 nm, well below bulk, correlating with a transition from modified-Wulff shapes to pentagonal bipyramids, consistent with finite element predictions. This quantitative framework for mapping size-dependent strain and pseudosymmetry enables precise design and control of functional crystalline solids and phase transformation for catalysis, photonics, electronics, and energy storage.

Lin, Oliver [University of Illinois at Urbana-Cham↗
DOE OSTI2025

Structural Distortions and Uniaxial Negative Thermal Expansion in the Polar Dion–Jacobson Oxide RbNdTa 2 O 7

We provide deeper insight into the crystal structures, sequential structural phase transitions (I2cm → Cmce → I4/mcm → P4/mmm), thermal expansion, and electronic properties of the n = 2 Dion–Jacobson polar oxide RbNdTa 2 O 7 , through X-ray powder diffraction, neutron powder diffraction, Raman studies, and density functional theory calculations. We observed a uniaxial negative thermal expansion (NTE) across the first-order transition, I2cm → Cmce, where the unit cell contracts along the c-axis, which is driven by a contraction of the NdTa 2 O 6 layer. Here, this NTE occurs within the temperature range of the first-order phase transition and contrasts with the corkscrew mechanism typically observed in Ruddlesden–Popper phases. In RbNdTa 2 O 7 , the I2cm (hybrid improper ferroelectric) → Cmce (antipolar) transition involves crucial changes in the bond lengths of Nd and Ta polyhedra, coupled with polar to antipolar displacement of the Nd ions, leading to a net contraction in the NdTa 2 O 6 layer along the c-axis, while preserving the overall octahedral tilting magnitude. This transition highlights the intricate interplay between the Nd and Ta coordination and the associated TaO 6 distortions. Temperature-dependent Raman spectra analysis further confirms the first-order structural transition and associated NTE, providing evidence for increased bond stiffness across this transition. Additionally, using neutron powder diffraction, we have determined that the transition I4/mcm → P4/mmm occurs at approximately 1150 K. Finally, we have calculated from DFT + U, the partial density of states, the energy bandgaps, and effective masses of the charge carriers of the polar ground structure.

Chemical structure↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2023

Nanoprecipitates to Enhance Radiation Tolerance in High-Entropy Alloys

The growth of advanced energy technologies for power generation is enabled by the design, development, and integration of structural materials that can withstand extreme environments, such as high temperatures, radiation damage, and corrosion. High-entropy alloys (HEAs) are a class of structural materials in which suitable chemical elements in four or more numbers are mixed to typically produce single-phase concentrated solid solution alloys (CSAs). Many of these alloys exhibit good radiation tolerance like limited void swelling and hardening up to relatively medium radiation doses (tens of displacements per atom (dpa)); however, at higher radiation damage levels (>50 dpa), some HEAs suffer from considerable void swelling limiting their near-term acceptance for advanced nuclear reactor concepts. In this study, we developed a HEA containing a high density of Cu-rich nanoprecipitates distributed in the HEA matrix. The Cu-added HEA, NiCoFeCrCu 0.12 , shows excellent void swelling resistance and negligible radiation-induced hardening upon irradiation up to high radiation doses (i.e., higher than 100 dpa). The void swelling resistance of the alloy is measured to be significantly better than NiCoFeCr CSA and austenitic stainless steels. Density functional theory simulations predict lower vacancy and interstitial formation energies at the coherent interfaces between Cu-rich nanoprecipitates and the HEA matrix. The alloy maintained a high sink strength achieved via nanoprecipitates and the coherent interface with the matrix at a high radiation dose (~50 dpa). From our experiments and simulations, the effective recombination of radiation-produced vacancies and interstitials at the coherent interfaces of the nanoprecipitates is suggested to be the critical mechanism responsible for the radiation tolerance of the alloy. Furthermore, the materials design strategy based on incorporating a high density of interfaces can be applied to high-entropy alloy systems to improve their radiation tolerance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electric-field-induced crossover of polarization reversal mechanisms in Al 1–x Sc x N ferroelectrics

Scandium-doped aluminum nitride, Al 1–x Sc x N, represents a new class of displacive ferroelectric materials with high polarization and sharp hysteresis along with high-temperature resilience, facile synthesizability and compatibility with standard CMOS fabrication techniques. The fundamental physics behind the transformation of unswitchable piezoelectric AlN into switchable Al–Sc–N ferroelectrics depends upon important atomic properties such as local structure, dopant distributions and the presence of competing mechanism of polarization switching in the presence of an applied electric-field that have not been understood. We computationally synthesize Al 1–x Sc x N to quantify the inhomogeneity of Sc distribution and phase segregation, and characterize its crystal and electronic structure as a function of Sc-doping. Nudged elastic band calculations of the potential energy surface and quantum molecular dynamics simulations of direct electric-field-driven ferroelectric switching reveal a crossover between two polarization reversal mechanisms—inhomogeneous nucleation-and-growth mechanism originating near Sc-rich regions in the limit of low applied fields and nucleation-limited-switching in the high-field regime. Finally, understanding polarization reversal pathways for these two mechanisms as well as the role of local Sc concentration on activation barriers provides design rules to identify other combinations of dopant elements, such as Zr, Mg etc. to synthesize superior AlN-based ferroelectric materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Failure investigation of nuclear grade POCO graphite target in high energy neutrino physics through numerical simulation

In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.

43 PARTICLE ACCELERATORS↗