Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Cobalt‐Catalyzed Electroreductive Alkylation of Unactivated Alkyl Chlorides with Conjugated Olefins

Abstract Reactions of unactivated alkyl chlorides under mild and sustainable conditions are rare compared to those of alkyl bromides or iodides. As a result, synthetic methods capable of modifying the vast chemical space of chloroalkane reagents, wastes, and materials are limited. We report the cobalt‐catalyzed reductive addition of unactivated alkyl chlorides to conjugated alkenes. Co‐catalyzed activation of alkyl chlorides is performed under electroreductive conditions, and the resulting reactions constitute formal alkyl‐alkyl bond formation. In addition to developing an operationally simple methodology, detailed mechanistic studies provide insights into the elementary steps of a proposed catalytic cycle. In particular, we propose a switch in the mechanism of C−Cl bond activation from nucleophilic substitution to halogen atom abstraction, which is critical for efficiently generating alkyl radicals. These mechanistic insights were leveraged in designing ligands that enable couplings of primary, secondary, and tertiary alkyl chlorides.

Al Zubaydi, Samir↗
DOE OSTI2020

Spin-phonon relaxation from a universal ab initio density-matrix approach

Designing new quantum materials with long-lived electron spin states urgently requires a general theoretical formalism and computational technique to reliably predict intrinsic spin relaxation times. We present a new, accurate and universal first-principles methodology based on Lindbladian dynamics of density matrices to calculate spin-phonon relaxation time of solids with arbitrary spin mixing and crystal symmetry. This method describes contributions of Elliott-Yafet and D’yakonov-Perel’ mechanisms to spin relaxation for systems with and without inversion symmetry on an equal footing. We show that intrinsic spin and momentum relaxation times both decrease with increasing temperature; however, for the D’yakonov-Perel’ mechanism, spin relaxation time varies inversely with extrinsic scattering time. We predict large anisotropy of spin lifetime in transition metal dichalcogenides. The excellent agreement with experiments for a broad range of materials underscores the predictive capability of our method for properties critical to quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

An accurate perturbative approach to redshift space clustering of biased tracers in modified gravity

We extend the scale-dependent Gaussian Streaming Model (GSM) to produce analytical predictions for the anisotropic redshift-space correlation function for biased tracers in modified gravity models. Employing the Convolution Lagrangian Perturbation Theory (CLPT) re-summation scheme, with a local Lagrangian bias schema provided by the peak-background split formalism, we predict the necessary ingredients that enter the GSM, the real-space halo pairwise velocity and the pairwise velocity dispersion. We further consider effective field theory contributions to the pairwise velocity dispersion in order to model correctly its large scale behavior. We apply our method on two widely-considered modified gravity models, the chameleon-screened ƒ(R) Hu-Sawicki model and the nDGP Vainshtein model and compare our predictions against state-of-the-art N-body simulations for these models. We demonstrate that the GSM approach to predict the monopole and the quadrupole of the redshift-space correlation function for halos, gives very good agreement with the simulation data, for a wide range of screening mechanisms, levels of screening and halo masses at z = 0.5 and z = 1. Our work shows the applicability of the GSM, for cosmologies beyond GR, demonstrating that it can serve as a powerful predictive tool for the next stage of cosmological surveys like DESI, Euclid, LSST and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Angular-momentum projection in coupled-cluster theory: Structure of 34 Mg

Single- reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for 8 Be and 20 Ne using a two-nucleon interaction from chiral effective field theory and for pf-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus 34 Mg and find agreement with data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

A stochastic model of cell survival for high-Z nanoparticle radiotherapy

Purpose: The authors present a stochastic framework for the assessment of cell survival in gold nanoparticle radiotherapy. Methods: The authors derive the equations for the effective macroscopic dose enhancement for a population of cells with nonideal distribution of gold nanoparticles (GNP), allowing different number of GNP per cell and different distances with respect to the cellular target. They use the mixed Poisson distribution formalism to model the impact of the aforementioned physical factors on the effective dose enhancement. Results: The authors show relatively large differences in the estimation of cell survival arising from using approximated formulae. They predict degeneration of the cell killing capacity due to different number of GNP per cell and different distances with respect to the cellular target. Conclusions: The presented stochastic framework can be used in interpretation of experimental cell survival or tumor control probability studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Updated Virophage Taxonomy and Distinction from Polinton-like Viruses

