Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlocality”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Shortwave feedbacks and El Nino-Southern Oscillation: Forced ocean coupled ocean-atmosphere experiments

Changes in tropical sea surface temperature (SST) can produce changes in cloudiness that modify incoming solar shortwave (SW) radiation, which in turn affects SST. The effects of this negative feedback on Pacific interannual variability are examined in forced ocean model and hybrid coupled ocean-atmosphere model simulations. Two empirical schemes are used to model the large-scale, low-frequency response of surface SW to SST anomalies. The first scheme attempts to account for the nonlocal nature of the atmospheric response to SST based patterns of covariability analyzed through singular value decomposition. In the observations the primary coupled mode of variability is composed of a SW anomaly in the central Pacific that covaries with anomalous SST in the eastern Pacific. This is applied in the model as a nonlocal feedback. The second scheme examines the effects of a purely local feedback with a spatially varying coefficient of magnitude chosen similar to the first scheme. In almost all cases the second scheme behaved similarly to the first, presumably because the correlation scale of SST is large enough for El Nino-Southern Oscillation (ENSO) dynamics that there is little sensitivity to the local approximation in the SW feedback. In simulations forced by time series of observed wind stress the SW feedback induced very minor SST damping. Results for a simpified heat budget analysis showed that while the SW feedback increased the local heat flux damping on SST, it also induced a mean shallowing of the mixed layer. The resulting changes in both the local mean vertical temperature gradient and the zonal velocity response to the wind stress acted to oppose the local heat flux damping effects. When the observed SW anomalies were applied to forced simulations, the simulated SST anomalies were modified as expected, and agreement with observed SST improved. In coupled simulations the SW feedbacks had greater impact than in the case of specified stress. The main effects were to decrease the magnitude of the warm and cold SST anomalies in the central Pacific, while leaving the pattern and evolution of ENSO anomalies essentially unchanged elsewhere. The SW feedbacks thus produce a modest improvement of the model ENSO SST pattern compared with observations, although they tended to shorten the period of the model ENSO cycle. Overall the results suggest that large-scale SW feedbacks are of quantitative importance to simulating some aspects of the ENSO cycle but are not critical to the overall occurrence of the phenomenon.

Waliser, Duane E.↗
NASA NTRS2013

Cloud Simulations in Response to Turbulence Parameterizations in the GISS Model E GCM

The response of cloud simulations to turbulence parameterizations is studied systematically using the GISS general circulation model (GCM) E2 employed in the Intergovernmental Panel on Climate Change's (IPCC) Fifth Assessment Report (AR5).Without the turbulence parameterization, the relative humidity (RH) and the low cloud cover peak unrealistically close to the surface; with the dry convection or with only the local turbulence parameterization, these two quantities improve their vertical structures, but the vertical transport of water vapor is still weak in the planetary boundary layers (PBLs); with both local and nonlocal turbulence parameterizations, the RH and low cloud cover have better vertical structures in all latitudes due to more significant vertical transport of water vapor in the PBL. The study also compares the cloud and radiation climatologies obtained from an experiment using a newer version of turbulence parameterization being developed at GISS with those obtained from the AR5 version. This newer scheme differs from the AR5 version in computing nonlocal transports, turbulent length scale, and PBL height and shows significant improvements in cloud and radiation simulations, especially over the subtropical eastern oceans and the southern oceans. The diagnosed PBL heights appear to correlate well with the low cloud distribution over oceans. This suggests that a cloud-producing scheme needs to be constructed in a framework that also takes the turbulence into consideration.

Yao, Mao-Sung↗
NASA NTRS2016

Impact of Source Region on the delta18 O signal in snow: A case study from Mount Wrangell Alaska

The stable isotopic composition of water in ice cores is an important source of information on past climate variability. At its simplest level, the underlying assumption is that there is an empirical relationship between the normalized difference in the concentration for these stable isotopes and a specified local temperature at the ice core site. There are, however, nonlocal processes, such as a change in source region or a change in the atmospheric pathway, which can impact the stable isotope signal, thereby complicating its use as a proxy for temperature. In this paper, the importance of these nonlocal processes are investigated through the analysis of the synoptic-scale circulation during a snowfall event at the summit of Mount Wrangell (62 deg N, 144 deg W; 4300 m MSL) in south-central Alaska. During this event there was, over a 1-day period in which the local temperature was approximately constant, a change in delta18 O that exceeded half that normally seen to occur in the region between summer and winter. As shall be shown, this arose from a change in the source region, from the subtropical eastern Pacific to northeastern Asia, for the snow that fell on Mount Wrangell during the event.

