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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Deep learning on the 2-dimensional Ising model to extract the crossover region with a variational autoencoder

The 2-dimensional Ising model on a square lattice is investigated with a variational autoencoder in the non-vanishing field case for the purpose of extracting the crossover region between the ferromagnetic and paramagnetic phases. The encoded latent variable space is found to provide suitable metrics for tracking the order and disorder in the Ising configurations that extends to the extraction of a crossover region in a way that is consistent with expectations. The extracted results achieve an exceptional prediction for the critical point as well as agreement with previously published results on the configurational magnetizations of the model. The performance of this method provides encouragement for the use of machine learning to extract meaningful structural information from complex physical systems where little a priori data is available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Anyon superconductivity from topological criticality in a Hofstadter–Hubbard model

We argue that the combination of strong repulsive interactions and high magnetic fields can generate electron pairing and superconductivity. Inspired by the large lattice constants of moiré materials, which make large flux per unit cell accessible at laboratory fields, we study the triangular lattice Hofstadter-Hubbard model at one-quarter flux quantum per plaquette, where previous literature has argued that a chiral spin liquid separates a weak-coupling integer quantum Hall phase and a strong-coupling topologically trivial antiferromagnetic insulator at a density of one electron per site. We argue that topological superconductivity emerges upon doping in the vicinity of the integer quantum Hall to chiral spin liquid transition. We employ exact diagonalization and density matrix renormalization group methods to examine this theoretical scenario and find that electronic pairing indeed occurs on both sides of criticality over a remarkably broad range of interaction strengths. On the chiral spin liquid side, our results provide a concrete model realization of the long-hypothesized mechanism of anyon superconductivity. Our study thus establishes a beyond-Bardeen-Cooper-Schrieffer route to electron pairing in a well-controlled limit, relying crucially on the interplay between electron correlations and band topology.

Hubbard model↗
NASA NTRS1972

Feasibility of wake vortex monitoring systems for air terminals

Wake vortex monitoring systems, especially those using laser Doppler sensors, were investigated. The initial phases of the effort involved talking with potential users (air traffic controllers, pilots, etc.) of a wake vortex monitoring system to determine system requirements from the user's viewpoint. These discussions involved the volumes of airspace to be monitored for vortices, and potential methods of using the monitored vortex data once the data are available. A subsequent task led to determining a suitable mathematical model of the vortex phenomena and developing a mathematical model of the laser Doppler sensor for monitoring the vortex flow field. The mathematical models were used in combination to help evaluate the capability of laser Doppler instrumentation in monitoring vortex flow fields both in the near vicinity of the sensor (within 1 kilometer and at long ranges(10 kilometers).

Wilson, D. J.↗
DOE OSTI2021

Summer Internship Project: Set Topology [Slides]

Why a Set Topology? 1. Expand the array of tools given to physicists, educators, and students. 2. Through usages of PIC, MPM, and Monte Carlo. 3. My work at LANL was specified to handle requirements for the usage of MPM with FleCSI. 4. What we want it to do - phase change, simulating plasma, and other Multiphysics applications. 5. In FleCSI we already have mesh structures, we needed a way to represent those particles and certain behavior's of these particles that clients desire.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Symmetry breaking forms split-off flat bands in quantum oxides controlling metal versus insulator phases

The crystal structure used as input to electronic structure calculation of conventionally bonded solids consists primarily of the standard crystallographic degrees of freedom. Quantum solids sometimes have additional microscopic degrees of freedom (m-DOF) nested within the crystallographic structure. These might include local motifs including positional (Peierls dimers, deformed octahedra, Jahn-Teller distortions), magnetic (local moment configurations) and dipolar (local ferroelectric configurations). Such motifs can be observed experimentally via local probes that avoid averaging, and theoretically as distinct total energy lowering features relative to more simplified “average crystallographic structures”. Here we examine the ability of electronic structure methods independent of strong correlation physics to explain the broad phenomenology of metal-insulator selectivity in quantum oxides by energy-lowering symmetry breaking. We do this by avoiding the restriction of considering only (i) electron-electron interactions in a fixed unresponsive lattice—as done in Mott strong correlation explanations—allowing, however, (ii) local positional and magnetic motifs that coexist within a crystallographic landscape. We find in a broad range of quantum oxides with different symmetry breaking modes the formation of split-off flat bands that can also control metallic versus insulating characteristics. The split-off band effect is common to magnetic and non-magnetic materials, including both binary and ternary oxides. One finds that when such local symmetry breaking motifs are considered in mean-field-like (e.g. density functional theory) electronic structure calculations of quantum oxides, they explain many of the observed trends in metal vs insulator phases discussed otherwise in prior literature via strong correlation in symmetry unbroken view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Spatial and temporal evaluations of the liquid argon purity in ProtoDUNE-SP

