SEARCH · Engineering Papers

Results for “time discretization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Utah FORGE: 2024 Discrete Fracture Network Model Data

The Utah FORGE 2024 Discrete Fracture Network (DFN) Model dataset provides a set of files representing discrete fracture network modeling for the FORGE site near Milford, Utah. The dataset includes four distinct DFN model file sets, each corresponding to different time frames and modeling approaches in 2024. These models characterize both natural and induced fractures in the geothermal reservoir, which consists of crystalline granitic and metamorphic rock approximately 8,000 feet below the ground surface. The dataset includes a reference DFN model from February 2024 that incorporates planar fractures and well trajectories, as well as upscaled permeability, porosity, compressibility, and storage values on specified grids. Additionally, there are models based on new microseismic (MEQ) data from May and July 2024, including fracture planes fitted to the latest MEQ catalog datasets, tensile fractures from hydraulic stimulation, and an alternative connected DFN for modeling purposes. Coordinate data is provided in both global and local frames, with detailed instructions on the transformations used to align with principal stress orientations. The dataset also includes notes and calculation files for estimating fracture sizes and differences between various fracture sets. There are subfolders for Global Coordinates and Local Coordinates. To move from the global to the local coordinate frame, fractures and wells were a) rotated 20 degrees counterclockwise looking down about the global point (335376.400482041, 4263189.99998761, 250.093546450195) to better align with the principal stresses; and b) translated by (-335408.68, -4263010.9, 1150). Upscaled permeability values using the _XYZ suffix show directions with respect to the global XYZ coordinate frame, while those using the _IJK suffix are aligned with local coordinate frame.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Real-time well integrity monitoring in underground gas storage wells using distributed temperature and strain sensing: a field demonstration

Here, this article presents the first successful field demonstration of a combined distributed temperature and strain sensing (DTSS) system installed directly on newly replaced tubing in a 5400-ft-deep operational underground gas storage well. The DTSS system uses a single optical fiber to monitor temperature and strain in real-time, providing a cost-effective solution for long-term well integrity assessment. In this study, the strain–stress correlation of the tubing—representative of material behavior analysis—is investigated as a potential method for monitoring tubing integrity throughout its lifetime. Moreover, the DTSS system’s capability to support both continuous and discrete monitoring is evaluated by comparing future data with historical records, enabling the early detection of issues such as material fatigue, corrosion, or deformation. Overall, the work examines the effectiveness and scalability of the DTSS system for real-time monitoring of well operations and integrity in a newly replaced well.

Distributed Strain Sensing↗

DOE OSTI2022

Sub Coulomb barrier d+ 208 Pb scattering in the time-dependent basis function approach

We employ the non-perturbative time-dependent basis function (tBF) approach to study the scattering of the deuteron on 208 Pb below the Coulomb barrier. We obtain the bound and discretized scattering states of the projectile, which form the basis representation of the tBF approach, by diagonalizing a realistic Hamiltonian in a large harmonic oscillator basis. We find that the higher-order inelastic scattering effects are noticeable for sub barrier scatterings with the tBF method. We have successfully reproduced experimental sub Coulomb barrier elastic cross section ratios with the tBF approach by considering only the electric dipole (E1) component of the Coulomb interaction between the projectile and the target during scatterings. Here, we find that the correction of the polarization potential to the Rutherford trajectory is dominant in reproducing the data at very low bombarding energies, whereas the role of internal transitions of the deuteron projectile induced by the E1 interaction during the scattering becomes increasingly significant at higher bombarding energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Optimizing multigrid reduction-in-time and Parareal coarse-grid operators for linear advection

Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent partial differential equations (PDEs) in modern high-performance computing environments. While these techniques have been very successful for parabolic equations, it has often been observed that their performance suffers dramatically when applied to advection-dominated problems or purely hyperbolic PDEs using standard rediscretization approaches on coarse grids. In this paper, we apply MGRIT or Parareal to the constant-coefficient linear advection equation, appealing to existing convergence theory to provide insight into the typically nonscalable or even divergent behavior of these solvers for this problem. To overcome these failings, we replace rediscretization on coarse grids with improved coarse-grid operators that are computed by applying optimization techniques to approximately minimize error estimates from the convergence theory. Therefore, one of our main findings is that, in order to obtain fast convergence as for parabolic problems, coarse-grid operators should take into account the behavior of the hyperbolic problem by tracking the characteristic curves. Our approach is tested for schemes of various orders using explicit or implicit Runge–Kutta methods combined with upwind-finite-difference spatial discretizations. In all cases, we obtain scalable convergence in just a handful of iterations, with parallel tests also showing significant speed-ups over sequential time-stepping.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Effects of Dead‐End Fractures on Non‐Fickian Transport in Three‐Dimensional Discrete Fracture Networks

Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.

