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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Spontaneous magnetization of collisionless plasma

Here we study within a fully kinetic framework the generation of “seed” magnetic fields through the Weibel instability, driven in an initially unmagnetized plasma by a large-scale shear force. We develop an analytical model that describes the development of thermal pressure anisotropy via phase mixing, the ensuing exponential growth of magnetic fields in the linear Weibel stage, and the saturation of the Weibel instability when the seed magnetic fields become strong enough to instigate gyromotion of particles and thereby inhibit their free-streaming. The predicted scaling dependencies of the saturated fields on key parameters (e.g., ratio of system scale to electron skin depth and forcing amplitude) are confirmed by two-dimensional and three-dimensional particle-in-cell simulations of an electron–positron plasma. This work demonstrates the spontaneous magnetization of a collisionless plasma through large-scale motions as simple as a shear flow and therefore has important implications for magnetogenesis in dilute astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Uncertainty Quantification in Atomistic Modeling of Metals and Its Effect on Mesoscale and Continuum Modeling: A Review

The design of next-generation alloys through the integrated computational materials engineering (ICME) approach relies on multiscale computer simulations to provide thermodynamic properties when experiments are difficult to conduct. Atomistic methods such as density functional theory (DFT) and molecular dynamics (MD) have been successful in predicting properties of never before studied compounds or phases. However, uncertainty quantification (UQ) of DFT and MD results is rarely reported due to computational and UQ methodology challenges. Over the past decade, studies that mitigate this gap have emerged. These advances are reviewed in the context of thermodynamic modeling and information exchange with mesoscale methods such as the phase-field method (PFM) and calculation of phase diagrams (CALPHAD). The importance of UQ is illustrated using properties of metals, with aluminum as an example, and highlighting deterministic, frequentist, and Bayesian methodologies. Finally, challenges facing routine uncertainty quantification and an outlook on addressing them are also presented.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Metallothionein expression in Eastern mosquitofish ( Gambusia holbrooki ) exposed to 137 Cs and lead: implications for using metallothionein as a metal exposure biomarker

Given the range of threats facing aquatic ecosystems, it is critical we have monitoring tools to quickly and efficiently identify stressors within ecosystems. Metallothionein’s (MTs) are a group of proteins widely used as biomarkers of metal exposure; however, their suitability has been questioned due to their influenceable nature and ability to be induced by various stressors. Through our study, we sought to (1) evaluate induction of MT in eastern mosquitofish (Gambusia holbrooki) exposed to a radionuclide ( 137 Cs) and metal (Pb), (2) assess the influence of in situ 137 Cs pre-exposure on Pb-induced MT expression, and (3) quantify changes in MT concentrations during Pb depuration. We conducted our experiment in three phases using mosquitofish as a model species: Phase I—pre-exposure to 137 Cs in the field, Phase II—exposure to Pb in the lab immediately following 137 Cs exposure, and Phase III—depuration in the laboratory. We quantified MT concentrations over time and assessed differences among treatments. We did not find in situ exposure of mosquitofish to 137 Cs to induce MT, nor did pre-exposure to 137 Cs appear to influence induction of MT following Pb exposure. Metallothionein was rapidly induced in mosquitofish following Pb exposure; however, MT induction was also observed in control fish. Additionally, there was no apparent relationship between MT and Pb body burdens. During depuration, we found MT concentrations to be highly variable among fish demonstrating no clear response. Furthermore, our study demonstrates the high variability exhibited by MT, the need for further validation of MT as a metal exposure biomarker, and the value in additional research on these proteins under varying exposure scenarios, particularly environmentally realistic concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Transport studies in a gate-tunable three-terminal Josephson junction

Josephson junctions with three or more superconducting leads have been predicted to exhibit topological effects in the presence of few conducting modes within the interstitial normal material. Such behavior, of relevance for topologically protected quantum bits, would lead to specific transport features measured between terminals, with topological phase transitions occurring as a function of phase and voltage bias. Although conventional, two-terminal Josephson junctions have been studied extensively, multiterminal devices have received relatively little attention to date. Motivated in part by the possibility to ultimately observe topological phenomena in multiterminal Josephson devices, as well as their potential for coupling gatemon qubits, here we describe the superconducting features of a top-gated mesoscopic three-terminal Josephson device. The device is based on an InAs two-dimensional electron gas proximitized by epitaxial aluminum. We map out the transport properties of the device as a function of bias currents, top gate voltage, and magnetic field. Finally, we find a very good agreement between the zero-field experimental phase diagram and a resistively and capacitively shunted junction computational model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Deploying a Model Predictive Traffic Signal Control Algorithm - A Field Deployment Experiment Case Study