Virophages are small dsDNA viruses that hijack the machinery of giant viruses during the co-infection of a protist (i.e., microeukaryotic) host and represent an exceptional case of “hyperparasitism” in the viral world. While only a handful of virophages have been isolated, a vast diversity of virophage-like sequences have been uncovered from diverse metagenomes. Their wide ecological distribution, idiosyncratic infection and replication strategy, ability to integrate into protist and giant virus genomes and potential role in antiviral defense have made virophages a topic of broad interest. However, one limitation for further studies is the lack of clarity regarding the nomenclature and taxonomy of this group of viruses. Specifically, virophages have been linked in the literature to other “virophage-like” mobile genetic elements and viruses, including polinton-like viruses (PLVs), but there are no formal demarcation criteria and proper nomenclature for either group, i.e., virophage or PLVs. Here, as part of the ICTV Virophage Study Group, we leverage a large set of genomes gathered from published datasets as well as newly generated protist genomes to propose delineation criteria and classification methods at multiple taxonomic ranks for virophages ‘sensu stricto’, i.e., genomes related to the prototype isolates Sputnik and mavirus. Based on a combination of comparative genomics and phylogenetic analyses, we show that this group of virophages forms a cohesive taxon that we propose to establish at the class level and suggest a subdivision into four orders and seven families with distinctive ecogenomic features. Finally, to illustrate how the proposed delineation criteria and classification method would be used, we apply these to two recently published datasets, which we show include both virophages and other virophage-related elements. Overall, we see this proposed classification as a necessary first step to provide a robust taxonomic framework in this area of the virosphere, which will need to be expanded in the future to cover other virophage-related viruses such as PLVs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Numerical evaluation of the bispectrum in multiple field inflation—the transport approach with code

We present a complete framework for numerical calculation of the power spectrum and bispectrum in canonical inflation with an arbitrary number of light or heavy fields. Our method includes all relevant effects at tree-level in the loop expansion, including (i) interference between growing and decaying modes near horizon exit; (ii) correlation and coupling between species near horizon exit and on superhorizon scales; (iii) contributions from mass terms; and (iv) all contributions from coupling to gravity. We track the evolution of each correlation function from the vacuum state through horizon exit and the superhorizon regime, with no need to match quantum and classical parts of the calculation; when integrated, our approach corresponds exactly with the tree-level Schwinger or 'in-in' formulation of quantum field theory. In this paper we give the equations necessary to evolve all two- and three-point correlation functions together with suitable initial conditions. The final formalism is suitable to compute the amplitude, shape, and scale dependence of the bispectrum in models with | f {sub NL}| of order unity or less, which are a target for future galaxy surveys such as Euclid, DESI and LSST. As an illustration we apply our framework to a number of examples, obtaining quantitatively accurate predictions for their bispectra for the first time. Two accompanying reports describe publicly-available software packages that implement the method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

Here, we investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core–valence separated (FC-CVS) EOM-CCSD approach – a successful variant of EOM-CCSD tailored towards core excitations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Relativistic Douglas–Kroll–Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas–Kroll–Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin–orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin–orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Chemistry↗
DOE OSTI2016

SU-G-BRC-16: Theory and Clinical Implications of the Constant Dosimetric Leaf Gap (DLG) Approximation

Purpose: Commercial dose calculation algorithms incorporate a single DLG value for a given beam energy that is applied across an entire treatment field. However, the physical processes associated with beam generation and dose delivery suggest that the DLG is not constant. The aim of this study is to evaluate the variation of DLG among all leaf pairs, to quantify how this variation impacts delivered dose, and to establish a novel method to correct dose distributions calculated using the approximation of constant DLG. Methods: A 2D diode array was used to measure the DLG for all 60 leaf pairs at several points along each leaf pair travel direction. This approach was validated by comparison to DLG values measured at select points using a 0.6 cc ion chamber with the standard formalism. In-house software was developed to enable incorporation of position dependent DLG values into dose distribution optimization and calculation. The accuracy of beam delivery of both the corrected and uncorrected treatment plans was studied through gamma pass rate evaluation. A comparison of DVH statistics in corrected and uncorrected treatment plans was made. Results: The outer 20 MLC leaf pairs (1.0 cm width) have DLG values that are 0.32 mm (mean) to 0.65 mm (maximum) lower than the central leaf-pair. VMAT plans using a large number of 1 cm wide leaves were more accurately delivered (gamma pass rate increased by 5%) and dose coverage was higher (D100 increased by 3%) when the 2D DLG was modeled. Conclusion: Using a constant DLG value for a given beam energy will result in dose optimization, dose calculation and treatment delivery inaccuracies that become significant for treatment plans with high modulation complexity scores delivered with 1 cm wide leaves.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

