Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

20% of US electricity from wind will have limited impacts on system efficiency and regional climate

Abstract Impacts from current and future wind turbine (WT) deployments necessary to achieve 20% electricity from wind are analyzed using high resolution numerical simulations over the eastern USA. Theoretical scenarios for future deployments are based on repowering (i.e. replacing with higher capacity WTs) thus avoiding competition for land. Simulations for the contemporary climate and current WT deployments exhibit good agreement with observed electricity generation efficiency (gross capacity factors (CF) from simulations = 45–48%, while net CF for WT installed in 2016 = 42.5%). Under the scenario of quadrupled installed capacity there is a small decrease in system-wide efficiency as indicated by annual mean CF. This difference is approximately equal to that from the two simulation years and may reflect saturation of the wind resource in some areas. WT modify the local near-surface climate in the grid cells where they are deployed. The simulated impact on near-surface climate properties at both the regional and local scales does not increase with increasing WT installed capacity. Climate impacts from WT are modest compared to regional changes induced by historical changes in land cover and to the global temperature perturbation induced by use of coal to generate an equivalent amount of electricity.

17 WIND ENERGY↗
DOE OSTI2022

Testing the conservative character of particle simulations: I. Canonical and noncanonical guiding center model in Boozer coordinates

The guiding center (GC) Lagrangian in Boozer coordinates for toroidally confined plasmas can be cast into canonical form by eliminating terms containing the covariant component B Ψ P of the magnetic field vector with respect to the poloidal flux function Ψ P . In an unperturbed plasma, B Ψ P can be eliminated via exact coordinate transformations, but, in general, one relies on approximations, assuming that the effect of B Ψ P is small. In this work, we are interested in the question whether Hamiltonian conservation laws are still satisfied when B Ψ P is retained in the presence of fluctuations. Considering fast ions in the presence of a shear Alfvén wave field with fixed amplitude, fixed frequency, and a single toroidal mode number n, we show that simulations using the code ORBIT with and without B Ψ P yield practically the same resonant and nonresonant GC orbits. The numerical results are consistent with theoretical analyses (presented in the appendix), which show that the unabridged GC Lagrangian with B Ψ P retained yields equations of motion that possess two key properties of Hamiltonian flows: (i) phase space conservation and (ii) energy conservation. As counter-examples, we also show cases where energy conservation (ii) or both conservation laws (i) and (ii) are broken by omitting certain small terms. When testing the conservative character of the simulation code, it is found to be beneficial to apply perturbations that do not resemble normal (eigen) modes of the plasma. The deviations are enhanced and, thus, more easily spotted when one inspects wave-particle interactions using nonnormal modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Evaluation of Electrical Resistivity Tomography to Monitor the Transport of Past Releases Beneath Tank Farms

Underground storage tanks at the Hanford Site, in southeastern Washington State, hold radioactive waste generated from four decades of plutonium production. The 149 single-shell tanks and the 28 double-shell tanks have all exceeded their initial design life of approximately 25 years. At least 67 tanks are assumed to have leaked in the past, resulting in radioactive releases into the vadose zone. Gamma ray logging within dry monitoring wells is currently the primary method for tracking the migration of leaked tank waste through the vadose zone. While this approach provides an accurate assessment of radioactive contamination, that information is only provided near (within ~1m) the borehole, leaving most of the vadose zone unmonitored, particularly the important region directly beneath the tank. This report describes a numerical study that investigates the feasibility and performance of time-lapse 3D electrical resistivity tomography (ERT) for long-term monitoring of a hypothetical tank waste location and migration through the vadose zone. ERT is a method of remotely imaging the bulk electrical conductivity (EC) of the subsurface, which is significantly impacted by the presence of conductive solid and liquid tank waste. The release of liquid tank wastes increases subsurface fluid conductivity and saturation over time, creating a target to use time-lapse ERT for long-term monitoring. Although the presence of metallic infrastructure can cause ERT interference, recent advancements in ERT data processing enable the deleterious effects of buried metallic infrastructure (e.g. pipes, wellbore casings, tanks) to be removed to better determine the liquid tank waste migration over time. Three hypothetical realistic scenarios were simulated in the ERT evaluation. The first two scenarios assume the same leak amount and rate (i.e., between 1/1/1951 and 12/31/1951 at the rate of 347 m 3 per year) but different leaky tanks. Scenario 1 assumes leaks under tank B-102, which is located on the edge of the B-tank farm and surrounded by a few metallic infrastructure including cased pipes/wells/tanks. Scenario 2 assumes leaks under tank B-108, which is located near the center of the B-tank farm and surrounded by larger amount of metallic infrastructure than B-102. Scenario 3 assumes the same metallic infrastructure as B-102, with a more recent contaminant leak that was simulated to have occurred between 1/1/2018 and 12/31/2023 at a rate of 1.89 m 3 per year. The leak time in Scenarios 1 and 2 corresponds to a historical overfill event in 1951 and Scenario 3 corresponds to a recent found tank leak in 2019. In each scenario, a “true” bulk EC model vs. time reflecting contaminant migration was generated. ERT data was simulated from these “true” bulk EC models and a time-lapse ERT inversion produced “imaged” bulk EC vs. time. Three electrode configurations in two, four and eight boreholes surrounding the leak tank were used in the ERT simulations in each scenario. These borehole configurations were considered logistically feasible and cost-effective for monitoring. The hypothetical ERT boreholes are assumed to have non-metallic casing. By comparing the “imaged” bulk EC with the “true” bulk EC, it was demonstrated that the three configurations of wells used (two, four, and eight wells) were able to successfully monitor the migration of tank leaks through the vadose zone, with bulk EC resolution increasing with the number of down borehole ERT arrays for all scenarios. Therefore, the use of eight boreholes to perform ERT monitoring beneath the tanks provided the best spatiotemporal information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition

This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Higher-order topological superconductivity in twisted bilayer graphene

In this work, we show that introducing spin-singlet or spin-triplet superconductivity into twisted bilayer graphene induces higher-order topological superconductivity. Multiple copies of C 2z ⁢T-protected Majorana Kramers pairs appear at corners on pairing domain walls. The topology originates from the anomaly analyzed in Song et al.—the absence of a lattice support—of the single-valley band structure of twisted bilayer graphene, which is protected by C 2z ⁢T and approximate particle-hole symmetry $\mathcal{P}$. We prove that any pairing (spin-singlet or spin-triplet) term preserving valley-U(1), spin-SU(2), time-reversal, C 2z ⁢T, and $\mathcal{P}$ must drive the system into a higher-order topological superconductor. Here spin-SU(2) is the global spin-SU(2) for the singlet pairing and is broken to U(1) for the triplet pairing. Using a Dirac Hamiltonian, we derive the corner modes and confirm with numerics. These corner states are stable even if $\mathcal{P}$ is weakly broken, which is true in experimental setups. Finally, we suggest experimental detection via the fractional Josephson effect in a TBG-TSC Josephson junction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Polaritonic Chern Insulators in Monolayer Semiconductors

Systems with strong light-matter interaction open up new avenues for studying topological phases of matter. Examples include exciton polaritons, mixed light-matter quasiparticles, where the topology of the polaritonic band structure arises from the collective coupling between matter wave and optical fields strongly confined in periodic dielectric structures. Distinct from light-matter interaction in a uniform environment, the spatially varying nature of the optical fields leads to a fundamental modification of the well-known optical selection rules, which were derived under the plane wave approximation. Here, in this study, we identify polaritonic Chern insulators by coupling valley excitons in transition metal dichalcogenides to photonic Bloch modes in a dielectric photonic crystal slab. We show that polaritonic Dirac points, which are markers for topological phase transition points, can be constructed from the collective coupling between valley excitons and photonic Dirac cones in the presence of both time-reversal and inversion symmetries. Lifting exciton valley degeneracy by breaking time-reversal symmetry leads to gapped polaritonic bands with nonzero Chern numbers. Through numerical simulations, we predict polaritonic chiral edge states residing inside the topological gaps. Our Letter paves the way for the further study of strong exciton-photon interaction in nanophotonic structures and for exploring polaritonic topological phases and their practical applications in polaritonic devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Determination of the nonlinear thermo-optic coefficient of silicon nitride and oxide using an effective index method

There is little literature characterizing the temperature-dependent thermo-optic coefficient (TOC) for low pressure chemical vapor deposition (LPCVD) silicon nitride or plasma enhanced chemical vapor deposition (PECVD) silicon dioxide at temperatures above 300 K. In this study, we characterize these material TOC’s from approximately 300-460 K, yielding values of (2.51 ± 0.08) · 10 −5 K −1 for silicon nitride and (5.67 ± 0.53) · 10 −6 K −1 for silicon oxide at room temperature (300 K). We use a simplified experimental setup and apply an analytical technique to account for thermal expansion during the extraction process. We also show that the waveguide geometry and method used to determine the resonant wavelength have a substantial impact on the precision of our results, a fact which can be used to improve the precision of numerous ring resonator index sensing experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Filtered Rayleigh-Ritz is all you need

