Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient structure”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

The onset of the new solar modulation cycle in 1987-1988 as a function of heliocentric radius and latitude

The onset of a new solar modulation cycle were studied in 1987 and 1988, on the basis of data on counting rates of particles with E greater than 60 MeV, observed at the IMP, Voyager (V), and Pioneer (P) satellites. It was found that the decrease at earth was rapid after the intensity maximum in early 1987 and was closely correlated with the increase in the average tilt of the heliospheric current sheet. The initial rapid intensity decrease was found to be related to three small Forbush decreases which are superimposed on a more gradual decrease. By the end of 1988, the intensity decreased by about 40 percent at earth, 30 percent at V2, and 18 percent at P10. This overall decrease was accompanied by an increase in the average integral radial gradient as well as a decrease in the radial dependence of the gradient. The onset of the new solar modulation cycle seems to be related to the complete altering of the solar magnetic structure as observed on the surface of the sun.

Webber, W. R.↗
NASA NTRS1991

Modelling the pressure-strain correlation of turbulence - An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗
NASA NTRS1996

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗
NASA NTRS1998

Effect of track structure and radioprotectors on the induction of oncogenic transformation in murine fibroblasts by heavy ions

The oncogenic potential of high-energy 56Fe particles (1 GeV/nucleon) accelerated with the Alternating Gradient Synchrotron at the Brookhaven National Laboratory was examined utilizing the mouse C3H 10T1/2 cell model. The dose-averaged LET for high-energy 56Fe is estimated to be 143 keV/micrometer with the exposure conditions used in this study. For 56Fe ions, the maximum relative biological effectiveness (RBEmax) values for cell survival and oncogenic transformation were 7.71 and 16.5 respectively. Compared to 150 keV/micrometer 4He nuclei, high-energy 56Fe nuclei were significantly less effective in cell killing and oncogenic induction. The prostaglandin E1 analog misoprostol, an effective oncoprotector of C3H 10T1/2 cells exposed to X rays, was evaluated for its potential as a radioprotector of oncogenic transformation with high-energy 56Fe. Exposure of cells to misoprostol did not alter 56Fe cytotoxicity or the rate of 56Fe-induced oncogenic transformation.

NASA Discipline Radiation Health↗
NASA NTRS2006

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗
NASA NTRS2015

Analysis of Advanced Thermoelectric Materials and Their Functional Limits

The world's demand for energy is increasing dramatically, but the best energy conversion systems operate at approximately 30% efficiency. One way to decrease energy loss is in the recovery of waste heat using thermoelectric (TE) generators. A TE generator is device that generates electricity by exploiting heat flow across a thermal gradient. The efficiency of a TE material for power generation and cooling is determined by the dimensionless Figure of Merit (ZT): ZT = S(exp. 2)sigmaT/κ: where S is the Seebeck coefficient, sigma is the electrical conductivity, T is the absolute temperature, and κ is the thermal conductivity. The parameters are not physically independent, but intrinsically coupled since they are a function of the transport properties of electrons. Traditional research on TE materials has focused on synthesizing bulk semiconductor-type materials that have low thermal conductivity and high electrical conductivity affording ZT values of 1. The optimization of the σ/κ ratio is difficult to achieve using current material formats, as these material constants are complementary. Recent areas of research are focusing on using nanostructural artifacts that introduce specific dislocations and boundary conditions that scatter the phonons. This disrupts the physical link between thermal (phonon) and electrical (electron) transport. The result is that κ is decreased without decreasing σ. These material formats give ZT values of up to 2 which represent approximately 18% energy gain from waste heat recovery. The next challenge in developing the next generation of TE materials with superior performance is to tailor the interconnected thermoelectric physical parameters of the material system. In order to approach this problem, the fundamental physics of each parameter S, sigma, and κ need to be physically understood in their context of electron/phonon interaction for the construction of new high ZT thermoelectric devices. Is it possible to overcome the physical limit imposed by of the effect of phonon lattice oscillation and energetic electrons towards thermal conductivity? Is the Seebeck coefficient, based on the difference in voltage over temperature gradient ( deltaV/deltaT), an intrinsic parameter of each material? All these parameters were manipulated using nano-bridge and twin-lattice structural concepts at the NASA Langley Research Center. This talk will review the current trend of TE research to optimize the ZT and discuss about new approaches on increasing ZT within functional limits of each parameter.

