Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear processes”

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 199 records · Page 11

Project ORION: Orbital Debris Removal Using Ground-Based Sensors and Lasers

About 100,000 pieces of 1 to 10-cm debris in low-Earth orbit are too small to track reliably but large enough to cripple or destroy spacecraft. The ORION team studied the feasibility of removing the debris with ground-based laser impulses. Photoablation experiments were surveyed and applied to likely debris materials. Laser intensities needed for debris orbit modification call for pulses on the order of lOkJ or continuous wave lasers on the order of 1 MW. Adaptive optics are necessary to correct for atmospheric turbulence. Wavelength and pulse duration windows were found that limit beam degradation due to nonlinear atmospheric processes. Debris can be detected and located to within about 10 microrads with existing radar and passive optical technology. Fine targeting would be accomplished with laser illumination, which might also be used for detection. Bistatic detection with communications satellites may also be possible. We recommend that existing technology be used to demonstrate the concept at a loss of about $20 million. We calculate that an installation to clear altitudes up to 800 km of 1 to 10-cm debris over 2 years of operation would cost about $80 million. Clearing altitudes up to 1,500 km would take about 3 years and cost about $160 million.

Campbell, J. W.↗

The Universe at Moderate Redshift

The report covers the work done in the past year and a wide range of fields including properties of clusters of galaxies; topological properties of galaxy distributions in terms of galaxy types; patterns of gravitational nonlinear clustering process; development of a ray tracing algorithm to study the gravitational lensing phenomenon by galaxies, clusters and large-scale structure, one of whose applications being the effects of weak gravitational lensing by large-scale structure on the determination of q(0); the origin of magnetic fields on the galactic and cluster scales; the topological properties of Ly(alpha) clouds the Ly(alpha) optical depth distribution; clustering properties of Ly(alpha) clouds; and a determination (lower bound) of Omega(b) based on the observed Ly(alpha) forest flux distribution. In the coming year, we plan to continue the investigation of Ly(alpha) clouds using larger dynamic range (about a factor of two) and better simulations (with more input physics included) than what we have now. We will study the properties of galaxies on 1 - 100h(sup -1) Mpc scales using our state-of-the-art large scale galaxy formation simulations of various cosmological models, which will have a resolution about a factor of 5 (in each dimension) better than our current, best simulations. We will plan to study the properties of X-ray clusters using unprecedented, very high dynamic range (20,000) simulations which will enable us to resolve the cores of clusters while keeping the simulation volume sufficiently large to ensure a statistically fair sample of the objects of interest. The details of the last year's works are now described.

Cen, Renyue↗

Partial Wave Analysis of Coupled Photonic Structures

The very high quality factors sustained by microcavity optical resonators are relevant to applications in wavelength filtering, routing, switching, modulation, and multiplexing/demultiplexing. Increases in the density of photonic elements require that attention be paid to how electromagnetic (EM) coupling modifies their optical properties. This is especially true when cavity resonances are involved, in which case, their characteristics may be fundamentally altered. Understanding the optical properties of microcavities that are near or in contact with photonic elements---such as other microcavities, nanostructures, couplers, and substrates---can be expected to advance our understanding of the roles that these structures may play in VLSI photonics, biosensors and similar device technologies. Wc present results from recent theoretical studies of the effects of inter- and intracavity coupling on optical resonances in compound spherical particles. Concentrically stratified spheres and bispheres constituted from homogeneous and stratified spheres are subjects of this investigation. A new formulation is introduced for the absorption of light in an arbitrary layer of a multilayered sphere, which is based on multiple reflections of the spherical partial waves of the Lorenz-Mie solution for scattering by a sphere. Absorption efficiencies, which can be used to profile cavity resonances and to infer fluorescence yields or the onset of nonlinear optical processes in the microcavities, are presented. Splitting of resonances in these multisphere systems is paid particular attention, and consequences for photonic device development and possible performance enhancements through carefully designed architectures that exploit EM coupling are considered.

