Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Phase Field Model”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Hybrid models for the simulation of microstructural evolution influenced by coupled, multiple physical processes

Most materials microstructural evolution processes progress with multiple processes occurring simultaneously. In this work, we have concentrated on the processes that are active in nuclear materials, in particular, nuclear fuels. These processes are coarsening, nucleation, differential diffusion, phase transformation, radiation-induced defect formation and swelling, often with temperature gradients present. All these couple and contribute to evolution that is unique to nuclear fuels and materials. Hybrid model that combines elements from the Potts Monte Carlo, phase-field models and others have been developed to address these multiple physical processes. These models are described and applied to several processes in this report. An important feature of the models developed are that they are coded as applications within SPPARKS, a Sandiadeveloped framework for simulation at the mesoscale of microstructural evolution processes by kinetic Monte Carlo methods. This makes these codes readily accessible and adaptable for future applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Multi-scale simulation of high burnup UO2 nuclear fuel during loss-of-coolant accident conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized that pulverization is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted, and the results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Diffuse interface approach to modeling crystal plasticity with accommodation of grain boundary sliding

In diffuse-interface (phase field) models of polycrystalline solids, the grain boundaries naturally possess a narrow region of finite thickness. A Ginzburg-Landau type kinetic equation for plastic deformation of a solid is derived that combines crystal plasticity and J 2 plasticity. Within such a framework for a polycrystal the grain interior can be defined by crystal plasticity, and the grain boundary region can be assigned J 2 type plasticity. The two are connected smoothly. This enables the modeling framework to accommodate grain boundary sliding (GBS), an important deformation mechanism for creep. The relaxation of elastic modulus of a polycrystalline solid and the power-law creep (both accommodated by GBS) are reported on by 2D and 3D simulations, respectively, and the results are compared to the 2D finite-element simulations of Ghahremani and Crossman-Ashby. A strong grain shape dependence and orientation dependence (for non-equiaxed grains) of the GBS effect are predicted by this model.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Structural deformation and transformation of θ'-Al 2 Cu precipitate in Al matrix via interfacial diffusion

Metastable θ'-Al 2 Cu precipitates enhance mechanical properties of Al alloy, but gradually transform to θ phase due to diffusion related coarsening at higher temperatures (>200°C). To improve understanding of the phase instability and transformation of θ'-Al 2 Cu, we investigated the interfacial atomic mobility of θ'-Al2Cu precipitate in Al matrix via molecular dynamics (MD). To characterize the interfacial atomic mobility, diffusion properties (activation energy, diffusivity, and jump attempt frequency) with respect to atomic species, interfacial structure, and temperature are calculated using atomic trajectories from MD. For the enhanced accuracy of this analysis, especially at low temperatures, temporal scale of molecular dynamics is significantly extended by employing the parallel replica dynamics. This study, as the first MD investigation on atomic diffusion of θ'-Al 2 Cu in Al matrix, suggests and discusses (i) structural deformation activated at the semi-coherent θ'-Al 2 Cu/Al interfaces, (ii) transformation initiated at the edge of the coherent θ'-Al 2 Cu/Al interface, and (iii) directional dependence of diffusion properties. Overall, these findings are expected to be employed in a larger scale phase-field modeling of precipitate growth, and to contribute to the understanding of phase transformation process and the development of alloys with θ'-phase retained to higher temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A thermodynamically consistent finite strain phase field approach to ductile fracture considering multi-axial stress states