This paper presents a field deployment experiment of a real-time traffic signal control algorithm. We implemented the model predictive control (MPC) algorithm based on the virtual phase-link (VPL) model. We selected the deployment locations and times based on an energy saving potential concept. We developed a set of experiment systems, which included sensing, processing, and actuating components, to enable field deployment. We tested the systems rigorously before the experiment days. We reported the key procedures on the experiment days, including the steps taken, the real-time control procedure, and the monitoring of the experiment. We evaluated the impact of the deployment by looking at the changes in delay and energy consumption.

deployment↗
DOE OSTI2024

Anisotropic Thermal Conductivity Oscillations in Relation to the Putative Kitaev Spin Liquid Phase of α-RuCl 3

In the presence of an external magnetic field, the Kitaev model could host either gapped topological anyons or gapless Majorana fermions. In α-RuCl 3 , the gapped and gapless cases are only separated by a 30° rotation of the in-plane magnetic field vector. The presence or absence of the spectral gap is key for understanding the thermal transport behavior in α-RuCl 3 . Here, we study the anisotropy of the oscillatory features of thermal conductivity in α-RuCl 3 . We examine the oscillatory features of thermal conductivities (κ / / a, κ / / b) with fixed external fields and find distinct behavior for the gapped (B / / a) and gapless (B / / b) scenarios. Furthermore, we track the evolution of thermal resistivity (λ α ) and its oscillatory features with the rotation of in-plane magnetic fields from B / / b to B / / a. The thermal resistivity λ(B, φ) displays distinct rotational symmetries before and after the emergence of the field-induced Kitaev spin liquid phase. These results suggest that oscillatory features of thermal conductivity in α-RuCl 3 are closely linked to the putative Kitaev spin liquid phase and its excitations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Solving nuclear structure problems with the adaptive variational quantum algorithm

Here we use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the sd and pf shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

The cosmological moduli problem and naturalness

Nowadays, the cosmological moduli problem (CMP) comes in three parts: 1. potential violation of Big-Bang nucleosynthesis (BBN) constraints from late decaying moduli fields, 2. the moduli-induced gravitino problem wherein gravitinos are overproduced and their decays violate BBN or dark matter overproduction bounds and 3. the moduli-induced lightest SUSY particle (LSP) overproduction problem. Also, the CMP may be regarded as either a problem or else a solution to scenarios with dark matter over- or under-production. We examine the cosmological moduli problem and its connection to electroweak naturalness. We calculate the various two-body decay widths of a light modulus field into MSSM particles and gravitinos within general supersymmetric models. We include both phase space and mixing effects. We examine cases without and with helicity suppression of modulus decays to gravitinos (cases 1 & 2) and/or gauginos (cases A & B). For case B1, we evaluate regions of gravitino mass m 3/2 vs. modulus mass mΦ parameter space constrained by BBN, by overproduction of gravitinos and by overproduction of neutralino dark matter, along with connections to naturalness. For this case, essentially all of parameter space is excluded unless mΦ ≳ 2.5 × 10 3 TeV with mΦ< 2m 3/2 For a potentially most propitious case B2 with Φ decay to Higgs and matter turned off, then modulus branching fractions to SUSY and to gravitinos become highly suppressed at large m Φ . But since the modulus number density increases faster than the branching fractions decrease, there is still gross overproduction of neutralino dark matter. We also show that in this scenario the thermally produced gravitino problem is fixed by huge entropy dilution, but non-thermal gravitino production from moduli decay remains a huge problem unless it is kinematically suppressed with m Φ < 2m 3/2 . In a pedagogical appendix, we present detailed calculations of modulus field two-body decay widths.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

General Approach to Solving Spin Moiré Superstructures