Fuller, Kirk A.↗

Phase Matching of Diverse Modes in a WGM Resonator

Phase matching of diverse electromagnetic modes (specifically, coexisting optical and microwave modes) in a whispering-gallery-mode (WGM) resonator has been predicted theoretically and verified experimentally. Such phase matching is necessary for storage of microwave/terahertz and optical electromagnetic energy in the same resonator, as needed for exploitation of nonlinear optical phenomena. WGM resonators are used in research on nonlinear optical phenomena at low optical intensities and as a basis for design and fabrication of novel optical devices. Examples of nonlinear optical phenomena recently demonstrated in WGM resonators include low-threshold Raman lasing, optomechanical oscillations, frequency doubling, and hyperparametric oscillations. The present findings regarding phase matching were made in research on low-threshold, strongly nondegenerate parametric oscillations in lithium niobate WGM resonators. The principle of operation of such an oscillator is rooted in two previously observed phenomena: (1) stimulated Raman scattering by polaritons in lithium niobate and (2) phase matching of nonlinear optical processes via geometrical confinement of light. The oscillator is partly similar to terahertz oscillators based on lithium niobate crystals, the key difference being that a novel geometrical configuration of this oscillator supports oscillation in the regime. The high resonance quality factors (Q values) typical of WGM resonators make it possible to achieve oscillation at a threshold signal level much lower than that in a non-WGM-resonator lithium niobate crystal.

Savchenkov, Anatoliy↗

Multi-Step Deep Reactive Ion Etching Fabrication Process for Silicon-Based Terahertz Components

A multi-step silicon etching process has been developed to fabricate silicon-based terahertz (THz) waveguide components. This technique provides precise dimensional control across multiple etch depths with batch processing capabilities. Nonlinear and passive components such as mixers and multipliers waveguides, hybrids, OMTs and twists have been fabricated and integrated into a small silicon package. This fabrication technique enables a wafer-stacking architecture to provide ultra-compact multi-pixel receiver front-ends in the THz range.

Jung-Kubiak, Cecile↗

Nonlinear Görtler Vortices and Their Secondary Instability in a Hypersonic Boundary Layer

Nonlinear development of the Görtler instability over a concave surface gives rise to a highly distorted inflectional flow field in the boundary layer that leads to both wall-normal and spanwise gradients in the flow. Such nonlinear structures are susceptible to strong, high-frequency secondary instabilities that may lead to the onset of laminar-turbulent transition. The present numerical study uses direct numerical simulations and linear secondary instability theory to investigate finite amplitude Görtler vortices and their secondary instability characteristics, respectively, in the hypersonic flow over an axisymmetric cone with a concave aft body. To complement previous studies in the literature wherein the Görtler instability was usually studied for a flat plate and initiated at some upstream location by imposing an eigenfunction as the inflow condition or by blowing and suction at the wall, the present investigation is focused on fully realizable Görtler instability that is excited by an azimuthally periodic array of surface protuberances. Furthermore, while the previous work had mostly focused on the secondary instability of Görtler vortices with cross-plane velocity contours that resembled bell-shaped structures, the present results confirm that fully developed mushroom structures also exist in the hypersonic regime when the Görtler vortex amplitude is sufficiently large. Computations further reveal that the dominant modes of secondary instability in these mushroom-shaped structures correspond to an antisymmetic (i.e., sinuous) “stem” mode that concentrates within the strong, nearly wall-normal internal shear layers surrounding the stem regions underneath the caps of the mushroom structures. Additionally, there exist a multitude of other significantly unstable secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes originating from the Mack mode instability play an important role during the nonlinear breakdown process.

Laminar turbulent transition↗

Revisiting the Effectiveness Of HCHO/NO2 Ratios for Inferring Ozone Sensitivity to Its Precursors Using High Resolution Airborne Remote Sensing Observations in A High Ozone Episode During the KORUS-AQ Campaign

