Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time minimization”

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 559 records · Page 31

Optimal gas cooling with minimal aero-opticaldistortion for next generation high-average-power lasers (Progress report, Year 2)

Next-generation lasers that combine high-peak-power ultrashort pulses with high repetition rates have the potential to efficiently deliver the high-average-power outputs required for inertial fusion energy generation, laser-based particle acceleration, and to generate compact secondary x- and 𝛾-ray radiation sources, among others. The combination of long gain length with high surface-to-volume ratio to remove residual heat makes gas-cooled multi-slab amplifier heads an inherently aperture- and average power scalable technology that enables this new class of solid-state lasers, and their thermal management becomes the most pressing technical challenge. As the index-of refraction depends on the gas density, heat removal from the slab surfaces unavoidably entails scattering of propagated light from regions of high temperature fluctuations. The requirement to maximize convective heat extraction from the gain medium under the constraint of minimum optical aberration of the beam renders amplifier slab cooling a technological challenge at the intersection of aero-optics, fluid-mechanics and heat transfer. The goal of this research is to leverage computational fluid dynamics and aero optics simulations of different levels of fidelity to devise a numerical modeling and optimization framework for gas-cooled multi-slab amplifiers. The project will be conducted in close collaboration with the Advanced Photon Technologies Group at Lawrence Livermore National Laboratory, where a gas-cooled test stand for laser materials and thermal management structures will be developed and installed. The final goal of the collaborative effort is to manufacture and benchmark a physical prototype of an numerically optimized amplifier head design. Measurement data on solid-state and thermo-optical distortions will be used for the validation of the numerical models. The key technical challenge in this multi-objective design problem is that of balancing optimal heat removal with minimal optical distortion. To facilitate computationally efficient design optimization and, at the same time, deepen our basic understanding of the interaction of light with turbulence, a complementary study that combines high and low-fidelity flow simulations and aero-optical models is proposed. Along the lines of these two objectives and solution strategies, the project is organized into two thrusts.

42 ENGINEERING↗

A Quantitative Approach for Demonstrating Plume Stabilization Under CCS Policy Frameworks

Conference presentation at Carbon Capture, Utilization, and Storage (CCUS) Conference 2024, Houston, Texas, March 11–13, 2024. The study highlighted the operator’s responsibilities for demonstrating plume stabilization in partial fulfillment of postinjection monitoring and site closure requirements. The metric dA/dt was chosen because of rules and regulations that call out the need to measure the change in area of the plume and plume extents over time. As dA/dt approached dA/dtcritical, the CO 2 plume’s growth was both minimal and predictable. Providing operators with a quantitative approach to determining plume stabilization is important for moving carbon capture and storage (CCS) projects beyond feasibility phases of development.

01 COAL, LIGNITE, AND PEAT↗

The BTSbot-nearby Discovery of SN 2024jlf: Rapid, Autonomous Follow-up Probes Interaction in an 18.5 Mpc Type IIP Supernova

We present observations of the Type IIP supernova (SN) SN 2024jlf, including spectroscopy beginning just 0.7 days (∼17 hr) after first light. Rapid follow-up was enabled by the new BTSbot-nearby program, which involves autonomously triggering target-of-opportunity requests for new transients in Zwicky Transient Facility data that are coincident with nearby (D < 60 Mpc) galaxies and identified by the BTSbot machine learning model. Early photometry and nondetections shortly prior to first light show that SN 2024jlf initially brightened by >4 mag day −1 , quicker than ∼90% of Type II SNe. Early spectra reveal weak flash ionization features: narrow, short-lived (1.3 < τ[days] < 1.8) emission lines of Hα, He II , and C IV . Assuming a wind velocity of v w = 50 km s −1 , these properties indicate that the red supergiant progenitor exhibited enhanced mass loss in the last year before explosion. We constrain the mass-loss rate to $1{0}^{-4}\lt \dot{M}\,[{M}_{\odot }\,{\mathrm{yr}}^{-1}]\lt 1{0}^{-3}$ by matching observations to model grids from two independent radiative hydrodynamics codes. BTSbot-nearby automation minimizes spectroscopic follow-up latency, enabling the observation of ephemeral early-time phenomena exhibited by transients.

core-collapse supernovae↗

Safe Operations at Roadway Junctions: Intelligent Roadway Infrastructure as Functional Interlocking