Kim, Hyun Jung↗
NASA NTRSDate not supplied

Three-Dimensional In Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions, which shed new light on the selection of dendritic array structures, and on roles of macroscopic curvatures and subgrain boundaries between grains with a small misorientation with respect to the temperature gradient in the spatiotemporal evolution of the primary cellular spacing.

Kaihua Ji↗
DOE OSTI2023

HELEN: Traveling Wave SRF Linear Collider Higgs Factory

Traveling wave SRF accelerating structures offer several advantages over the traditional standing wave structures: substantially lower $H_pk/E_acc$ and lower $E_pk/E_acc$, ratios of peak magnetic field and peak electric field to the accelerating gradient, respectively, together with substantially higher $R/Q$. In this paper we discuss how a linear collider Higgs Factory HELEN can be built using TW-based SRF linacs. We cover a plan to address technological challenges and describe ways to upgrade the collider luminosity and energy.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Chemical cartography with LAMOST and Gaia reveal azimuthal and spiral structure in the Galactic disc

Chemical Cartography, or mapping, of our Galaxy has the potential to fully transform our view of its structure and formation. In this work, we use chemical cartography to explore the metallicity distribution of OBAF-type disc stars from the LAMOST survey and a complementary sample of disc giant stars from Gaia DR3. We use these samples to constrain the radial and vertical metallicity gradients across the Galactic disc. We also explore whether there are detectable azimuthal variations in the metallicity distribution on top of the radial gradient. For the OBAF-type star sample from LAMOST, we find a radial metallicity gradient of Δ[Fe/H]/ΔR ~–0.078 ± 0.001 dex kpc –1 in the plane of the disc and a vertical metallicity gradient of Δ[Fe/H]/ΔZ ~–0.15 ± 0.01 dex kpc –1 in the solar neighbourhood. The radial gradient becomes shallower with increasing vertical height, while the vertical gradient becomes shallower with increasing Galactocentric radius, consistent with other studies. We also find detectable spatially dependent azimuthal variations on top of the radial metallicity gradient at the level of ~0.10 dex. Interestingly, the azimuthal variations appear be close to the Galactic spiral arms in one data set (Gaia DR3) but not the other (LAMOST). Furthermore, these results suggest that there is azimuthal structure in the Galactic metallicity distribution and that in some cases it is co-located with spiral arms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

MOOSE-Based Fast Reactor Core Bowing Capabilities: Coupled Structural Mechanics – Thermal Fluids Demonstration and Related Verification Efforts

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. Verification and assessment efforts continued on the Multiphysics Object Oriented Simulation Environment (MOOSE) capabilities relevant for modeling thermo-mechanical core bowing behavior. IAEA Verification Problem 4, which was started in FY23, was further examined with MOOSE capabilities to rectify discrepancies observed in previous years when compared to IAEA benchmark participant data. Meshing and postprocessing capabilities in MOOSE were also advanced by other teams and utilized this year. A thermal fluids-structural mechanical coupling demonstration has performed on 7-assemblyand 19-assembly fast reactor assembly configurations using MOOSE. Subchannel capabilities are used to calculate coolant temperature, and heat conduction capabilities calculate duct wall temperature as well as heat transfer through the inter-assembly gap. Structural mechanical capabilities then deform the mesh, accounting for contact between assemblies, according to the temperature gradients calculated by the thermal solvers. Power distributions are imposed rather than calculated to demonstrate different deformations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Analysis of the HAE activity in the TJ-II stellarator using a Landau closure model

The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.

Plasma confinement↗
DOE OSTI2026

Structural uniformity and compositional homogeneity of solid-phase alloyed rod