MRI-Based Radiotherapy Planning to Reduce Rectal Dose in Excess of Tolerance

Background and Purpose. Chronic rectal toxicity significantly decreases the quality of life for men who receive radiotherapy for prostate cancer. The most significant predictor of rectal toxicity is rectal dose-volume exceeding tolerance. To minimize the volume of rectum in the high dose field, it is essential to accurately define the prostate-rectum interface. This can be challenging to do by computed tomography (CT) imaging alone. The current study was undertaken to formally demonstrate in a clinical trial setting that image-guided intensity-modulated radiation therapy (IG-IMRT) planning using magnetic resonance imaging (MRI) can reduce the volume of rectum exceeding 70 Gy, a validated metric that predicts the risk of late rectal toxicity. Materials and Methods. This prospective single-arm study enrolled 15 men treated with IG-IMRT for localized prostate cancer. All participants received a dedicated 3 Tesla MRI examination of the prostate in addition to a pelvic CT examination for treatment planning. Two volumetric modulated arc therapy (VMAT) plans with a prescription dose of 79.2 Gy were designed using identical constraints based on CT- and MRI-defined consensus volumes. The volume of rectum exposed to 70 Gy or more was compared using the Wilcoxon paired signed rank test. Results. For CT-based treatment plans, the median volume of rectum receiving 70 Gy or more was 9.3 cubic centimeters (cc) (IQR 7.0 to 10.2) compared with 4.9 cc (IQR 4.1 to 7.8) for MRI-based plans. This resulted in a median volume reduction of 2.1 cc (IQR 0.5 to 5.3, P < .001). Conclusions. Using MRI to plan prostate IG-IMRT to a dose of 79.2 Gy reduces the volume of rectum receiving radiation dose in excess of tolerance (70 Gy or more) and should be considered in men who are at high risk for late rectal toxicity and are not good candidates for other rectal sparing techniques such as hydrogel spacer. This trial is registered with NCT02470910.

Schmidt, Daniel R.↗
DOE OSTI2024

The Transactive Energy Network Template Metamodel