The nonlinear chemical processes involved in ozone production (P(O3)) have necessitated using proxy indicators to convey information about the primary dependence of P(O3) on volatile organic compounds (VOCs) or nitrogen oxides (NOx). In particular, the ratio of remotely sensed columns of formaldehyde (HCHO) to nitrogen dioxide (NO2) has been widely used for studying O3 sensitivity. Previous studies found that the errors in retrievals and the incoherent relationship between the column and the near-surface concentrations are a barrier in applying the ratio in a robust way. In addition to these obstacles, we provide calculational-observational evidence, using an ensemble of 0-D photochemical box models constrained by DC-8 aircraft measurements on an ozone event during the Korea-United States Air Quality (KORUS-AQ) campaign over Seoul, to demonstrate the chemical feedback of NO2 on the formation of HCHO is a controlling factor for the transition line between NOx-sensitive and NOx-saturated regimes. A fixed value (~2.7) of the ratio of the chemical loss of NOx (LNOx) to the chemical loss of HO2+RO2 (LROx) perceptibly differentiates the regimes. Following this value, data points with a ratio of HCHO/NO2 less than 1 can be safely classified as NOx-saturated regime, whereas points with ratios between 1 and 4 fall into one or the other regime. We attribute this mainly to the HCHO-NO2 chemical relationship causing the transition line to occur at larger (smaller) HCHO/NO2 ratios in VOC-rich (VOC-poor) environments. We then redefine the transition line to LNOx/LROx~2.7 that accounts for the HCHO-NO2 chemical relationship leading to HCHO = 3.7 × (NO2 – 1.14 × 1016 molec.cm-2). Although the revised formula is locally calibrated (i.e., requires for readjustment for other regions), its mathematical format removes the need for having a wide range of thresholds used in HCHO/NO2 ratios that is a result of the chemical feedback. Therefore, to be able to properly take the chemical feedback into consideration, the use of HCHO = a × (NO2 – b) formula should be preferred to the ratio in future works. We then use the Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument to study O3 sensitivity in Seoul. The unprecedented spatial (250 × 250 m2) and temporal (~every 2 h) resolutions of HCHO and NO2 observations form the sensor enhance our understanding of P(O3) in Seoul; rather than providing a crude label for the entire city, more in-depth variabilities in chemical regimes are observed that should be able to inform mitigation strategies correspondingly.

KORUS-AQ; formaldehyde; NOx; VOC; OMI; GeoTASO; oz↗

Exploring the Capabilities of a Machine Learning Algorithm to Detect Space Weather-Significant Emerging Active Regions

Active regions are a source of various phenomena responsible for Space Weather disturbances; therefore, developing a technology for early warning about upcoming magnetic activity is crucial to mitigate its impact. However, observational limitations and the high nonlinearity of processes associated with the accumulation of magnetic flux and its interaction with the surrounding plasma during the emergence through the convection zone make early activity detection a challenging problem. To address these challenges, we developed a physics-driven machine learning model that allows us to detect active regions (ARs) before they become visible on the solar surface by analyzing the power spectra of acoustic oscillations observed by the SDO/HMI instrument. This study is based on a time series of Doppler shift maps of 31x31-degree areas tracked with the Carrington rotation rate for four days before and after the emergence. The Doppler shift time series are processed into the oscillation power maps for four frequency ranges and accompanied by line-of-sight magnetograms and the continuum intensity maps from SDO/HMI. The resulting data are converted into a 1D time series representing the mean temporal variations of these quantities. The redacted time series are used as input to predict AR emergence using the Long Short Term Memory (LSTM) method. The training of the LSTM model is based on 40 ARs, which includes an independent analysis for each sub region that exhibits AR emergence or remains quiet. The emergence of magnetic flux (defined as a decrease of the continuum intensity) was detected with the developed LSTM algorithm from 5 to 48 hours before the reported time by NOAA. The developed model is capable of pointing to the time and location of active region formation. In this presentation, we discuss reasons that impact how early in advance the model can identify the upcoming activity and the possibility of improving the current predictive skills and steps to transition to the operational forecast.