Automated vehicle (AV) technology is quickly maturing, and the corresponding infrastructure systems that evaluate traffic and communicate to vehicles requires sophisticated sensing and perception technologies, referred to as intelligent roadway infrastructure (IRI), to complement emerging AV capabilities. IRI provides signals to vehicles, indicating right-of-way for vehicles and communicating to approaching AVs that no other vehicle is failing to yield. This capability, denoted as safety-affirmative signaling, provides a green light or a green arrow as appropriate and affirms through communication links to connected vehicles when it is safe to proceed. About 36% of collisions occur at intersections, with most occurring upon left turns (22.2%) or crossing over (12.6%), and only a small percentage (1.2%) while turning right at an intersection. Of all intersection crashes about half (52.5%) of those vehicles were traveling through a signalized intersection 2. Safety-affirmative signaling would guarantee safety of AV fleet vehicles, by providing the interlocking principle, a term from automated train control that only allows progression through a railway intersection after affirming no opportunity for a crash exists. IRI through safety-affirmative signaling would bring performance and safety to complex roadway intersections where AV transit fleet service is most needed, as well as safety benefits to traditional, non-automated vehicles and vulnerable road users. The implementation of IRI has functional, programmatic, and technical challenges. Research work performed at the National Renewable Energy Laboratory (NREL) in an integrative approach encapsulating these themes, and termed infrastructure perception and control (IPC) is motivated by improved performance (travel time), improved safety (reduced collisions), and improved energy efficiency (less fuel burned and minimized production of greenhouse gases). IPC is intended not only for roadway and intersection applications but also in extension to inform complementary buildings and grid systems to enable better co-management, as vehicles and their charging needs become increasingly integrated into the built environment. The NREL IPC project presents an open-source framework, architecture, and supporting technology to implement IRI, addressing critical issues such as fusion of data, reliability, standardization of data interfaces, and confidence of detection. The framework is informed by previous experience in U.S. Department of Defense research technology, specifically in the use of radar to detect, identify, and track aerial threats. These principles combined with multi-sensor fusion provides for a complete digital twin with known and measurable confidence and accuracy from which safety-affirmative signaling can be developed and deployed.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗

Spent-fuel recycling molecule degradation product standards can be produced with preparative liquid chromatography

Solvent systems used for spent nuclear fuel separation schemes are subject to intense radiation fields from the spent fuel, resulting in a degradation of separation performance over time from both the radiolytic destruction of active compounds (extractants, holdback agents, etc.) and the production of degradation compounds (compounds produced from the destruction of the active compounds). Development of a quantitative understanding of the effects of radiolysis on separations is impeded by a near-complete lack of availability of degradation products from chemical manufacturers. Synthesis of these degradation products through traditional means is both cost and time prohibitive, inhibiting advancement of new fuel separation processes. We hypothesize that comprehensive sets of degradation products can be synthesized using preparative radiolysis, and purified using high performance liquid chromatography in quantities sufficient for unequivocal identification and subsequent use as quantitative analytical studies. Achieving this objective will facilitate development of a much more accurate understanding of radiolysis reactions and their effects on solvent extraction processes, positioning Idaho National Laboratory (INL) in a leadership role in the science of fuel cycle separations. We propose producing these degradation compounds by irradiating concentrated samples of solvent ligands in the INL gamma irradiator, separating the resulting mixture with preparative chromatography, and removing the mobile phase in the separated fractions with a vacuum concentrator. This will yield milligram to gram quantities of purified degradation products, sufficient for fundamental research, at a fraction of the time and cost of synthesis. Facile production of separation molecule degradation products at minimal cost will enable the fundamental research that will accelerate the development of fuel cycles that will ensure the competitiveness of nuclear energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Error-controlled, progressive, and adaptable retrieval of scientific data with multilevel decomposition

Extreme-scale simulations and high-resolution instruments have been generating an increasing amount of data, which poses significant challenges to not only data storage during the run, but also post-processing where data will be repeatedly retrieved and analyzed for a long period of time. The challenges in satisfying a wide range of post-hoc analysis needs while minimizing the I/O overhead caused by inappropriate and/or excessive data retrieval should never be left unmanaged. In this paper, we propose a data refactoring, compressing, and retrieval framework capable of 1) fine-grained data refactoring with regard to precision; 2) incrementally retrieving and recomposing the data in terms of various error bounds; and 3) adaptively retrieving data in multi-precision and multi-resolution with respect to different analysis. With the progressive data re-composition and the adaptable retrieval algorithms, our framework significantly reduces the amount of data retrieved when multiple incremental precision are requested and/or the downstream analysis time when coarse resolution is used. Experiments show that the amount of data retrieved under the same progressively requested error bound using our framework is 64% less than that using state-of-the-art single-error-bounded approaches. Parallel experiments with up to 1, 024 cores and ~ 600 GB data in total show that our approach yields 1.36× and 2.52× performance over existing approaches in writing to and reading from persistent storage systems, respectively.

