Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

Measurements and Model Improvement: Insight into NWP Model Error Using Doppler Lidar and Other WFIP2 Measurement Systems

Abstract Doppler-lidar wind-profile measurements at three sites were used to evaluate NWP model errors from two versions of NOAA’s 3-km-grid HRRR model, to see whether updates in the latest version 4 reduced errors when compared against the original version 1. Nested (750-m grid) versions of each were also tested to see how grid spacing affected forecast skill. The measurements were part of the field phase of the Second Wind Forecasting Improvement Project (WFIP2), an 18-month deployment into central Oregon–Washington, a major wind-energy-producing region. This study focuses on errors in simulating marine intrusions, a summertime, 600–800-m-deep, regional sea-breeze flow found to generate large errors. HRRR errors proved to be complex and site dependent. The most prominent error resulted from a premature drop in modeled marine-intrusion wind speeds after local midnight, when lidar-measured winds of greater than 8 m s −1 persisted through the next morning. These large negative errors were offset at low levels by positive errors due to excessive mixing, complicating the interpretation of model “improvement,” such that the updates to the full-scale versions produced mixed results, sometimes enhancing but sometimes degrading model skill. Nesting consistently improved model performance, with version 1’s nest producing the smallest errors overall. HRRR’s ability to represent the stages of sea-breeze forcing was evaluated using radiation budget, surface-energy balance, and near-surface temperature measurements available during WFIP2. The significant site-to-site differences in model error and the complex nature of these errors mean that field-measurement campaigns having dense arrays of profiling sensors are necessary to properly diagnose and characterize model errors, as part of a systematic approach to NWP model improvement. Significance Statement Dramatic increases in NWP model skill will be required over the coming decades. This paper describes the role of major deployments of accurate profiling sensors in achieving that goal and presents an example from the Second Wind Forecast Improvement Program (WFIP2). Wind-profile data from scanning Doppler lidars were used to evaluate two versions of HRRR, the original and an updated version, and nested versions of each. This study focuses on the ability of updated HRRR versions to improve upon predicting a regional sea-breeze flow, which was found to generate large errors by the original HRRR. Updates to the full-scale HRRR versions produced mixed results, but the finer-mesh versions consistently reduced model errors.

Meteorology & Atmospheric Sciences↗

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

Robust and ultrafast state preparation by ramping artificial gauge potentials

Abstract The implementation of static artificial magnetic fields in ultracold atomic systems has become a powerful tool, e.g. for simulating quantum-Hall physics with charge-neutral atoms. Taking an interacting bosonic flux ladder as a minimal model, we investigate protocols for adiabatic state preparation via magnetic flux ramps. Considering the fact that it is actually the artificial vector potential (in the form of Peierls phases) that can be experimentally engineered in optical lattices, rather than the magnetic field, we find that the time required for adiabatic state preparation dramatically depends on which pattern of Peierls phases is used. This can be understood intuitively by noting that different patterns of time-dependent Peierls phases that all give rise to the same magnetic field ramp, generally lead to different artificial electric fields during the ramp. As an intriguing result, we find that an optimal choice allows for preparing the ground state almost instantaneously in the non-interacting system, which can be related to the concept of counterdiabatic driving. Remarkably, we find extremely short preparation times also in the strongly-interacting regime. Our findings open new possibilities for robust state preparation in atomic quantum simulators.

Wang, Botao (ORCID:0000000282202452)↗

A Landau–Devonshire analysis of strain effects on ferroelectric Al 1-x Sc x N