Liquid argon time projection chambers (LArTPCs) rely on highly pure argon to ensure that ionization electrons produced by charged particles reach readout arrays. ProtoDUNE Single-Phase (ProtoDUNE-SP) was an approximately 700-ton liquid argon detector intended to prototype the Deep Underground Neutrino Experiment (DUNE) Far Detector Horizontal Drift module. It contains two drift volumes bisected by the cathode plane assembly, which is biased to create an almost uniform electric field in both volumes. The DUNE Far Detector modules must have robust cryogenic systems capable of filtering argon and supplying the TPC with clean liquid. This paper will explore comparisons of the argon purity measured by the purity monitors with those measured using muons in the TPC from October 2018 to November 2018. A new method is introduced to measure the liquid argon purity in the TPC using muons crossing both drift volumes of ProtoDUNE-SP. For extended periods on the timescale of weeks, the drift electron lifetime was measured to be above 30 ms using both systems. A particular focus will be placed on the measured purity of argon as a function of position in the detector.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Emergence of high piezoelectricity from competing local polar order-disorder in relaxor ferroelectrics

Relaxor ferroelectrics are known for outstanding piezoelectric properties, finding a broad range of applications in advanced electromechanical devices. Decoding the origins of the enhanced properties, however, have long been complicated by the heterogeneous local structures. Here, we employ the advanced big-box refinement method by fitting neutron-, X-ray-based total scattering, and X-ray absorption spectrum simultaneously, to extract local atomic polar displacements and construct 3D polar configurations in the classical relaxor ferroelectric Pb(Mg 1/3 Nb 2/3 )O 3 –PbTiO 3 . Our results demonstrate that prevailing order-disorder character accompanied by the continuous rotation of local polar displacements commands the composition-driven global structure evolution. The omnidirectional local polar disordering appears as an indication of macroscopic relaxor characteristics. Combined with phase-field simulations, it demonstrates that the competing local polar order-disorder between different states with balanced local polar length and direction randomness leads to a flattening free-energy profile over a wide polar length, thus giving rise to high piezoelectricity. Our work clarifies that the critical structural feature required for high piezoelectricity is the competition states of local polar rather than relaxor.

42 ENGINEERING↗
DOE OSTI2022

Dynamical topological phase realized in a trapped-ion quantum simulator

Nascent platforms for programmable quantum simulation offer unprecedented access to new regimes of far-from-equilibrium quantum many-body dynamics in almost isolated systems. Here achieving precise control over quantum many-body entanglement is an essential task for quantum sensing and computation. Extensive theoretical work indicates that these capabilities can enable dynamical phases and critical phenomena that show topologically robust methods to create, protect and manipulate quantum entanglement that self-correct against large classes of errors. However, so far, experimental realizations have been confined to classical (non-entangled) symmetry-breaking orders. In this work, we demonstrate an emergent dynamical symmetry-protected topological phase, in a quasiperiodically driven array of ten 171 Yb + hyperfine qubits in Quantinuum’s System Model H1 trapped-ion quantum processor. This phase shows edge qubits that are dynamically protected from control errors, cross-talk and stray fields. Crucially, this edge protection relies purely on emergent dynamical symmetries that are absolutely stable to generic coherent perturbations. Here, this property is special to quasiperiodically driven systems: as we demonstrate, the analogous edge states of a periodically driven qubit array are vulnerable to symmetry-breaking errors and quickly decohere. Our work paves the way for implementation of more complex dynamical topological orders that would enable error-resilient manipulation of quantum information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Chemomechanics in alloy phase stability