Liang, Xin↗

Image Deconvolution and Point-spread Function Reconstruction with STARRED: A Wavelet-based Two-channel Method Optimized for Light-curve Extraction

We present starred, a point-spread function (PSF) reconstruction, two-channel deconvolution, and light-curve extraction method designed for high-precision photometric measurements in imaging time series. An improved resolution of the data is targeted rather than an infinite one, thereby minimizing deconvolution artifacts. In addition, starred performs a joint deconvolution of all available data, accounting for epoch-to-epoch variations of the PSF and decomposing the resulting deconvolved image into a point source and an extended source channel. The output is a high-signal-to-noise-ratio, high-resolution frame combining all data and the photometry of all point sources in the field of view as a function of time. Of note, starred also provides exquisite PSF models for each data frame. We showcase three applications of starred in the context of the imminent LSST survey and of JWST imaging: (i) the extraction of supernovae light curves and the scene representation of their host galaxy; (ii) the extraction of lensed quasar light curves for time-delay cosmography; and (iii) the measurement of the spectral energy distribution of globular clusters in the "Sparkler," a galaxy at redshift z = 1.378 strongly lensed by the galaxy cluster SMACS J0723.3-7327. starred is implemented in jax, leveraging automatic differentiation and graphics processing unit acceleration. This enables the rapid processing of large time-domain data sets, positioning the method as a powerful tool for extracting light curves from the multitude of lensed or unlensed variable and transient objects in the Rubin-LSST data, even when blended with intervening objects.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing Dynamic Structure in Battery Electrodes by Time-Resolved Cryo-TEM

In recent years, cryogenic transmission electron microscopy (cryo-TEM) has enabled high-resolution characterization of sensitive battery materials by minimizing electron beam-induced artifacts and damage. Success of this technique relies on the preparation of thin, rapidly frozen samples, generally by disassembling batteries under inert atmosphere, transferring materials of interest to a TEM grid, and finally plunge freezing into a cryogen. In material degradation studies, this extensive time between electrochemical cycling and cryo-TEM characterization leaves room for structural relaxation, diffusion, and other dynamic processes that make it difficult to precisely correlate the imaged structure with the native structure that evolves during battery cycling or aging. Here, we present a method to integrate battery cycling with fast preparation of electrode samples for cryo-TEM. This enables higher fidelity between the structures characterized and the electrochemical state of interest, which we use to study deformation in silicon nanoparticle anodes for lithium-ion batteries.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗

Just-in-Time Compilation and Link-Time Optimization for OpenMP Target Offloading

Following the mass adoption of external accelerators for high performance computing, the overall performance of many applications has become increasingly dependent on relatively small accelerated kernels. As static analysis is fundamentally limited by dynamic values and external definitions, standard ahead-of-time compilation is not always sufficient to achieve the best performance. Furthermore, many users looking to port an existing application to run on an external accelerator will not want to fundamentally restructure their programs. These and other problems can be addressed through both link-time optimization (LTO) and just-in-time (JIT) compilation, but until now had sparse and inconsistent support from the compiler. In this work, we present a new compilation method that enables device-side LTO as well as a transparent JIT compilation tool-chain for OpenMP target offloading. Our contributions include an entirely new device linking and embedding scheme to enable LTO as well as a novel JIT engine to efficiently optimize OpenMP offloading regions at run-time. We also introduce a persistent caching system to improve end-to-end runtime using the JIT engine and minimize kernel launching overheads. We measure the performance of our LTO and JIT implementation via several real-world scientific applications. With our optimizations we observe significant improvements through LTO on large applications as well as significant end-to-end execution time improvement using JIT.

Tian, Shilei↗

Stochastic equation of motion approach to fermionic dissipative dynamics. I. Formalism

In this work, we establish formally exact stochastic equation of motion (SEOM) theory to describe the dissipative dynamics of fermionic open systems. The construction of the SEOM is based on a stochastic decoupling of the dissipative interaction between the system and fermionic environment, and the influence of environmental fluctuations on the reduced system dynamics is characterized by stochastic Grassmann fields. Meanwhile, numerical realization of the time-dependent Grassmann fields has remained a long-standing challenge. To solve this problem, we propose a minimal auxiliary space (MAS) mapping scheme with which the stochastic Grassmann fields are represented by conventional c-number fields along with a set of pseudo-levels. This eventually leads to a numerically feasible MAS-SEOM method. The important properties of the MAS-SEOM are analyzed by making connection to the well-established time-dependent perturbation theory and the hierarchical equations of motion theory. The MAS-SEOM method provides a potentially promising approach for the accurate and efficient simulation of fermionic open systems at ultra-low temperatures.