Heliophysics↗

Response of convectively controlled burning to nonlinear disturbances.

A process where energy release rate is proportional to the square root of the Reynolds number of a burning particle or surface is used to demonstrate some distinctive dynamic properties of a disturbed convectively controlled burning process. The instantaneous energy release rates caused by assumed periodic flow field disturbances are numerically evaluated and examined. Correlation coefficients which express the energy released in-phase with the pressure disturbance are evaluated and analytical solutions for these coefficients are derived. The response of the process is shown to be highly sensitive to the harmonic distortion of the disturbance. For some disturbances the energy released in-phase with the pressure disturbance is an order of magnitude larger than that for linear (sinusoidal) disturbances. The results show that harmonic distortion increases the coupling between the burning process and the flow field. This amplifying effect of harmonic distortion is suppressed when a disturbance becomes steep-fronted; therefore, steepening can act to limit the equilibrium amplitudes in disturbed systems. Amplification is also suppressed when the velocity disturbance is phase shifted with respect to the pressure disturbance or when the process is exposed to a high steady velocity.

Heidmann, M. F.↗

Effects of collisions on the nonlinear particle dynamics in the magnetotail

The effects of collisional processes on the nonlinear particle dynamics in the magnetotail are considered.A simple collision operator is developed to model the effects of pitch-angle and energy scattering. It is found that the phase space partition persists for up to moderate scattering amplitudes in pitch-angle and energy, and that certain distribution function features are robust even in the presence of large amplitude collisions. It is shown that if the collisions are due to short scale length electrostatic fields, excessively large field amplitudes are required to signifiantly alter the phase space structures and the resulting distribution function features.

Holland, Daniel L.↗

Comparative Study of Three High Order Schemes for LES of Temporally Evolving Mixing Layers

Three high order shock-capturing schemes are compared for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) ranging from the quasi-incompressible regime to highly compressible supersonic regime. The considered high order schemes are fifth-order WENO (WENO5), seventh-order WENO (WENO7) and the associated eighth-order central spatial base scheme with the dissipative portion of WENO7 as a nonlinear post-processing filter step (WENO7fi). This high order nonlinear filter method (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The LES results by WENO7fi using the same scheme parameter agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002), whereas results by WENO5 and WENO7 compare poorly with experimental data and DNS computations.

Yee, Helen M. C.↗

Spectral line shapes. Volume 4 - Proceedings of the Eighth International Conference, College of William and Mary, Williamsburg, VA, June 9-13, 1986

The present conference discusses topics in line broadening and shift phenomena in low-to-moderately dense plasmas, resonances in the continuum, hot and dense laser-produced plasmas, collisional broadening and intermolecular potentials, molecular broadening and collision-induced spectroscopy, collisional redistribution and laser-assisted collision processes, coherent and nonlinear spectroscopy, astrophysical line-shapes, and short-wavelength lasers. Attention is given to the widths and shifts of Be ion lines, the Stark broadening of Xe II and Fe I lines, perturbed and unperturbed series in the line-shapes of autoionizing Rydberg series, the line-broadening of hydrogenic ions in dense plasmas, Doppler-free line-shapes, the Franck-Condon principle, the dynamics of laser-assisted reactive collisions, the gas-kinetic effects of light, collision-induced processes in nonlinear optics, and line-shaping phenomena in cool stars.

Exton, Reginald J.↗

On the Application of Hilbert Spectral Analysis for Climate Studies