While transactive energy, which is defined as an allocation of electricity based on dynamically discovered values or prices, has been extensively studied, its uptake and use has been slow. This report describes a tool, the transactive network template, which should hasten the creation and uptake of transactive energy networks. Some basic principles of transactive energy are familiar from existing wholesale electricity markets. Locational prices are calculated today for zones within bulk electric transmission systems. Locational prices differ while accounting for the locational costs of electricity generation and the losses and constraints incurred when electricity is transmitted from generators and distributed to consumers. A transactive energy network might include these transmission zones. However, current research strives to apply transactive energy also in electricity distribution circuits, buildings, and even for individual generating and consuming devices. At the same time, researchers explore how to apply transactive energy in real time during increasingly shorter time intervals. Automated computational agents become necessary as transactive energy becomes applied to smaller circuit zones and at faster dynamic timescales. A transactive energy network is an example of a multi-agent system. Each zone in the network is represented by its transactive agent, which makes decisions for and acts on behalf of a business entity that is responsible for and manages one of the circuit regions. A transactive energy network is also an example of a decentralized, distributed control system. Control decisions and responsibilities are distributed among the network’s transactive agents. The transactive agents are independent; that is, there typically is no centralized authority or oversight function. Instead, transactive agents exchange transactive signals and thereby negotiate the prices and quantities of electricity that they will exchange. Initially, the circuit regions and responsibilities of transactive agents appear to be very dissimilar. Each circuit region may comprise transmission, distribution, or building-level circuits. Each has a unique position and electrical connectivity within the transactive energy network. Each possesses unique assets that either generate or consume electricity, and these (e.g., renewable energy generator, diesel generator, aggregate utility load, building load, space conditioning, refrigerator, etc.) may further differ in their price flexibility and in their strategies for responding to dynamic electricity prices. Given such diversity, an implementer’s first inclination might be to start from scratch to define all these devices and to engineer their seemingly unique interactions. Given that each implementer’s perspective may be narrow within a transactive energy network, it is unlikely that uniquely engineered systems would interact well. This is where the transactive network template is applicable. The transactive network template is a metamodel that has been developed to guide implementers as they configure their own transactive agent within a network of such agents. The object-oriented design of the transactive network template provides basic code object types that may be used and extended by implementers to represent each of the assets in their circuit region. These objects further facilitate the transactive agent’s necessary computations, which are divided among responsibilities to schedule power usage, balance electric supply and demand, and coordinate the exchange of electricity with the other transactive agents. This report addresses the conceptual transactive network template design. Implementers are directed to more formal design documents and reference implementations. A Python™-based1 reference implementation of the transactive network template has been coded, and three implementations have been configured to represent a national laboratory and two university campuses. Version 2 of the transactive node template generalizes the market class and its methods to facilitate multiple, and more diverse market coordination mechanisms than were facilitated by and demonstrated using Version 1. Version 3 includes new Appendix B, which addresses the designs of methods that would make dynamic prices track approved electricity rates. In the future, the author wishes to make the transactive network template more generally applicable to networks that require more accurate power flow. Development of the transactive network template is jointly funded by the U.S. Department of Energy (DOE) Energy Efficiency and Renewable Energy and the DOE Office of Electricity. In late 2015, one of the first projects to be funded by the DOE Grid Laboratory Modernization Laboratory Consortium was the Clean Energy and Transactive Campus project, led by Pacific Northwest National Laboratory. DOE funds were matched by an investment by the Washington Department of Commerce through its Clean Energy Fund. The transactive network template was developed to guide the implementation of transactive energy networks within this project’s scope.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

A Quantum-mechanical Treatment of Electron Broadening in Strong Magnetic Fields

Abstract Spectral line-shape fitting is an extremely useful tool in determining the gravity of white dwarf stars. This method is so far limited to nonmagnetic white dwarfs largely because the theory of line broadening in high magnetic fields is not as complete as in the nonmagnetic case. Current Stark+Zeeman models treat plasma particles classically and ignore the motion of the nucleus. We develop the formalism for a quantum-mechanical treatment of the perturbing electrons and include the nuclear motion as part of the broadening and explore their relative importance. The conditions we explore are those found in white dwarf and neutron star atmospheres. We find that, contrary to previous studies, the quantized perturbing electrons create more broadening than perturbers on a straight-path trajectory. Additionally, the quantization of the plasma electrons gives rise to resonances away from the line center. The nuclear motion creates an additional electric field, which also leads to an increase in line broadening; however, this effect in neutron star atmospheres is not as large as previously estimated. This suggests that neutron star spectral lines are sensitive to density and that their mass and radius can be obtained from spectral line fitting, which would help constrain the neutron star equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Novel Radiometry for In-Pipe Robotic Inspection of Holdup Deposits in Gaseous Diffusion Piping - 19503

This paper details the first complete, DOE approved, deployed and operational radiometry for robotic nondestructive assay (NDA) of holdup deposits in gaseous diffusion piping. Some features, like peak-finding, Compton correction, detector efficiency, and material properties are characteristic of typical manual methods. Most features are specific to in-pipe robotic deployment. These include in-motion radiometry, self-determined location, geometric modeling of thick deposits, bounding of self-attenuation thickness and auto-checking of replicate measurements. Significant are the means for determining all the properties, equations, and constants for automatically computing and displaying all the quantities, uncertainties and information needed for NDA reporting, analysis, and review. Earlier work by this team demonstrated rudimentary but convincing robotic in-pipe NDA of holdup deposits in pipes. Cold testing with high-enrichment mat sources and hot tests with low U-235 loadings together succeeded to exhibit convincing proof-of-principle. This research-grade radiometry implementation lacked accurate odometry, efficiency calibration, high-fidelity detector modeling, attenuation modeling, forward geometric modeling, auto-segmenting, uncertainty modeling, replicate checking, deconvolution and much more. This work explains the technical basis behind robotic in-pipe NDA. It includes discussion of a new method for modeling of self-attenuation using in-pipe deposit geometry information, detail of the method's total measurement uncertainty, and explanation of calibration methods and results. Beyond development of automation-specific radiometric elements, the paper highlights analyses, testing, and technical bases specific to and supportive of the methodology. Examples include collimation characterization, detector efficiency determination, auto-calibration methodology, and odometric replicate-checking. Results of this methodology are formalized in a DOE EM Technical Basis Document, implemented in software, formally verified in review and acceptance testing and proceduralized as a means for D and D NDA at the Portsmouth enrichment facility. Measurement of U-235 quantity is a well-known cost and schedule driver in D and D of every defunct gaseous diffusion enrichment facility in the world. The radiometry methodologies and results detailed in this paper were the high hurdles. Having established the technical basis for this automated NDA, the opportunity is to apply this to many sizes and robot forms for addressing holdup deposit NDA throughout the D and D community. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Development of Human Factors Awareness for Criticality Safety Specialists at Westinghouse Nuclear Fuels - 14223

