Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data Intensive Science”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Discrimination between individual dust and bioparticles using aerosol time-of-flight mass spectrometry

Ice nucleating particles (INPs) impact cloud properties and precipitation processes through their ability to trigger cloud glaciation. Dust and bioparticles are two important sources of INPs that have markedly different atmospheric loadings and ice nucleating efficiencies. In situ identification of the sources of INPs in clouds has been accomplished using single particle mass spectrometry (SPMS). However, external mixtures of dust and bioparticles present a unique challenge as they have overlapping mass spectral ion signatures, complicating their unambiguous identification. This study presents a detailed discussion of dust and bioparticle SPMS signatures, uniting data from a broad array of studies. As emphasized, the ion signals from both dust and bioparticles are highly sensitive to ionization conditions. To understand the observed variations, we characterize the mass spectral dependence of distinct dust and bioparticle samples using total positive ion intensity (TPII) as an indicator of the laser pulse energy each particle encountered. Through this analysis, a broad range of characteristic biogenic low intensity ion peaks that may be useful to distinguish bioparticles from dust became apparent and are highlighted. Insights informed by this analysis were utilized to identify bioparticles in ambient SPMS data. Ambient particles exhibiting both dust and characteristic biogenic spectral fingerprints were excluded from the bioparticle classification. Although bioparticles only made up 0.2% of all sampled particles, their abundance was moderately correlated with INP concentrations measured at 15 °C.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Atomic structures determined from digitally defined nanocrystalline regions

Nanocrystallography has transformed our ability to interrogate the atomic structures of proteins, peptides, organic molecules and materials. By probing atomic level details in ordered sub-10 nm regions of nanocrystals, scanning nanobeam electron diffraction extends the reach of nanocrystallography and in principle obviates the need for diffraction from large portions of one or more crystals. Scanning nanobeam electron diffraction is now applied to determine atomic structures from digitally defined regions of beam-sensitive peptide nanocrystals. Using a direct electron detector, thousands of sparse diffraction patterns over multiple orientations of a given crystal are recorded. Each pattern is assigned to a specific location on a single nanocrystal with axial, lateral and angular coordinates. This approach yields a collection of patterns that represent a tilt series across an angular wedge of reciprocal space: a scanning nanobeam diffraction tomogram. Using this diffraction tomogram, intensities can be digitally extracted from any desired region of a scan in real or diffraction space, exclusive of all other scanned points. Intensities from multiple regions of a crystal or from multiple crystals can be merged to increase data completeness and mitigate missing wedges. It is demonstrated that merged intensities from digitally defined regions of two crystals of a segment from the OsPYL/RCAR5 protein produce fragment-based ab initio solutions that can be refined to atomic resolution, analogous to structures determined by selected-area electron diffraction. In allowing atomic structures to now be determined from digitally outlined regions of a nanocrystal, scanning nanobeam diffraction tomography breaks new ground in nanocrystallography.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Structural Studies of Amorphous Materials by Fluctuation Electron Microscopy

Fluctuation Electron Microscopy (FEM) is a technique that examines the fluctuations in electron scattering across a uniformly thin amorphous sample. The statistics of the intensity fluctuations, mean and variance, reveal any underlying medium-range order present in the structure. The goals of this project were: (1) To determine the fundamentals of the scattering physics that gives rise to the variance signal in fluctuation electron microscopy (FEM); (2) To use these discoveries to find ways to quantify FEM; (3) To apply the FEM method to interesting and technologically important families of amorphous materials, particularly those with important applications in energy-related processes. Excellent progress was made in items (1) and (2). In stage (3) we did not examine the metamict zircons, as proposed. Instead, we examined films of polycrystalline and amorphous semi-conducting diamond. Significant accomplishments are: (1) A Reverse Monte Carlo procedure was successfully implemented to invert FEM data into a structural model. This is computer-intensive, but it demonstrated that diffraction and FEM data from amorphous silicon are most consistent with a paracrystallite model. This means that there is more diamond-like topology present in amorphous silicon than is predicted by the continuous random network model. (2) There is significant displacement decoherence arising in diffraction from amorphous silicon and carbon. The samples are being bombarded by the electron beam and atoms do not stay still while being irradiated – much more than was formerly understood. The atom motions cause the destructive and constructive interferences in the diffraction pattern to fluctuate with time, and it is the time-averaged speckle that is being measured. The variance is reduced by a factor m, 4 ≤ m ≤ 1000, relative to that predicted by kinematical scattering theory. (3) Speckle intensity obeys a gamma distribution, where the mean intensity $ \overline{I}\ $ and m are the two parameters governing the shape of the gamma distribution profile. m is determined by the illumination spatial coherence, which is normally very high, and mostly by the displacement decoherence within the sample. (4) Amorphous materials are more affected by the electron beam than are crystalline materials. Different samples exhibit different disruptibility, as measured by the effective values of m that fit the data. (5) Understanding the origin of the displacement decoherence better should lead to efficient methods for computing the observed variance from amorphous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The interaction of reactants, intermediates and products with Cu ions in Cu-SSZ-13 NH3 SCR catalysts: An energetic and ab initio X-ray absorption modeling study