The Hilbert spectral analysis (Huang et al, 1998, Proceedings of the Royal Society of London, A 454, pp 903-995) consisted of two steps: First, the data has to be reduced into a finite number of Intrinsic Mode Function by the Empirical Mode Decomposition method, then the resulting Intrinsic Mode Functions are converted to time-frequency-energy distribution through Hilbert transform. In this approach, the Empirical Mode Functions served as the basis functions with which the data is expanded. This basis function is adaptive, and the decomposition is nonlinear. Furthermore, as the Hilbert transform is a singular transform, it retains a high degree of local information. The instantaneous frequency is determined by differentiation of the phase function; therefore, there is no restriction of the 'uncertainty principle' for all the time-frequency analysis resulting from a priori basis approach. With the adaptive basis and the instantaneous frequency, the Hilbert Spectral analysis can represent data from nonlinear and nonstationary processes without resorting to the harmonics. Another advantage of using instantaneous frequency is the ability to find out frequency from limited length of data, which is a critical problem in climate studies. As the processes driving the climate changes could be both nonlinear and nonstationary, the Hilbert Spectral Analysis could be of great use in examining the underlying mechanisms. A preliminary study based on the length of day data will be presented as example for the application of the Hilbert Spectral Analysis for climate study.

Huang, Norden E.↗

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. IV. Responses after spectacle-induced adaptation

The horizontal angular vestibuloocular reflex (VOR) evoked by sinusoidal rotations from 0.5 to 15 Hz and acceleration steps up to 3,000 degrees /s(2) to 150 degrees /s was studied in six squirrel monkeys following adaptation with x2.2 magnifying and x0.45 minimizing spectacles. For sinusoidal rotations with peak velocities of 20 degrees /s, there were significant changes in gain at all frequencies; however, the greatest gain changes occurred at the lower frequencies. The frequency- and velocity-dependent gain enhancement seen in normal monkeys was accentuated following adaptation to magnifying spectacles and diminished with adaptation to minimizing spectacles. A differential increase in gain for the steps of acceleration was noted after adaptation to the magnifying spectacles. The gain during the acceleration portion, G(A), of a step of acceleration (3,000 degrees /s(2) to 150 degrees /s) increased from preadaptation values of 1.05 +/- 0.08 to 1.96 +/- 0.16, while the gain during the velocity plateau, G(V), only increased from 0.93 +/- 0.04 to 1.36 +/- 0.08. Polynomial fits to the trajectory of the response during the acceleration step revealed a greater increase in the cubic than the linear term following adaptation with the magnifying lenses. Following adaptation to the minimizing lenses, the value of G(A) decreased to 0.61 +/- 0.08, and the value of G(V) decreased to 0.59 +/- 0.09 for the 3,000 degrees /s(2) steps of acceleration. Polynomial fits to the trajectory of the response during the acceleration step revealed that there was a significantly greater reduction in the cubic term than in the linear term following adaptation with the minimizing lenses. These findings indicate that there is greater modification of the nonlinear as compared with the linear component of the VOR with spectacle-induced adaptation. In addition, the latency to the onset of the adapted response varied with the dynamics of the stimulus. The findings were modeled with a bilateral model of the VOR containing linear and nonlinear pathways that describe the normal behavior and adaptive processes. Adaptation for the linear pathway is described by a transfer function that shows the dependence of adaptation on the frequency of the head movement. The adaptive process for the nonlinear pathway is a gain enhancement element that provides for the accentuated gain with rising head velocity and the increased cubic component of the responses to steps of acceleration. While this model is substantially different from earlier models of VOR adaptation, it accounts for the data in the present experiments and also predicts the findings observed in the earlier studies.

NASA Discipline Neuroscience↗

Use of networked workstations for parallel nonlinear structural dynamic simulations of rotating bladed-disk assemblies