Phase field models for ductile fracture have gained significant attention in the last two decades due to their ability in implicitly tracking the nucleation and propagation of cracks. However, most crack phase field formulations for elastoplastic solids focus only on the effects of plastic deformation, and do not consider the different multi-axial stress states that may arise in practical designs. In this work, a thermodynamically consistent phase field approach coupled with finite strain plasticity, considering multi-axial stress states is presented. In order to account for the coupling between plasticity and stress states, the Stress-Weighted Ductile Fracture Model (SWDFM) is utilized. The SWDFM represents a criterion for predicting ductile crack initiation under both monotonic and cyclic loadings based on histories of an internal plastic variable, stress triaxiality, and the Lode angle parameter. The excellent performance of the SWDFM for predicting ductile crack initiation motivates for its incorporation into a phase field approach for predicting both crack initiation and propagation through degradation of the fracture toughness. Moreover, based on the second law of thermodynamics, exact requirements are imposed on the rate at which the fracture toughness can evolve. Here, a novel function for degrading the plastic yield surface during the evolution of damage is introduced. This function, in line with experimental observations, leads to an accumulation of plastic deformation in damaged regions of a solid, and avoids numerical instabilities arising from concentrations of large plastic deformations in severely damaged regions. For validating the proposed model, results of computational simulations are compared to data from selected tests considering different multi-axial stress states. Comparisons of the numerical results with data from laboratory experiments demonstrate the capabilities of the proposed framework.

42 ENGINEERING↗
DOE OSTI2016

Exploring Polarization Rotation Phase Instabilities in Super-Tetragonal BiFeO3 Epitaxial Thin Films and Their Technological Implications

Many functional properties of ferroelectrics are underlain by structural instabilities, which render these materials very susceptible to small alternating applied fields (electric, mechanical, etc.) through certain constitutive coupling relations, e.g., elastic compliance and piezoelectric response, and often such instabilities can be shifted by static applied fields thus meaning tunable dynamic properties. Structural instabilities are naturally accommodated on the brink of morphotropic phase boundaries (MPB s) where multiple phases of small energy difference coexist in different crystallographic forms. Canonical MPB is realized through compositional mixture, as is typically exemplified by Pb(Zr1-xTix)O3 solid solutions and relaxor ferroelectrics of (1-x)PbMg1/3Nb2/3O3-xPbTiO3. More recently, a strain-driven MPB has been discovered in BiFeO3 (BFO) thin films epitaxially grown on LaAlO3 (LAO) crystal substrates (which imposes about -4.5% in-plane strains). Such an MPB is in between a rhombohedral (R) phase that bulk BFO exhibits and a so-called super-tetragonal (T) phase, which name hints at its giant lattice axial ratio (c/a ~ 1.23) and accordingly high electric polarization (~1.5 C m-2). The discovery of an MPB in BFO has revealed another facet of this multiferroic system, further adding opportunities to its many exotic functionalities such as domain wall conduction, magnetoelectric and photovoltaic effects As with other MPB s, large electric-field induced strains as well as more underlying lattice softening effects are observed near this MPB promising piezoelectric-based applications. In addition, T-phase BFO itself shows distinct properties, e.g., electronic band gap and optical absorption, from the R-phase and the resultant switching effects between them may also be exploitable. However, unlike conventional ferroelectric oxides where the phases across an MPB usually have subtle difference caused by rotations of an ion off-centering polarization, the BFO system bears multiple structural degrees of freedom, in particular antiferrodistortive modes of oxygen octahedral tilt, and a multitude of structural transition paths are thereby facilitated. Moreover, since the MPB of BFO is driven by epitaxial strain, it is sensitive to the strain relaxation related to film thickness and growth conditions, and above some critical thickness the films appear in the form of a hierarchical mixed-phase microstructure involving several coexistent polymorphs with distinct lattice distortion (monoclinic phases) and tilts. Elastic and electrical heterogeneities are necessary consequences of such complex microstructure, which couple to the intrinsic order parameters and expectedly have profound influence on the structural dynamics and material properties. All these make it a demanding task to obtain a deep understanding of this MPB system on par with its application prospect. The thin-film material form also brings in experimental restrictions to the MPB phase transition studies of BFO since many pertinent techniques fail to operate at the nanoscale or suffer from formidably weak signals. Therefore, static structural characterizations using e.g. X-ray diffraction, electron microscopy and Raman scattering have prevailed thus far. Recently, we adapted band-excitation piezoresponse spectroscopy (BEPS) to probe the bias-induced R/T phase transition dynamics of BFO/SrTiO3 (STO) thin films, and revealed the soft-mode elastic behavior of the transition. The efficacy of our method, along with standard piezoresponse force microscopy (PFM), can be significantly leveraged by integrating a variety of local and/or global excitations, e.g., tip pressure (10 s GPa attainable), heating, photo-irradiation, available to modern scanning probe platforms. This can thus afford unique opportunities to survey the structural dynamics of ferroelectric materials coupled to those field variables, enabling rapid discovery or validation of their functional properties. In this study, we focus on the T-phase BFO/LAO system; we have examined its phase transition behavior not only due to local coaction of tip bias and loading force, but also under device-level global electric fields in a coplanar capacitor. The intrinsic elastic softening phenomena near the structural transitions are comprehensively revealed by BEPS and corroborated by phase-field modeling. Our findings may open a new pathway for technological utilization of the MPB phase instabilities of BFO.