Recently, a host of exciting magnetic textures such as topologically protected skyrmion lattices has been discovered in several bulk metallic lanthanide compounds. In addition to hosting skyrmion phases, a hallmark of this class of materials is the appearance of numerous spin textures characterized by superposition of multiple magnetic modulations: spin moiré superlattices. In order to understand the multitude of complex phases often present in these materials, we require a general-purpose experimental and theoretical framework. Here, we demonstrate such an approach in EuAg4⁢Sb2 by comprehensively characterizing and modeling its three complex zero-field magnetic textures. Systematic symmetry-breaking experiments using uniaxial strain determine that the ground-state incommensurate magnetic phase (ICM1) is single 𝑞, meaning the magnetic moments modulate along one magnetic propagation vector. In contrast, ICM2 and ICM3 are both double 𝑞, meaning they are formed from the superposition of two sinusoidal spin modulations, i.e., spin moiré superlattices. Further, through application of polarized small-angle neutron scattering and spherical neutron polarimetry, we demonstrate that ICM1 is a single-𝑞 cycloid and ICM2 and ICM3 are double-𝑞 vortex lattices. Despite the quasi-two-dimensional nature of EuAg4⁢Sb2, the modulations propagate out of the ab plane, leading to a shift of the spin texture between triangular lattice planes. Further, the ICM3 to ICM2 transition includes an unusual 45° rotation of the magnetic vortex lattice. Motivated by the coexistence of such drastically different phases in this compound, we conclude by developing a phenomenological model that sheds light on the energetic origins of these varied phases. Our experimental probes and theoretical modeling definitively characterize three different and tunable phases in one material and provide insight for the design of topological spin-texture materials.

Neves, Paul [Massachusetts Institute of Technology↗
DOE OSTI2021

Integrable model of topological SO(5) superfluidity

Assisted by general symmetry arguments and a many-body invariant, we introduce a phase of matter that constitutes a topological SO(5) superfluid. Key to this finding is the realization of an exactly solvable model that displays some similarities with a minimal model of superfluid 3 He. We study its quantum phase diagram and correlations, and find exotic superfluid as well as metallic phases in the repulsive sector. At the critical point separating trivial and nontrivial superfluid phases, our Hamiltonian reduces to the globally SO(5)-symmetric Gaudin model with a degenerate ground manifold that includes quartet states. Most importantly, the exact solution permits uncovering of an interesting non-pair-breaking mechanism for superfluids subject to external magnetic fields. Furthermore, nonintegrable modifications of our model lead to a strong-coupling limit of our metallic phase with a ground-state manifold that shows an extensive entropy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

All Optical Neural Networks for Low Power Edge Computing

We developed a simplistic physics-based model of an all-optical neural network that mimics the encoder part of an autoencoder neural network for image compression. Our approach relies on the generation of a MATLAB-based model for both data compression and decompression and utilizes MATLAB's built-in autoencoder networks in combination with simple propagation of optical fields between layers constituting phase elements via Fourier transform. We optimize the phase elements using the particle swarm optimization technique and using our model, we demonstrate a compression ratio of 25% for 2828-pixel input images containing numeric digits from 0 to 9.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Machine learning for the identification of phase transitions in interacting agent-based systems: A Desai-Zwanzig example

Deriving closed-form analytical expressions for reduced-order models, and judiciously choosing the closures leading to them, has long been the strategy of choice for studying phase- and noise-induced transitions for agent-based models (ABMs). In this paper, we propose a data-driven framework that pinpoints phase transitions for an ABM—the Desai-Zwanzig model—in its mean-field limit, using a smaller number of variables than traditional closed-form models. To this end, we use the manifold learning algorithm Diffusion Maps to identify a parsimonious set of data-driven latent variables, and we show that they are in one-to-one correspondence with the expected theoretical order parameter of the ABM. We then utilize a deep learning framework to obtain a conformal reparametrization of the data-driven coordinates that facilitates, in our example, the identification of a single parameter-dependent ordinary differential equation (ODE) in these coordinates. Additionally, we identify this ODE through a residual neural network inspired by a numerical integration scheme (forward Euler). We then use the identified ODE—enabled through an odd symmetry transformation—to construct the bifurcation diagram exhibiting the phase transition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Phase diagram of the anisotropic triangular lattice Hubbard model

In a recent study [Phys. Rev. X 10, 021042 (2020)], we showed using large-scale density matrix renormalization group (DMRG) simulations on infinite cylinders that the triangular lattice Hubbard model has a chiral spin liquid phase. In this work, we introduce hopping anisotropy in the model, making one of the three distinct bonds on the lattice stronger or weaker compared with the other two. We implement the anisotropy in two inequivalent ways, one which respects the mirror symmetry of the cylinder and one which breaks this symmetry. In the full range of anisotropy, from the square lattice to weakly coupled one-dimensional chains, we find a variety of phases. Near the isotropic limit we find the three phases identified in our previous work: metal, chiral spin liquid, and 120$^\circ$ spiral order; we note that a recent paper suggests the apparently metallic phase may actually be a Luther-Emery liquid, which would also be in agreement with our results. Furthermore, when one bond is weakened by a relatively small amount, the ground state quickly becomes the square lattice N\'{e}el order. When one bond is strengthened, the story is much less clear, with the phases that we find depending on the orientation of the anisotropy and on the cylinder circumference. While our work is to our knowledge the first DMRG study of the anisotropic triangular lattice Hubbard model, the overall phase diagram we find is broadly consistent with that found previously using other methods, such as variational Monte Carlo and dynamical mean field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

On the use of charged particles to characterize precipitation in irradiated reactor pressure vessel steels with a wide range of compositions

Nuclear reactor lifetimes may be limited by nano-scale Cu-Mn-Ni-Si precipitates (CRPs and MNSPs) that form under neutron irradiation (NI) of pressure vessel (RPV) steels, resulting in hardening and ductile to brittle transition temperature increases (embrittlement). Physical models of embrittlement must be based on characterization of precipitation as a function of the combination of metallurgical and irradiation variables. Here we focus on rapid and convenient charged particle irradiations (CPI) to both: a) compare to precipitates formed in NI; and, b) use CPI to efficiently explore precipitation in steels with a very wide range of compositions. Atom probe tomography (APT) comparisons show NI and CPI for similar bulk steel solute contents yield nearly the same precipitate compositions, albeit with some differences in their number density, size and volume fraction (f) dose (dpa) dependence. Yet, the overall precipitate evolutions are very similar. Advanced high Ni (>3 wt%) RPV steels, with superior unirradiated properties, were also investigated at high CPI dpa. For typical Mn contents, MNSPs have Ni 16 Mn 6 Si 7 or Ni 3 Mn 2 Si phase type compositions, with f values that are close to the equilibrium phase separated values. In contrast, in steels with very low Mn and high Ni, Ni 2-3 Si silicide phase type precipitate compositions are observed; and when Ni is low, the precipitate compositions are close to the MnSi phase field. Low Mn significantly reduces, but does not eliminate, precipitation in high Ni steels. A comparison of dispersed barrier model predictions with measured hardening data suggests that the Ni-Si dominated precipitates are weaker dislocation obstacles than the G phase type MNSPs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗
DOE OSTI2021

