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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantum paramagnetism in a non-Kramers rare-earth oxide: Monoclinic Pr 2 Ti 2 O 7

Little is so far known about the magnetism of the A 2 B 2 O 7 monoclinic layered perovskites that replace the spin-ice supporting pyrochlore structure for r A /r B > 1.78. We show that high quality monoclinic Pr 2 Ti 2 O 7 single crystals with a three-dimensional network of non-Kramers Pr 3+ ions that interact through edge-sharing superexchange interactions, form a singlet ground-state quantum paramagnet that does not undergo any magnetic phase transitions down to, at least, 1.8 K. The chemical phase stability, structure, and magnetic properties of the layered perovskite Pr 2 Ti 2 O 7 were investigated using x-ray diffraction, transmission electron microscopy, and magnetization measurements. Synthesis of polycrystalline samples with the nominal compositions of Pr 2 Ti 2+x O 7 (–0.16 ≤ x ≤ 0.16 ) showed that deviations from the Pr 2 Ti 2 O 7 stoichiometry lead to secondary phases of related structures including the perovskite phase Pr 2/3 TiO 3 and the orthorhombic phases Pr 4 Ti 9 O 24 and Pr 2 TiO 5 . No indications of site disordering (stuffing and antistuffing) or vacancy defects were observed in the Pr 2 Ti 2 O 7 majority phase. A procedure for growth of high-structural-quality stoichiometric single crystals of Pr 2 Ti 2 O 7 by the traveling solvent floating zone method is reported. Thermomagnetic measurements of single-crystalline Pr 2 Ti 2 O 7 reveal an isolated singlet ground state that we associate with the low-symmetry crystal electric-field environments that split the (2J + 1 = 9)-fold degenerate spin-orbital multiplets of the four differently coordinated Pr 3+ ions into 36 isolated singlets resulting in an anisotropic temperature-independent van Vleck susceptibility at low T. Here, a small isotropic Curie term is associated with 0.96(2)% noninteracting Pr 4+ impurities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Manipulating Surface Termination of Perovskite Manganate for Oxygen Activation

For ABO 3 perovskite oxides, one of the key issues limiting their utilization in heterogeneous catalysis is the dominant presence of catalytically inactive A-site cations at the surface. The engineering of B-site terminated perovskites is considered as an effective method to address this issue, especially when dealing with Mn/Co-based perovskite catalysts. However, to date, such a strategy has not been fully successful and remains a major challenge in the field. Herein, a Mn-terminated La 0.45 Sr 0.45 MnO 3 (B-LSM) is successfully synthesized via a one-pot hydrothermal method, in which low-valence Mn ions partially occupy the A site to form the active Mn-excess phase. Additionally, experimental results and theoretical calculations reveal that the presence of the surface Mn termination in B-LSM optimizes the hybrid orbitals of Mn 3d-O 2p and promotes the activation of surface lattice oxygen, where the pristine inert lattice O 2- is evolved into active and stable lattice O 2-x . Such structural optimization significantly reduces the activation energy barriers on going from O 2 ⁻ species to important intermediate O - species during O 2 activation. Moreover, this results in good stability and Pt-like activity for the B-LSM during CO oxidation. This work offers a new chemical route for the design of advanced perovskite-type oxides possessing novel functions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Novel Phases and Dynamics of Multi-orbital Mott Insulators

The goal of this project is to gain theoretical insight into magnetism in the presence of spin-orbit coupling (SOC) for multi-orbital Mott insulators. Inspired by the exact solution of the Kitaev model that harbors a quantum spin liquid with novel excitations and struck by its relevance for materials with anisotropic orbital interactions, the PI will explore various new classes of 4d and 5d transition metal oxides. Starting with all electron Hamiltonians, the PI will derive minimal magnetic models. The aim will be to understand the role played by orbital frustration, even in the absence of any geometric frustration, in creating orbitally ordered as well as spin-orbital liquid phases. Orbital frustration arises primarily from the directional or anisotropic nature of d-orbitals in contrast to the isotropic nature of the spin degree of freedom. These models will be investigated by a variety of theoretical and numerical methods, including exact diagonalization, mean field theories, density matrix renormalization group and quantum Monte Carlo methods. Testable predictions for three experiments: nuclear magnetic resonance (NMR) and resonant x-ray scattering (RXS) to probe orbital ordering and pump-probe experiments to probe quantum dynamics will test the validity of these models for materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Predictions of the glass transition temperature and viscosity of organic aerosols from volatility distributions

