Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Landau equation”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Phase field dislocation dynamics formulation coupled with Fourier based micromechanics solver and its application to grain boundary–dislocation interactions

A new phase field dislocation dynamics (PFDD) formulation for homogeneous and heterogeneous materials is presented, which couples micromechanical solvers and the time-dependent Ginzburg–Landau equation. The strain fields are obtained from the micromechanical solver by solving the Lippmann–Schwinger equation and then used to define energy terms to model the evolution of the dislocations. Grain boundary (GB)–dislocation interactions are studied using the coupled PFDD formulation and by describing GBs as inclusions. GB energy and stiffness tensors are computed from molecular statics simulations, and a newly proposed lattice energy term that is dependent on the GB energy is considered in the calculations. Interaction of a screw dislocation with minimum energy and metastable states of low and high angle ⟨110⟩ symmetric tilt grain boundaries are studied. We show good agreement between predictions from our PFDD formulation and molecular dynamics simulations of grain boundary–dislocation interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Instabilities of the normal state in current-biased narrow superconducting strips

We study the current-voltage characteristic of narrow superconducting strips in the gapless regime near the critical temperature in the framework of the Ginzburg-Landau model. Our focus is on its instabilities occurring at high current biases. The latter are consequences of dynamical states with periodic phase-slip events in space and time. We analyze their structure and derive the value of the reentrance current at the onset of the instability of the normal state. It is expressed in terms of the kinetic coefficient of the time-dependent Ginzburg-Landau equation and calculated numerically.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Superconducting contact and quantum interference between two-dimensional van der Waals and three-dimensional conventional superconductors

Atomically thin two-dimensional (2D) transition-metal dichalcogenide (TMD) superconductors enable uniform, flat, and clean van der Waals tunneling interfaces, motivating their integration into conventional superconducting circuits. However, fully superconducting contact must be made between the 2D material and three-dimensional (3D) superconductors to employ the standard microwave drive and readout of qubits in such circuits. We present a method for creating zero-resistance contacts between 2D NbSe 2 and 3D aluminum that behave as Josephson junctions (JJs) with large effective areas compared to 3D-3D JJs. The devices formed from 2D TMD superconductors are strongly influenced by the geometry of the flakes themselves as well as the placement of the contacts to bulk 3D superconducting leads. We present a model for the supercurrent flow in a 2D-3D superconducting structure by a numerical solution of the Ginzburg-Landau equations and find good agreement with experiment. Finally, these results demonstrate a crucial step towards a new generation of hybrid superconducting quantum circuits.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗
DOE OSTI2021

NUMERICAL CALCULATION OF LOSSES OF TRAPPED VORTICES UNDER STRONG RF MEISSNER CURRENT AND DC SUPERHEATING FIELD IN TYPE II SUPERCONDUCTORS

Research on the vortex dynamics and enhancing of superheating field in superconductors has attracted much attention in accelerator physics community to develop next-generation high-performance accelerator cavities. However, the extreme dynamics of curvilinear elastic vortices driven by very strong currents close to the depairing limit or superheating field of a superconductor with a nanostructured surface has not been well understood. We calculated the superheating field Hsh and critical momentum kc characterizing the wavelength of the instability ?m of the Meissner state to flux penetration by solving numerically the Ginzburg?Landau equations. A bulk superconductor, superconductor with the inhomogeneous surface disorder (S-S), and multilayered surface (S-I-S) have been thoroughly investigated in this work. Our result showed that S-S and S-I-S structures can enhance the superheating field well above their clean limit. In this work extensive numerical simulation of the power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was performed. Our simulations take into account the nonlinear vortex line tension, vortex mass, Bardeen?Stephen viscous vortex drag applicable at low fields, and nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient ?(v) at high fields and pinning force. We showed that the LO decrease of ?(v) with the vortex velocity v could radically change the field dependence of the surface resistance Ri(H) caused by trapped vortices. At low frequencies Ri(H) exhibits a conventional increase with H. However, as frequency increases, the surface resistance becomes a nonmonotonic function of H which decreases with H at higher fields irrespective of the pinning distribution. Overheating can mask the descending field dependence of Ri(H) as frequency increases. Our numerical simulations also show that the LO effect can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along with the vortex when a vortex is pinned strongly to one end. Nonlinear losses of trapped vortices in thick films under high-amplitude RF fields as functions of frequency, mean free path and pinning characteristics have been calculated.

