Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlocal models”

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

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators

Here, we conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C 2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C 2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C 2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Machine Learning Global Simulation of Nonlocal Gravity Wave Propagation

Global climate models typically operate at a grid resolution of hundreds of kilometers and fail to resolve atmospheric mesoscale processes, e.g., clouds, precipitation, and gravity waves (GWs).Model representation of these processes and their sources is essential to the global circulation and planetary energy budget, but subgrid scale contributions from these processes are often only approximately represented in models using parameterizations. These parameterizations are subject to approximations and idealizations, which limit their capability and accuracy. The most drastic of these approximations is the “single-column approximation” which completely neglects the horizontal evolution of these processes, resulting in key biases in current climate models. With a focus on atmospheric GWs, we present the first-ever global simulation of atmospheric GW fluxes using machine learning (ML) models trained on the WINDSET dataset to emulate global GW emulation in the atmosphere, as an alternative to traditional single-column parameterizations. Using an Attention U-Net-based architecture trained on globally resolved GW momentum fluxes, we illustrate the importance and effectiveness of global nonlocality, when simulating GWs using data-driven schemes.

Aman Gupta↗
DOE OSTI2026

Instabilities and phase transitions in architected metamaterials: a gradient-enhanced continuum approach

Architected metamaterials such as foams and lattices exhibit a wide range of properties governed by microstructural instabilities and emerging phase transitions. Their macroscopic response–including energy dissipation during impact, large recoverable deformations, morphing between configurations, and auxetic behavior–remains difficult to capture with conventional continuum models, which often rely on discrete approaches that limit scalability. In this work, we propose a nonlocal continuum formulation that captures both stable and unstable responses of elastic architected metamaterials. The framework extends anisotropic hyperelasticity by introducing nonlocal variables and internal length scales reflective of microstructural features. Local polyconvex free-energy models are systematically augmented with two families of non-(poly)convex energies, enabling both metastable and bistable responses. Implementation in a finite element framework enables solution using a hybrid monolithic–staggered strategy. Simulations capture densification fronts, forward and reverse transitions, hysteresis loops, imperfection sensitivity, and globally coordinated auxetic modes. Overall, this framework provides a robust foundation for accelerated modeling of instability-driven phenomena in architected metamaterials, while enabling extensions to anisotropic, dissipative, and active systems as well as integration with data-driven and machine learning approaches.

42 ENGINEERING↗
DOE OSTI2022

Development of an x-ray radiography platform to study laser-direct-drive energy coupling at the National Ignition Facility

A platform has been developed to study laser-direct-drive energy coupling at the National Ignition Facility (NIF) using a plastic sphere target irradiated in a polar-direct-drive geometry to launch a spherically converging shock wave. To diagnose this system evolution, eight NIF laser beams are directed onto a curved Cu foil to generate Heα line emission at a photon energy of 8.4 keV. These x rays are collected by a 100-ps gated x-ray imager in the opposing port to produce temporally gated radiographs. The platform is capable of acquiring images during and after the laser drive launches the shock wave. A backlighter profile is fit to the radiographs, and the resulting transmission images are Abel inverted to infer radial density profiles of the shock front and to track its temporal evolution. Here, the measurements provide experimental shock trajectories and radial density profiles that are compared to 2D radiation-hydrodynamic simulations using cross-beam energy transfer and nonlocal heat-transport models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Eddy diffusivity operator in homogeneous isotropic turbulence

Here, we use the recently developed macroscopic forcing method to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However, for small scales eddy diffusivity is found to vanish inversely proportional to the wave number. Behavior at all scales is reasonably captured by a nonlocal eddy diffusivity operator modeled as D/$\sqrt{I - l^2∇^2}$, where D is the eddy diffusivity in the Boussinesq limit, and l is a constant on the order of the large-eddy length. These results present a direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Determining Exact RANS Operators with the Macroscopic Forcing Method (Final Report)

