Search NASA⌕ Search

SEARCH · Search NASA

Results for “space interfaces”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Simulating quantum-classical interfaces via the Lindblad master equation

In hybrid quantum systems, the interface between quantum and classical domains is essential for the generation, control, and measurement of quantum states. Quantum-classical interfaces (QCIs) are ubiquitous in devices such as optical modulators, quantum sensors, and signal processors, where classical signals influence quantum dynamics. In this paper, we employ the Lindblad master equation to simulate the evolution of a quantum system interacting with a classical control system. Our model captures both linear and nonlinear interactions by incorporating first- and second-order susceptibilities, and it quantifies the influence of externally applied control parameters on decoherence and state evolution. As an illustrative example, we analyze an optical modulator and demonstrate how variations in material response and drive conditions affect photon statistics, coherence, and phase-space distributions. In conclusion, the findings offer a path to an all-encompassing model for understanding and optimizing QCIs, with wide-ranging implications for the performance, design, and robustness of next-generation quantum devices.

Quantum engineering↗
DOE OSTI2024

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The Electron Spectro-Microscopy (ESM) Beamline at NSLS-II

Photoelectron spectroscopy is a primary tool for the study of the electronic structure of materials and the chemical composition of surfaces. High-resolution angle-resolved photoemission spectroscopy (ARPES) has the unique ability to map the energy bands in momentum space. Furthermore, going beyond the single particle picture, the self-energy corrections caused by correlations in solids can be extracted from the analysis of the emission line shape. The current level of refinement, in terms of energy and angular resolution (ΔE < 1 meV, Δθ < 0.1°), makes the technique sensitive to the lowest energy excitations and the dynamics of electrons, which in turn virtually determine all the macroscopic properties of any system and govern the chemical, electrical, magnetic, and physical processes. Similarly important, X-ray photoelectron microscopy (XPEEM), combined with the low-energy electron microscopy (LEEM), is indispensable in probing the complexity of chemical, structural, electronic and magnetic properties of surfaces and shallow interfaces, with the spatial resolution of few tens of nanometer (nm). The Electron-Spectro-Microscopy beamline (ESM) has been recently commissioned at NSLS-II and is now in operation. The primary spectroscopic technique is photoemission, performed over a wide energy range with control of light polarization and in a variety of flux/resolution conditions. The beamline has two experimental end stations that allow to perform ARPES and XPEEM/LEEM, separately. The ARPES end station focuses on high energy-resolution work, with spot-size of a few microns. The XPEEM/LEEM end station is a full-field microscope (XPEEM) operating either with the synchrotron generated X-rays (XPEEM), or with an internal electron gun (LEEM). Spatial resolution is crucial in studies of newly synthesized complex materials since they are often initially available only as small specimens (typically micron size). Furthermore, chemical inhomogeneities on surfaces are often an integral part of surface chemical processes. Finally, the ESM beamline with X-ray spots of few microns is optimized to study the electronic structure of novel materials with microscopy capabilities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Laboratory Demonstration of Collisionless Blob Formation via Laser-Produced Plasma Self-Focusing

Strongly localized, propagating plasma density structures that are capable of crossing magnetic field lines are known as “blobs.” Here we demonstrate a novel mechanism for the formation and propagation of an ion gyroradius-scale blob-cavity structure at the interface between a super-Alfvénic laser-produced plasma (LPP) and an ambient magnetized plasma. The LPP self-focuses along the edge of the diamagnetic cavity which results in a dense, jetlike structure as compared to ballistic motion. This collimated flow couples momentum to the ambient plasma through a collisionless process known as Larmor coupling. The Larmor electric fields locally displace the ambient ions forming a blob above the LPP flow. In the region between a gyrating blob and collimated LPP flow, a secondary cavity of expelled magnetic field forms. In conclusion, these findings are supported by particle-in-cell simulations that replicate the blob formation mechanism and provide insight to similar processes in space, astrophysical, and laboratory settings characterized by ion kinetic scales.

Fluorescence spectroscopy↗
DOE OSTI2025

Advanced Manufactured Strain Sensors for Extreme Environments