Pathirana, Walive↗
DOE OSTI2022

A Phase-Field Model for In-Space Manufacturing of Binary Alloys

The integrity of the final printed components is mostly dictated by the adhesion between the particles and phases that form upon solidification, which is a major problem in printing metallic parts using available In-Space Manufacturing (ISM) technologies based on the Fused Deposition Modeling (FDM) methodology. Understanding the melting/solidification process helps increase particle adherence and allows to produce components with greater mechanical integrity. We developed a phase-field model of solidification for binary alloys. The phase-field approach is unique in capturing the microstructure with computationally tractable costs. The developed phase-field model of solidification of binary alloys satisfies the stability conditions at all temperatures. The suggested model is tuned for Ni-Cu alloy feedstocks. We derived the Ginzburg-Landau equations governing the phase transformation kinetics and solved them analytically for the dilute solution. We calculated the concentration profile as a function of interface velocity for a one-dimensional steady-state diffuse interface neglecting elasticity and obtained the partition coefficient, k, as a function of interface velocity. Numerical simulations for the diluted solution are used to study the interface velocity as a function of undercooling for the classic sharp interface model, partitionless solidification, and thin interface.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Global weak solutions for a nonlocal multispecies Fokker–Planck–Landau system

The global-in-time existence of weak solutions to a spatially homogeneous multispecies Fokker–Planck–Landau system for plasmas in the three-dimensional whole space is shown. The Fokker–Planck–Landau system is a simplification of the Landau equations assuming a linearized, velocity-independent, and isotropic kernel. The resulting equations depend nonlocally and nonlinearly on the moments of the distribution functions via the multispecies local Maxwellians. Furthermore, the existence proof is based on a three-level approximation scheme, energy and entropy estimates, as well as compactness results, and it holds for both soft and hard potentials.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Landau-Lifshitz-Bloch equation for ferrimagnets with higher-order interaction

In this work, we present a micromagnetic formulation for modeling the magnetization dynamics and thermal equilibrium in ferrimagnetic materials at low and elevated temperatures. The formulation is based on a mean-field approximation (MFA). In this formulation, the ferrimagnet is described micromagnetically by two coupled sublattices with corresponding interactions, including inter- and intrasublattice micromagnetic exchange as well as four-spin interactions described as an intersublattice molecular field with a cubic dependence of the magnetization. The MFA is used to derive a Landau-Lifshitz-Bloch type equation for ferrimagnetic material, including cases with a ferromagnetic-antiferromagnetic phase transitions. For validation, the results obtained via the presented model are compared with recent experimental data for phase transitions in FeRh.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Coherent Spin Pumping in a Strongly Coupled Magnon-Magnon Hybrid System

Here, we experimentally identify coherent spin pumping in the magnon-magnon hybrid modes of yttrium iron garnet/permalloy (YIG/Py) bilayers. By reducing the YIG and Py thicknesses, the strong interfacial exchange coupling leads to large avoided crossings between the uniform mode of Py and the spin wave modes of YIG enabling accurate determination of modification of the linewidths due to the dampinglike torque. We identify additional linewidth suppression and enhancement for the in-phase and out-of-phase hybrid modes, respectively, which can be interpreted as concerted dampinglike torque from spin pumping. Furthermore, varying the Py thickness shows that both the fieldlike and dampinglike couplings vary like 1/$\sqrt{𝑡_{Py}}$, verifying the prediction by the coupled Landau-Lifshitz equations.

ferromagnetic resonance↗
DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Novel signatures of radiation reaction in electron–laser sidescattering

In this article we investigate novel signatures of radiation reaction via the angular deflection of an electron beam colliding at 90 degrees with an intense laser pulse. Due to the radiation reaction effect, the electrons can be deflected towards the beam axis for plane wave backgrounds, which is not possible in the absence of radiation reaction effects. The magnitude and size of the deflection angle can be controlled by tailoring the laser pulse shapes. The effect is first derived analytically using the Landau–Lifshitz equation, which allows to determine the important scaling behavior with laser intensity and particle energy. We then move on to full scale 3D Monte Carlo simulations to verify the effect is observable with present day laser technology. We investigate the opportunities for an indirect observation of laser depletion in such side scattering scenarios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Structure preservation using discrete gradients in the Vlasov-Poisson-Landau system