Han, Lu (ORCID:0000000339002225)↗

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗

A high power, low temperature molten sodium battery

Low temperature molten sodium batteries promise low-cost, grid-scale energy storage using earth-abundant materials. To be economically viable, they must demonstrate high current (power) at relevant discharge times. Here a low temperature (135 °C) molten sodium battery with a NaI–AlCl 3 molten salt catholyte and NaSICON separator is explored, minimizing cell ohmic resistances and avoiding precipitation reactions in the molten salt catholyte. Operating currents were increased by as much as 100× compared to the energy-dense baseline design (197 Wh kg −1 and 248 mAh cm −2 ). Cells were cycled 100 times at 50 mA cm −2 (150 mW cm −2 ), averaging 99.8% coulombic efficiency and 79.7% energy efficiency for 30% of the theoretical capacity. Charging currents up to 250 mA cm −2 (1070 mW cm −2 ) are demonstrated for 10% of the theoretical capacity. Toward extending the accessible capacity in these high current systems, the volume ratio of catholyte molten salt to a graphite felt current collector was varied, altering (1) the local current density in the graphite felt and (2) the discharge time at a fixed current density. Optimizing the catholyte : felt ratio dramatically increased capacity utilization to 60% (102 mAh cm −2 ) at 50 mA cm −2 (150 mW cm −2 ) charge and 10 mA cm −2 discharge. This catholyte : felt ratio of 1.4 significantly decreased the cell impedance, minimizing electrode blocking effects seen previously in NaI–AlCl 3 catholytes. Together, these high areal loadings and impressive power outputs demonstrate how low temperature molten sodium batteries can compete with commercialized sodium systems operating at more than twice the temperature.

NaSICON↗

Optimal sensor network design for multi-scale, time-varying differential algebraic equation systems: Application to an entrained-flow gasifier refractory brick

An algorithm for optimal sensor network design for multi-scale, time-varying differential algebraic equation systems with non-separable dynamics is presented. As the process is time-varying, an integral normalized posterior error covariance of a multi-scale filter is minimized to obtain the optimal sensor locations. For reducing the computational cost, an adaptive sampling rate approach is considered for the slowly-varying variables. The algorithm is applied to a smart refractory brick with embedded sensors as part of an entrained-flow gasifier. Thermistors and interdigital capacitors are considered as candidate measurement technologies for estimating temperature and slag penetration profile along the gasifier wall. Finally, when the optimal set of sensors obtained from the algorithm is used for estimating temperature and slag penetration profiles in a multi-scale Kalman filter framework, satisfactory estimates are obtained despite high measurement noise and model mismatch.

42 ENGINEERING↗

A Workflow for Accelerating Multimodal Data Collection for Electrodeposited Films

Abstract Future machine learning strategies for materials process optimization will likely replace human capital-intensive artisan research with autonomous and/or accelerated approaches. Such automation enables accelerated multimodal characterization that simultaneously minimizes human errors, lowers costs, enhances statistical sampling, and allows scientists to allocate their time to critical thinking instead of repetitive manual tasks. Previous acceleration efforts to synthesize and evaluate materials have often employed elaborate robotic self-driving laboratories or used specialized strategies that are difficult to generalize. Herein we describe an implemented workflow for accelerating the multimodal characterization of a combinatorial set of 915 electroplated Ni and Ni–Fe thin films resulting in a data cube with over 160,000 individual data files. Our acceleration strategies do not require manufacturing-scale resources and are thus amenable to typical materials research facilities in academic, government, or commercial laboratories. The workflow demonstrated the acceleration of six characterization modalities: optical microscopy, laser profilometry, X-ray diffraction, X-ray fluorescence, nanoindentation, and tribological (friction and wear) testing, each with speedup factors ranging from 13–46x. In addition, automated data upload to a repository using FAIR data principles was accelerated by 64x.

36 MATERIALS SCIENCE↗

The superior hydrothermal stability of Pd/SSZ-39 in low temperature passive NOx adsorption (PNA) and methane combustion

