Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Axion detection via superfluid 3 He ferromagnetic phase and quantum measurement techniques

We propose to use the nuclear spin excitation in the ferromagnetic A 1 phase of the superfluid 3 He for the axion dark matter detection. This approach is striking in that it is sensitive to the axion-nucleon coupling, one of the most important features of the QCD axion introduced to solve the strong CP problem. We review a quantum mechanical description of the nuclear spin excitation and apply it to the estimation of the axion-induced spin excitation rate. We also describe a possible detection method of the spin excitation in detail and show that the combination of the squeezing of the final state with the Josephson parametric amplifier and the homodyne measurement can enhance the sensitivity. It turns out that this approach gives good sensitivity to the axion dark matter with the mass of $\mathcal{O}$(1) μeV depending on the size of the external magnetic field. We estimate the parameters of experimental setups, e.g., the detector volume and the amplitude of squeezing, required to reach the QCD axion parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Moment method and continued fraction expansion in Floquet operator Krylov space

Recursion methods such as Krylov techniques map complex dynamics to an effective noninteracting problem in one dimension. For example, the operator Krylov space for Floquet dynamics can be mapped to the dynamics of an edge operator of the one-dimensional Floquet inhomogeneous transverse field Ising model (ITFIM), where the latter, after a Jordan-Wigner transformation, is a Floquet model of noninteracting Majorana fermions and the couplings correspond to Krylov angles. We present an application of this showing that a moment method exists where given an autocorrelation function, one can construct the corresponding Krylov angles and from that the corresponding Floquet ITFIM. Consequently, when no solutions for the Krylov angles are obtained, it indicates that the autocorrelation is not generated by unitary dynamics. We highlight this by studying certain special cases: stable m-period dynamics derived using the method of continued fractions, exponentially decaying, and power-law decaying stroboscopic dynamics. Remarkably, our examples of stable m-period dynamics correspond to m-period edge modes for the Floquet ITFIM where, deep in the chain, the couplings correspond to a critical phase. Furthermore our results pave the way to engineer Floquet systems with desired properties of edge modes and also provide examples of persistent edge modes in gapless Floquet systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Generalized boost transformations in finite volumes and application to Hamiltonian methods

The investigation of hadron interactions within lattice QCD has been facilitated by the well-known quantisation condition, linking scattering phase shifts to finite-volume energies. Additionally, the ability to utilise systems at finite total boosts has been pivotal in smoothly charting the energy-dependent behaviour of these phase shifts. The existing implementations of the quantization condition at finite boosts rely on momentum transformations between rest and moving frames, defined directly in terms of the energy eigenvalues. This energy dependence is unsuitable in the formulation of a Hamiltonian. In this work, we introduce a novel approach to generalise the three-momentum boost prescription, enabling the incorporation of energy-independent finite-volume Hamiltonians within moving frames. We demonstrate the application of our method through numerical comparisons, employing a phenomenological ππ scattering example.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Critical fluid dynamics in two and three dimensions

We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Bridging length scales in hard materials with ultra-small angle X-ray scattering – a critical review

Owing to their exceptional properties, hard materials such as advanced ceramics, metals and composites have enormous economic and societal value, with applications across numerous industries. Understanding their microstructural characteristics is crucial for enhancing their performance, materials development and unleashing their potential for future innovative applications. However, their microstructures are unambiguously hierarchical and typically span several length scales, from sub-ångstrom to micrometres, posing demanding challenges for their characterization, especially for in situ characterization which is critical to understanding the kinetic processes controlling microstructure formation. This review provides a comprehensive description of the rapidly developing technique of ultra-small angle X-ray scattering (USAXS), a nondestructive method for probing the nano-to-micrometre scale features of hard materials. USAXS and its complementary techniques, when developed for and applied to hard materials, offer valuable insights into their porosity, grain size, phase composition and inhomogeneities. We discuss the fundamental principles, instrumentation, advantages, challenges and global status of USAXS for hard materials. Using selected examples, we demonstrate the potential of this technique for unveiling the microstructural characteristics of hard materials and its relevance to advanced materials development and manufacturing process optimization. We also provide our perspective on the opportunities and challenges for the continued development of USAXS, including multimodal characterization, coherent scattering, time-resolved studies, machine learning and autonomous experiments. Our goal is to stimulate further implementation and exploration of USAXS techniques and inspire their broader adoption across various domains of hard materials science, thereby driving the field toward discoveries and further developments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Al–W gradient density materials—Processing and dynamic ramp compression

