Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 433 records · Page 24

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dimerization and spin decoupling in a two-leg Heisenberg ladder with frustrated trimer rungs

We study the antiferromagnetic spin-half Heisenberg ladder in the presence of an additional frustrating rung spin which is motivated and relevant also for the description of real two-dimensional materials such as the two-dimensional trimer magnet Ba 4 Ir 3 O 10 . We study the zero-temperature phase diagram, where we combine numerical and analytical methods into an overall consistent description. All numerical simulations are also accompanied by studies of the dynamical spin structure factor obtained via the density matrix renormalization group. Overall, we find in the regime of strong rung coupling a gapped dimerized phase related to competing symmetry sectors in Hilbert space that ultimately results in frustration-driven spin-Peierls transition. In the weak rung-coupling regime, the system is uniform, yet shows a gapped spinon continuum together with a sharp coherent low-energy branch which renders the system critical overall. In either case, the additional rung spin quickly get sidelined and nearly decouple in the regime when their bare coupling to the ladder is weakened

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High-order ( N = 4–6) multiphoton absorption and mid-infrared Kerr nonlinearity in GaP, ZnSe, GaSe, and ZGP crystals

We report a study of high-order multiphoton absorption, nonlinear refraction, and their anisotropy in four notable mid-infrared χ (2) crystals: GaP, ZnSe, GaSe, and ZGP using femtosecond pulses at 2.35 μm with peak intensity in excess of 200 GW/cm 2 . We found that the nonlinear absorption obeys a perturbation model with multiphoton absorption (MPA) orders from N = 4–6, in agreement with the bandgaps of the crystals. A study of the role of free carrier absorption, performed by changing the pulse duration from 30 to 70 fs while maintaining a constant peak intensity, showed that, at our intensity levels, free carriers generated during the MPA process absorb much more strongly than would be expected from their known linear absorption cross section. Possible mechanisms for this anomalous behavior are high-field effects, such as intravalley scattering in the conduction band and absorption to higher lying bands. Nonlinear refractive indices were measured using (i) closed aperture Z-scan and (ii) spectral broadening due to self-phase modulation, both methods agreeing with each other.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Learning heterogeneous reaction kinetics from X-ray movies pixel-by-pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries 1 and electrocatalysts 2. Experimental characterizations of such materials by operando microscopy produce rich image datasets 3–6, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation 7. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

Chueh, William↗

helioLuna: Prospect for cislunar spacecraft and near-earth asteroid detection using heliostat fields at night

We experimentally and computationally investigate a proposed frequency-domain method for detecting and tracking cislunar spacecraft and near-earth asteroids using heliostat fields at night. Unlike imaging, which detects spacecraft and asteroids by their streak in sidereally-fixed long-exposure photographs, our proposed detection method oscillates the orientation of heliostats concentrating light from the stellar field and measures the light’s photocurrent power spectrum at sub-milliHertz resolution. If heliostat oscillation traces out an ellipse fixed in the galactic coordinate system, spacecraft or asteroids produce a peak in photocurrent power spectrum at a frequency slightly shifted from the starlight peak. The frequency shift is on the scale of milliHertz and proportional to apparent angular rate relative to sidereal. Relative phase corresponds to relative angular position, enabling tracking. A potential advantage of this frequency-domain method over imaging is that detectivity improves with apparent angular rate and number of heliostats. Since heliostats are inexpensive compared to an astronomical observatory and otherwise unused at night, the proposed method may cost-effectively augment observatory systems such as NASA’s Asteroid Terrestrial-impact Last Alert System (ATLAS).

79 ASTRONOMY AND ASTROPHYSICS↗

Strong-field nonsequential double photoionization using virtual-detector theory with path summation

We present an ab initio study of the nonsequential strong-field ionization dynamics of a model two-electron atom with helium character. Single- and double-ionization events are characterized and displayed using detector signals extracted at different points in the two-electron two-dimensional space. The double photoelectron momentum distribution is calculated via coherent path-summation over virtual-particle trajectories. A comparison is made between the momentum distributions obtained with the virtual-detector method, the Schrödinger wave function, and the time-dependent surface flux method developed by Tao and Scrinzi [New J. Phys. 14, 013021 (2012)]. Insights into different ionization and electron recollision pathways are gained from detailed virtual-particle tracking and energy-time readouts. Furthermore, this study demonstrates the extension of virtual-detector theory to strong-field multielectron quantum dynamics and highlights the importance of the evolving quantum phase in quasiclassical electron propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Regulation of Cathode Mass and Charge Transfer by Structural 3D Engineering for Protonic Ceramic Fuel Cell at 400 °C