We successfully synthesized uniform SSZ-39 with an average crystal size of about a micron. Pd (0.7 - 3 wt%) was supported on SSZ-39 with Si/Al ratio ~12. The as-synthesized materials were characterized by FTIR, XRD, Helium Ion Microscopy, HAADF-STEM imaging, 27 Al, 29 Si and H solid state NMR spectroscopic techniques. FTIR studies with CO and NO probe molecules reveal that the 0.7 wt% Pd/SSZ-39 material with Si/Al ~12 has the majority of Pd dispersed atomically as isolated Pd(II) and Pd(II)-OH centers, and thus can be used as a low-temperature passive NOx adsorber. Pd(II)-NO, Pd(II)(OH)(NO) and Pd(II)(CO)(NO) complexes form during PNA in this material. We compare this PNA material directly with the Pd/SSZ-13 system (with Si/Al ratio ~12) and show its superior hydrothermal stability. Remarkably, Pd/SSZ-39 with Si/Al ratio ~12 survives hydrothermal aging up to 815 ºC in 10% H 2 O/Air vapor for 16 hours without significant loss in activity. The SSZ-39 crystal structure remains intact during hydrothermal aging up to 1,000 ºC as we elucidate it with XRD and HAADF-STEM imaging/EDS mapping. However, changes to the framework during such harsh hydrothermal treatment significantly change the NOx release profiles during PNA as evidenced by high-field 27 Al NMR on fresh and aged Pd/SSZ-39 samples as well as PNA performance measurements. Besides PNA application, these hydrothermally very stable materials (3 wt% Pd on SSZ-39 with Si/Al ratio ~12) can be used as a robust methane combustion catalyst under industrially relevant conditions (GHSV~600,000hr -1 ). This catalyst shows minimal deactivation after both harsh hydrothermal aging at 750 and 800 ºC, and prolonged time on stream (105 hrs) at 425 °C. In contrast, both 3wt% Pd/alumina and 3wt% SSZ-13 supported samples lose a significant portion of their activity. Our study opens new avenues to prepare the most hydrothermally stable known Pd/zeolite materials with applications for adsorption and catalytic hydrocarbon combustion under industrially relevant conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

System risk quantification and decision making support using functional modeling and dynamic Bayesian network

Risk-informed decision-making requires a probabilistic assessment of the likelihood of success of control action, given the system status. This paper presents a systematic state transition modeling approach integrating dynamic probabilistic risk assessment with a decision-making process using a dynamic Bayesian network (DBN) coupled with functional modeling. A functional model designed with multilevel flow modeling (MFM) technique was used to build a system state structure inferred by energy, mass, and information flow so that one can verify the developed model with respect to system functionality. The MFM model represents the causal relationship among the nodes, which captures the structure of process parameters and control units. Each node may have multiple possible states, and the DBN structured by the MFM model represents the time-domain transitions among the defined states. Furthermore, the MFM-DBN integrated state transition modeling is a white-box approach that allows one to draw the system's risk profile by updating the system states and supports the decisions probabilistically with physical inference. An example of a simple heating system has been used to illustrate this process, including decision-making support based on quantitative risk profile. For demonstrating its applicability to a complex system operational decision making, a case study of station blackout accident scenario leading to the seal loss of coolant accident in a nuclear power plant is presented. The proposed approach effectively provided the risk profile along time for each option so that the operators can make the best decision, which minimizes the plant risk.

42 ENGINEERING↗

Single pair of Weyl nodes in the spin-canted structure of EuCd 2 As 2

We report time reversal symmetry breaking Weyl semimetals are unique among Weyl materials in allowing the minimal number of Weyl points, thus offering the clearest signatures of the associated physics. Here we present neutron diffraction, density-functional theory, and transport measurement results which indicate that EuCd 2 As 2 , under ambient field, strain, and pressure, is such a material with a single pair of Weyl points. Our work reveals a magnetic structure (magnetic space group C2'/m') with Eu moments pointing along the [210] direction in plane and canted ~ 30° out of plane. Density functional theory calculations using this structure show that the observed canting drastically alters the relevant electronic bands, relative to the in-plane order, leading to a single set of well-defined Weyl points. Furthermore, we find the canting angle can tune the distance of the Weyl points above the Fermi level, with the smallest distance at low canting angles. Finally, transport measurements of the anomalous Hall effect and longitudinal magnetoresistance exhibit properties indicative of a chiral anomaly, thus supporting the neutron scattering and DFT results suggesting EuCd 2 As 2 is close to the ideal situation of the Weyl hydrogen atom.

36 MATERIALS SCIENCE↗

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