The principal objective of this research is to investigate, develop and demonstrate coarse-grained, parallel-processing strategies for nonlinear dynamic simulations for rotating bladed-disk assemblies. The parallel -processing strategies addressed include numerical algorithms for parallel nonlinear solutions and techniques to effect load balancing among processors. The parallel environment employed is a distributed-memory, coarse-grained one consisting of networked workstations. A parallel explicit time integration method has been implemented for transient nonlinear solutions of rotationg bladed-disk assemblies. Automatic domain partitioning techniques have been investigated for load balancing among processors. Advanced computing environments, data structures and interactive computer graphics all contribute to an integrated parallel finite element analysis system to facilitate more efficient and powerful dynamic simulations.

Hsieh, Shang-Hsien↗

EFQPSK Versus CERN: A Comparative Study

This report presents a comparative study on Enhanced Feher's Quadrature Phase Shift Keying (EFQPSK) and Constrained Envelope Root Nyquist (CERN) techniques. These two techniques have been developed in recent times to provide high spectral and power efficiencies under nonlinear amplifier environment. The purpose of this study is to gain insights into these techniques and to help system planners and designers with an appropriate set of guidelines for using these techniques. The comparative study presented in this report relies on effective simulation models and procedures. Therefore, a significant part of this report is devoted to understanding the mathematical and simulation models of the techniques and their set-up procedures. In particular, mathematical models of EFQPSK and CERN, effects of the sampling rate in discrete time signal representation, and modeling of nonlinear amplifiers and predistorters have been considered in detail. The results of this study show that both EFQPSK and CERN signals provide spectrally efficient communications compared to filtered conventional linear modulation techniques when a nonlinear power amplifier is used. However, there are important differences. The spectral efficiency of CERN signals, with a small amount of input backoff, is significantly better than that of EFQPSK signals if the nonlinear amplifier is an ideal clipper. However, to achieve such spectral efficiencies with a practical nonlinear amplifier, CERN processing requires a predistorter which effectively translates the amplifier's characteristics close to those of an ideal clipper. Thus, the spectral performance of CERN signals strongly depends on the predistorter. EFQPSK signals, on the other hand, do not need such predistorters since their spectra are almost unaffected by the nonlinear amplifier, Ibis report discusses several receiver structures for EFQPSK signals. It is observed that optimal receiver structures can be realized for both coded and uncoded EFQPSK signals with not too much increase in computational complexity. When a nonlinear amplifier is used, the bit error rate (BER) performance of the CERN signals with a matched filter receiver is found to be more than one decibel (dB) worse compared to the bit error performance of EFQPSK signals. Although channel coding is found to provide BER performance improvement for both EFQPSK and CERN signals, the performance of EFQPSK signals remains better than that of CERN. Optimal receiver structures for CERN signals with nonlinear equalization is left as a possible future work. Based on the numerical results, it is concluded that, in nonlinear channels, CERN processing leads towards better bandwidth efficiency with a compromise in power efficiency. Hence for bandwidth efficient communications needs, CERN is a good solution provided effective adaptive predistorters can be realized. On the other hand, EFQPSK signals provide a good power efficient solution with a compromise in band width efficiency.

Borah, Deva K.↗

Bilinear modeling and nonlinear estimation

New methods are illustrated for online nonlinear estimation applied to the lateral deflection of an elastic beam on board measurements of angular rates and angular accelerations. The development of the filter equations, together with practical issues of their numerical solution as developed from global linearization by nonlinear output injection are contrasted with the usual method of the extended Kalman filter (EKF). It is shown how nonlinear estimation due to gyroscopic coupling can be implemented as an adaptive covariance filter using off-the-shelf Kalman filter algorithms. The effect of the global linearization by nonlinear output injection is to introduce a change of coordinates in which only the process noise covariance is to be updated in online implementation. This is in contrast to the computational approach which arises in EKF methods arising by local linearization with respect to the current conditional mean. Processing refinements for nonlinear estimation based on optimal, nonlinear interpolation between observations are also highlighted. In these methods the extrapolation of the process dynamics between measurement updates is obtained by replacing a transition matrix with an operator spline that is optimized off-line from responses to selected test inputs.

Dwyer, Thomas A. W., III↗