Experimentation at irradiation test facilities are essential for reducing the innovation time of developmental fuels, fuel cladding, and structural materials employed in next-generation nuclear reactors. However, due to the harsh conditions generated in such reactors and limited instrumentation space, the evaluation of a material’s mechanical properties is often limited to characterization after the materials have been removed from reactor conditions, or post-irradiation examination. These experiments are costly, time-consuming, and fail to capture the critical time-evolving phenomena that occur during the irradiation experiments. Advanced manufactured digital image correlation patterns and strain sensor devices serve as two promising technologies that can be deployed in the confined and challenging orientations of these irradiation experiments while also providing key insight on salient materials phenomena (i.e., mechanical properties). To guide the development of the printed strain sensors prior to their deployment in critical experiments, the adhesion strength between the substrate and printed film interface is measured via tensile testing and a non-contact laser-induced spallation technique. The establishment of these process control steps helped guide the successful fabrication and testing of direct-write strain sensing devices discussed in this work. The fabrication process controls are necessary for enabling the sustained operation of these strain sensing device through experimentation and minimize the potential for premature failure.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Advanced Manufactured Strain Sensors for Extreme Environments

Experimentation at irradiation test facilities are essential for reducing the innovation time of developmental fuels, fuel cladding, and structural materials employed in next-generation nuclear reactors. However, due to the harsh conditions generated in such reactors and limited instrumentation space, the evaluation of a material’s mechanical properties is often limited to characterization after the materials have been removed from reactor conditions, or post-irradiation examination. These experiments are costly, time-consuming, and fail to capture the critical time-evolving phenomena that occur during the irradiation experiments. Advanced manufactured digital image correlation patterns and strain sensor devices serve as two promising technologies that can be deployed in the confined and challenging orientations of these irradiation experiments while also providing key insight on salient materials phenomena (i.e., mechanical properties).To guide the development of the printed strain sensors prior to their deployment in critical experiments, the adhesion strength between the substrate and printed film interface is measured via tensile testing and a non-contact laser-induced spallation technique. The establishment of these process control steps helped guide the successful fabrication and testing of direct-write strain sensing devices discussed in this work. The fabrication process controls are necessary for enabling the sustained operation of these strain sensing device through experimentation and minimize the potential for premature failure.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

A Computational Tool Compatible with NEAMS Code Packages for Optimizing the Shape of Nuclear Reactor Components and of Whole Core Performance

We designed and implemented a shape optimization tool that functions with NEAMS codes, and that nuclear scientists and engineers can employ to optimize the shape of individual components and the whole core under the applicable single- or multi-physics model comprising the employed code(s). The shape-optimization tool enables varying the geometric shape itself as well as its dimensions to yield, potentially, new component designs that are not limited by the designer’s intuition and previous experience. In cases where the optimal-shape object is an individual component, we provide the capability for additional verification that the whole-core performance using the optimized component performs better, under the prescribed optimization criteria, than the initial design. Our shape-optimization tool couples to NEAMS codes via a flexible input- composer interface and enables the user to constrain the shape’s evolution to ensure the component’s manufacturability. Finally, we demonstrate our shape-optimization tool with single- and multi-physics NEAMS codes. This objective is motivated by the recent advances in manufacturing technology that, combined with rising interest in novel reactor concepts, are creating new opportunities for innovation in the design of individual components that affect the performance of the full reactor system. In particular, Additive Manufacturing (AM) enables mass production of highly precise, intricate and complex component shapes that are not feasible with traditional manufacturing techniques. To accomplish this goal we developed and implemented in MOOSE: (1) discrete shape optimization capability based on a state-space search that uses Artificial Intelligence strategies to find the optimal state/shape; (2) smooth shape optimization tool that employs PETSc’s toolkit for advanced optimization (TAO) to optimize node-displacement of the components’ model sidesets; (3) hierarchical core optimization workflow that recognizes the repeating patterns typical in a nuclear reactor and performs the optimization one level at a time with increasing length scale. Each of these tools is equipped with user-specified constraints to avoid optimal shapes that are not manufacturable. The developed shape optimization tool is verified and demonstrated on various nuclear reactor core components and models. The optimization process accounts for tightly coupled physics that govern the behavior of these target reactors, and exercises several NEAMS codes in a coupled multiphysics fashion. The impact of the delivered shape optimization tool will materialize in the optimal design, from the outset, of advanced reactors currently contemplated to regain the US’s leadership in nuclear energy R&D. Novel reactor concepts, e.g. Molten Salt Reactors, and sizes/capacities, e.g. micro- reactors, provide a unique opportunity to optimize performance from the early stages of development, before the investment in components’ production lines, validation experiments, and licensing regimes make future improvements in performance prohibitively expensive and force sub-optimal performance on the affected reactor concept in perpetuity. This benefit will be realized by the delivered shape optimization tool regardless of the applicable manufacturing process whether traditional or AM, thereby broadening the impact of this project on current and future reactor concepts and technologies

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Phosphonic-Acid-Reinforced Polymer Hole Transport Layers for Deployable p-i-n Perovskite Photovoltaics