Snow↗
DOE OSTI2020

Dirichlet absorbing boundary conditions for classical and peridynamic diffusion-type models

Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of the PD diffusion models proposed so far are applied to bounded domains only. In this study, we propose an effective way to handle unbounded domains both with PD and classical diffusion models. For the former, we employ a meshfree discretization, whereas for the latter the finite element method (FEM) is employed. The proposed ABCs are time-dependent and Dirichlet-type, making the approach easy to implement in the available models. The performance of the approach, in terms of accuracy and stability, is illustrated by numerical examples in 1D, 2D, and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Optimal design of chemoepitaxial guideposts for the directed self-assembly of block copolymer systems using an inexact Newton algorithm

Directed self-assembly (DSA) of block copolymers (BCPs) is one of the most promising developments in the cost-effective production of nanoscale devices. The process makes use of the natural tendency for BCP melts to form nanoscale structures upon phase separation. The phase separation can be directed through the use of chemically patterned substrates to promote the formation of morphologies that are essential to the production of semiconductor devices. Moreover, the design of substrate pattern can be formulated as an optimization problem for which we seek optimal substrate designs that effectively produce given target morphologies. In this paper, we adopt a phase field model given by a nonlocal Cahn–Hilliard partial differential equation (PDE) based on the minimization of the Ohta–Kawasaki free energy, and present an efficient PDE-constrained optimization framework for the optimal design problem. The design variables are the locations of circular- or strip-shaped guiding posts that are used to model the substrate chemical pattern. To solve the ensuing optimization problem, we propose a variant of an inexact Newton conjugate gradient algorithm tailored to this problem. Additionally, we demonstrate the effectiveness of our computational strategy on numerical examples that span a range of target morphologies. Owing to our second-order optimizer and fast state solver, the numerical results demonstrate five orders of magnitude reduction in computational cost over previous work. The efficiency of our framework and the fast convergence of our optimization algorithm enable us to rapidly solve the optimal design problem in not only two, but also three spatial dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Many-Body Scars as a Group Invariant Sector of Hilbert Space

We present a class of Hamiltonians H for which a sector of the Hilbert space invariant under a Lie group G, which is not a symmetry of H, possesses the essential properties of many-body scar states. These include the absence of thermalization and the “revivals” of special initial states in time evolution. Here, a particular class of examples concerns interacting spin-1/2 fermions on a lattice consisting of N sites (it includes deformations of the Fermi-Hubbard model as special cases), and we show that it contains two families of N + 1 scar states. One of these families, which was found in recent literature, comprises the well-known η-pairing states. We find another family of scar states that is U(N) invariant. Both families and most of the group-invariant scar states produced by our construction, in general, give rise to the off-diagonal long range order, which survives at high temperatures and is insensitive to the details of the dynamics. Such states could be used for reliable quantum information processing because the information is stored nonlocally and thus cannot be easily erased by local perturbations. In contrast, other scar states we find are product states, which could be easily prepared experimentally. The dimension of scar subspace is directly controlled by the choice of group G and can be made exponentially large.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Generative Ensemble Regression: Learning Particle Dynamics from Observations of Ensembles with Physics-informed Deep Generative Models