Although human factors assessment typically involves specialist expertise, it is considered that enhanced local knowledge of this area for criticality safety specialists supports the incorporation of human factors consideration into tasks such as safety assessment, design and hazard reviews, development of operational and engineered criticality controls, training, incident investigation and emergency response. At the Westinghouse Springfields Fuel Manufacturing Site in the United Kingdom, specialist human factor support has been engaged to identify key competencies and training plans for criticality safety specialists. In addition, discussions have been held with the UK's Office for Nuclear Regulation (ONR) to outline the proposal and to discuss the approach. Following completion of an internationally recognised human factors training course the learning is now being applied by criticality safety specialists. Recent examples of learning applied include: 1) Review of past events and incidents to support the consideration of operator/equipment interfaces, leading to improved designs and processes. 2) Identification of critical tasks and the level of human factor assessment that is likely to be required. To support this, guidance has been produced by the criticality section for staff and contractors performing criticality safety assessments. 3) Assessment of the 'aims and operability' of engineered and procedural controls derived as part of the criticality safety assessment process. This involves participation of plant management and operators. 4) Consideration of human factors in emergency response activities such as infrequent tasks and tasks being undertaken under time pressures etc. This has been examined in recent emergency exercises and in the development of a recent criticality exercise scenario. 5) Criticality training has been reviewed with Plant Management, focussing on the methods utilised to deliver training (e.g. classroom based, computer training, on plant walk-downs); the aims and scope of training and how understanding is validated. 6) The link between criticality safety assessment and how the operations are / will be conducted has been reviewed as part of formal and informal criticality meetings and inspections. A key aspect of this is the criticality specialist's interaction with both plant management and operational staff. The development of human factors awareness into the criticality safety specialist's competency framework continues to be developed at Westinghouse but real benefits to criticality safety are already being demonstrated by this approach. (author)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Covariant phase space with boundaries

The covariant phase space method of Iyer, Lee, Wald, and Zoupas gives an elegant way to understand the Hamiltonian dynamics of Lagrangian field theories without breaking covariance. The original literature however does not systematically treat total derivatives and boundary terms, which has led to some confusion about how exactly to apply the formalism in the presence of boundaries. In particular the original construction of the canonical Hamiltonian relies on the assumed existence of a certain boundary quantity “B”, whose physical interpretation has not been clear. We here give an algorithmic procedure for applying the covariant phase space formalism to field theories with spatial boundaries, from which the term in the Hamiltonian involving B emerges naturally. Our procedure also produces an additional boundary term, which was not present in the original literature and which so far has only appeared implicitly in specific examples, and which is already nonvanishing even in general relativity with sufficiently permissive boundary conditions. The only requirement we impose is that at solutions of the equations of motion the action is stationary modulo future/past boundary terms under arbitrary variations obeying the spatial boundary conditions; from this the symplectic structure and the Hamiltonian for any diffeomorphism that preserves the theory are unambiguously constructed. We show in examples that the Hamiltonian so constructed agrees with previous results. We also show that the Poisson bracket on covariant phase space directly coincides with the Peierls bracket, without any need for non-covariant intermediate steps, and we discuss possible implications for the entropy of dynamical black hole horizons.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