The long-term durability prospects of halide perovskite solar cells are rapidly improving; however, the interface between the hole transport layer (HTL) and the perovskite remains a source of degradation. Alone, polymer- or carbazole-based HTLs suffer from incomplete coverage of the underlying indium tin oxide glass, leading to degradation and compromised performance. Here, we show a multi-HTL approach whereby a polymer HTL is reinforced using a phosphonic acid modification leading to better protection of the buried perovskite interface and more columnar growth of perovskite film, resulting in an ~40-mV open-circuit voltage (VOC) improvement indicative of suppressed interfacial recombination across multiple p-i-n device architectures. Solar cells with this reinforced HTL show higher tolerance to several accelerated stress tests. We report, among the best durabilities for unencapsulated cells, a T90 ~3,000 h (T80 ~5,900 h) at 65 degrees C under continuous 1.2 sun AM 1.5G illumination and maximum power point tracking, representing a nearly 4-fold increase compared with [2-(9H-carbazol-9-yl)ethyl]phosphonic acid (2PACz)-only devices. Furthermore, we deployed a device with this reinforced HTL on a cube satellite, with long-duration operational space testing results exceeding T80 for the complete mission duration of ~100 days in low Earth orbit.

14 SOLAR ENERGY↗
DOE OSTI2026

Urea-to-Ammonia Conversion at Proteus mirabilis Modified Pt–Ni/BDD Electrodes

Efficient wastewater recycling technologies are essential for long-duration space missions and sustainable water management on Earth. Here, a bioelectrochemical system integrating Proteus mirabilis with a platinum–nickel-modified boron-doped diamond electrode (Pt–Ni/BDDE) for urea-to-ammonia conversion in synthetic urine is presented. Immobilized P. mirabilis catalyzes enzymatic ureolysis, converting urea into ammonia, which is subsequently oxidized electrochemically, and no direct electrochemical urea oxidation is observed. Cyclic voltammetry (CV) of P. mirabilis on Pt–Ni/BDDE in 0.1 M urea and synthetic urine revealed a broad anodic oxidation peak at approximately 0.55–0.75 V vs Ag/AgCl (sat. KCl), corresponding to ammonia oxidation. Control experiments using bare BDDE and Pt–Ni/BDDE in synthetic urine showed no oxidation peak, establishing that the bioelectrocatalytic response originates exclusively from the bioelectrode interface. Chronoamperometry studies revealed that immobilization potential and time critically influenced bacterial adhesion and electrochemical response, with optimal conditions yielding a maximum current density of 0.0055 mA·cm –2 . These quantitative results establish that microbial ureolysis can be efficiently coupled with advanced electrode materials for urea-to-ammonia conversion, offering a promising self-sustaining strategy for urine processing and water recovery in closed-loop life support systems.

Ammonia↗
DOE OSTI2025

PubChemLite Plus Collision Cross Section (CCS) Values for Enhanced Interpretation of Nontarget Environmental Data