Here, we propose a new method for inferring the governing stochastic ordinary differential equations (SODEs) by observing particle ensembles at discrete and sparse time instants, i.e., multiple “snapshots.” Particle coordinates at a single time instant, possibly noisy or truncated, are recorded in each snapshot but are unpaired across the snapshots. By training a physics-informed generative model that generates “fake” sample paths, we aim to fit the observed particle ensemble distributions with a curve in the probability measure space, which is induced from the inferred particle dynamics. We employ different metrics to quantify the differences between distributions, e.g., the sliced Wasserstein distances and the adversarial losses in generative adversarial networks. We refer to this method as generative “ensemble-regression” (GER), in analogy to the classic “point-regression,” where we infer the dynamics by performing regression in the Euclidean space. We illustrate the GER by learning the drift and diffusion terms of particle ensembles governed by SODEs with Brownian motions and Lévy processes up to 100 dimensions. We also discuss how to treat cases with noisy or truncated observations. Apart from systems consisting of independent particles, we also tackle nonlocal interacting particle systems with unknown interaction potential parameters by constructing a physics-informed loss function. Finally, we investigate scenarios of paired observations and discuss how to reduce the dimensionality in such cases by proving a convergence theorem that provides theoretical support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A data-driven peridynamic continuum model for upscaling molecular dynamics

Nonlocal models, including peridynamics, often use integral operators that embed lengthscales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials’ behavior, we propose a learning framework to extract, from MD data, an optimal Linear Peridynamic Solid (LPS) model as a surrogate for MD displacements. To maximize the accuracy of the learnt model we allow the peridynamic influence function to be partially negative, while preserving the well-posedness of the resulting model. To achieve this, we provide sufficient well-posedness conditions for discretized LPS models with sign-changing influence functions and develop a constrained optimization algorithm that minimizes the equation residual while enforcing such solvability conditions. This framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. We illustrate the efficacy of the proposed approach with several numerical tests for single layer graphene. Our two-dimensional tests show the robustness of the proposed algorithm on validation data sets that include thermal noise, different domain shapes and external loadings, and discretizations substantially different from the ones used for training.

homogenization↗
DOE OSTI2024

Unraveling Electronic and Vibrational Coherences Following a Charge Transfer Process in a Photosystem II Reaction Center

A reaction center is a unique biological system that performs the initial charge separation within a Photosystem II (PSII) multiunit enzyme, which eventually drives the catalytic water-splitting in plants and algae. The possible role of quantum coherences coinciding with the energy and charge transfer processes in PSII reaction center is one of the active areas of research. Here, we study these quantum coherences by using a numerically exact method on an excitonic dimer model, including linear vibronic coupling and employing optimal parameters from experimental two-dimensional coherent spectroscopic measurements. This enables us to precisely capture the excitonic interaction between pigments and the dissipation of the energy from electronic and charge-transfer (CT) states to the protein environment. We employ the time nonlocal (TNL) quantum master equation to calculate the population dynamics, which yields numerically reliable results. The calculated results show that, due to the strong dissipation, the lifetime of electronic coherence is too short to have direct participation in the charge transfer processes. However, there are long-lived vibrational coherences present in the system at frequencies close to the excitionic energy gap. These are strongly coupled with the electronic coherences, which makes the detection of the electronic coherences with conventional techniques very challenging. Additionally, we unravel the strong excitonic interaction of radical pair (PD1 and PD2) in the reaction center, which results in a long-lived electronic coherence of >100 fs, even at room temperature. Our work provide important physical insight to the charge separation process in PSII reaction center, which may be helpful for better understanding of photophysical processes in other natural and artificial light-harvesting systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

High‐Order Nonlinear Frequency Conversion in Transparent Conducting Oxide Thin Films

The study of conductive oxides has gained momentum within the photonics community due to their unique linear and nonlinear optical properties. Despite recent experiments reporting on high harmonic generation from thin films, the optical/electronic behavior of these compounds at the nanoscale is still not fully understood due to the lack of a suitable theoretical model. In the present work, aluminum zinc oxide is excited near its epsilon-near-zero crossing point using incident femtosecond pulses having peak power densities in the 1 TW cm −2 range. A relatively efficient frequency up-conversion including even and odd harmonics up to the seventh order is observed. A hydrodynamic-Maxwell theoretical approach is adopted, capable of simultaneously taking into account linear and nonlinear dispersions, nonlocal effects, surface, magnetic, and bulk nonlinearities in a spectral region that spans over two and a half octaves from the UV to the NIR region. The study enables a deeper understanding of the fundamental material parameters regulating optical nonlinearities, providing important insights to engineer this class of materials for applications in sensing, ultra-fast physics, and spectroscopy.

material modeling↗
DOE OSTI2020

Almost-Rigidity of Frameworks