Abstract Lowering the operating temperature (ideally below 400 °C) for solid oxide fuel cell (SOFC) technology deployment has been an important transition that introduces the benefit of reduced operational costs and system durability. However, the key technical issue limiting the transition is the sluggish cathodic performance, namely the oxygen reduction reaction (ORR) rate of the conventional sponge‐like cathode dramatically drops as the temperature reduces. In this paper, 3D engineering of a cathode is conducted on a protonic ceramic fuel cell to obtain an enhanced ORR between 400 and 600 °C. Compared with a cell using a conventional sponge‐like cathode, 3D engineering improves the cathode ORR by 41% at 400 °C with a peak power density of 0.410 W cm −2 . A phase field simulation is applied to assist the engineering by understanding the competition between the cathode mass and charge transfer with different cathode porosities. The results show that structural engineering of existing well‐developed cathodes is a simple and effective method to promote cathode ORR for low temperature SOFC by regulating the mass and charge transfer.

Bian, Wenjuan↗

Submonomer synthesis of sequence defined peptoids with diverse side-chains.

Peptoids are a diverse family of sequence-defined oligomers of N-substituted glycine monomers, that can be readily accessed by the solid-phase submonomer synthesis method. Due to the versatility and efficiency of this chemistry, and the easy access to hundreds of potential monomers, there is an enormous potential sequence space that can be explored. This has enabled researchers from many different fields to custom-design peptoid sequences tailored to a wide variety of problems in biomedicine, nanoscience and polymer science. Here we provide detailed protocols for the synthesis of peptoids, using optimized protocols that can be performed by non-chemists. The submonomer method is fully compatible with Fmoc-peptide synthesis conditions, so the method is readily automated on existing automated peptide synthesizers using protocols provided here. Although the submonomer synthesis for peptoids is well established, there are special considerations required in order to access many of the most useful and desirable sidechains. Here we provide methods to include most of the amino-acid-like side chains, some of the most important non-natural monomer classes, as well as the creation of peptoid conjugates and peptide-peptoid hybrids.

Connolly, Michael D↗

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗

Plasma Assisted Combustion and Chemical Processing: Chapter 9 - Plasma Diagnostics

Plasma dynamics and chemistry have a broad range of timescales from picoseconds to milliseconds. In addition, it involves nonequilibrium energy transfer between electrons, ions, electronically and vibrationally excited states, radicals, intermediate species, and reactants and products as well as surface charges and chemistry. To understand plasma physics and chemistry, it is essential to conduct time and space resolved, quantitative detection of nonequilibrium temperature distributions, electron energy and number density, electric field, and species concentrations. There are enormous publications and review articles on this subject. The focus of this chapter is to be placed on the most recent progress in gas phase plasma properties and chemistry, especially on optical emission spectroscopy, laser absorption spectroscopy, Faraday rotational spectroscopy, Raman and Thompson scattering, femtosecond and picosecond (fs/ps) coherent anti-Stokes Raman scattering (CARS) spectroscopy, and electric field-induced second harmonic generation methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Designing open quantum systems with known steady states: Davies generators and beyond

We provide a systematic framework for constructing generic models of nonequilibrium quantum dynamics with a target stationary (mixed) state. Our framework identifies (almost) all combinations of Hamiltonian and dissipative dynamics that relax to a steady state of interest, generalizing the Davies’ generator for dissipative relaxation at finite temperature to nonequilibrium dynamics targeting arbitrary stationary states. We focus on Gibbs states of stabilizer Hamiltonians, identifying local Lindbladians compatible therewith by constraining the rates of dissipative and unitary processes. Moreover, given terms in the Lindbladian not compatible with the target state, our formalism identifies the operations – including syndrome measurements and local feedback – one must apply to correct these errors. Our methods also reveal new models of quantum dynamics: for example, we provide a “measurement-induced phase transition” in which measurable two-point functions exhibit critical (power-law) scaling with distance at a critical ratio of the transverse field and rate of measurement and feedback. Time-reversal symmetry – defined naturally within our formalism – can be broken both in effectively classical and intrinsically quantum ways. Our framework provides a systematic starting point for exploring the landscape of dynamical universality classes in open quantum systems, as well as identifying new protocols for quantum error correction.