We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials using the classic Landau–Devonshire approach. Electrostrictive and dielectric stiffness coefficients of Al 1-x Sc x N with a wurtzite structure (6 mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first-order phase transition based on the dielectric stiffness relationships. The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics. These analysis results reported here explain experimentally reported sensitivity of the coercive field to elastic strain/stress in Al 1-x Sc x N films and would enable further thermodynamic analysis via phase field simulation and related methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Methodology of the Extraction of Multi-Differential Cross Sections of Charged-Current $\nu_\mu$-Argon Interactions in MicroBooNE Using Wire-Cell $\nu_\mu$CC Selection

Neutrino physics experiments rely on accurate nuclear-interaction models that are in part guided by experimental observations. Multi-differential cross section measurements are particularly valuable for exploring the underlying physics described in these models, making them a significant step forward in the field. This note presents methodology in preparation to unfold neutrino flux-averaged double and triple-differential cross sections of the inclusive muon neutrino charged-current interaction on Argon. This work builds upon the existing framework and methodology presented in recent energy dependent cross section measurements at MicroBooNE using the Wire-Cell tomographic event reconstruction. The signal definition, choice of binning, and handling of estimated detector uncertainties are discussed, and the MicroBooNE simulation model is validated over the multi-dimensional phase space of muon momentum, muon polar angle, and visible hadronic energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fractons from Vector Gauge Theory

Motivated by the prediction of fractonic topological defects in a quantum crystal, we utilize a reformulated elasticity duality to derive a description of a fracton phase in terms of coupled vector U(1) gauge theories. The fracton order and restricted mobility emerge as a result of an unusual Gauss law where electric field lines of one gauge field act as sources of charge for others. At low energies this vector gauge theory reduces to the previously studied fractonic symmetric tensor gauge theory. Here, we construct the corresponding lattice model and a number of generalizations, which realize fracton phases via a condensation of stringlike excitations built out of charged particles, analogous to the p-string condensation mechanism of the gapped X-cube fracton phase.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Photoelectron spectra in circularly and elliptically polarized laser pulses

We present results of numerical simulations of the time-dependent Schrödinger equation and theoretical analysis concerning the interaction of a rare gas atom with circularly and elliptically polarized laser pulses. In agreement with recent observation in circularly polarized fields, the photoelectron energy spectra for counterrotating electrons are peaked at lower kinetic energy than those for corotating electrons. We show that this difference can be interpreted as being due to the additional pathways to ionization that are available for the counterrotating electrons only. Furthermore, our results show, in agreement with earlier work, that the offset angle by which the emission of electrons is rotated in an elliptically polarized field increases with each successive above-threshold ionization (ATI) order and is larger for the emission from counterrotating states as compared to that from corotating states. As a result, a simple model based on the interference and relative phase difference between just three continuum states provides remarkable agreement and once again emphasizes the importance of the additional ionization pathways available for the counterrotating electrons.

74 ATOMIC AND MOLECULAR PHYSICS↗

Complex magnetic phases in the polar tetragonal intermetallic NdCoGe 3

Polar materials can host a variety of topologically significant magnetic phases, which often emerge from a modulated magnetic ground state. Relatively few noncentrosymmetric tetragonal materials have been shown to host topological spin textures and new candidate materials are necessary to expand the current theoretical models. This manuscript reports on the anisotropic magnetism in the polar, tetragonal material NdCoGe 3 via thermodynamic and neutron diffraction measurements. The previously reported H - T phase diagram is updated to include several additional phases, which exist for both H = 0 and with an applied field H ⊥ c . Neutron diffraction data reveal that the magnetic structures below T N 1 = 3.70 K and T N 2 = 3.50 K are incommensurate, with a ground-state magnetic order that is incommensurate in all directions with the propagation vector $\vec{k}$ = (0.494, 0.0044, 0.385) at 1.8 K. A unique magnetic structure solution is not achievable, but the possible single-$\vec{k}$ and $\vec{k}$ spin models are discussed. These results demonstrate that NdCoGe 3 hosts complicated magnetic order derived from modulated magnetic moments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗

Topological or rotational non-Abelian gauge fields from Einstein-Skyrme holography

We report analytically known states at non-zero temperature which may serve as a powerful tool to reveal common topological and thermodynamic properties of systems ranging from the QCD phase diagram to topological phase transitions in condensed matter materials. In the holographically dual gravity theory, these are analytic solutions to a five-dimensional non-linear-sigma (Skyrme) model dynamically coupled to Einstein gravity. This theory is shown to be holographically dual to \( \mathcal{N} \) = 4 Super-Yang-Mills theory coupled to an SU(2)-current. All solutions are fully backreacted asymptotically Anti-de Sitter (AdS) black branes or holes. One family of global AdS black hole solutions contains non-Abelian gauge field configurations with positive integer Chern numbers and finite energy density. Larger Chern numbers increase the Hawking-Page transition temperature. In the holographically dual field theory this indicates a significant effect on the deconfinement phase transition. Black holes with one Hawking temperature can have distinct Chern numbers, potentially enabling topological transitions. A second family of analytic solutions, rotating black branes, is also provided. These rotating solutions induce states with propagating charge density waves in the dual field theory. We compute the Hawking temperature, entropy density, angular velocity and free energy for these black holes/branes. These correspond to thermodynamic data in the dual field theory. For these states the energy-momentum tensor, (non-)conserved current, and topological charge are interpreted.

79 ASTRONOMY AND ASTROPHYSICS↗

Solid–liquid transition in a skyrmion matter

We report Monte-Carlo studies of the orientational order and melting of a 2D skyrmion lattice containing more than one million spins. Two models have been investigated, a microscopic model of lattice spins with Dzyaloshinskii–Moryia interaction that possesses skyrmions, and the model in which skyrmions are treated as point particles with repulsive interaction derived from a spin model. They produce similar results. The skyrmion lattice exhibits a sharp one-step transition between solid and liquid phases on temperature and the magnetic field. This solid–liquid transition is characterized by the kink in the magnetization. Here, the field-temperature phase diagram is computed. We show that the application of the field gradient to a 2D system of skyrmions produces a solid–liquid interface that must be possible to observe in experiments.

2D melting↗

Modeling the Simultaneous Dropout of Energetic Electrons and Protons by Magnetopause Shadowing

Abstract Magnetopause shadowing (MPS) effect could drive a concurrent dropout of radiation belt electrons and ring current protons. However, its relative role in the dropout of both plasma populations has not been well quantified. In this work, we study the simultaneous dropout of MeV electrons and 100s keV protons during an intense geomagnetic storm in May 2017. A radial diffusion model with an event‐specific last closed drift shell is used to simulate the MPS loss of both populations. The model well captures the fast shadowing loss of both populations at L * > 4.6, while the loss at L * < 4.6, possibly due to the electromagnetic ion cyclotron wave scattering, is not captured. The observed butterfly pitch angle distributions of electron fluxes in the initial loss phase are well reproduced by the model. The initial proton losses at low pitch angles are underestimated, potentially also contributed by other mechanisms such as field line curvature scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Theoretical description of atomtronic Josephson junctions in an optical lattice

Experimental realizations of “atomtronic” Josephson junctions have recently been created in annular traps in relative rotation with respect to potential barriers that generate the weak links. If these devices are additionally subjected to an optical lattice potential, then they can incorporate strong-coupling Mott physics within the design, which can modify the behavior and can allow for interesting new configurations of barriers and of superfluid flow patterns. We examine theoretically the behavior of a Bose superfluid in an optical lattice in the presence of an annular trap and a barrier across the annular region which acts as a Josephson junction. As the superfluid is rotated, circulating super-currents appear. Beyond a threshold superfluid velocity, phase slips develop, which generate vortices. We use a finite temperature strong-coupling expansion about the mean-field solution of the Bose Hubbard model to calculate various properties of such devices. Additionally, we discuss some of the rich behavior that can result when there are Mott regions within the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental evidence for the two-path description of neutron spin echo

Here, we describe an experiment that strongly supports a two-path interferometric model in which the spin-up and spin-down components of each neutron propagate coherently along spatially separated parallel paths in a typical neutron spin-echo small-angle scattering (SESANS) experiment. Specifically, we show that the usual semi-classical, single-path treatment of Larmor precession of a polarized neutron in an external magnetic field predicts a damping as a function of the spin-echo length of the SESANS signal obtained with a periodic phase grating when the transverse width of the neutron wave packet is finite. However, no such damping is observed experimentally, implying either that the Larmor model is incorrect or that the transverse extent of the wave packet is very large. In contrast, we demonstrate theoretically that a quantum-mechanical interferometric model in which the two mode-entangled (i.e., intraparticle entangled) spin states of a single neutron are separated in space when they interact with the grating accurately predicts the measured SESANS signal, which is independent of the wave packet width.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Field Induced Density Wave in a Kagome Superconductor

On the kagome lattice, electrons benefit from the simultaneous presence of band topology, flat electronic bands, and Van Hove singularities, forming competing or cooperating orders. Understanding the interrelation between these distinct order parameters remains a significant challenge, leaving much of the associated physics unexplored. In the kagome superconductor KV 3 Sb 5 , which exhibits a charge density wave (CDW) state below T ≃ 78 K, we uncover an unpredicted field-induced phase transition below 6 K. The observed transition is marked by a hysteretic anomaly in the resistivity, nonlinear electrical transport, and a change in the symmetry of the electronic response as probed via the angular dependence of the magnetoresistivity. These observations surprisingly suggest the emergence of an unanticipated broken symmetry state coexisting with the original CDW. To understand this experimental observation, we developed a theoretical minimal model for the normal state inside the high-temperature parent CDW phase, where an incommensurate CDW order emerges as an instability subleading to superconductivity. Here, the incommensurate CDW emerges when superconducting fluctuations become fully suppressed by large magnetic fields. Our results suggest that, in kagome superconductors, quantum states can either coexist or are nearly degenerate in energy, indicating that these are rich platforms to expose new correlated phenomena.

Charge density waves↗

Development and testing of a solar flux heating freeze recovery system for molten salt parabolic troughs

A major obstacle preventing commercialization of molten salt parabolic trough plants is freeze protection and recovery of the solar field. One of the most common nitrate salt compositions under consideration allows for solar field outlet temperatures up to 565°C and begins solidifying above 240°C. Using molten salt directly in the solar field can significantly reduce thermal energy storage costs and the higher temperature compared to oil heat transfer fluids can power a more efficient power cycle, but the very high freezing point is cause for concern. This work builds on a previous modeling effort investigating different methods for melting Solar Salt frozen in the solar field, which suggested the viability of a novel solar heating method to lower costs of the freeze protection system. In this work higher fidelity models are developed to better investigate the melting process with a more detailed 3D geometry which includes the insulated bellows. The thermal-fluid model of the melting process is updated to include temperature dependent density and data for solid phase properties. In addition, a finite element model is developed to resolve thermal stresses at the point with highest thermal gradients. Results with the new model confirm the viability of using controllable solar flux heating to thaw salt frozen in the solar field without damaging the receiver; however, the presence of non-illuminated sections significantly slows down the melting process. Solar heating simulations with pauses off-sun of 120 and 60 s require 18.5 and 11.9 h to melt salt from a night-time temperature of 10°C with 1000 W m2 DNI, which results in maximum thermal stresses of 48.2 and 53.5 MPa, respectively. Adding 150 W m−1 to the solar heating simulation with a 60 s pause off-sun decreases both the melt time and thermal stresses due to the more uniform heating, suggesting a combined heating method may be the best option.

14 SOLAR ENERGY↗

Circuit complexity near critical points

Here, we consider the Bose–Hubbard model in two and three spatial dimensions and numerically compute the quantum circuit complexity of the ground state in the Mott insulator and superfluid phases using a mean field approximation with additional quadratic fluctuations. After mapping to a qubit system, the result is given by the complexity associated with a Bogoliubov transformation applied to the reference state taken to be the mean field ground state. In particular, the complexity has peaks at the O(2) critical points where the system can be described by a relativistic quantum field theory. Given that we use a Gaussian approximation, near criticality the numerical results agree with a free field theory calculation. To go beyond the Gaussian approximation we use general scaling arguments that imply that, as we approach the critical point t → t c , there is a non-analytic behavior in the complexity c 2 (t) of the form |c 2 (t) – c 2 (t c )| ~ |t – t c | νd , up to possible logarithmic corrections. Here d is the number of spatial dimensions and ν is the usual critical exponent for the correlation length ξ ~ |t – t c | –ν . As a check, for d = 2 this agrees with the numerical computation if we use the Gaussian critical exponent $v=\frac{1}{2}$. Finally, using AdS/CFT methods, we study higher dimensional examples and confirm this scaling argument with non-Gaussian exponent ν for strongly interacting theories that have a gravity dual.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