Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlocal interactions”

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 OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

nPINNs: Nonlocal physics-informed neural networks for a parametrized nonlocal universal Laplacian operator. Algorithms and applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius $\delta$ goes to zero, and to the fractional Laplacian as $\delta$ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We also provide theoretical convergence rates with respect to $\delta$ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, $\delta$ and $\alpha$, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order $\alpha(y)$, which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as $\delta$. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

nPINNs: nonlocal Physics-Informed Neural Networks for a parametrized nonlocal universal Laplacian operator. Algorithms and Applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius δ goes to zero, and to the fractional Laplacian as δ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We provide theoretical convergence rates with respect to δ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, δ and α, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order α(γ), which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as δ. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

COLLISIONLESS DAMPING AT ELECTRON SCALES IN SOLAR WIND TURBULENCE

The dissipation of turbulence in the weakly collisional solar wind plasma is governed by unknown kinetic mechanisms. Two candidates have been suggested to play an important role in the dissipation, collisionless damping via wave-particle interactions and dissipation in small-scale current sheets. High resolution spacecraft measurements of the turbulent magnetic energy spectrum provide important constraints on the dissipation mechanism. The limitations of popular fluid and hybrid numerical schemes for simulation of the dissipation of solar wind turbulence are discussed, and instead a three-dimensional kinetic approach is recommended. We present a three-dimensional nonlinear gyrokinetic simulation of solar wind turbulence at electron scales that quantitatively reproduces the exponential form of the turbulent magnetic energy spectrum measured in the solar wind. A weakened cascade model that accounts for nonlocal interactions and collisionless Landau damping also quantitatively agrees with the observed exponential form. These results establish that a turbulent cascade of kinetic Alfven waves that is terminated by collisionless Landau damping is sufficient to explain the observed magnetic energy spectrum in the dissipation range of solar wind turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Large‐Area, High‐Numerical‐Aperture, Freeform Metasurfaces

Abstract Nanophotonic devices are optical platforms capable of unprecedented wavefront control. To push the limits of experimental device performance, scalable design methodologies that combine the simplicity and fabricability of conventional design paradigms with the extended capabilities of freeform optimization are required. A novel gradient‐based design framework for large‐area freeform metasurfaces is introduced in which nonlocal interactions between simply shaped nanostructures, placed on an irregular lattice, are tailored to produce high‐order hybridized modes that support customizable large‐angle scattering profiles. Utilizing this approach, multifunctional super‐dispersive metalenses are designed and experimentally demonstrated. The approach to high‐numerical‐aperture radial metalenses capable of diffraction limited focusing and the generation of donut‐shaped point spread functions is also extended. It is anticipated that these concepts will have utility in super‐resolution microscopy, particle trapping, additive manufacturing, and metrology applications that require ultra‐high numerical apertures.

Optics↗
DOE OSTI2013

Global current profile effects on the evolution and saturation of magnetic islands

The nonlinear evolution of magnetic islands is investigated by means of numerical simulations. The growth and saturation of the island are shown to depend not only on asymptotic tearing mode stability parameter {Delta}{sup Prime} but also on the initial global current profile. Presence of the external current sheets leads to the formation of different island state for the same value of {Delta}{sup Prime }. It is found that the flow vorticity generated by the interacting current sheets is an important feature of nonlocal interactions and nonlocal effects in magnetic islands dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Local interactions in renormalization methods for Navier-Stokes turbulence

Two distinct renormalization-group (RG) approaches are applied to Navier-Stokes turbulence: epsilon-RG and recursive RG. Epsilon-RG takes into account only nonlocal interactions and utilizes an infinitesimal subgrid (unresolvable scale) shell limit. Recursive RG takes into account both nonlocal and local interactions and does not require an infinitesimal subgrid shell limit to be taken. The role of local interactions and the introduction of RG-induced nonlinearities are discussed and clarified.

Zhou, YE↗
DOE OSTI2019

Quantum Monte Carlo Methods in Nuclear Physics: Recent Advances

In recent years, the combination of precise quantum Monte Carlo (QMC) methods with realistic nuclear interactions and consistent electroweak currents, in particular those constructed within effective field theories (EFTs), has led to new insights in light and medium-mass nuclei, neutron matter, and electroweak reactions. For example, with the same chiral interactions, QMC calculations can reproduce binding energies and radii for light nuclei, n –α scattering phase shifts, and the neutron matter equation of state. This compelling new body of work has been made possible both by advances in QMC methods for nuclear physics, which push the bounds of applicability to heavier nuclei and to asymmetric nuclear matter, and by the development of local chiral EFT interactions up to next-to-next-to-leading order and minimally nonlocal interactions including Δ degrees of freedom. In this review, we discuss these recent developments and give an overview of the exciting results for nuclei, neutron matter and neutron stars, and electroweak reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Accelerating magnonic simulations with the pseudospectral Landau-Lifshitz equation