BiFeO3↗
DOE OSTI2016

High-Burnup-Structure (HBS): Model Development in MARMOT for HBS Formation and Stability Under Radiation and High Temperature

A detailed phase field model for the formation of High Burnup Structure (HBS) was developed and implemented in MARMOT. The model treats the HBS formation as an irradiation-induced recrystallization. The model takes into consideration the stored energy associated with dislocations formed under irradiation. The accumulation of radiation damage, hence, increases the system free energy and triggers recrystallization. The increase in the free energy due to the formation of new grain boundaries is offset by the reduction in the free energy by creating dislocation-free grains at the expense of the deformed grains. The model was first used to study the growth of recrystallized flat and circular grains. The model reults were shown to agree well with theorrtical predictions. The case of HBS formation in UO2 was then investigated. It was found that a threshold dislocation density of (or equivalently a threshold burn-up of 33-40 GWd/t) is required for HBS formation at 1200K, which is in good agrrement with theory and experiments. In future studies, the presence of gas bubbles and their effect on the formation and evolution of HBS will be considered.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Creating physically-based three-dimensional microstructures: Bridging phase-field and crystal plasticity models.

In order to better incorporate microstructures in continuum scale models, we use a novel finite element (FE) meshing technique to generate three-dimensional polycrystalline aggregates from a phase field grain growth model of grain microstructures. The proposed meshing technique creates hexahedral FE meshes that capture smooth interfaces between adjacent grains. Three dimensional realizations of grain microstructures from the phase field model are used in crystal plasticity-finite element (CP-FE) simulations of polycrystalline a -iron. We show that the interface conformal meshes significantly reduce artificial stress localizations in voxelated meshes that exhibit the so-called "wedding cake" interfaces. This framework provides a direct link between two mesoscale models - phase field and crystal plasticity - and for the first time allows mechanics simulations of polycrystalline materials using three-dimensional hexahedral finite element meshes with realistic topological features.

36 MATERIALS SCIENCE↗
DOE OSTI2016

On the kinetics of dendritic sidebranching: A three dimensional phase field study

The underlying mechanism for dendritic sidebranching was studied using 3-D phase field modeling. Results showed that in 3-D the requirement of applying the random thermal noise to induce dendritic sidebranching (i.e., normally the case for 2-D phase field simulations) was fully relaxed. The stretching of the secondary or higher order arms occurred spontaneously and symmetrically as the growth of the dendrite. With periodic external perturbation and if the stimulating frequency was lower than a critical value, both tip velocity and sidebranching would get completely synchronized with the perturbation. Whereas if the perturbation frequency was higher than the critical value, rather than increasing, the sidebranching frequency would become stable and maintain at the same magnitude as that of the natural sidebranching, i.e., when no external perturbation was applied. It was shown that the underlying mechanism for sidebranching was deterministic rather than stochastic, and anisotropy tendency and curvature effect were shown to be the most important influence factors. Moreover, the difference of the anisotropy tendency would lead to an uneven distribution of curvature on the solid/liquid interface, i.e., formation of concave and convex geometries. The growth of these geometries would subsequently break the initial spherical structure of solid seed and lead to further sidebranching.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2020

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing the Critical Nucleus Size in the Metal-Insulator Phase Transition of VO 2