Abstract. Volatility and viscosity are important properties of organic aerosols (OA), affecting aerosol processes such as formation, evolution, and partitioning of OA. Volatility distributions of ambient OA particles have often been measured, while viscosity measurements are scarce. We have previously developed a method to estimate the glass transition temperature (Tg) of an organic compound containing carbon, hydrogen, and oxygen. Based on analysis of over 2400 organic compounds including oxygenated organic compounds, as well as nitrogen- and sulfur-containing organic compounds, we extend this method to include nitrogen- and sulfur-containing compounds based on elemental composition. In addition, parameterizations are developed to predict Tg as a function of volatility and the atomic oxygen-to-carbon ratio based on a negative correlation between Tg and volatility. This prediction method of Tg is applied to ambient observations of volatility distributions at 11 field sites. The predicted Tg values of OA under dry conditions vary mainly from 290 to 339 K and the predicted viscosities are consistent with the results of ambient particle-phase-state measurements in the southeastern US and the Amazonian rain forest. Reducing the uncertainties in measured volatility distributions would improve predictions of viscosity, especially at low relative humidity. We also predict the Tg of OA components identified via positive matrix factorization of aerosol mass spectrometer (AMS) data. The predicted viscosity of oxidized OA is consistent with previously reported viscosity of secondary organic aerosols (SOA) derived from α-pinene, toluene, isoprene epoxydiol (IEPOX), and diesel fuel. Comparison of the predicted viscosity based on the observed volatility distributions with the viscosity simulated by a chemical transport model implies that missing low volatility compounds in a global model can lead to underestimation of OA viscosity at some sites. The relation between volatility and viscosity can be applied in the molecular corridor or volatility basis set approaches to improve OA simulations in chemical transport models by consideration of effects of particle viscosity in OA formation and evolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Investigating the impact of donor atoms of single-source precursors on f-element nanomaterial formation and mechanisms of initial bond cleavage

Nanotechnology is advancing exponentially, with novel discoveries and applications emerging every day across various fields, including environmental remediation, biomedicine, electronics, and textiles. Nanoparticles (NPs) are defined as particles with sizes ranging from 1 to 100 nm, and their properties are influenced by factors such as composition, phase, size, and shape. Single-source precursors (SSPs) offer a promising method for controlling the synthesis of phase-pure nanoparticles with the desired stoichiometry and have been successfully applied to prepare transition metal and lanthanide NPs. SSPs are molecular compounds that contain all the necessary components to form the final nanomaterial and decompose under mild conditions, such as heating, to create the solid-state product. The objective of this work was to garner a mechanistic understanding of cerium oxide NP formation from cerium SSPs. To achieve these goals, first, ligands with oxygen and sulfur donor atoms were synthesized; among the array of ligands synthesized, some were novel. These ligands were then chelated to cerium ions to synthesize cerium precursors. Initially, known cerium alkoxide precursors were synthesized. Subsequently, these basic cerium complexes were further manipulated through acid-base protonolysis reactions to install a diverse set of ligands. Within this body of work, 10 novel cerium precursors were successfully synthesized. The ligands and the SSPs were analyzed using nuclear magnetic resonance (NMR) and thermogravimetric analysis coupled with differential scanning calorimetry (TGA-DSC). These SSPs were then introduced to a collision-induced dissociation mass spectrometer (CID-MS), and their dissociation pathways were studied.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2022

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
DOE OSTI2021

Colloidal quantum dot based infrared detectors: extending to the mid-infrared and moving from the lab to the field