Guo, Jinkang [Department of Physics and Center for↗

Industrial PLC Network Modeling and Parameter Identification Using Sensitivity Analysis and Mean Field Variational Inference

A multiple input multiple output (MIMO) power line communication (PLC) model for industrial facilities was developed that uses the physics of a bottom-up model but can be calibrated like top-down models. The PLC model considers 4-conductor cables (three-phase conductors and a ground conductor) and has several load types, including motor loads. The model is calibrated to data using mean field variational inference with a sensitivity analysis to reduce the parameter space. The results show that the inference method can accurately identify many of the model parameters, and the model is accurate even when the network is modified.

47 OTHER INSTRUMENTATION↗

Kinematic Variables and Feature Engineering for Particle Phenomenology

Kinematic variables have been playing an important role in collider phenomenology, as they expedite discoveries of new particles by separating signal events from unwanted background events and allow for measurements of particle properties such as masses, couplings, spins, etc. For the past 10 years, an enormous number of kinematic variables have been designed and proposed, primarily for the experiments at the Large Hadron Collider, allowing for a drastic reduction of high-dimensional experimental data to lower-dimensional observables, from which one can readily extract underlying features of phase space and develop better-optimized data-analysis strategies. We review these recent developments in the area of phase space kinematics, summarizing the new kinematic variables with important phenomenological implications and physics applications. We also review recently proposed analysis methods and techniques specifically designed to leverage the new kinematic variables. As machine learning is nowadays percolating through many fields of particle physics including collider phenomenology, we discuss the interconnection and mutual complementarity of kinematic variables and machine learning techniques. We finally discuss how the utilization of kinematic variables originally developed for colliders can be extended to other high-energy physics experiments including neutrino experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Understanding and control of Zener pinning via phase field and ensemble learning

Zener pinning refers to the dispersion of fine particles which influences grain size distribution via movement of grain boundaries in a polycrystalline material. Grain size distribution in polycrystals has a significant impact on their properties including physical, chemical, mechanical, and optical to name a few. We explore the use of Phase-field modeling and machine-learning techniques to understand and improve the control of grain size distribution via Zener pinning in polycrystalline materials. We develop a machine learning model that determines the relative importance of various parameters to exercise microstructure control via Zener pinning. Our workflow combines high-throughput phase-field simulations and machine learning to address the computational bottlenecks associated with large-scale simulations as well as identify features necessary for microstructure control in polycrystals. A random forest (RF) regression model was developed to predict grain sizes based on five Phase-field model parameters, achieving an average prediction error of 0.72 nm for the training data and 1.44 nm for the test data. The importance of the input parameters is analyzed using the SHapley Additive exPlanations (SHAP) approach which reveals that diffusivity, volume fraction, and particle diameter are the most important parameters in determining the final grain size. These findings will allow us to select the best second-phase particles, optimize grain size distributions and thus design microstructures with the desired properties. The developed method is a highly versatile and generalizable approach that can be used to assess the combined effects of individual features in the presence of multiple variables.

36 MATERIALS SCIENCE↗

Development and testing of coarse-grained models for ultrasonic simulations of cast austenitic stainless steel

Ultrasonic inspection of cast austenitic stainless steel (CASS) in the nuclear industry is particularly challenging because of sound field scatter and attenuation caused by the coarse-grained microstructure. Modeling and simulation are important tools in ultrasonic testing, as they can be used to help address key aspects of inspections, such as developing new probe designs, predicting inspection reliability, and testing phased-array focal laws. However, developing a useful and reliable CASS model is challenging due to the many grain interfaces and crystalline orientations that must be captured. We demonstrate a method of creating a realistic CASS model that is usable in CIVA, a commercially available modeling and simulation software platform. Using polished and chemically etched sections, we generate models of a coarse-grained equiaxed specimen and a columnar specimen. We also test an alternative method of generating a coarse-grained model using Voronoi regions. We qualitatively compare sound field scatter and quantitatively compare sound field attenuation and beam partitioning in simulated sound fields to those of laboratory-measured sound fields. Results show that the Voronoi models perform as well as or better than the models based on actual grain morphology. Here we also show that model-to-model randomness in Voronoi grain structure can impact the magnitude of a simulated echo response by a factor of two or more. Although CASS models are potentially a good depiction of reality for a given scenario, they should not be considered representative since CASS morphology can change significantly from specimen to specimen or within the same specimen.

36 MATERIALS SCIENCE↗

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