Materials with high-density gradients are desired for controlling loading paths in dynamic compression, important for studying material properties in extreme conditions and inertial confinement fusion. The large density difference between Al and W makes them ideal choices for producing gradient density materials, but their extremely different melting temperatures make them challenging to fabricate simultaneously. We report a method for producing Al–W porosity-free materials with a fourfold increase in density (2.7–11 g/cm 3 ) across the composition range, from Al-rich to W-rich, without intermetallic phase formation. This was achieved by understanding the aluminum-dominated densification behavior and examining the influence of pressure and temperature on the densification of Al–W composites. Dynamic compression experiments conducted with the Al–W gradient density material produced shock ramp compressions as expected based on the designed composition, and the performed hydrodynamics simulations showed excellent agreement with experimental results. The results demonstrate that current activated pressure-assisted densification allows for the easy and rapid fabrication of gradient density materials with significant density gradients and tailored compositions, facilitating precise control of the loading paths. These materials have the potential to create customized pressure drives for advancing the fields of material science in extreme environments and dynamic compression.

Alloys↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quantum Thermodynamics of Nonequilibrium Processes in Lattice Gauge Theories

A key objective in nuclear and high-energy physics is to describe nonequilibrium dynamics of matter, e.g., in the early Universe and in particle colliders, starting from the standard model of particle physics. Classical computing methods, via the framework of lattice gauge theory, have experienced limited success in this mission. Quantum simulation of lattice gauge theories holds promise for overcoming computational limitations. Because of local constraints (Gauss’s laws), lattice gauge theories have an intricate Hilbert-space structure. This structure complicates the definition of thermodynamic properties of systems coupled to reservoirs during equilibrium and nonequilibrium processes. We show how to define thermodynamic quantities such as work and heat using strong-coupling thermodynamics, a framework that has recently burgeoned within the field of quantum thermodynamics. Our definitions suit instantaneous quenches, simple nonequilibrium processes undertaken in quantum simulators. To illustrate our framework, we compute the work and heat exchanged during a quench in a Z 2 lattice gauge theory coupled to matter in 1+1 dimensions. Here, the thermodynamic quantities, as functions of the quench parameter, evidence a phase transition. For general thermal states, we derive a simple relation between a quantum many-body system’s entanglement Hamiltonian, measurable with quantum-information-processing tools, and the Hamiltonian of mean force, used to define strong-coupling thermodynamic quantities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗
DOE OSTI2024

Enabling resonant ultrasound spectroscopy in high magnetic fields

Resonant ultrasound spectroscopy (RUS) is a powerful method to determine elastic constants with high accuracy and precision from a single measurement of the mechanical resonances of a sample. Conventionally, the quantitative extraction of elastic moduli with RUS assumes free boundary conditions which can often lead to the adoption of unstable sample positioning between ultrasonic transducers that is incompatible with extreme environments like high magnetic fields. We show that, under specific conditions, introducing a small amount of adhesive between a RUS sample and ultrasonic transducers introduces a perturbation to the free resonance condition which can be accounted for by a simple model. This means elastic constants can be determined to within the uncertainty of conventional RUS, but with significant improvements including sample stability and control of sample orientation. We demonstrate the efficacy of this approach with measurements on a range of materials including room temperature measurements on polycrystalline metals, temperature-dependent measurements of the structural phase transition in strontium titanate single crystals, and magnetic field-dependent measurements of magnetic phase transitions in gadolinium polycrystals up to 14 T.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Improving generative model-based unfolding with Schrödinger bridges