Here, we extend the mathematical theory of rigidity of frameworks (graphs embedded in d-dimensional space) to consider nonlocal rigidity and flexibility properties. We provide conditions on a framework under which (I) as the framework flexes continuously it must remain inside a small ball, a property we call “almost-rigidity”; (II) any other framework with the same edge lengths must lie outside a much larger ball; (III) if the framework deforms by some given amount, its edge lengths change by a minimum amount; (IV) there is a nearby framework that is prestress stable, and thus rigid. The conditions can be tested efficiently using semidefinite programming. The test is a slight extension of the test for prestress stability of a framework, and gives analytic expressions for the radii of the balls and the edge length changes. Examples illustrate how the theory may be applied in practice, and we provide an algorithm to test for rigidity or almost-rigidity. We briefly discuss how the theory may be applied to tensegrities. © 2020 Wiley Periodicals LLC.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗
DOE OSTI2024

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗
DOE OSTI2022

On the differential representation and color-kinematics duality of AdS boundary correlators

The AdS boundary correlators and their dual correlation functions of boundary operators have been the main dynamic observables of the holographic duality relating a bulk AdS theory and a boundary conformal field theory. We show that tree-level AdS boundary correlators for generic states can be expressed as nonlocal differential operators of a certain structure acting on contact Witten diagrams. We further write the boundary correlators in a form that is very similar to flat space amplitudes, with Mandelstam variables replaced by certain combinations of single-state conformal generators, prove that all tree-level AdS boundary correlators have a differential representation, and detail the conversion of such differential expressions to position space. We illustrate the construction through the computation of the boundary correlators of scalars coupled to gluons and gravitons; when converted to position space, they reproduce known results. Color-kinematics duality and BCJ relations can be defined in analogy with their flat space counterparts, and are respected by the scalar correlators with a gluon exchange. We also discuss potential approaches to the double copy and find that its direct generalization may require nontrivial extensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Measuring quantum discord at the LHC

There has been an increasing interest in exploring quantities associated with quantum information at colliders. We perform a detailed analysis describing how to measure the quantum discord in the top anti-top quantum state at the Large Hadron Collider (LHC). While for pure states, quantum discord, entanglement, and Bell nonlocality all probe the same correlations, for mixed states they probe different aspects of quantum correlations. The quantum discord, in particular, is interesting because it aims to encapsulate all correlations between systems that cannot have a classical origin. We employ two complementary approaches for the study of the top anti-top system, namely the decay method and the kinematic method. We highlight subtleties associated with measuring discord for reconstructed quantum states at colliders. Usually quantum discord is difficult to compute due to an extremization that must be performed. We show, however, that for the $t\bar{t}$ system this extremization can be performed analytically and we provide closed-form formulas for the quantum discord. We demonstrate that with current LHC datasets, quantum discord can be observed at 3.6 – 5.7σ, depending on the signal region, with the decay method and can be measured at a precision of 0.1 – 0.2% with the kinematic method. At the high luminosity LHC, the observation of quantum discord is expected to be > 5σ using both the decay and kinematic methods and can be measured with a precision of 5% with the decay method and 0.05% with the kinematic method. Additionally, we identify the kinematic cuts at the LHC to isolate the $t\bar{t}$ state that is separable but has non-zero discord. By systematically investigating quantum discord for the first time through a detailed collider analysis, this work expands the toolkit for quantum information studies in particle physics and lays the groundwork for deeper insights into the quantum properties in high-energy collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Holographic scattering requires a connected entanglement wedge

In AdS/CFT, there can exist local 2-to-2 bulk scattering processes even when local scattering is not possible on the boundary; these have previously been studied in con- nection with boundary correlation functions. We show that boundary regions associated with these scattering configurations must have O(1/G N ) mutual information, and hence a connected entanglement wedge. One of us previously argued for this statement from the boundary theory using operational tools in quantum information theory. We improve that argument to make it robust to small errors and provide a proof in the bulk using focusing arguments in general relativity. We also provide a direct link to entanglement wedge reconstruction by showing that the bulk scattering region must lie inside the con- nected entanglement wedge. Our construction implies the existence of nonlocal quantum computation protocols that are exponentially more efficient than the optimal protocols currently known.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