Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential geometry”

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

Fluoropolymer Aging Phenomena

This project aims to investigate structural and morphological changes in fluoropolymers induced by both processing and aging phenomena such as heat, time, and irradiation. Fluoropolymers are an important class of thermoplastics that are broadly used in industry as o-rings and seals when chemical resistance and thermal performance are important. They are also often used in thin film geometries as binders in batteries, insulation layers, water barriers, and anti-fouling coatings. For all of these applications, mechanical integrity and aging are important aspects of their use. Changes in the crystalline morphology are known to impact the mechanical properties of the polymer and can lead to failure. The broad goals of this project are to correlate compositional changes, morphological changes, and mechanical changes. report focuses on the morphological changes that are correlated with nuclear magnetic resonance (NMR), calorimetry and mechanical testing. We plan to compare structure and morphology measurements across complementary bulk and surface techniques. We use in situ atomic force microscopy (AFM) to measure domain shape and crystallization kinetics; in situ NMR to measure crystalline composition; and differential scanning calorimetry (DSC) to measure melting temperatures.

36 MATERIALS SCIENCE↗
NASA NTRS2005

Vibration Based Crack Detection in a Rotating Disk: Experimental Results - Part 2

This paper describes the experimental results concerning the detection of a crack in a rotating disk. The goal was to utilize blade tip clearance and shaft vibration measurements to monitor changes in the system's center of mass and/or blade deformation behaviors. The concept of the approach is based on the fact that the development of a disk crack results in a distorted strain field within the component. As a result, a minute deformation in the disk's geometry as well as a change in the system's center of mass occurs. Here, a notch was used to simulate an actual crack. The vibration based experimental results failed to identify the existence of a notch when utilizing the approach described above, even with a rather large, circumferential notch (l.2 in.) located approximately mid-span on the disk (disk radius = 4.63 in. with notch at r = 2.12 in.). This was somewhat expected, since the finite element based results in Part 1 of this study predicted changes in blade tip clearance as well as center of mass shifts due to a notch to be less than 0.001 in. Therefore, the small changes incurred by the notch could not be differentiated from the mechanical and electrical noise of the rotor system. Although the crack detection technique of interest failed to identify the existence ofthe notch, the vibration data produced and captured here will be utilized in upcoming studies that will focus on different data mining techniques concerning damage detection in a disk.

Gyekenyesi, Andrew L.↗
DOE OSTI2023

Theory of an Automatic Seepage Meter and Ramifications for Applications

A new approach for measuring fluxes across surface water—groundwater interfaces was recently proposed. The Automatic Seepage Meter (ASM) is equipped with a precise water level sensor and digital memory that analyzes water level time series in a vertical tube inserted into a streambed. The ability to infer flux values with high temporal resolution relies on an accurate interpretation of water level dynamics inside the tube. Here, we reduce the three-dimensional hydrodynamic problem that describes the ASM water level in a variety of field conditions to a single ordinary differential equation. This novel general analytical solution for estimating ASM responses is more comprehensive and flexible than previous approaches and is applicable to the entire range of field conditions, including steady or transient stream stages, evaporation, rainfall, and noise. For example, our analysis determines the timing of the nonmonotonic ASM response to a monotonic linear stream stage variation and explains previously used empirical parabolic approximation for estimating fluxes. We present algorithms for simultaneous inference of vertical interface flux and hydraulic conductivity values together with an example code. We quantify how the accuracy of parameter estimation depends on test duration and noise amplitude and propose how our analysis can be used to optimize field test protocols. On this basis, changing the ASM geometry by increasing the radius and decreasing tube insertion depth may enable ASM field test protocols that estimate interface flux and hydraulic conductivity faster while maintaining desired accuracy. Potential applications of joint parameter estimation are suggested.

47 OTHER INSTRUMENTATION↗
NASA NTRS1989

Contrasting landform perception with varied radar illumination geometries and at simulated resolutions of Venera and Magellan