Copper-exchanged SSZ-13 (Cu-SSZ-13) provides excellent catalytic activity and hydrothermal stability in the selective catalytic reduction (SCR) of NOx by using NH3 as a reductant. In this contribution, we compare the most likely positions for Cu with a single charge compensating Al atom (ZCu), including the effect of the adsorption of reactants, intermediates and products of a NH3 SCR reaction by using first-principles calculations based on density functional theory. We find that the 6-membered ring (6MR) site is the most energetically favorable while the 8-membered ring (8MR) sites are less favorable with energy differences about 0.5 eV with respect to the 6MR site. Upon molecular adsorption, the energy differences between Cu in the 8MR and 6MR sites decreases and almost disappears, which means that the thermodynamic stability of ZCu with molecularly adsorbed species is nearly identical for Cu in different sites. Considering more complicated scenarios of NO (or CO) adsorption, we find that the co-adsorption of 2 NO (or 2 CO) molecules weakens their interaction with Cu, as well as their co-adsorption with OH and H2O. The X-ray absorption near edge structure (XANES) of Cu in Cu-SSZ-13 under different conditions was also modeled from first principles. We find that the pre-edge feature in the K-edge XANES of Cu in clean Cu-SSZ-13 with Cu in a 8MR site is induced by a 1s to 4s transition while there is no pre-edge peak when Cu is in a 6MR site. Molecular adsorption onto Cu in a 6MR or a 8MR site results in the same K-edge XANES. When NO is adsorbed onto Cu in SSZ-13, we find that the origin of the small pre-edge peak in the K-edge XANES lies in a 1s to 3d electron transition. On the other hand, when CO or N2 is adsorbed onto an isolated Cu atom in SSZ-13, a transition from the 1s state to the hybridized state of 3d and 4p results in a pre-edge peak in the K-edge of the XANES spectrum. It is also found that the edge position and intensity of the maximum peak of Cu K-edge XANES depend on the Cu locations and the adsorbed molecules. By studying the Cu K-edge XANES of ZCu with adsorption of species of different oxidizing power, we conclude that a higher oxidation state of Cu will result in a higher energy edge position in the Cu K-edge XANES. These results are expected to greatly aid our understanding of Cu-SSZ-13 SCR catalysts. This work was supported by institutional funds provided to JSM from the Voiland School of Chemical Engineering and Bioengineering. Financial support was also provided by the National Science Foundation GOALI program under contract Nos. CBET-1258717, CBET-1258715 and CBET-1258690. We thank Prof. Fabio Ribeiro, Mr. Atish Parekh, Prof. W. F. Schneider, Mr. Christopher Paolucci, Mr. Trunjoyo Anggara and Prof. Jeff Miller for stimulating discussions on the modeling of the XANES spectrum and Ms. Alyssa Hensley for her comments on the manuscript. J.Sz. acknowledges the financial support of this work by the U.S. Department of Energy (DOE), Office of Energy Efficiency and Renewable Energy, Vehicle Technologies Program. A portion of the computer time for the computational work was performed using EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at PNNL. PNNL is a multi-program national laboratory operated for the US DOE by Battelle.

Zhang, Renqin↗
DOE OSTI2019

Elucidating the nature of crystallization kinetics in Zr 46 Cu 46 Al 8 metallic glass through simultaneous WAXS/SAXS measurements