Solid-phase processes have emerged as an alternative to fusion-based alloying to avoid coarse microstructures, undesirable phase formation, and high energy consumption. However, achieving uniform distribution of alloying elements during friction-based processing remains challenging due to highly heterogeneous thermomechanical conditions. This work evaluates the structural uniformity and compositional homogeneity of Al–Cu–Zn alloyed rods produced by friction extrusion (FE) and establishes the role of the rotational speed to feed rate ratio (N/V) on alloying effectiveness. A systematic matrix of FE experiments was conducted at constant extrusion ratio with N/V values ranging from 3.7 to 300. Compositional uniformity was assessed along the rod length (ICP-OES), in three dimensions (X-ray computed tomography), and at the microscale (SEM–EDS), supported by a gray-level co-occurrence matrix (GLCM)–based homogeneity metric. Smoothed particle hydrodynamics (SPH) simulations were used to reveal material flow and thermomechanical fields. Results show that N/V = 100 produces a high-shear mixing zone that eliminates the unmixed core and enables near-full dissolution and dispersion of Cu and Zn. At lower N/V, a laminar flow region persists at the rod center, causing segregation and large composition gradients. The combined experimental–computational analysis provides mechanistic insight into the transition from fragmented particle dispersion to thermomechanically assisted metallurgical mixing. This study establishes processing–structure relationships for solid-phase alloying and provides guidance for achieving homogenized compositions comparable to wrought alloys via rapid, scalable FE processing.

Aluminum↗
NASA NTRS2018

Mangrove Canopy Height Globally Related to Precipitation and Cyclone Frequency

Mangrove wetlands are among the most productive and carbon-dense ecosystems in the world. Their structural attributes vary considerably across spatial scales, yielding large uncertainties in regional and global estimates of carbon stocks. Here, we present a global analysis of mangrove canopy height gradients and aboveground carbon stocks based on remotely sensed measurements and field data. Our study highlights that precipitation, temperature and cyclone frequency explain 74% of the global trends in maximum canopy height, with other geophysical factors influencing the observed variability at local and regional scales. We find the tallest mangrove forests in Gabon, equatorial Africa, where stands attain 62.8 m. The total global mangrove carbon stock (above- and belowground biomass, and soil) is estimated at 5.03 Pg, with a quarter of this value stored in Indonesia. Our analysis implies sensitivity of mangrove structure to climate change, and offers a baseline to monitor national and regional trends in mangrove carbon stocks.

Carbon cycle↗
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Dawn-dusk asymmetries in the low-latitude boundary layer arising from the Kelvin-Helmholtz instability: A particle simulation

Along the low-latitude boundary layer (LLBL), the Kelvin-Helmholtz (K-H) instability can provide a means for injection of solar wind plasma across closed field lines of the magnetosphere. A fully electromagnetic, two-dimensional (three-velocity) particle code is used to investigate dawn-dusk asymmetries that can arise from velocity differences due to gradient drifts and electric field gradients, and ion acceleration across the narrow field transition layers in the dawn and dusk flanks. The model includes the dawn and dusk sides of the LLBL simultaneously in a slab geometry, incorporating seperate populations of ions and electrons for the magnetosphere and the magnetosheath plasmas. We report several effects: (1) asymmetries in the observed morphology of the turbulent structures, with familiar fluidlike vortex formation on the duskside, but tongues of magnetosheath plasma penetrating into the magnetosphere on the dawnside; (2) the formation of discrete current layers, characterized by strong currents and sharp gradients in the magnetic field, and discrete charge layers, having net charges and constant, weaker currents; (3) increasing asymmetry in dawn/dusk behavior with a decrease of initial currents; (3) increasing asymmetry in dawn/dusk behavior with a decrease of initial boundary layer width; and (4) enhancement of the dawn-to-dusk electric field as magnetosheath particles and momenta are transported across the magnetopause.

Wilber, M.↗
DOE OSTI2026

Solid state 59 Co NMR study of a Np(VII) compound