The pseudospectral Landau-Lifshitz (PS-LL) model can describe atomic-scale magnetic exchange interactions within a continuum framework. This is achieved by employing a convolution kernel that models the nonlocal interaction in a grid-independent manner. Even though the PS-LL was originally introduced to address atomic exchange, any nonlocal kernel can be modeled. In the field of magnonics, the dipole field is fundamental to describe the dispersion relation of magnons, the quasiparticle representation of angular momentum. Because dipole-dipole interactions are long-range, numerical approaches typically rely on convolutions. Here, we demonstrate that the PS-LL model can be used to perform magnonic simulations with a single convolution kernel derived from analytical solutions. We demonstrate a twofold increase in computational speed compared with the full dipole calculation. This approach is valid insofar as the excitations are linear, which is typically the case for magnons. Our results have the potential to accelerate magnonic research, particularly for the inverse design method, where several simulations must be performed to achieve the desired outcome.

Mathematics and computing↗
DOE OSTI2016

Interaction trajectory of solitons in nonlinear media with an arbitrary degree of nonlocality

The interaction trajectory of solitons in nonlocal nonlinear media is investigated. A simple differential equation describing the interaction trajectories is derived based on the light ray equation. Numerical calculations are carried out to illustrate the interaction trajectories with different parameters. The results show that the degree of nonlocality greatly affects the interaction of solitons. For a strongly nonlocal case, the interaction trajectory can be described by a cosine function. Analytical expressions describing the trajectory and the oscillation period are obtained. For generally and weakly nonlocal cases, the interaction trajectories still oscillate periodically, however it is no longer sinusoidal and the oscillation period increases with the nonlocal degree decreasing. In addition, the trajectory of two solitons launched with a relative angle at the entrance plane is investigated. It is found that there exists a critical angle. When the initial relative angle is larger than the critical angle, the two solitons do not collide on propagation. The influence of the degree of nonlocality on the critical angle is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Neutron matter from local chiral effective field theory interactions at large cutoffs

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts compared to nonlocal interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV ≤ Λ 𝑐 ≤ 700MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton 𝛽-decay half-life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Origin and effect of nonlocality in a layered composite.

A simple demonstration of nonlocality in a heterogeneous material is presented. By analysis of the microscale deformation of a two-component layered medium, it is shown that nonlocal interactions necessarily appear in a homogenized model of the system. Explicit expressions for the nonlocal forces are determined. The way these nonlocal forces appear in various nonlocal elasticity theories is derived. The length scales that emerge involve the constituent material properties as well as their geometrical dimen- sions. A peridynamic material model for the smoothed displacement eld is derived. It is demonstrated by comparison with experimental data that the incorporation of non- locality in modeling dramatically improves the prediction of the stress concentration in an open hole tension test on a composite plate.

Silling, Stewart Andrew↗
DOE OSTI2021

Ground-state properties of light 4 n self-conjugate nuclei in ab initio no-core Monte Carlo shell model calculations with nonlocal N N interactions

We report J π = 0 + ground-state energies and point-proton radii of 4 He, 8 Be, 12 C, 16 O, and 20 Ne nuclei calculated by the ab initio no-core Monte Carlo shell model with the JISP16 and Daejeon16 nonlocal NN interactions. Ground-state energies are obtained in the basis spaces up to seven oscillator shells (N shell =7) with several oscillator energies (ℏω) around the optimal oscillator energy for the convergence of ground-state energies. These energy eigenvalues are extrapolated to obtain estimates of converged ground-state energies in each basis space using energy variances of computed energy eigenvalues. We further extrapolate these energy-variance-extrapolated energies obtained in the finite basis spaces to infinite basis-space results with an empirical exponential form. This form features a dependence on the basis-space size but is independent of the value of ℏω used for the harmonic-oscillator basis functions. Point-proton radii for these states of atomic nuclei are also calculated following techniques employed for the energies. From these results, it is found that the Daejeon16 NN interaction provides good agreement with experimental data up to approximately 16 O, while the JISP16 NN interaction provides good agreement with experimental data up to approximately 12 C. Beyond these nuclei, the interactions produce overbinding accompanied by radii that are too small. These findings suggest and encourage further revisions of nonlocal NN interactions towards the investigation of nuclear structure in heavier-mass regions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Deterministic Fast Scrambling with Neutral Atom Arrays