This report contains a compilation of key results and findings from the ACT project “Determining Exact RANS Operators with the Macroscopic Forcing Method.” The Macroscopic Forcing Method (MFM), a numerical tool for determining closure operators, is used to measure eddy diffusivity moments in Rayleigh-Taylor (RT) instability. It is first applied to low-Atwood 2D RT instability; that work is then extended to 3D RT at different finite Atwood numbers. It is found that nonlocality is important for modeling the mean scalar transport closure operator in RT mixing. Additionally, work is done to improve the statistical convergence of MFM for chaotic problems like RT mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Emission lines from X-ray-heated accretion disks in low-mass X-ray binaries

We investigate the structure of accretion disks illuminated by X-rays from a central compact object in a binary system. X-rays can photoionize the upper atmosphere of the disk and form an accretion disk corona (ADC) where emission lines can form. We construct a model to calculate the vertical structure and the emission spectrum of the ADC with parameters appropriate to low-mass X-ray binaries. These models are made by nonlocal thermodynamic equilibrium calculations of ion and level populations and include a large number of atomic processes for 10 cosmically abundant elements. Transfer of radiation is treated by using the escape probability formalism. The vertical temperature profile of the ADC consists of a Compton-heated region and a mid-T zone where the temperature is approximately 10(exp 6) K. A thermal instability occurs close to the disk photosphere and causes the temperature of the ADC to drop abruptly from 10(exp 6) K to several times 10(exp 4) K. The emission spectrum in the optical, ultraviolet, extreme ultraviolet, and X-ray range is discussed and compared with the observations.

Ko, Yuan-Kuen↗
DOE OSTI2020

Geometrical quasi-ballistic effects on thermal transport in nanostructured devices

In this work, we employ thermoreflectance thermal imaging to directly measure the steady-state two-dimensional (2D) temperature field generated by nanostructured heat sources deposited on silicon substrate with different geometrical configurations and characteristic sizes down to 400nm. The analysis of the results using Fourier’s law not only breaks down as size scales down, but it also fails to capture the impact of the geometry of the heat source. The substrate effective Fourier thermal conductivities fitted to wire-shaped and circular-shaped structures with identical characteristic lengths are found to display up to 40% mismatch. Remarkably, a hydrodynamic heat transport model reproduces the observed temperature fields for all device sizes and shapes using just intrinsic Si parameters, i.e., a geometry and size-independent thermal conductivity and nonlocal length scale. The hydrodynamic model provides insight into the observed thermal response and of the contradictory Fourier predictions. We discuss the substantial Silicon hydrodynamic behavior at room temperature and contrast it to InGaAs, which shows less hydrodynamic effects due to dominant phonon-impurity scattering.

42 ENGINEERING↗
DOE OSTI2024

DiscoFluxM version 1.x

DiscoFluxM is software for modeling the deformation behavior of crystalline materials within the third-party finite-element based software, MOOSE. The software enables full-field crystal plasticity finite element simulations of deformation processes of metals and other crystalline materials. DiscoFluxM is an implementation of the DiscoFlux theoretical and algorithmic framework for coupling single-crystal plasticity with dislocation transport. The continuum mesoscale modeling framework is founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise spatially-local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations.

Luscher, Darby↗
DOE OSTI2023

Exact and Model Exchange-Correlation Potentials for Open-Shell Systems