Np(VII) compounds with [Co(NH 3 ) 6 ] 3+ cations were synthesized and structurally examined using powder X-ray diffraction. Multiple phases were observed, consisting of octahedral [Co(NH 3 ) 6 ] 3+ cations, discrete tetragonal bipyramidal [NpO 4 (OH) 2 ] 3− anions, and waters of hydration. Electric field gradient tensors at Co sites were measured by solid state 59 Co nuclear magnetic resonance (NMR) spectroscopy and compared with theoretical calculations. The relative contributions of the chemical shift and electric field tensors to the NMR lineshape were determined by recording spectra at field strengths of 7.04 and 11.74 Tesla. Further, the evolution of structure and morphology as a function of sample age and the effects on NMR spectral parameters has also been investigated. These results demonstrate the use of NMR at multiple fields to expand understanding of the stability and electronic structure of high valent neptunium compounds.

heptavalent neptunium↗
DOE OSTI2023

Predictive understanding of the surface tension and velocity of sound in ionic liquids using machine learning

Knowledge of the physical properties of ionic liquids (ILs), such as the surface tension and speed of sound, is important for both industrial and research applications. Unfortunately, technical challenges and costs limit exhaustive experimental screening efforts of ILs for these critical properties. Previous work has demonstrated that the use of quantum-mechanics-based thermochemical property prediction tools, such as the conductor-like screening model for real solvents, when combined with machine learning (ML) approaches, may provide an alternative pathway to guide the rapid screening and design of ILs for desired physiochemical properties. However, the question of which machine-learning approaches are most appropriate remains. In the present study, we examine how different ML architectures, ranging from tree-based approaches to feed-forward artificial neural networks, perform in generating nonlinear multivariate quantitative structure–property relationship models for the prediction of the temperature- and pressure-dependent surface tension of and speed of sound in ILs over a wide range of surface tensions (16.9–76.2 mN/m) and speeds of sound (1009.7–1992 m/s). The ML models are further interrogated using the powerful interpretation method, shapley additive explanations. We find that several different ML models provide high accuracy, according to traditional statistical metrics. The decision tree-based approaches appear to be the most accurate and precise, with extreme gradient-boosting trees and gradient-boosting trees being the best performers. However, our results also indicate that the promise of using machine-learning to gain deep insights into the underlying physics driving structure–property relationships in ILs may still be somewhat premature.

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

Direct laser writing of volumetric gradient index lenses and waveguides

Abstract Direct laser writing (DLW) has been shown to render 3D polymeric optical components, including lenses, beam expanders, and mirrors, with submicrometer precision. However, these printed structures are limited to the refractive index and dispersive properties of the photopolymer. Here, we present the subsurface controllable refractive index via beam exposure (SCRIBE) method, a lithographic approach that enables the tuning of the refractive index over a range of greater than 0.3 by performing DLW inside photoresist-filled nanoporous silicon and silica scaffolds. Adjusting the laser exposure during printing enables 3D submicron control of the polymer infilling and thus the refractive index and chromatic dispersion. Combining SCRIBE’s unprecedented index range and 3D writing accuracy has realized the world’s smallest (15 µm diameter) spherical Luneburg lens operating at visible wavelengths. SCRIBE’s ability to tune the chromatic dispersion alongside the refractive index was leveraged to render achromatic doublets in a single printing step, eliminating the need for multiple photoresins and writing sequences. SCRIBE also has the potential to form multicomponent optics by cascading optical elements within a scaffold. As a demonstration, stacked focusing structures that generate photonic nanojets were fabricated inside porous silicon. Finally, an all-pass ring resonator was coupled to a subsurface 3D waveguide. The measured quality factor of 4600 at 1550 nm suggests the possibility of compact photonic systems with optical interconnects that traverse multiple planes. SCRIBE is uniquely suited for constructing such photonic integrated circuits due to its ability to integrate multiple optical components, including lenses and waveguides, without additional printed supports.

36 MATERIALS SCIENCE↗