Recent work has shown that the (block) Lanczos algorithm can be used to extract approximate energy spectra and matrix elements from (matrices of) correlation functions in quantum field theory, and identified exact coincidences between Lanczos analysis methods and others. In this work, we note another coincidence: the Lanczos algorithm is equivalent to the well-known Rayleigh-Ritz method applied to Krylov subspaces. Rayleigh-Ritz provides optimal eigenvalue approximations within subspaces; we find that spurious-state filtering allows these optimality guarantees to be retained in the presence of statistical noise. We explore the relation between Lanczos and Prony's method, their block generalizations, generalized pencil of functions (GPOF), and methods based on the generalized eigenvalue problem (GEVP), and find they all fall into a larger "Prony-Ritz equivalence class", identified as all methods which solve a finite-dimensional spectrum exactly given sufficient correlation function (matrix) data. This equivalence allows simpler and more numerically stable implementations of (block) Lanczos analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Temporal evolution of the light emitted by a thin, laser-ionized plasma source

We present an experimental and simulation-based investigation of the temporal evolution of light emission from a thin, laser-ionized helium plasma source. We demonstrate an analytic model to calculate the approximate scaling of the time-integrated, on-axis light emission with the initial plasma density and temperature, supported by the experiment, which enhances the understanding of plasma light measurement for plasma wakefield accelerator (PWFA) plasma sources. Our model simulates the plasma density and temperature using a split-step Fourier code and a particle-in-cell code. A fluid simulation is then used to model the plasma and neutral density, and the electron temperature as a function of time and position. We then show the numerical results of the space-and-time-resolved light emission and that collisional excitation is the dominant source of light emission. We validate our model by measuring the light emitted by a laser-ionized plasma using a novel statistical method capable of resolving the nanosecond-scale temporal dynamics of the plasma light using a cost-effective camera with microsecond-scale timing jitter. This method is ideal for deployment in the high radiation environment of a particle accelerator that precludes the use of expensive nanosecond-gated cameras. Our results show that our models can effectively simulate the dynamics of a thin, laser-ionized plasma source. In addition, this work provides a detailed understanding of the plasma light measurement, which is one of the few diagnostic signals available for the direct measurement of PWFA plasma sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Toward Verification of RANS Simulations of the T-Tube Modular Divertor Using Large Eddy Simulations of Impinging Turbulent Plane Jets

Turbulent impinging jets have been proposed to cool high heat flux plasma-facing components such as the solid tungsten target plates of the divertor in long-pulse magnetic fusion energy reactors. In particular, the T-tube modular divertor, originally developed by the ARIES Team, consists of two concentric cylindrical tubes where helium flows through a slot in the inner tube, forming an approximately planar jet that impinges upon and cools the inner surface of the pressure boundary (namely, the outer tube) and the ~15-cm 2 plasma-facing W target. The objective of this work is to demonstrate that large eddy simulations (LESs) accurately simulate the thermal transport in canonical flows that comprise the cooling flow in the T-tube, as well as validate temperatures from LES with experimental measurements in a simplified T-tube geometry. Wall‑resolved LESs, validated by experimental data and verified by direct numerical simulations (DNSs), provide benchmark data for two canonical flows in the T‑tube, namely, planar impinging and wall jets, for Reynolds numbers Re B = 4 × 10 3 to 2 × 10 4 . Our LES results are within 4% to 12% root-mean-square error (RMSE) of surface Nusselt number distributions (Nu) from experiments and DNSs. The validated LES results are then used as the ground truth to evaluate four Reynolds‑averaged Navier-Stokes (RANS) turbulence closures, namely, the k‑ω SST, realizable k‑ε, GEKO, and γ‑SST models. The k‑ω SST model has the best overall performance in terms of heat transfer, giving surface Nu within 12% RMSE of the LES results for high‑ReB impinging jets and reduced overprediction in the wall‑jet region. The GEKO model with default constants has the next best performance, providing slightly better Nu predictions for low ReB impinging jets (versus k-ω SST) but worse overall performance over the full range of ReB studied here. The realizable k‑ε turbulence model significantly overestimates turbulence near the stagnation point, while the γ‑SST model suppresses near‑wall production, biasing the simulations toward simulating laminar surface heat transfer. Simulations of the simplified T‑tube show that LES and RANS simulations with the k‑ω SST model give nearly identical average heat transfer coefficients (HTCs) over the impingement surface. The realizable k‑ε model predicts significantly lower wall temperatures due to overestimation of HTC in the outlet flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Single-ancilla ground state preparation via Lindbladians