Simultaneous wide- and small-angle X-ray scattering (WAXS/SAXS) measurements were performed to investigate the low-temperature crystallization kinetics of Zr 46 Cu 46 Al 8 bulk metallic glass during isothermal annealing. Quantitative analysis of the WAXS intensity data indicated that the crystallization process was dominated by steady-state homogeneous nucleation and three-dimensional diffusion-controlled growth. This observation was corroborated by quantitative analysis of the SAXS data showing a saturation of the growth of nanoscale crystallites. Comparison of simultaneously measured WAXS and SAXS data showed evidence of composition fluctuation prior to crystallization, suggesting that the crystallization was facilitated by nanoscale phase separation which established regions of favorable compositions catalyzing crystallization. The temperature dependence of crystallization mechanisms and the resulting microstructures in metallic glass are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Climatology of Clear-Air Turbulence and Mountain Wave Turbulence Throughout the Stratosphere Using MERRA-2 Data

Forty years of MERRA-2 reanalysis data is used to calculate twenty-nine unique indices for both clear-air turbulence (CAT) and mountain wave turbulence (MWT) throughout the stratosphere. The indices are combined into a single estimate of turbulence intensity for both CAT and MWT for five latitudinal domains. These indices have previously been applied at commercial flight levels, but this work is the first time such indices have been evaluated for estimating turbulence throughout the entire stratosphere. The results indicate that turbulence throughout the stratosphere is generally around four-times more intense in the winter season and increases in intensity with altitude. The data suggests the strongest source of CAT in each hemisphere is the polar vortex in the winter stratosphere above 50 hPa. The strongest sources of MWT are the Himalayas and the Andes in the lower stratosphere. The relationship between CAT and the quasi-biennial oscillation is also examined and indicates that easterly tropical winds are associated with the most intense turbulence. The trends in the climatology of CAT over forty years were compared indicating a significant increase in turbulence intensity in the upper stratosphere, and a weaker but still significant decrease in turbulence intensity in the lower stratosphere.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

Characteristics of Mini-Magnetospheres Formed by Paleo-Magnetic Fields of Mars

The intensely and non-uniformly magnetized crustal sources generate an effective large-scale magnetic field. In the Southern hemisphere the strongest crustal fields lead to the formation of large-scale mini-magnetospheres. In the Northern hemisphere, the crustal fields are rather weak and there are only isolated mini-magnetospheres. Re-connection with the interplanetary magnetic field (IMF) occurs in many localized regions. This may occur not only in cusp-like structures above nearly vertical field anomalies but also in halos extending several hundreds of kilometers from these sources. Re-connection will permit solar wind (SW) and more energetic particles to precipitate into and heat the neutral atmosphere. Electron density profiles of the ionosphere of Mars derived from radio occultation data obtained by the Radio Science Mars Global Surveyor (MGS) experiment are concentrated in the near polar regions. The effective scale-height of the neutral atmosphere density in the vicinity of the ionization peak has been derived for each of the profiles studied. The effective scale-heights have been compared with the crustal magnetic fields measured by the MGS Magnetometer/Electron Reflectometer (MAG/ER) experiment. A significant difference between the large-scale mini-magnetospheres and regions outside of them has been found. The neutral atmosphere is cooler inside the large-scale mini-magnetospheres. It appears that outside of the cusps the strong crustal magnetic fields prevent additional heating of the neutral atmosphere by direct interaction of the SW. The scale-height of the neutral atmosphere density derived from the experiment with the MGS Accelerometer has been compared with MAG/ER data. The scale-height was found to be usually larger than mean value near the boundaries of potential mini-magnetospheres and around cusps . It may indicate that the paleo-magnetic/IMF field re-connection is characteristic of the mini-magnetospheres at Mars.

N F Ness↗
DOE OSTI2010

An imaging proton spectrometer for short-pulse laser plasma experiments

The ultraintense short pulse laser pulses incident on solid targets can generate energetic protons. In addition to their potentially important applications such as in cancer treatments and proton fast ignition, these protons are essential to understand the complex physics of intense laser plasma interaction. To better characterize these laser-produced protons, we designed and constructed a novel spectrometer that will not only measure proton energy distribution with high resolution but also provide its angular characteristics. The information obtained from this spectrometer compliments those from commonly used diagnostics including radiochromic film packs, CR39 nuclear track detectors, and nonimaging magnetic spectrometers. The basic characterizations and sample data from this instrument are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Sensitivity of the KM3NeT detector to a neutrino flux from the Fermi bubbles

The Fermi Large Area Telescope data has provided evidence for a high-intensity emission of high-energy gamma rays with a E{sup −2} spectrum from two large bubbles above and below the Galactic Center. Hadronic mechanisms were proposed for this gamma-ray emission making the Fermi bubbles promising source candidates of high-energy neutrino emission. In this work preliminary simulation results regarding the detectability of high-energy neutrinos from the Fermi bubbles with the future multi-km{sup 3} neutrino telescope KM3NeT are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Calibration of scintillation-light filters for neutron time-of-flight spectrometers at the National Ignition Facility