Quantum dots (QDs) that absorb in the mid-wave infrared (MWIR) regime (3–5 μm) have recently generated significant interest as possible detector materials for MWIR cameras, with promises to reduce materials and device fabrication costs and potentially increase device operating temperatures. However, these materials have been primarily explored in single pixel devices – not multi-pixel cameras – and their current performance lags behind that of commercially available MWIR cameras based on single crystalline materials. To realize the potential of QD MWIR photodetectors with imaging capabilities, improvements are needed in the quality of the MWIR-absorbing QD materials, ligands dictating inter-QD charge transport, ordering of QD films, device architectures, and scalability of these methods to larger areas. With three families of QDs being researched – lead-, silver-, and mercury-based chalcogenides – as well as a multitude of possible capping ligands and film deposition techniques, the experimental phase space for MWIR QDs is vast. Here, we provide a roadmap, considering the pros and cons of various film deposition and ligand exchange techniques, as well as reintroduce lessons learned over years of research on QD film formation. We also extracted and created a database of reported quantum dot photodetectors’ performance and fabrication methods and have developed an interactive data visualization dashboard for this database, which provides researchers in the field a quick snapshot of the existing state of the art. Finally, we outline figures of merit and information that should be presented in papers moving forward which would help to clarify the reported results, as well as offer ideas for future steps. We provide a database visualization of the QD photodetector literature at https://public.tableau.com/app/profile/tom.nakotte/viz/IRQDphotodetectors_16384709473480/Dashboard1#1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Niobium hydride formation in superconducting qubit thin films

The formation of nonsuperconducting hydrides in 160–170-nm-thick films of niobium is examined. We identify six elastically distinct orientation relationships ɛ−Nb 4 ⁢H 3 takes within the matrix, solid-solution 𝛼−NbH 𝑥 phase. We employ a phase field model to assess the impact of elastic energy induced by the strain of phase transformation on the morphology and transformation dynamics of ɛ precipitates within a thin film. We consider the dimensions of the thin film, crystallographic growth direction, and diffusion rates to predict the timescale of hydride evolution. Leveraging the finite element method, we predict two-dimensional and three-dimensional equilibrium shapes of ɛ-hydrides within a bulk sample and in a thin film that has a traction-free surface. Our results suggest that niobium hydrides migrate to the free surface of the film. Precipitates which reach the free surface coarsen, while precipitates within the film dissolve. Precipitates in both two dimensions and three dimensions experience a repulsive interaction force at the free surface, that is attractive in the bulk, shown in experiment and theory of previous studies.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Increased phase precision of spatial light modulators using irrational slopes: application to attosecond metrology

The ability of spatial light modulators (SLMs) to modify the amplitude and phase of light has proved them invaluable to the optics and photonics community. In many applications, the bit-depth of SLMs is a major limiting factor dictated by a digital processor. As a result, there is usually a compromise between refresh speed and bit-depth. Here, we present a method to increase the effective bit-depth of SLMs, which utilizes a linear slope, as is commonly applied to deal with the zeroth-order effect. This technique was tested using two interferometric transient absorption spectroscopy setups. Through the high harmonic generation in gases producing a train of attosecond pulses and harmonics from solids in the ultraviolet, two pulses are generated that interfere in the far field providing a measurement of the optical phase. An increase in the precision far beyond the limit dictated by the digital processor in the bit-depth was found.

Harrison, Geoffrey R. (ORCID:0000000152141866)↗
DOE OSTI2020

A Case Study of Nonresonant Mode 3-s ULF Waves Observed by MMS

The nature of the 3-s ultralow frequency (ULF) wave in the Earth's foreshock region and the associated wave-particle interaction are not yet well understood. In this work, we investigate the 3-s ULF waves using Magnetospheric Multiscale (MMS) observations. By combining the plasma rest frame wave properties obtained from multiple methods with the instability analysis based on the velocity distribution in the linear wave stage, the ULF wave is determined to be due to the ion/ion nonresonant mode instability. The interaction between the wave and ions is analyzed using the phase relationship between the transverse wave fields and ion velocities and using the longitudinal momentum equation. During the stage when ULF waves have sinusoidal waveforms up to |dB|/|B 0 | ~ 3, where dB is the wave magnetic field and B0 is the background magnetic field, the wave electric fields perpendicular to B 0 do negative work to solar wind ions; along B 0 , a longitudinal electric field develops, but the V × B force is stronger and leads to solar wind ion deceleration. During the same wave stage, the backstreaming beam ions gain energy from the transverse wave fields and get deceleration along B 0 by the longitudinal electric field. The ULF wave leads to electron heating, preferentially in the direction perpendicular to the local magnetic field. Secondary waves are generated within the ULF waveforms, including whistler waves near half of the electron cyclotron frequency, high-frequency electrostatic waves, and magnetosonic whistler waves. The work improves the understanding of the nature of 3-s ULF waves and the associated wave-particle interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Recent Advances in Small Angle X-ray Scattering for Superlattice Study