The high sensitivity of imaging radars to slope at moderate to low incidence angles enhances the perception of linear topography on images. It reveals broad spatial patterns that are essential to landform mapping and interpretation. As radar responses are strongly directional, the ability to discriminate linear features on images varies with their orientation. Landforms that appear prominent on images where they are transverse to the illumination may be obscure to indistinguishable on images where they are parallel to it. Landform detection is also influenced by the spatial resolution in radar images. Seasat radar images of the Gran Desierto Dunes complex, Sonora, Mexico; the Appalachian Valley and Ridge Province; and accreted terranes in eastern interior Alaska were processed to simulate both Venera 15 and 16 images (1000 to 3000 km resolution) and image data expected from the Magellan mission (120 to 300 m resolution. The Gran Desierto Dunes are not discernable in the Venera simulation, whereas the higher resolution Magellan simulation shows dominant dune patterns produced from differential erosion of the rocks. The Magellan simulation also shows that fluvial processes have dominated erosion and exposure of the folds.

Ford, J. P.↗
NASA NTRS2018

Instability Wave–Streak Interactions in a High Mach Number Boundary Layer at Flight Conditions

The interaction of stationary streaks undergoing nonmodal growth with modally unstable instability waves in a hypersonic boundary-layer flow is studied using numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment; namely, a 7 degree half-angle, circular cone with 2.5 mm nose radius, freestream Mach number equal to 5.30, unit Reynolds number equal to 13.42 m -1 , and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the vehicle. This paper investigates the nonlinear evolution of initially linear optimal disturbances that evolve into finite-amplitude streaks, followed by an analysis of the modal instability characteristics of the perturbed, streaky boundary-layer flow. The investigation is performed with stationary direct numerical simulations (DNS) and plane-marching parabolized stability equations (PSE), in conjunction with partial-differential-equation-based planar eigenvalue analysis. The overall effect of streaks is to reduce the peak amplification factors of instability waves, indicating a possible downstream shift in the onset of laminar-turbulent transition. The present study confirms previous findings that the mean flow distortion of the nonlinear streak perturbation reduces the amplification rates of the Mack-mode instability. More importantly, however, the present results demonstrate that the spanwise varying component of the streak can produce a larger effect on the Mack-mode amplification. The study with selected azimuthal wavenumbers for the stationary streaks reveals that a wavenumber of approximately 1.4 times larger than the optimal wavenumber is more effective in stabilizing the planar Mack-mode instabilities. In the absence of unstable first-mode waves for the present cold-wall condition, transition onset is expected to be delayed until the peak streak amplitude increases to nearly 35 percent of the freestream velocity, when intrinsic instabilities of the boundary-layer streaks begin to dominate the transition process. For streak amplitudes below that limit a significant net stabilization is achieved, yielding a potential transition delay that can exceed 100 percent of the length of the laminar region in the uncontrolled case.

Boundary layer stability↗
DOE OSTI2026

The Structure-Properties Relationship of Alternative Bismaleimide Variants for Candidacy for Additive Manufacturing

Modernizing the manufacturing of high-performance polymer foams such as amino-poly(oxadiazole) bismaleimide (APO-BMI), a bismaleimide resin with superior thermal and compressive strength that incorporates additive manufacturing (AM) techniques, is crucial for its applications, but the parameters for AM can be challenging based on the physical properties of the associated monomer. For our applications, selective laser sintering (SLS) is typically used. SLS is a 3D printing technique that allows for complex shapes and geometries without structural supports while also providing high resolution material. However, printing thermosets like APO-BMI with SLS is challenging due to the complex melting and curing considerations required when selecting parameters. Additionally, the temperature difference between melting and curing of APO-BMI is over a hundred º C, which makes selecting a sintering window especially difficult. This work explores structural modifications of APO-BMI that may be more amenable for selective laser sintering. The effects of how different structural changes such as substitution pattern, heteroatom identity in the bridge, and bridge length affect the thermal properties of the material were also compared. All APO variants were found to have a smaller temperature window between the melting and curing peaks based on differential scanning calorimetry (DSC) which is advantageous for SLS. Small structural changes significantly altered the melting and curing properties of APO. Additionally, DSC revealed significant polymorphisms in APO-BMI and other APO variants which could be attributed to differences in thermal history and would need to be considered when adapting for SLS.

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

CCD photometry of 2060 Chiron, 1991 January

Observations of 2060 Chiron was performed on 7 to 8 Jan. 1991 with the Mt. Palomar 1.52 m telescope in the Gunn-R passband. On-chip field stars were used to perform differential reductions. The repeatability of the 5.9 hour light curve was excellent, both within a night and from night to night. No evidence for short-term secular variations similar to those seen last year by both Luu and Jewitt (1990) and Buratti and Dunbar (1991) is seen in the new light curve. Chiron's rotational light curve appears strikingly similar to that obtained a year earlier by Luu and Jewitt (1990), both in amplitude and shape. Both light curves show strongly correlated changes over a timescale of perhaps 15 minutes. These same features were marginally visible in the 1986 light curve. Such behavior is believed to be evidence that Chiron may be more aspherical than the 4 percent intensity variation might otherwise indicate, and favors a viewing geometry where the subearth latitude is rather low. Chiron was much fainter in 1985, when a partial light curve was obtained by Marcialis. Due to the lower sampling rate of these early data, no conclusions can be made regarding the high-frequency light curve structure back then. All three of these light curves differ significantly from that obtained by Buratti and Dunbar (1991), one week before the observations of Luu and Jewitt. The Chiron field was calibrated using Landolt standards on Ut 15 Mar. 1991. A mean R-magnitude of 15.6 + or - 0.1 was found. Variability of 2060 Chiron was demonstrated over timescales of minutes, hours, and years. An intense campaign was urged to monitor the photometric behavior of Chiron throughout the 1990s.

Buratti, B. J.↗
NASA NTRS1999

Early root cap development and graviresponse in white clover (Trifolium repens) grown in space and on a two-axis clinostat

White clover (Trifolium repens) was germinated and grown in microgravity aboard the Space Shuttle (STS-60, 1994; STS-63, 1995), on Earth in stationary racks and in a slow-rotating two-axis clinostat. The objective of this study was to determine if normal root cap development and early plant gravity responses were dependent on gravitational cues. Seedlings were germinated in space and chemically fixed in orbit after 21, 40, and 72 h. Seedlings 96 h old were returned viable to earth. Germination and total seedling length were not dependent on gravity treatment. In space-flown seedlings, the number of cell stories in the root cap and the geometry of central columella cells did not differ from those of the Earth-grown seedlings. The root cap structure of clinorotated plants appeared similar to that of seedlings from microgravity, with the exception of three-day rotated plants, which displayed significant cellular damage in the columella region. Nuclear polarity did not depend on gravity; however, the positions of amyloplasts in the central columella cells were dependent on both the gravity treatment and the age of the seedlings. Seedlings from space, returned viable to earth, responded to horizontal stimulation as did 1 g controls, but seedlings rotated on the clinostat for the same duration had a reduced curvature response. This study demonstrates that initial root cap development is insensitive to either chronic clinorotation or microgravity. Soon after differentiation, however, clinorotation leads to loss of normal root cap structure and plant graviresponse while microgravity does not.

manned↗
DOE OSTI2021

Bayesian sparse learning with preconditioned stochastic gradient MCMC and its applications

Deep neural networks have been successfully employed in an extensive variety of research areas, including solving partial differential equations. Despite its significant success, there are some challenges in effectively training DNN, such as avoiding overfitting in over-parameterized DNNs and accelerating the optimization in DNNs with pathological curvature. Here, we propose a Bayesian type sparse deep learning algorithm. The algorithm utilizes a set of spike-and-slab priors for the parameters in the deep neural network. The hierarchical Bayesian mixture will be trained using an adaptive empirical method. That is, one will alternatively sample from the posterior using preconditioned stochastic gradient Langevin Dynamics (PSGLD), and optimize the latent variables via stochastic approximation. The sparsity of the network is achieved while optimizing the hyperparameters with adaptive searching and penalizing. A popular SG-MCMC approach is Stochastic gradient Langevin dynamics (SGLD). However, considering the complex geometry in the model parameter space in nonconvex learning, updating parameters using a universal step size in each component as in SGLD may cause slow mixing. To address this issue, we apply a computationally manageable preconditioner in the updating rule, which provides a step-size parameter to adapt to local geometric properties. Moreover, by smoothly optimizing the hyperparameter in the preconditioning matrix, our proposed algorithm ensures a decreasing bias, which is introduced by ignoring the correction term in the preconditioned SGLD. According to the existing theoretical framework, we show that the proposed algorithm can asymptotically converge to the correct distribution with a controllable bias under mild conditions. Numerical tests are performed on both synthetic regression problems and learning solutions of elliptic PDE, which demonstrate the accuracy and efficiency of the present work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A simple wire-coil resistive heater for high temperature radial x-ray diffraction in a diamond anvil cell

Diamond anvil cells are commonly used at synchrotron x-ray diffraction beamlines to study structural and thermoelastic properties of materials at high pressures. In a radial geometry, where the x-ray probe is oriented perpendicular to the axis of force, the deformation and strength of materials can be measured in situ. Because the anelastic and failure properties of materials depend strongly on temperature, many applications would benefit from the ability to measure high pressure radial diffraction in elevated and accurately controlled thermal environments. Previous work to introduce high temperature to radial diamond anvil cells has been largely limited to laser heating, with relatively scant efforts to resistively heat the sample. Here, we report a relatively straightforward adaptation of a simple wire coil heater, with in situ high-temperature radial diffraction performed on tungsten carbide up to 573 K at beamline 12.2.2 of the Advanced Light Source. In conclusion, the results demonstrate that the differential stress supported by WC decreases with increasing temperature: the differential stress on the basal (001) and pyramidal (101) planes decreased 6.6% and 5.5%, respectively, while the (100) plane only saw a 2.7% decrease, in agreement with previous studies.

Diamond anvil cells↗
NASA NTRS2004

Nonlinear Oscillations and Flow of Gas Within Closed and Open Conical Resonators

A dissonant acoustic resonator with a conical shaped cavity was tested in four configurations: (A) baseline resonator with closed ends and no blockage; (B) closed resonator with internal blockage; (C) ventilated resonator with no blockage; and (D) ventilated resonator with an applied pressure differential. These tests were conducted to investigate the effects of blockage and ventilation holes on dynamic pressurization. Additionally, the investigation was to determine the ability of acoustic pressurization to impede flow through the resonator. In each of the configurations studied, the entire resonator was oscillated at the gas resonant frequency while dynamic pressure, static pressure, and temperature of the fluid were measured. In the final configuration, flow through the resonator was recorded for three oscillation conditions. Ambient condition air was used as the working fluid. The baseline results showed a marked reduction in the amplitude of the dynamic pressure waveforms over previously published studies due to the use of air instead of refrigerant as the working fluid. A change in the resonant frequency was recorded when blockages of differing geometries were used in the closed resonator, while acoustic pressure amplitudes were reduced from baseline measurements. A sharp reduction in the amplitude of the acoustic pressure waves was expected and recorded when ventilation ports were added. With elevated pressure applied to one end of the resonator, flow was reduced by oscillating the cavity at the fluid fundamental resonant frequency compared to cases without oscillation and oscillation off-resonance.

Daniels, Christopher↗
DOE OSTI2026

Directional dependence of equilibrated TLD-400 chips in multiple radiation fields

Fielding of thermoluminescent dosimeters (TLDs) for measurement of photon radiation dose in experiments is the standard practice for γ-irradiation facilities, pulsed power x-ray facilities, and reactor facilities at Sandia National Laboratories (SNL). Due to the high-dose experimental conditions and the mixed ( 1 n, γ) fields in these facilities, SNL radiation metrologists have historically used CaF 2 :Mn TLDs (also known as TLD-400). Recent inquiries to the radiation metrology staff have raised concerns that the aluminum-equilibrated TLD-400 chips may exhibit a directionally dependent response. The metrologists were asked whether the dose measured by the chip may be impacted by the angle of incidence on the equilibrated TLD. To provide a thorough answer to this query, a set of adjoint Monte Carlo radiation transport calculations was performed for three different equilibrated TLDs as well as bare TLD chips using the Integrated Tiger Series (ITS) code. The key feature of the adjoint calculations performed in this study is that each photon escaping the modeled geometry was tallied into angular bins to provide information on both the energy and angular dependence of the equilibrated chip. After the adjoint Monte Carlo calculations were completed, the resulting energy-dependent response function for each angular bin was convolved with multiple photon energy spectra representing various radiation facilities at SNL. The radiation facilities selected for analysis span a range of photon energies from approximately 1 keV up to approximately 20 MeV. Thus, the presented results are applicable to a wide variety of radiation facilities around the world. Although the bare TLD-400 chip was expected to display the largest variation due to the photon angle of impact on the dosimeter, the dosimeter with the thinnest aluminum equilibrator (SNL thin equilibrated TLD-400) was determined to have the biggest differential between the impact angle with the maximum dose ( D max ) and the impact angle with the minimum dose ( D min ). However, the SNL normal equilibrated TLD-400 chip demonstrated a dramatic reduction in that differential between maximum and minimum dose angles. The reduction in this differential is one of the dominant factors in experimenters’ choice to field these dosimeters at SNL radiation facilities. The results from the PNNL (Hanford) energy-flattening field capsule are consistent with the previous publications. The differential dose responses due to the photon impact angle indicate that experimenters should strive to field their TLD-400 dosimeters in a consistent manner to avoid additional uncertainty in the measurements based on dosimeter orientation.

Adjoint radiation transport↗
NASA NTRS1994

Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. These algorithms are now used in the practical engineering design of devices such as cars and airplanes as well as theoretical studies of complex phenomena such as fluid turbulence. In past years, limited computer resources have motivated the development of efficient numerical methods in computational fluid dynamics (CFD) utilizing structured meshes. These meshes are comprised of systematic arrays of quadrilateral or hexahedral cells. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Structured meshes also permit the use of highly efficient solution techniques such as alternating direction implicit (ADI) iteration schemes or multigrid. Following the dramatic improvement in computing speed in recent years, emphasis has shifted towards the design of algorithms capable of treating complex geometries. The automatic generation of structured grids about complex geometries is problematic. Unstructured grids offer one promising alternative technique for treating these general geometries. Unstructured meshes have irregular connectivity and usually contain triangles and/or quadrilaterals in two dimensions and tetrahedra and/or hexahedra in three dimensions. The generation and use of unstructured grids poses new challenges in computational fluid dynamics. This is true for both grid generation as well as for the design of algorithms for flow solution. The purpose of these notes is to present recent developments in the unstructured grid generation and flow solution technology.

Barth, T. J.↗
DOE OSTI2025

One-shot learning for solution operators of partial differential equations

Learning and solving governing equations of a physical system, represented by partial differential equations (PDEs), from data is a central challenge in many areas of science and engineering. Traditional numerical methods can be computationally expensive for complex systems and require complete governing equations. Existing data-driven machine learning methods require large datasets to learn a surrogate solution operator, which could be impractical. Here, we propose a solution operator learning method that requires only one PDE solution, i.e., one-shot learning, along with suitable initial and boundary conditions. Leveraging the locality of derivatives, we define a local solution operator in small local domains, train it using a neural network, and use it to predict solutions of new input functions via mesh-based fixed-point iteration or meshfree neural-network based approaches. We test our method on various PDEs, complex geometries, and a practical spatial infection spread application, demonstrating its effectiveness and generalization capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Design of experiments to spectroscopically characterize radiation flow in stochastic media

Precise characterization of experimental radiation flow is required to validate the high energy density physics models, numerical methods, and codes that are used to simulate radiation-hydrodynamics phenomena such as thermal radiation transport in stochastic media. The Cassio code is used to simulate thermal radiation flow through inhomogeneous, stochastic-media-foam configurations containing optically thick clumps dispersed within an optically thin background aerogel. Cassio can model small inhomogeneous problems directly, but most problems require approximations to meet computer limitations on run-times and memory usage. Various examples of these approximations are methods that produce, in one calculation, an ensemble-averaged solution and associated standard deviation; reduced spatial dimensionality with approximate geometries; and full material homogenization with no geometric detail. Cassio simulations are used to design experiments at the OMEGA-60 Laser Facility that can measure the radiation flow using the spatially resolved COAX absorption spectroscopy diagnostic. The experimental platforms flow radiation through foam targets ranging from a background-only aerogel, to a single configuration of a specified stochastic medium, to a fully homogenized foam of the background and clump materials. Under constant total clump mass, larger clumps (here, larger than 10 μm diameter) will mix more slowly with the background such that the bulk radiation flow is faster than it would be in a fully homogenized material. The COAX platform can be used to infer temperature and density profiles in both the background material and clumps, simultaneously, and therefore to differentiate radiation flow in a range of stochastic and homogeneous media.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

An improved method for the prediction of completely three-dimensional aerodynamic load distributions of configurations with leading edge vortex separation

The application of a higher-order subsonic potential flow panel method to the solution of three-dimensional flow about wing and wing-body combinations with leading-edge vortex separation is presented. The governing equations are the linear flow differential equation and nonlinear boundary conditions which require that the flow be parallel to the wing and body surfaces and that the free vortex sheet, springing from the leading and trailing edges, be aligned with the local flow and support no pressure jump. The vortex core is modeled as a simple line vortex which receives vorticity from the free sheet through a connecting sheet. The Kutta condition is imposed on all appropriate edges of the wing. This set of nonlinear equations is solved by an iterative procedure. The Goethert rule accounts for compressibility. The method has been programmed for the CDC 6600. Delta wings, gothic wings, arrow wings, cambered wings, and wing with body have been analyzed. Initial studies involving variations of panel density, vortex sheet sizing, Jacobian update, and initial geometry demonstrate that the present method generally exhibits good convergence characteristics.

Rubbert, p. E.↗
DOE OSTI2021

Searching water megamasers by using mid-infrared spectroscopy (I): Possible mid-infrared indicators

ABSTRACT Water megamasers at 22 GHz with a gas disc configuration in galaxies provide the most precise measurements of supermassive black hole masses, as well as independent constraints on the Hubble constant in the nearby universe. The existence of other maser types, such as jet or outflow masers, represents another tracer for AGN science. However, the detection rate of water megamasers in galaxies is extremely low. Over 40 yr, only ∼ 160 galaxies are found to harbour maser emission, and ∼30 per cent of them show features in their maser emission that indicate a disc-like geometry. Therefore, increasing the detection rate of masers is a crucial task to allow expanding on maser studies. We present a comparison of mid-infrared spectroscopic data between a maser galaxy sample and a Seyfert 2 control sample. We find that maser galaxies show significant peculiarities in their mid-infrared spectra: (1) Maser galaxies tend to present stronger silicate absorption at τ9.7 $\mu\mathrm{ m}$ than the control sample, (2) PAH 11.3 $\mu$m emission in maser galaxies is much weaker than in the control sample, (3) spectral indices at 20–30 $\mu$m are steeper in maser galaxies than in the control sample and tend to be mid-infrared enhanced population. We conclude that there may be good indicators in mid-infrared and far-infrared which could differentiate maser and non-maser Seyfert 2 galaxies. Upcoming infrared facilities, such as the James Webb Space Telescope, may be able to exploit these and other useful criteria and tracers for water megamaser observations.

Lam, Man I.↗
DOE OSTI2024

Inverse design for waveguide dispersion with a differentiable mode solver

Inverse design of optical components based on adjoint sensitivity analysis has the potential to address the most challenging photonic engineering problems. However, existing inverse design tools based on finite-difference-time-domain (FDTD) models are poorly suited for optimizing waveguide modes for adiabatic transformation or perturbative coupling, which lies at the heart of many important photonic devices. Among these, dispersion engineering of optical waveguides is especially challenging in ultrafast and nonlinear optical applications involving broad optical bandwidths and frequency-dependent anisotropic dielectric material response. In this work, we develop gradient back-propagation through a general-purpose electromagnetic eigenmode solver and use it to demonstrate waveguide dispersion optimization for second harmonic generation with maximized phase-matching bandwidth. This optimization of three design parameters converges in eight steps, reducing the computational cost of optimization by ∼100x compared to exhaustive search and identifying new designs for broadband optical frequency doubling of laser sources in the 1.3–1.4 µm wavelength range. Furthermore, we demonstrate that the computational cost of gradient back-propagation is independent of the number of parameters, as required for optimization of complex geometries. This technique enables practical inverse design for a broad range of previously intractable photonic devices.

Gray, Dodd (ORCID:000000030469599X)↗