We describe a first-principles statistical mechanics method to calculate the free energies of crystalline alloys that depend on temperature, composition, and strain. The approach relies on an extension of the alloy cluster expansion to include an explicit dependence on homogeneous strain in addition to site occupation variables that track the degree of chemical ordering. The method is applied to the Si-Ge binary alloy and is used to calculate free energies that describe phase stability under arbitrary epitaxial constraints. We find that while the incoherent phase diagram (in which coexisting phases are not affected by coherency constraints) hosts a miscibility gap, coherent phase equilibrium predicts ordering and negative enthalpies of mixing. Instead of chemical instability, the chemomechanical free energy exhibits instabilities along directions that couple the composition of the alloy with a volumetric strain order parameter. Furthermore, this has fundamental implications for phase field models of spinodal decomposition as it indicates the importance of gradient energy coefficients that couple gradients in composition with gradients in strain.

Materials Science↗
DOE OSTI2021

Isotopic fractionation as in-situ sensor of subsurface reactive flow and precursor for rock failure

Greater utilization of subsurface reservoirs perturbs in-situ chemical-mechanical conditions with wide ranging consequences from decreased performance to project failure. Understanding the chemical precursors to rock deformation is critical to reducing the risks of these activities. To address this need, we investigated the coupled flow-dissolution- precipitation-adsorption reactions involving calcite and environmentally-relevant solid phases. Experimentally, we quantified (1) stable isotope fractionation processes for strontium during calcite nucleation and growth, and during reactive fluid flow; (2) consolidation behavior of calcite assemblages in the common brines. Numerically, we quantified water weakening of calcite using molecular dynamics simulations; and quantified the impact of calcite dissolution rate on macroscopic fracturing using finite element models. With microfluidic experiments and modeling, we show the effect of local flow fields on the dissolution kinetics of calcite. Taken together across a wide range of scales and methods, our studies allow us to separate the effects of reaction, flow, and transport, on calcite fracturing and the evolution of strontium isotopic signatures in the reactive fluids.

42 ENGINEERING↗
NASA NTRS1980

Optical phase distortion due to compressible flow over laser turrets

Analytical models for optical phase distortion due to compressible flow over a laser turret are developed. Phase distortion is calculated for both blunt and small perturbation turrets. For the blunt turret, the Janzen-Rayleigh technique is used to determine the flow field. Phase distortions of 2.2 wavelengths at 3.8 microns are calculated for the blunt turret. For small perturbation turrets, a versatile analytical model is developed for a turret on a fuselage with circular cross section. With a two dimensional Fourier series representation of the turret, any shape can be considered. Both subsonic and supersonic flows can be calculated. Phase distortions of 1.2 wavelengths at 3.8 microns are calculated for one turret at high subsonic Mach number. In addition to being of value for laser turrets, the methods are applicable to reconnaissance aircraft using photographic equipment and cruise missiles using celestial navigation.

Fuhs, A. E.↗
NASA NTRS1988

Interior impedance wedge diffraction with surface waves

The exact impedance wedge solution is evaluated asymptotically using the method of steepest descents for plane wave illumination at normal incidence. Uniform but different impedances on each face are considered for both soft and hard polarizations. The asymptotic solution isolates the incident, singly reflected, multiply reflected, diffracted, and surface wave fields. Multiply reflected fields of any order are permitted. The multiply reflected fields from the exact solution are written as ratios of auxiliary Maliuzhinets functions, whereas a geometrical analysis gives the reflected fields as products of reflection coefficients. These two representations are shown to be identical in magnitude, phase and the angular range over which they exist. The diffracted field includes four Fresnel transition functions as in the perfect conductor case, and the expressions for the appropriate discontinuities at the shadow boundaries are presented. The surface wave exists over a finite angular range and only for certain surface impedances. A surface wave transition field is included to retain continuity. Computations are presented for interior wedge diffractions although the formulation is valid for both exterior and interior wedges.

Balanis, Constantine A.↗
NASA NTRS1994

The Flux-integral Method for Multidimensional Convection and Diffusion