Finding relevant chemicals in the vast (known) chemical space is a major challenge for environmental and exposomics studies leveraging nontarget high resolution mass spectrometry (NT-HRMS) methods. Chemical databases now contain hundreds of millions of chemicals, yet many are not relevant. This article details an extensive collaborative, open science effort to provide a dynamic collection of chemicals for environmental, metabolomics, and exposomics research, along with supporting information about their relevance to assist researchers in the interpretation of candidate hits. The PubChemLite for Exposomics collection is compiled from ten annotation categories within PubChem, enhanced with patent, literature and annotation counts, predicted partition coefficient (logP) values, as well as predicted collision cross section (CCS) values using CCSbase. Monthly versions are archived on Zenodo under a CC-BY license, supporting reproducible research, and a new interface has been developed, including historical trends of patent and literature data, for researchers to browse the collection. This article details how PubChemLite can support researchers in environmental and exposomics studies, describes efforts to increase the availability of experimental CCS values, and explores known limitations and potential for future developments. The data and code behind these efforts are openly available.

PubChem↗
DOE OSTI2026

Tuning Methanol Transformation Pathways for Sustainable Steam Reforming: Na-Promotion Effects on Ag/m-ZrO 2 Catalysts

This work investigates the influence of sodium promotion on Ag/m-ZrO 2 catalysts for methanol steam reforming (MSR), focusing on activity, selectivity, surface chemistry, and mechanistic pathways. Temperature programmed reduction (TPR), XANES/EXAFS, CO 2 TPD, DRIFTS, and temperature programmed surface reaction methods were combined with fixed bed MSR testing to develop an integrated structure–function understanding of Na-modified Ag-ZrO 2 interfaces. Na addition systematically increases surface basicity, stabilizes strongly basic O 2− sites, and weakens the ν(CH) vibrational mode of surface formate, thereby facilitating C–H bond scission and accelerating decarboxylation to CO 2 . At moderate promoter levels (0.5–1.0 wt.% Na), the catalysts show significantly enhanced CO 2 selectivity and increased conversion relative to unpromoted Ag/m-ZrO 2 , while CH 4 formation remains negligible. Excessive Na (≥1.8 wt.%) leads to slower formate decomposition, greater carbonate stabilization, and suppressed conversion, revealing a narrow optimum around 1 wt.% Na. Short-term stability testing demonstrates steady conversion and product selectivity for both unpromoted and Na-promoted catalysts, with the latter maintaining markedly higher CO 2 selectivity. Although Pt/YSZ retains far superior intrinsic activity at ~10× higher space velocity, Ag offers a cost-advantaged alternative where lower cost metals are desirable. Collectively, these findings show that Na promotion enables tunable MSR selectivity on Ag/m-ZrO 2 by directing formate decomposition toward the CO 2 -forming pathway.

CO2 selectivity↗
DOE OSTI2024

Differentiable modeling and optimization of non-aqueous Li-based battery electrolyte solutions using geometric deep learning

Electrolytes play a critical role in designing next-generation battery systems, by allowing efficient ion transfer, preventing charge transfer, and stabilizing electrode-electrolyte interfaces. In this work, we develop a differentiable geometric deep learning (GDL) model for chemical mixtures, DiffMix, which is applied in guiding robotic experimentation and optimization towards fast charging battery electrolytes. In particular, we extend mixture thermodynamic and transport laws by creating GDL-learnable physical coefficients. We evaluate our model with mixture thermodynamics and ion transport properties, where we show improved prediction accuracy and model robustness of Diff-Mix than its purely data-driven variants. Furthermore, with a robotic experimentation setup, Clio, we improve ionic conductivity of electrolytes by over 18.8% within 10 experimental steps, via differentiable optimization built on DiffMix gradients. By combining GDL, mixture physics laws, and robotic experimentation, DiffMix expands the predictive modeling methods for chemical mixtures and enables efficient optimization in large chemical spaces.

25 - ENERGY STORAGE↗
DOE OSTI2025

Particle Acceleration in Collisionless Magnetically Arrested Disks