Fast scramblers are dynamical quantum systems that produce many-body entanglement on a timescale that grows logarithmically with the system size N. Here, we propose and investigate a family of deterministic, fast scrambling quantum circuits realizable in near-term experiments with arrays of neutral atoms. We show that three experimental tools—nearest-neighbor Rydberg interactions, global single-qubit rotations, and shuffling operations facilitated by an auxiliary tweezer array—are sufficient to generate nonlocal interaction graphs capable of scrambling quantum information using only O(logN) parallel applications of nearest-neighbor gates. These tools enable direct experimental access to fast scrambling dynamics in a highly controlled and programmable way and can be harnessed to produce highly entangled states with varied applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Dynamics of metastable contact soliton dissipative exchange flows in one-dimensional ferromagnetic channels

Dissipative exchange flows (DEFs) are large-amplitude boundary value solutions of ferromagnetic channels. In their low-injection limit, DEFs reduce to spin superfluids. However, in the strong injection limit, nonlinearities dominate close to the injection site and a soliton is formed; this solution has been termed a contact soliton dissipative exchange flow (CS-DEF). Here, in this work, we numerically investigate CS-DEF solutions in a moderate injection regime and a finite injection width. We find a solution where two metastable solitons coexist in the injection region. This solution is metastable in the sense that any perturbation to the system will eject one of the solitons out of the injection region. Moreover, soliton dynamics can be excited when two injection regions are separated by a certain distance. We find that the ensuing DEF between the solitons induces a steady-state dynamics in which metastable solitons are continually ejected and nucleated. Furthermore, and depending on the relative signs of the spin injections, the soliton dynamics possess a particular handedness and frequency related to the spin transfer torque delivered by the DEF. Our results provide insights into the transport of spin current by DEFs - where the interaction between DEFs and solitons suggests a mechanism for detaching contact-solitons from the injection boundary. Although this study focuses on the "nonlocal" interaction between solitons, it may lead to the investigation of new mechanisms for inserting solitons in a DEF, e.g., for discrete motion and transport of information over long distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nematic superconductivity in twisted bilayer graphene

Twisted bilayer graphene displays insulating and superconducting phases caused by exceptional flattening of its lowest energy bands. Superconductivity with highest T c appears near half filling of the valence band (n ≈ –2). The data show that in the large part of the superconducting dome near n = –2, threefold lattice rotation symmetry is broken in the superconducting phase, i.e., a superconductor is also a nematic. We perform a comprehensive analysis of superconductivity in twisted-bilayer graphene within an itinerant approach and present a mechanism for nematic superconductivity. We take as an input the fact that at dopings, where superconductivity has been observed, the Fermi energy lies in the vicinity of twist-induced Van Hove singularities in the density of states. We argue that the low-energy physics can be properly described by patch models with either six or twelve Van Hove points. We obtain pairing interactions for the patch models in terms of parameters of the microscopic model for the flat bands, which contains both local and twist-induced nonlocal interactions. Here, we show that the latter give rise to attraction in different superconducting channels. For six Van Hove points, there is just one attractive d-wave channel, and we find chiral d ± id superconducting order, which breaks time-reversal symmetry but leaves the lattice rotation symmetry intact. For twelve Van Hove points, we find two attractive channels, g and i waves, with almost equal coupling constants. We show that both order parameters are nonzero in the ground state and explicitly demonstrate that in this co-existence state the threefold lattice rotation symmetry is broken, i.e., the superconducting state is also a nematic. We find two possible nematic states, one is time-reversal symmetric, the other additionally breaks time-reversal symmetry. We apply our scenario to twisted bilayer graphene near n = –2 and argue that it is applicable also to other systems with two (or more) attractive channels with similar couplings as our reasoning for a nematic superconductivity is based on generic symmetry considerations.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Aspects of renormalization in finite-density field theory

We study the renormalization of the Fermi surface coupled to a massless boson near three spatial dimensions. For this, we set up a Wilsonian RG with independent decimation procedures for bosons and fermions, where the four-fermion interaction “Landau parameters” run already at tree level. Our explicit one-loop analysis resolves previously found obstacles in the renormalization of finite-density field theory, including logarithmic divergences in nonlocal interactions and the appearance of multilogarithms. The key aspects of the RG are the above tree-level running, and a UV-IR mixing between virtual bosons and fermions at the quantum level, which is responsible for the renormalization of the Fermi velocity. We apply this approach to the renormalization of 2 k F singularities, and to Fermi surface instabilities in a companion paper, showing how multilogarithms are properly renormalized. We end with some comments on the renormalization of finite-density field theory with the inclusion of Landau damping of the boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