We present a novel structure-preserving framework for solving the Vlasov-Poisson-Landau system of equations using a particle in cell (PIC) discretization combined with discrete gradient time integrators. The Vlasov-Poisson-Landau system is an accurate model for studying hot plasma dynamics at a kinetic scale where small-angle Coulomb collisions dominate. Our scheme guarantees conservation of mass, momentum and energy as well as preservation of the monotonicity of entropy production in both the time-continuous and discrete systems. We employ the conservative integrator for both the Hamiltonian Vlasov-Poisson equations and the dissipative Landau equation using the PETSc library (www.mcs.anl.gov/petsc) to showcase structure-preserving properties.

Discrete gradients↗
DOE OSTI2020

Dynamic pair-breaking current, critical superfluid velocity, and nonlinear electromagnetic response of nonequilibrium superconductors

Here, we report numerical calculations of a dynamic pair-breaking current density $J_d$ and a critical superfluid velocity $v_d$ in a nonequilibrium superconductor carrying a uniform, large-amplitude AC current density $J(t) = J_a$sinΩ$\textit{t}$ with Ω well below the gap frequency Ω $\ll Δ_0/ℏ$. The dependencies $J_d$(Ω, $\textit{T}$) and $v_d$(Ω, $\textit{T}$) near the critical temperature $T_c$ were calculated from either the full time-dependent nonequilibrium equations for a dirty $\textit{s}$-wave superconductor or the time-dependent Ginzburg-Landau (TDGL) equations for a gapped superconductor, taking into account the GL relaxation time of the order parameter $τ_{\text{GL}}$ and the inelastic electron-phonon relaxation time of quasiparticles $τ_E$. We show that both approaches give similar frequency dependencies of $J_d$(Ω) and $v_d$(Ω) which gradually increase from their static pair-breaking GL values $J_c$ and $v_c$ at Ω$τ_E \ll$1 to $\sqrt{2}J_c$ and $\sqrt{2}v_c$ at Ω$τ_E \gg$ 1. Here $J_d, v_d$ and a dynamic superheating field at which the Meissner state becomes unstable were calculated in two different regimes of a fixed AC current and a fixed AC superfluid velocity induced by the applied AC magnetic field $H = H_a$ sinΩ$\textit{t}$ in a thin superconducting filament or a type-II superconductor with a large GL parameter. We also calculated a nonlinear electromagnetic response of a nonequilibrium superconducting state, particularly a dynamic kinetic inductance and a dissipative quasiparticle conductivity, taking into account the oscillatory dynamics of superconducting condensate and the kinetics of quasiparticles driven by a strong AC current. It is shown that an AC current density produces multiple harmonics of the electric field, the amplitudes of the higher-order harmonics diminishing as $τ_E$ increases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture $n$ = 0 axisymmetric mode evolution in realistic tokamak configurations. As we report here, these improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Universal time-dependent Ginzburg-Landau theory

We study the hydrodynamics of superconductors within the framework of Schwinger-Keldysh effective field theory (EFT). We show that in the vicinity of the superconducting phase transition the most general leading-order EFT satisfying the local Kubo-Martin-Schwinger condition is described by a version of the time-dependent Ginzburg-Landau (TDGL) equations augmented with stochastic terms. This version of TDGL is applicable in the gapless regime independent of any microscopic details. Within this approach, it is possible to include systematically the effects of nonuniform temperature and heat conductivity, as well as explicit or spontaneous breaking of time reversal. Here, we also introduce a thermal version of the Josephson relation and use it to construct an exotic hydrodynamics describing a phase of matter where heat can flow without dissipation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Pseudospectral Landau-Lifshitz description of magnetization dynamics

Magnetic materials host a wealth of nonlinear dynamics, textures, and topological defects. This is possible due to the competition between strong nonlinearity and dispersion, which act at the atomic scale, as well as long-range interactions. However, these features are difficult to study analytically and numerically because of the vastly different temporal and spatial scales involved. Here, in this work, we present a pseudospectral approach for the Landau-Lifshitz equation that invokes energy and momentum conservation embodied in the magnon dispersion relation to accurately describe both atomic and continuum limits. Furthermore, this approach enables analytical study at every scale. We show the applicability of this model in both the continuum and atomic limit by investigating modulational instability and ultrafast evolution of magnetization due to transient grating, respectively, in a one-dimensional ferromagnetic chain with perpendicular magnetic anisotropy. This model provides the possibility of grid-independent multiscale numerical approaches that will enable the description of singularities within a single framework.

1-dimensional spin chains↗