Machine learning-based unfolding has enabled unbinned and high-dimensional differential cross section measurements. Two main approaches have emerged in this research area; one based on discriminative models and one based on generative models. The main advantage of discriminative models is that they learn a small correction to a starting simulation while generative models scale better to regions of phase space with little data. We propose to use Schrödinger bridges and diffusion models to create , an unfolding approach that combines the strengths of both discriminative and generative models. The key feature of is that its generative model maps one set of events into another without having to go through a known probability density as is the case for normalizing flows and standard diffusion models. We show that achieves excellent performance compared to state of the art methods on a synthetic Z + jets dataset. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A bi-channel aided stitching of atomic force microscopy images

Microscopy is an essential tool in scientific research, enabling the visualization of structures at micro- and nanoscale resolutions. However, the field of microscopy often encounters limitations in field-of-view (FOV), restricting the amount of sample that can be imaged in a single capture. To overcome this limitation, image stitching techniques have been developed to seamlessly merge multiple overlapping images into a single, high-resolution composite. The images collected from microscope need to be optimally stitched before accurate physical information can be extracted from post analysis. However, the existing stitching tools either struggle to stitch images together when the microscopy images are feature sparse or cannot address all the transformations of images when performing image stitching. To address these issues, we propose a bi-channel aided feature-based image stitching method and demonstrate its use on Atomic Force Microscopy (AFM) generated Pantoea sp. YR343 biofilm and PTO thin film sample images as experimental data. The topographical channel image of AFM data captures the morphological details of the sample, and a stitched topographical image is desired for researchers. We utilize the amplitude and phase channels of AFM data to maximize the matching features and to estimate the position of the original topographical images and show that the proposed bi-channel aided stitching method outperforms the traditional direct stitching approach in AFM topographical image stitching task. Here, we demonstrated the application on AFM, but similar approaches could be employed of optical microscopy with brightfield and fluorescence channels. We believe this proposed workflow can serve as a valuable augmentation strategy for microscopy image stitching tasks and will benefit the experimentalist to avoid erroneous analysis and discovery due to incorrect stitching.

Atomic force microscopy↗
DOE OSTI2025

Nontrivial critical behavior at magnetic transitions: A case study of Sm 7 ⁢Pd 3

We present a comprehensive analysis of the critical behavior of Sm 7 ⁢Pd 3 in the vicinity of its second-order magnetoelastic transition at 𝑇 c =173 K. The critical exponents (CEs) 𝛽 and 𝛾, determined using both the standard convergence procedure and the average normalized slope (ANS) method, diverge at 𝑇 c –-a characteristic typically associated with first-order transitions. Notably, none of the established universality classes satisfactorily describe the critical behavior of Sm 7 ⁢Pd 3 , and we discuss the possible origins of this deviation in the context of the strong spin-lattice coupling intrinsic to the sample. We emphasize the importance of accurately selecting the critical temperature and magnetic field ranges to ensure robust critical behavior analysis, and propose a quantitative approach to assess the reliability of the extracted CEs. Additionally, we demonstrate that in the ANS method, the critical exponents 𝛽 and 𝛾 should be calculated separately using data for 𝑇 ⩽ 𝑇 c and 𝑇 ⩾ 𝑇 c , respectively. In conclusion, our findings underscore the need for a revised theoretical framework to accurately describe second-order magnetoelastic transitions.

Ferrimagnets↗
DOE OSTI2025

Accelerating the discovery of low-energy structure configurations: A computational approach that integrates first-principles calculations, Monte Carlo sampling, and Machine Learning

Finding Minimum Energy Configurations (MECs) is essential in fields such as physics, chemistry, and materials science, as they represent the most stable states of the systems. In particular, identifying such MECs in multi-component alloys considered candidate PFMs is key because it determines the most stable arrangement of atoms within the alloy, directly influencing its phase stability, structural integrity, and thermo-mechanical properties. However, since the search space grows exponentially with the number of atoms considered, obtaining such MECs using computationally expensive first-principles DFT calculations often results in a cumbersome task. To escape the above compromise between physical fidelity and computational efficiency, we have developed a novel physics-based data-driven approach that combines Monte Carlo sampling, first-principles DFT calculations, and Machine Learning to accelerate the discovery of MECs in multi-component alloys. More specifically, we have leveraged well-established Cluster Expansion (CE) techniques with Local Outlier Factor models to establish strategies that enhance the reliability of the CE method. In this work, we demonstrated the capabilities of the proposed approach for the particular case of a tungsten-based quaternary high-entropy alloy. However, the method is applicable to other types of alloys and enables a wide range of applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Emulating ab initio computations of infinite nucleonic matter