A multiscale insight into the growth of h-BN: effect of the enclosure

Abstract There is a lack of knowledge on the fundamental growth mechanisms governing the characteristics of 2D materials synthesized by the chemical vapor deposition (CVD) technique and their correlation with experimentally controllable parameters, which hindered their wafer-scale synthesis. Here, we pursued an analytical and computational approach to access the system states that are not experimentally viable to address these critical needs. We developed a multiscale computational framework correlating the macroscale heat and mass flow with the mesoscale morphology of the as-grown 2D materials by solving the coupled system of heat/mass transfer and phase-field equations. We used hexagonal boron nitride (h-BN) as our model material and investigated the effect of substrate enclosure on its growth kinetics and final morphology. We revealed a lower concentration with a more uniform distribution on the substrate in an enclosed-growth than open-growth. It leads to a more uniform size distribution of the h-BN islands, consistent with existing experimental investigations.

Materials Science↗
DOE OSTI2025

XY-like incommensurate magnetic order in Ce 2 ⁢SnS 5

We report the synthesis of single crystals of Ce 2 ⁢SnS 5 through a two-stage chemical vapor transport method. The Ce 2 ⁢SnS 5 system is a member of the orthorhombic Pbam (No. 55) space group and realizes a distorted trigonal tricapped prism (TTP) crystal field around each cerium site. We characterized the sample through orientation-dependent magnetization and heat capacity measurements to probe the magnetic anisotropy in the system characteristic of XY-like anisotropic Heisenberg model behavior. Ce 2 ⁢SnS 5 furthermore enters a zero-field ordered phase under 𝑇 𝑁 =2.4 K; powder neutron diffraction measurements reveal incommensurate magnetic order near 𝑇 𝑁 . Furthermore, the system then locks into a commensurate, two-𝑞 magnetic structure below approximately 1.2 K. This commensurate structure belongs to the Shubnikov group 𝑃⁢𝑏′⁢𝑎′⁢𝑚′ ⁡(MSG 55.359) and realizes the propagation vectors $\overrightarrow{𝑞}$ = (1/3, 0, 0) and $\overrightarrow{𝑞}$ = (0, 0, 0).

Antiferromagnetism↗