The conventional approaches to the inverse density functional theory problem typically assume nondegeneracy of the Kohn–Sham (KS) eigenvalues, greatly hindering their use in open-shell systems. Here, we present a generalization of the inverse density functional theory problem that can seamlessly admit degenerate KS eigenvalues. Additionally, we allow for fractional occupancy of the Kohn–Sham orbitals to also handle noninteracting ensemble-v-representable densities, as opposed to just noninteracting pure-v-representable densities. We present the exact exchange-correlation (XC) potentials for six open-shell systems–four atoms (Li, C, N, and O) and two molecules (CN and CH 2 )–using accurate ground-state densities from configuration interaction calculations. We compare these exact XC potentials with model XC potentials obtained using nonlocal (B3LYP, SCAN0) and local/semilocal (SCAN, PBE, PW92) XC functionals. Although the relative errors in the densities obtained from these DFT functionals are of $O$(10 –3 to 10 –2 ), the relative errors in the model XC potentials remain substantially large–$O$(10 –1 to 10 0 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Comments on the present state of second-order closure models for incompressible flows

Second-order closure models account for history and nonlocal effects of the mean velocity gradients on the Reynolds stress tensor. Turbulent flows involving body forces or curvature, Reynolds stress relaxational effects, and counter-gradient transport are usually better described. The topics are presented in viewgraph form and include: (1) the Reynolds stress transport equation; (2) issues in second-order closure modeling; and (3) near wall models.

Speziale, Charles G.↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Learning stochastic dynamics and predicting emergent behavior using transformers

We show that a neural network originally designed for language processing can learn the dynamical rules of a stochastic system by observation of a single dynamical trajectory of the system, and can accurately predict its emergent behavior under conditions not observed during training. We consider a lattice model of active matter undergoing continuous-time Monte Carlo dynamics, simulated at a density at which its steady state comprises small, dispersed clusters. We train a neural network called a transformer on a single trajectory of the model. The transformer, which we show has the capacity to represent dynamical rules that are numerous and nonlocal, learns that the dynamics of this model consists of a small number of processes. Forward-propagated trajectories of the trained transformer, at densities not encountered during training, exhibit motility-induced phase separation and so predict the existence of a nonequilibrium phase transition. Transformers have the flexibility to learn dynamical rules from observation without explicit enumeration of rates or coarse-graining of configuration space, and so the procedure used here can be applied to a wide range of physical systems, including those with large and complex dynamical generators.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2023

Modeling the Chromosphere and Transition Region of Planet-hosting Star GJ 436

Ahead of upcoming space missions intending to conduct observations of low-mass stars in the ultraviolet (UV) spectral region it becomes imperative to simultaneously conduct atmospheric modeling from the UV to the visible (VIS) and near-infrared (NIR). Investigations on extended spectral regions will help to improve the overall understanding of the diversity of spectral lines arising from very different atmospheric temperature regions. Here we investigate atmosphere models with a chromosphere and transition region for the M2.5V star GJ 436, which hosts a close-in Hot Neptune. The atmosphere models are guided by observed spectral features from the UV to the VIS/NIR originating in the chromosphere and transition region of GJ 436. High-resolution observations from the Hubble Space Telescope and Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs (CARMENES) are used to obtain an appropriate model spectrum for the investigated M dwarf. We use a large set of atomic species considered in nonlocal thermodynamic equilibrium conditions within our PHOENIX model computations to approximate the physics within the low-density atmospheric regions. In order to obtain an overall match for the nonsimultaneous observations, it is necessary to apply a linear combination of two model spectra, where one of them better reproduces the UV lines while the other better represents the lines from the VIS/NIR range. This is needed to adequately handle different activity states across the observations.

stars↗
NASA NTRS1976

Solar astrophysics - Ghettosis from, or symbiosis with, stellar and galactic astrophysics

The purpose of the paper is to show how the solar-stellar symbiotic approach has led to the modeling of a star as a concentration of matter and energy. By 'solar-stellar symbiosis' is meant the philosophy of investigation according to which one asks what change in our general understanding of stellar structure and of stellar spectroscopic diagnostics is required to satisfy both the sun and an unusual star when, for example, some feature of an unusual star is discovered. The evolution of stellar models is traced, from walled, thermodynamic-equilibrium models to de-isolated models featuring transition zones and nonlocal thermodynamic equilibrium.

Pecker, J.-C.↗
NASA NTRS1981

The aluminum I autoionization doublet in the quiet solar spectrum

Observations are presented of the Al I autoionization doublet 1932 A and 1936 A in the quiet solar spectrum, obtained from the NRL slit spectrograph aboard Skylab and from the University of Hawaii Echelle Rocket Spectrograph. The observed profiles are compared with theoretical spectra computed for the Harvard Smithsonian Reference Atmosphere and the Vernazza, Avrett and Loeser (1976) solar models. It is found that nonlocal thermodynamic equilibrium effects are important in the line-formation problem and the synthetic spectra are in good agreeement with the data.

Heasley, J. N.↗