Sixty-four neutral density filters constructed of metal plates with 88 apertures of varying diameter have been radiographed with a soft x-ray source and CCD camera at National Security Technologies, Livermore. An analysis of the radiographs fits the radial dependence of the apertures’ image intensities to sigmoid functions, which can describe the rapidly decreasing intensity towards the apertures’ edges. The fitted image intensities determine the relative attenuation value of each filter. Absolute attenuation values of several imaged filters, measured in situ during calibration experiments, normalize the relative quantities which are now used in analyses of neutron spectrometer data at the National Ignition Facility.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Benchmarking SARS CoV-2 Infection in the Workplace to Support Continuity of Operations

Objective: The COVID-19 pandemic jeopardizes continuity of operations of workplaces and the health and safety of workers. Exemplar workplace-related SARS-CoV-2 benchmarks are described and illustrated with empirical data. Methods: Benchmarks were collected over a 9-month period on a large workplace (N = 5500+). These ranged from quantitative indices associated with RT-qPCR targeted testing and random surveillance screening, surveillance for new variants of SARS-CoV-2, intensive contact tracing, case management, return to work procedures, to monitoring of antibody seropositive status. Results: Data and analyses substantiated effectiveness of interventions. This was evidenced in suppressed infection rates, rapid case identification and isolation, acceptance of the program by employees, documentation of presumptive immunity, and working relationships with senior management. Conclusions: These SARS-CoV-2 exemplar benchmarks provided an evidence-base for practice and contributed strategically to organizational decisions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Multicode benchmark on simulated Ti K-edge x-ray absorption spectra of Ti-O compounds

X-ray absorption spectroscopy (XAS) is an element-specific materials characterization technique that is sensitive to structural and electronic properties. First-principles simulated XAS has been widely used as a powerful tool to interpret experimental spectra and draw physical insights. Recently, there has also been growing interest in building computational XAS databases to enable data analytics and machine learning applications. However, there are non-trivial differences among commonly used XAS simulation codes, both in underlying theoretical formalism and in technical implementation. Reliable and reproducible computational XAS databases require systematic benchmark studies. In this work, we benchmarked Ti K-edge XAS simulations of ten representative Ti-O binary compounds, which we refer to as the Ti-O-10 dataset, using three state-of-the-art codes: XSPECTRA, OCEAN and exciting. We systematically studied the convergence behavior with respect to the input parameters and developed a workflow to automate and standardize the calculations to ensure converged spectra. Our benchmark comparison considers a 35 eV spectral range starting from the K-edge onset, representative of widely used near-edge spectra. Quantitative comparison over this range is based on Spearman’s rank correlation score (r sp ). Our results show: (1) the two Bethe-Salpeter equation (BSE) codes (OCEAN and exciting) have excellent agreement with an average $r$ $sp$ of 0.998; (2) good agreement is obtained between the core-hole potential code (XSPECTRA) and BSE codes (OCEAN and exciting) with an average $r$ $sp$ of 0.990, and this smaller $r$ $sp$ reflects the noticeable differences in the main edge spectral shape that can be primarily attributed to the difference in the strength of the screened core-hole potential; (3) simulations from both methods overall reproduce well the main experimental spectral features of rutile and anatase, and the different treatments of the screened core-hole potential have visible impact on pre-edge intensities and the peak ratio of the main edge; (4) there exist moderate differences in the relative edge alignment of the three codes with a standard deviation of about 0.2 eV, which arise from multiple contributions including the frozen core approximation, final state effects, and different approximations used for the self-energy correction. In conclusion, our benchmark study provides important standards for first-principles XAS simulations with broad impact in data-driven XAS analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evidence for phonon hardening in laser-excited gold using x-ray diffraction at a hard x-ray free electron laser

Studies of laser-heated materials on femtosecond timescales have shown that the interatomic potential can be perturbed at sufficiently high laser intensities. For gold, it has been postulated to undergo a strong stiffening leading to an increase of the phonon energies, known as phonon hardening. Despite efforts to investigate this behavior, only measurements at low absorbed energy density have been performed, for which the interpretation of the experimental data remains ambiguous. By using in situ single-shot x-ray diffraction at a hard x-ray free-electron laser, the evolution of diffraction line intensities of laser-excited Au to a higher energy density provides evidence for phonon hardening.