Small-angle x-ray scattering is used for the structure determination of superlattice for its superior resolution, nondestructive nature, and high penetration power of x rays. With the advent of high brilliance x-ray sources and innovative computing algorithms, there have been notable advances in small angle x-ray scattering analysis of superlattices. High brilliance x-ray beams have made data analyses less model-dependent. Additionally, novel data acquisition systems are faster and more competitive than ever before, enabling a more accurate mapping of the superlattices' reciprocal space. Fast and high-throughput computing systems and algorithms also make possible advanced analysis methods, including iterative phasing algorithms, non-parameterized fitting of scattering data with molecular dynamics simulations, and the use of machine learning algorithms. As a result, solving nanoscale structures with high resolutions has become an attainable task. In this review, we highlight new developments in the field and introduce their applications for the analysis of nanoscale ordered structures, including nanoparticle supercrystals, nanoscale lithography patterns, and supramolecular self-assemblies. Particularly, we highlight the reciprocal space mapping techniques and the use of iterative phase retrieval algorithms. We also cover coherent-beam-based small angle x-ray scattering techniques such as ptychography and ptycho-tomography in view of the traditional small angle x-ray scattering perspective.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Quantifying particle-scale 3D granular dynamics during rapid compaction from time-resolved in situ 2D x-ray images

Understanding the particle-scale dynamics of granular materials during rapid compaction and flow is of fundamental importance for manufacturing, planetary science, geology, and defense applications. Time-resolved 2D radiography and static 3D x-ray tomography are powerful in situ tools for studying particle-scale dynamics but provide detail only in 2D or with significant time-scale limitations, respectively. Here, we introduce a new method that uses 2D in situ x-ray imaging for determining time-resolved 3D particle-scale dynamics in rapidly compressed granular materials. The method employs initial particle packing structures obtained from x-ray tomography, a 2D x-ray image generation algorithm, and an optimization algorithm. We first describe and validate the method using finite element simulations. Furthermore, we then apply the technique to x-ray phase-contrast images obtained during rapid compaction of granular materials with varying particle morphology and sample thickness. The depth-resolved particle-scale dynamics reveal complex velocity and porosity fields evolving heterogeneously along and perpendicular to the compaction direction. We characterize these features, their fluctuations near the compaction front, and the compaction front thickness. Our technique can be applied to understanding granular dynamics during rapid compaction events, and rearrangements during slower, but non-quasi-static, flows.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Photon-Instanton Collider Implemented by a Superconducting Circuit

Instantons, spacetime-localized quantum field tunneling events, are ubiquitous in correlated condensed matter and high-energy systems. However, their direct observation through collisions with conventional particles has not been considered possible. We show how recent advances in circuit quantum electrodynamics, specifically, the realization of galvanic coupling of a transmon qubit to a high-impedance transmission line, allows the observation of inelastic collisions of single microwave photons with instantons (phase slips). We develop a formalism for calculating the photon-instanton cross section, which should be useful in other quantum field theoretical contexts. In particular, we show that the inelastic scattering probability can significantly exceed the effect of conventional Josephson quartic anharmonicity and reach order-unity values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A windowed mean trajectory approximation for condensed phase dynamics

We propose a trajectory-based quasi-classical method for approximating dynamics in condensed phase systems. Building upon the previously developed optimized mean trajectory approximation that has been used to compute linear and nonlinear spectra, we borrow some ideas from filtering trajectory methods to obtain a novel semiclassical method for the dynamical propagation of density matrices. This new approximation is tested rigorously against standard multistate electronic models, spin-boson models, and models of the Fenna–Matthews–Olson complex. For dissipative systems, the current method is significantly better or as good as many other semiclassical methods available, especially at low temperatures and for off-diagonal density matrix elements, whereas for scattering models, the current method bears similar limitations as mean-field propagation schemes. All results are tested against the numerically exact hierarchical equations of motion method. In conclusion, the new method shows excellent agreement across various parameter regimes with numerically exact results, highlighting the robustness and accuracy of our approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