The flux-integral method is a procedure for constructing an explicit, single-step, forward-in-time, conservative, control volume update of the unsteady, multidimensional convection-diffusion equation. The convective plus diffusive flux at each face of a control-volume cell is estimated by integrating the transported variable and its face-normal derivative over the volume swept out by the convecting velocity field. This yields a unique description of the fluxes, whereas other conservative methods rely on nonunique, arbitrary pseudoflux-difference splitting procedures. The accuracy of the resulting scheme depends on the form of the subcell interpolation assumed, given cell-average data. Cellwise constant behavior results in a (very artificially diffusive) first-order convection scheme. Second-order convection-diffusion schemes correspond to cellwise linear (or bilinear) subcell interpolation. Cellwise quadratic subcell interpolants generate a highly accurate convection-diffusion scheme with excellent phase accuracy. Under constant-coefficient conditions, this is a uniformly third-order polynomial interpolation algorithm (UTOPIA).

Leonard, B. P.↗
NASA NTRS2003

Effects of Convective Transport of Solute and Impurities on Defect-Causing Kinetics Instabilities in Protein Crystallization

Insight into the crystallization processes of biological macromolecules into crystals or aggregates can provide valuable guidelines in many fundamental and applied fields. Such insight will prompt new means to regulate protein phase transitions in-vivo, e.g., polymerization of hemoglobin S in the red cells, crystallization of crystallins in the eye lens, etc. Understanding of protein crystal nucleation will help achieve narrow crystallite size distributions, needed for sustained release of pharmaceutical protein preparations such as insulin or interferon. Traditionally, protein crystallization studies have been related to the pursuit of crystal perfection needed to improve the structure details provided by x-ray, electron or neutron diffraction methods. Crystallization trials for the purposes of structural biology carried out in space have posed an intriguing question related to the inconsistency of the effects of the microgravity growth on the quality of the crystals.

Vekilov, Peter G.↗
NASA NTRS2016

Wave Telescope Technique for MMS Magnetometer

Multipoint measurements are a powerful method in studying wavefields in space plasmas.The wave telescope technique is tested against magnetic field fluctuations in the terrestrial magnetosheath measured by the four Magnetospheric Multiscale (MMS) spacecraft on a spatial scale of about 20 km.The dispersion relation diagram and the wave vector distribution are determined for the first time in the ion-kinetic range. Moreover, the dispersion relation diagram is determined in a proxy plasma restframe by regarding the low-frequency dispersion relation as a Doppler relation and compensating for the apparent phase velocity. Fluctuations are highly compressible, and the wave vectors have an angle of about 60 from the mean magnetic field. We interpret that the measured fluctuations represent akinetic-drift mirror mode in the magnetosheath which is dispersive and in a turbulent state accompanied by a sideband formation.

Magnetometer↗
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

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

Block Mott insulating state induced by next-nearest-neighbor hopping in the S = 3 2 zigzag chain BaCoTe 2 O 7

Quasi-one-dimensional correlated electronic multiorbital systems with either ladder or chain geometries continue attracting considerable interest due to their complex electronic phases arising from the interplay of the hopping matrix, the crystal-field splitting, the electronic correlations (Hubbard repulsion U and Hund coupling J H ), and strong quantum fluctuations. Recently, the intriguing cobalt zigzag chain system BaCoTe 2 ⁢ O 7 , with electronic density n=7, was prepared experimentally. Here, we systematically study the electronic and magnetic properties of this quasi-one-dimensional compound from the theoretical perspective. Based on first-principles density functional theory calculations, strongly anisotropic one-dimensional electronic Co 3⁢d bands were found near the Fermi level. By evaluating the relevant hopping amplitudes, we provide the magnitude and origin of the nearest-neighbor (NN) and next-nearest-neighbor (NNN) hopping matrices in BaCoTe 2 ⁢ O 7 . With this information, we constructed a three-orbital electronic Hubbard model for this zigzag chain system, and studied two cases: with only a NN hopping matrix, and with NN plus NNN hopping matrices. Introducing the Hubbard and Hund couplings and studying the model via the density matrix renormalization group method, we constructed the ground-state phase diagram. A robust staggered ↑ - ↓ - ↑ - ↓ antiferromagnetic (AFM) region was found when only the NN hopping matrix in the chain direction was employed. However, for the realistic case where the NNN hopping matrix is also included, the dominant state becomes instead a block AFM ↑ - ↑ - ↓ - ↓ order, in agreement with experiments. The system displays Mott insulator characteristics with three half-filled orbitals, when the block AFM order is stable. In conclusion, our results for BaCoTe 2 ⁢O 7 provide guidance to experimentalists and theorists working on this zigzag one-dimensional chain and related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