36 MATERIALS SCIENCE↗
DOE OSTI2017

High-Sensitivity Fast Neutron Detector KNK-2-8M

The design of the fast neutron detector KNK-2-8M is outlined. The results of he detector study in the pulse counting mode with pulses from {sup 238}U nuclei fission in the radiator of the neutron-sensitive section and in the current mode with separation of functional section currents are presented. The possibilities of determination of the effective number of {sup 238}U nuclei in the radiator of the neutron-sensitive section are considered. The diagnostic capabilities of the detector in the counting mode are demonstrated, as exemplified by the analysis of reference data on characteristics of neutron fields in the BR-1 reactor hall. The diagnostic capabilities of the detector in the current mode are demonstrated, as exemplified by the results of measurements of {sup 238}U fission intensity in the power startup of the BR-K1 reactor in the fission pulse generation mode with delayed neutrons and the detector placed in the reactor cavity in conditions of large-scale variation of the reactor radiation fields.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Reflectometry of X-ray Whispering Gallery Waves Propagating along Liquid Meniscuses

X-ray reflectometry of whispering galleries (WGs), which propagate along meniscuses of deionized water or silica hydrosols (Ludox{sup ®} SM) enriched by Cs{sup +} ions, was analyzed for the first time. The measurements were performed using the diffractometer with a moving tube-detector system. The X-ray beam rotation angle reached a maximum value of 4° on a silica hydrosol sample at CuK{sub α} radiation, while for H{sub 2}O it is only 1°. The intensities of WG modes propagating near the surface of a concave meniscus were found from the rigorous solution of the respective Helmholtz equation and applying the Monte Carlo technique for the random surface. For analysis of intensities of the X-ray scattering we used a hybrid model of the liquid surface with the upper transition layer having sinusoidal ripples in which the concentration of levitating Cs{sup +} ions has a maximum near the surface, derived from the experiment, while the concentration maximum for SiO{sub 2} particles with sizes of ~6–15 nm is positioned at the depth of ~15 nm from the hydrosol’s surface. Our numerical results for reflectances are in a good agreement with the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fast and accurate learned multiresolution dynamical downscaling for precipitation

Abstract. This study develops a neural-network-based approach for emulating high-resolution modeled precipitation data with comparable statistical properties but at greatly reduced computational cost. The key idea is to use combination of low- and high-resolution simulations (that differ not only in spatial resolution but also in geospatial patterns) to train a neural network to map from the former to the latter. Specifically, we define two types of CNNs, one that stacks variables directly and one that encodes each variable before stacking, and we train each CNN type both with a conventional loss function, such as mean square error (MSE), and with a conditional generative adversarial network (CGAN), for a total of four CNN variants. We compare the four new CNN-derived high-resolution precipitation results with precipitation generated from original high-resolution simulations, a bilinear interpolater and the state-of-the-art CNN-based super-resolution (SR) technique. Results show that the SR technique produces results similar to those of the bilinear interpolator with smoother spatial and temporal distributions and smaller data variabilities and extremes than the original high-resolution simulations. While the new CNNs trained by MSE generate better results over some regions than the interpolator and SR technique do, their predictions are still biased from the original high-resolution simulations. The CNNs trained by CGAN generate more realistic and physically reasonable results, better capturing not only data variability in time and space but also extremes such as intense and long-lasting storms. The new proposed CNN-based downscaling approach can downscale precipitation from 50 to 12 km in 14 min for 30 years once the network is trained (training takes 4 h using 1 GPU), while the conventional dynamical downscaling would take 1 month using 600 CPU cores to generate simulations at the resolution of 12 km over the contiguous United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Statistics of the charge spectrum of photo-multipliers and methods for absolute calibration

We derive the full expression for the shape of the charge spectrum that results from the illumination of a photo-multiplier tube. The derivation is for low intensity illumination with constant gain, a common condition for most nuclear and particle physics applications. Under these conditions, it is shown that an analytic expression can be formulated that describes details of the spectrum including the pedestal and dark noise with excellent fidelity to allow statistical fits to data. The derivation and full formula using either Gaussian or Poisson models for gain, and its limiting forms under various simplifying assumptions are presented with strategies on their use. The analytic description can be used to formulate data acquisition strategies to perform precise absolute calibration of photo-multipliers, the digitizers, and the data acquisition system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