We design a quantum algorithm for ground state preparation in the early fault tolerant regime. As a Monte Carlo style quantum algorithm, our method features a Lindbladian where the target state is stationary. The construction of this Lindbladian is algorithmic and should not be seen as a specific approximation to some weakly coupled system-bath dynamics in nature. Our algorithm can be implemented using just one ancilla qubit and efficiently simulated on a quantum computer. It can prepare the ground state even when the initial state has zero overlap with the ground state, bypassing the most significant limitation of methods like quantum phase estimation. As a variant, we also propose a discrete-time algorithm, demonstrating even better efficiency and providing a near-optimal simulation cost depending on the desired evolution time and precision. Numerical simulations using Ising and Hubbard models demonstrate the efficacy and applicability of our method. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗
DOE OSTI2025

Performance of wave function and Green's function methods for non-equilibrium many-body dynamics

Theoretical descriptions of the non-equilibrium dynamics of quantum many-body systems essentially employ either (i) explicit treatments, relying on the truncation of the expansion of the many-body wave function, (ii) compressed representations of the many-body wave function, or (iii) evolution of an effective (downfolded) representation through Green's functions. In this work, we select representative cases of each of the methods and address how these complementary approaches capture the dynamics driven by intense field perturbations to non-equilibrium states. Under strong driving, the systems are characterized by strong entanglement of the single-particle density matrix and natural populations approaching those of a strongly interacting equilibrium system. We generate a representative set of results that are numerically exact and form a basis for a critical comparison of the distinct families of methods. We demonstrate that the compressed formulation based on similarity-transformed Hamiltonians (coupled-cluster approach) is practically exact in weak fields and, hence, weakly or moderately correlated systems. Coupled cluster, however, struggles for strong driving fields, under which the system exhibits strongly correlated behavior, as measured by the von Neumann entropy of the single-particle density matrix. The dynamics predicted by Green's functions in the (widely popular) G W approximation are less accurate, but improve significantly upon the mean-field results in the strongly driven regime. Published by the American Physical Society 2025

Reeves, Cian C. (ORCID:0009000642581845)↗
DOE OSTI2023

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
DOE OSTI2021

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗
DOE OSTI2023

Universal linear intensity transformations using spatially incoherent diffractive processors

Abstract Under spatially coherent light, a diffractive optical network composed of structured surfaces can be designed to perform any arbitrary complex-valued linear transformation between its input and output fields-of-view (FOVs) if the total number ( N ) of optimizable phase-only diffractive features is ≥~2 N i N o , where N i and N o refer to the number of useful pixels at the input and the output FOVs, respectively. Here we report the design of a spatially incoherent diffractive optical processor that can approximate any arbitrary linear transformation in time-averaged intensity between its input and output FOVs. Under spatially incoherent monochromatic light, the spatially varying intensity point spread function ( H ) of a diffractive network, corresponding to a given, arbitrarily-selected linear intensity transformation, can be written as H ( m , n ; m ′, n ′) = | h ( m , n ; m ′, n ′)| 2 , where h is the spatially coherent point spread function of the same diffractive network, and ( m , n ) and ( m ′, n ′) define the coordinates of the output and input FOVs, respectively. Using numerical simulations and deep learning, supervised through examples of input-output profiles, we demonstrate that a spatially incoherent diffractive network can be trained to all-optically perform any arbitrary linear intensity transformation between its input and output if N ≥ ~2 N i N o . We also report the design of spatially incoherent diffractive networks for linear processing of intensity information at multiple illumination wavelengths, operating simultaneously. Finally, we numerically demonstrate a diffractive network design that performs all-optical classification of handwritten digits under spatially incoherent illumination, achieving a test accuracy of >95%. Spatially incoherent diffractive networks will be broadly useful for designing all-optical visual processors that can work under natural light.

36 MATERIALS SCIENCE↗