58 GEOSCIENCES↗

DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2026

Quantum Time-Space Tradeoffs for Matrix Problems

We consider the time and space required for quantum computers to solve a wide variety of problems involving matrices, many of which have only been analyzed classically in prior work. Our main results show that for a range of linear algebra problems—including matrix-vector product, matrix inversion, matrix multiplication and powering—existing classical time-space tradeoffs, several of which are tight for every space bound, also apply to quantum algorithms with at most a constant factor loss. For example, for almost all fixed matrices 𝐴, including the discrete Fourier transform matrix, we prove that quantum circuits with at most 𝑇 input queries and 𝑆 qubits of memory require 𝑇 = Ω⁢(𝑛 2 /𝑆) to compute matrix-vector product 𝐴⁢𝑥 for 𝑥 ∈{0,1 𝑛 . We similarly prove that matrix multiplication for 𝑛 ×𝑛 binary matrices requires 𝑇 = Ω⁢(𝑛 3 /$\sqrt{𝑆}$). Because many of our lower bounds are matched by deterministic algorithms with the same time and space complexity, our results show that quantum computers cannot provide any asymptotic advantage for these problems with any space bound. We obtain matching lower bounds for the stronger notion of quantum cumulative memory complexity—the sum of the space per layer of a circuit. We also consider Boolean (i.e., AND-OR) matrix multiplication and matrix-vector products, improving the previous quantum time-space tradeoff lower bounds for 𝑛 × 𝑛 Boolean matrix multiplication to 𝑇 = Ω⁢(𝑛 2.5 /𝑆 1/4 ) from 𝑇 = Ω⁢(𝑛 2.5 /𝑆 1/2 ). Our improved lower bound for Boolean matrix multiplication is based on a new coloring argument that extracts more from the strong direct product theorem that was the basis for prior work. To obtain our tight lower bounds for linear algebra problems, we require much stronger bounds than strong direct product theorems. We obtain these bounds by adding a new bucketing method to the quantum recording-query technique of Zhandry that lets us apply classical arguments to upper bound the success probability of quantum circuits.

lower bounds↗

DOE OSTI2022

Learning in continuous action space for developing high dimensional potential energy models

Reinforcement learning (RL) approaches that combine a tree search with deep learning have found remarkable success in searching exorbitantly large, albeit discrete action spaces, as in chess, Shogi and Go. Many real-world materials discovery and design applications, however, involve multi-dimensional search problems and learning domains that have continuous action spaces. Exploring high-dimensional potential energy models of materials is an example. Traditionally, these searches are time consuming (often several years for a single bulk system) and driven by human intuition and/or expertise and more recently by global/local optimization searches that have issues with convergence and/or do not scale well with the search dimensionality. Here, in a departure from discrete action and other gradient-based approaches, we introduce a RL strategy based on decision trees that incorporates modified rewards for improved exploration, efficient sampling during playouts and a “window scaling scheme" for enhanced exploitation, to enable efficient and scalable search for continuous action space problems. Using high-dimensional artificial landscapes and control RL problems, we successfully benchmark our approach against popular global optimization schemes and state of the art policy gradient methods, respectively. We demonstrate its efficacy to parameterize potential models (physics based and high-dimensional neural networks) for 54 different elemental systems across the periodic table as well as alloys. We analyze error trends across different elements in the latent space and trace their origin to elemental structural diversity and the smoothness of the element energy surface. Broadly, our RL strategy will be applicable to many other physical science problems involving search over continuous action spaces.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Transient disease dynamics across ecological scales

Analyses of transient dynamics are critical to understanding infectious disease transmission and persistence. Identifying and predicting transients across scales, from within-host to community-level patterns, plays an important role in combating ongoing epidemics and mitigating the risk of future outbreaks. Moreover, greater emphases on non-asymptotic processes will enable timely evaluations of wildlife and human diseases and lead to improved surveillance efforts, preventive responses, and intervention strategies. Here, we explore the contributions of transient analyses in recent models spanning the fields of epidemiology, movement ecology, and parasitology. In addition to their roles in predicting epidemic patterns and endemic outbreaks, we explore transients in the contexts of pathogen transmission, resistance, and avoidance at various scales of the ecological hierarchy. Examples illustrate how (i) transient movement dynamics at the individual host level can modify opportunities for transmission events over time; (ii) within-host energetic processes often lead to transient dynamics in immunity, pathogen load, and transmission potential; (iii) transient connectivity between discrete populations in response to environmental factors and outbreak dynamics can affect disease spread across spatial networks; and (iv) increasing species richness in a community can provide transient protection to individuals against infection. Ultimately, we suggest that transient analyses offer deeper insights and raise new, interdisciplinary questions for disease research, consequently broadening the applications of dynamical models for outbreak preparedness and management.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

On the initiation and evolution of dielectric breakdown in auto-magnetizing liner experiments

Auto-magnetizing (AutoMag) liners are cylindrical tubes composed of discrete metallic helices encapsulated in insulating material; when driven with a ~2 MA, ~100-ns prepulse on the 20 MA, 100-ns rise time Z accelerator, AutoMag targets produced >150 T internal axial magnetic fields. Once the current rise rate of the pulsed power driver reaches sufficient magnitude, the induced electric fields in the liner cause dielectric breakdown of the insulator material and, with sufficient current, the cylindrical target radially implodes. The dielectric breakdown process of the insulating material in AutoMag liners has been studied in experiments on the 500–900 kA, ~100-ns rise time Mykonos accelerator. Multi-frame gated imaging enabled the first time-resolved observations of photoemission from dynamically evolving plasma distributions during the breakdown process in AutoMag targets. Using magnetohydrodynamic simulations, we calculate the induced electric field distribution and provide a detailed comparison to the experimental data. We find that breakdown in AutoMag targets does not primarily depend on the induced electric field in the gaps between conductive helices as previously thought. Finally, to better control the dielectric breakdown time, a 12–32 mJ, 170 ps ultraviolet (λ = 266 nm) laser was implemented to irradiate the outer surface of AutoMag targets to promote breakdown in a controlled manner at a lower internal axial field. Here, the laser had an observable effect on the time of breakdown and subsequent plasma evolution, indicating that pulsed UV lasers can be used to control breakdown timing in AutoMag.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Critical and near-critical relaxation of holographic superfluids

We investigate the relaxation of holographic superfluids after quenches, when the end state is either tuned to be exactly at the critical point, or very close to it. By solving the bulk equations of motion numerically, we demonstrate that in the former case the system exhibits a power law falloff, as well as an emergent discrete scale invariance. The latter case is in the regime dominated by critical slowing down, and we show that there is an intermediate time range before the onset of late-time exponential falloff, where the system behaves similarly to the critical point with its power law falloff. We further postulate a phenomenological Gross-Pitaevskii-like equation (corresponding to model F of Hohenberg and Halperin) that is able to make quantitative predictions for the behavior of the holographic superfluid after near-critical quenches into the superfluid and normal phase. Intriguingly, all parameters of our phenomenological equation, which describes the nonlinear time evolution, may be fixed with information from the static equilibrium solutions and linear response theory. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2021

Pressure-Driven Fiber-Optic Sensor for Online Corrosion Monitoring

For many industrial applications, corrosion is a life-limiting phenomenon and therefore requires careful consideration and rigorous experimental testing before structural materials can be reliably deployed, particularly in harsh chemical environments. The traditional approach to measuring corrosion requires exposing many samples to the intended environment and then extracting them individually at discrete intervals for postexposure characterization. This approach does not provide a high degree of temporal resolution, nor does it provide any real-time information regarding dynamic changes in corrosion rates. Additionally, developing an online corrosion monitor capable of surviving harsh chemical environments could provide valuable information regarding the structural health of components and changing process conditions that could accelerate corrosion. To this end, a corrosion sensor was developed based on a pressure-driven Fabry-Pérot cavity (FPC). This sensor uses a pressure control system to internally pressurize the FPC formed between the sensor’s housing and a metal-embedded singlemode optical fiber. Simultaneous measurement of the change in FPC length using low-coherence interferometry and the applied pressure enables the calculation of the relative changes in the sensor’s diaphragm thickness due to corrosion on its outer surface. Measurements were made in situ while actively corroding the sensor and were validated against surveillance specimens that were corroded simultaneously. The uncertainties in the measured corrosion were analyzed using error propagation of the uncertainties in the measured pressures and displacements and were found to be <1% of the initial diaphragm thickness.

42 ENGINEERING↗

DOE OSTI2022

A massively parallel time-domain coupled electrodynamics–micromagnetics solver

We present a high-performance coupled electrodynamics–micromagnetics solver for full physical modeling of signals in microelectronic circuitry. The overall strategy couples a finite-difference time-domain approach for Maxwell’s equations to a magnetization model described by the Landau–Lifshitz–Gilbert equation. The algorithm is implemented in the Exascale Computing Project software framework, AMReX, which provides effective scalability on manycore and GPU-based supercomputing architectures. Furthermore, the code leverages ongoing developments of the Exascale Application Code, WarpX, which is primarily being developed for plasma wakefield accelerator modeling. Our temporal coupling scheme provides second-order accuracy in space and time by combining the integration steps for the magnetic field and magnetization into an iterative sub-step that includes a trapezoidal temporal discretization for the magnetization. The performance of the algorithm is demonstrated by the excellent scaling results on NERSC multicore and GPU systems, with a significant (59×) speedup on the GPU using a node-by-node comparison. We demonstrate the utility of our code by performing simulations of an electromagnetic waveguide and a magnetically tunable filter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Opportunities and Challenges to Capturing Distributed Battery Value via Retail Utility Rates and Programs

Distributed battery deployment is increasing with advanced metering, control, and communication technologies, leaving electric utilities with an under-utilized, flexible grid resource in aggregate. Rates can reflect locational and temporal prices while utility incentive-based programs allow DERs to provide direct grid services. However, utilities must balance accurately reflecting dynamic grid conditions versus simple and feasible design that encourages customer participation. Currently, most rates and incentive-based programs are simple, but as penetration of DER and advanced controls increase, dynamic designs could become prevalent. Utilities could encourage providing multiple services to optimize distributed battery dispatch and value streams, however, challenges persist when stacking services across distribution and bulk systems. A DER committed to multiple discrete services concurrently necessitates coordination between operators and a clear hierarchy of commitments. One way to address this is to separate commitments by time or capacity. For services that follow cyclic, predictable patterns, or those that are peak driven with predictability, an operator could ensure sufficient state of charge for participation, leaving time where a distributed battery could otherwise provide different services by segmenting participation temporally. To provide continuous or unexpected services, a battery operator may use state of charge management to reserve some percentage of the battery and segment participation by capacity. Macroeconomic trends, load patterns, generation profiles, and grid configurations drive variation in value and the subsequent implications for utility offerings and how a customer might participate. As distributed battery adoption increases, both regulators and utilities will need to ensure no adverse grid impacts and encourage provision of societal value beyond the customer domain.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗

DOE OSTI2022

3D, Real-time Label-free Visualization of Cellular Health

This work directly addresses a key question raised in TA7, Topic L1 of HDTRA1-16-24-FRCWMDCall, “Can quantum level changes be detected, mapped, and understood in the biological environment to allow faster diagnostic responses?” This work will use the inherently sensitive quantum property of a molecular excited state lifetime to visualize and map discrete, subtle nanometer-scale changes (e.g. NADPH vs NADH) that are key indicators of cellular health in response to chemical and biological threats and MCMs. It will also develop and exploit new, fast imaging light-sheet methods to allow faster diagnostic responses. This new imaging system will exploit time-correlated single photon counting and new advances in single-objective light-sheet microscopy (oblique plane microscopy) for fast, quantitative analysis of cellular redox state. We will add an entirely new dimension (excited state fluorescence lifetime) to oblique plane microscopy—a method that even without this added dimension was labeled by Nature Methods in 2021 as a ‘Method to Watch.’ Following instrument development and validation, we will explore cellular response to Burkholderia infections and antibiotic treatment in immortal cell lines followed by expansion to more realistic cellular environments (such as neuro-muscular junctions, NMJs) to test nerve agent simulants and their MCMs. As a microscopy platform, these methods hold the promise of being able to visualize a single infected cell (or a cellular compartment) and MCM treatment well before a larger and later organ/organism response.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Scalable freeform optimization of wide-aperture 3D metalenses by zoned discrete axisymmetry

We introduce a novel framework for design and optimization of 3D freeform metalenses that attains nearly linear scaling of computational cost with diameter, by breaking the lens into a sequence of radial “zones” with 𝑛-fold discrete axisymmetry, where 𝑛 increases with radius. This allows vastly more design freedom than imposing continuous axisymmetry, while avoiding the compromises of the locally periodic approximation (LPA) or scalar diffraction theory. Using a GPU-accelerated finite-difference time-domain (FDTD) solver in cylindrical coordinates, we perform full-wave simulation and topology optimization within each supra-wavelength zone. We validate our approach by designing millimeter and centimeter-scale, poly-achromatic, 3D freeform metalenses which outperform the state of the art. By demonstrating the scalability and resulting optical performance enabled by our “zoned discrete axisymmetry” (ZDA) and supra-wavelength domain decomposition, we highlight the potential of our framework to advance large-scale meta-optics and next-generation photonic technologies.

Sun, Mengdi [Wesleyan University]↗