In a first-order phase transition, critical nucleus size governs nucleation kinetics, but the direct experimental test of the theory and determination of the critical nucleation size have been achieved only recently in the case of ice formation in supercooled water. The widely known metal-insulator phase transition (MIT) in strongly correlated VO 2 is a first-order electronic phase transition coupled with a solid-solid structural transformation. It is unclear whether classical nucleation theory applies in such a complex case. Here we directly measure the critical nucleus size of the MIT by introducing size-controlled nanoscale nucleation seeds with focused ion irradiation at the surface of a deeply supercooled metal phase of VO 2 . The results compare favorably with classical nucleation theory and are further explained by phase-field modeling. This Letter validates the application of classical nucleation theory as a parametrizable model to describe phase transitions of strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Microstructure, Thermal, and Mechanical Properties Relationships in U and UZr Alloys (Final Report)

Uranium-zirconium (U-Zr) alloys are candidate fuel systems for transmutation based reactors that can be used to burn long-lived minor actinides and fission products in fast spectrum reactors. Metallic fuels have also been gaining more recent attention for applications as accident tolerant fuels, as well as for use in small modular reactors. This research focused on a “science-based” approach to capture the connections between U and U-Zr alloys’ three-dimensional (3-D) microstructure, thermal properties, and mechanical properties through closely coordinated experiments and modeling efforts from the unirradiated to the irradiated fuels. Advanced characterization and modeling techniques were used to understand irradiation-induced microstructural evolution and its direct impact on the thermal and mechanical properties of U and U-Zr fuel. Closely coordinated experiments and modeling were performed to provide crucial data that does not currently exist. Overall, this research spanned multiple length and time scales within the models and experiments. The scope of the research encompassed the understanding of the irradiation effects in U and various U-Zr alloys with particular attention paid to three task areas: (1) microstructural evolution, (2) in-situ/ex-situ thermal and mechanical properties, and (3) multiscale modeling connections to microstructure, thermal, and mechanical properties. This research resulted in (1) the 3-D characterization of neutron irradiated U-Zr fuel in multiple phase regions to better understand fission gas swelling and constituent redistribution, (2) development of a microstructural model linking thermal and mechanical properties via in situ Raman and nanoindentation, (3) and mesoscale phase field modeling was coupled with the AEH method in the MOOSE framework was used to calculate the effective thermal conductivities of U-Zr fuels consisting of α-U and δ-UZr 2 heterogeneous microstructures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

True random number generation using the spin crossover in LaCoO 3

While digital computers rely on software-generated pseudo-random number generators, hardware-based true random number generators (TRNGs), which employ the natural physics of the underlying hardware, provide true stochasticity, and power and area efficiency. Research into TRNGs has extensively relied on the unpredictability in phase transitions, but such phase transitions are difficult to control given their often abrupt and narrow parameter ranges (e.g., occurring in a small temperature window). Here we demonstrate a TRNG based on self-oscillations in LaCoO 3 that is electrically biased within its spin crossover regime. The LaCoO 3 TRNG passes all standard tests of true stochasticity and uses only half the number of components compared to prior TRNGs. Assisted by phase field modeling, we show how spin crossovers are fundamentally better in producing true stochasticity compared to traditional phase transitions. As a validation, by probabilistically solving the NP-hard max-cut problem in a memristor crossbar array using our TRNG as a source of the required stochasticity, we demonstrate solution quality exceeding that using software-generated randomness.

97 MATHEMATICS AND COMPUTING↗