We construct efficient emulators for the computation of the infinite nuclear matter equation of state. These emulators are based on the subspace-projected coupled-cluster method for which we here develop a new algorithm called small-batch voting to eliminate spurious states that might appear when emulating quantum many-body methods based on a non-Hermitian Hamiltonian. The efficiency and accuracy of these emulators facilitate a rigorous statistical analysis within which we explore nuclear matter predictions for > 10 6 different parametrizations of a chiral interaction model with explicit Δ -isobars at next-to-next-to leading order. Constrained by nucleon-nucleon scattering phase shifts and bound-state observables of light nuclei up to He 4 , we use history matching to identify nonimplausible domains for the low-energy coupling constants of the chiral interaction. Within these domains we perform a Bayesian analysis using sampling and importance resampling with different likelihood calibrations and study correlations between interaction parameters, calibration observables in light nuclei, and nuclear matter saturation properties. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Machine-Assisted Discovery of Integrable Symplectic Mappings

Integrable systems possess a hidden symmetry associated with the existence of conserved quantities known as integrals of motion. These systems play an important role in understanding general dynamics in accelerators and have potential for future designs. This work will cover two automated methods for finding integrable symplectic maps of the plane. The first algorithm is based on the observation that the evolution of an integrable system in phase space is confined to a lower-dimensional submanifold of a specific type. The second algorithm relies on an analysis of dynamical variables. Both methods rediscover some of the famous McMillan-Suris integrable mappings and ultra-discrete Painlev\'e equations. Over 100 new integrable families are presented and analyzed, some of which are isolated in the space of parameters, while others are families with one parameter (or the ratio of parameters) being either continuous or discrete. In addition, the newly discovered maps are related to a general 2D symplectic map through the use of discrete perturbation theory. A method is proposed for constructing smooth near-integrable dynamical systems based on mappings with polygon invariants.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Accelerating Resonance Searches via Signature-Oriented Pre-training

The search for heavy resonances beyond the Standard Model (BSM) is a key objective at the LHC. While the recent use of advanced deep neural networks for boosted-jet tagging significantly enhances the sensitivity of dedicated searches, it is limited to specific final states, leaving vast potential BSM phase space underexplored. We introduce a novel experimental method, Signature-Oriented Pre-training for Heavy-resonance ObservatioN (Sophon), which leverages deep learning to cover an extensive number of boosted final states. Pre-trained on the comprehensive JetClass-II dataset, the Sophon model learns intricate jet signatures, ensuring the optimal constructions of various jet tagging discriminates and enabling high-performance transfer learning capabilities. We show that the method can not only push widespread model-specific searches to their sensitivity frontier, but also greatly improve model-agnostic approaches, accelerating LHC resonance searches in a broad sense.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Disentangling High Harmonic Generation from Surface and Bulk States of a Topological Insulator

The discovery of topological phases has introduced a new dimension to materials science. Three-dimensional (3D) topological insulators (TIs) are a remarkable class of matter that is insulating in the bulk while hosting conductive topological surface states (TSSs) with unique charge and spin properties. High-order harmonic generation (HHG) has emerged as a powerful tool to probe condensed matter systems by providing insights into their electronic structure and dynamic behavior. Here, we investigate HHG in the prototype 3D-TI Bi$_2$Se$_3$. We demonstrate that the contributions of bulk and surface states to the harmonic emission can be controlled by tuning the thickness of thin film samples. An ultrathin (6 nm) film substantially enhances HHG from the surface states, while the bulk states dominate HHG in a thicker (50 nm) film. By applying a quasi-static terahertz perturbing field, we disentangle the bulk and surface responses and reveal the significant impact of the surface states' shift vector and Berry curvature on HHG. Our study provides effective methods for isolating the optical responses of TSSs from those of the bulk, which opens the door to resolving an ongoing debate regarding whether it is possible to reliably extract topological signatures in HHG.

Atomic Physics (physics.atom-ph)↗