We present the first collisionless realization of two-dimensional axisymmetric black hole accretion consistent with a persistent magnetically arrested disk state. The accretion flow, consisting of an ion-electron disk plasma combined with magnetospheric pair creation effects, is simulated using first-principles general-relativistic particle-in-cell methods. The simulation is evolved over significant dynamical timescales during which a quasisteady accretion state is reached with several magnetic flux eruption cycles. We include a realistic treatment of inverse Compton scattering and pair production, which allows for studying the interaction between the collisionless accretion flow and pair-loaded jet. Our findings indicate that magnetic flux eruptions associated with equatorial magnetic reconnection within the black hole magnetosphere and the formation of spark gaps are locations of maximal particle acceleration. Flux eruptions, starting near the central black hole, can trigger Kelvin-Helmholtz-like vortices at the jet-disk interface that facilitate efficient mixing between disk and jet plasma in this region. Transient periods of increased pair production following magnetic flux eruptions and reconnection events are responsible for most of the highly accelerated particles.

Accretion disk & black-hole plasma↗
DOE OSTI2024

Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation

Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.

58 GEOSCIENCES↗
DOE OSTI2024

Characterization of compositional variation and solidification condition on rapid solidification behavior of beta stabilized titanium

Rapid solidification (RS) is a field of study that has increased in importance due to the rise of additive manufacturing (AM) techniques. The breakdown of equilibrium at the solid-liquid interface results in an inability to use traditional solidification models that have been applied to great success in casting. As such, there is a current need for a high-throughput testing method to screen alloys at rapid solidification conditions without the significant cost for powder atomization. In this work, twin hammer splat quenching (SQ) is utilized and theoretically developed as a technique uniquely placed to test alloy systems for AM applications. Using analytical thermal spray models and numerical heat transfer simulations, the cooling rate and solidification conditions of splat quenching are compared across several parameters to evaluate the key parameters for generating cooling rates relevant to AM process. Once a framework for the cooling rates and relevancy of SQ are developed, SQ was used to process alloy compositions in the Ti5553 alloy space (AMS 7026) to study the formation of detrimental secondary phases occurring upon rapid quenching from the melt. A martensitic phase is found to occur at the lower end of the standard, but the composition listed in the standard were found to be safe. Finally, the rapid solidification behavior of Ti-5553 is evaluated by comparing splat quenching, laser welding, and plasma atomization that offer a wide window of solidification conditions. The segregation and solidification structure are analyzed and compared to models to provide a deeper understanding of the solidification behavior of beta stabilized titanium alloys.

Williamson, Christopher J.↗
DOE OSTI2025

ArborX 2.0

ArborX library tackles a problem of efficiently finding geometric objects that are close in space. Variations of this problem, such as finding the nearest neighbors of a point, or finding all objects within a certain distance, are inherent components of applications in many fields. The data may be large so that solving the problem efficiently may require significant computational resources, such as multiple processors or accelerators such as general purpose GPUs. ArborX' main advantage in its ability to solve large problems efficiently utilizing a combination of distributed and on-node parallelism. ArborX can be run efficiently on a wide variety of hardware, including GPUs from different vendors, which distinguishes it from other available libraries which typically choose only few of these. The other advantage is that it supports both types of user problems: spatial problems (useful for intersections and finding objects within certain distance), and nearest neighbor problems. ArborX also supports flexible interface in its interaction with a user. Particularly, it allows a user to call user's own function on a positive match, a functionality not rarely available in other libraries. ArborX implements construction and traversal algorithms using efficient tree structures, such as bounding volume hierarchy (BVH). At its core, ArborX uses linear BVH for its low construction cost and sufficient quality. ArborX implements both spatial and nearest-neighbor traversal algorithms. ArborX also provides several clustering algorithms (minimum spanning tree, DBSCAN, HDBSCAN*), interpolation using minimum least squares and ray tracing. ArborX is written using C++, and is parallelized using the message passing interface (MPI) for the distributed communication, and the Kokkos library for on-node parallelism. This approach allows ArborX to be run on a wide variety of hardware, from common laptops and desktops to supercomputers while using the same codebase.

Prokopenko, Andrey [Oak Ridge National Laboratory ↗
DOE OSTI2024

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