Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector field”

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,081 records · Page 60

Development of a Sunspot Tracking System

Large solar flares produce a significant amount of energetic particles which pose a hazard for human activity in space. In the hope of understanding flare mechanisms and thus better predicting solar flares, NASA's Marshall Space Flight Center (MSFC) developed an experimental vector magnetograph (EXVM) polarimeter to measure the Sun's magnetic field. The EXVM will be used to perform ground-based solar observations and will provide a proof of concept for the design of a similar instrument for the Japanese Solar-B space mission. The EXVM typically operates for a period of several minutes. During this time there is image motion due to atmospheric fluctuation and telescope wind loading. To optimize the EXVM performance an image motion compensation device (sunspot tracker) is needed. The sunspot tracker consists of two parts, an image motion determination system and an image deflection system. For image motion determination a CCD or CID camera is used to digitize an image, than an algorithm is applied to determine the motion. This motion or error signal is sent to the image deflection system which moves the image back to its original location. Both of these systems are under development. Two algorithms are available for sunspot tracking which require the use of only one row and one column of image data. To implement these algorithms, two identical independent systems are being developed, one system for each axis of motion. Two CID cameras have been purchased; the data from each camera will be used to determine image motion for each direction. The error signal generated by the tracking algorithm will be sent to an image deflection system consisting of an actuator and a mirror constrained to move about one axis. Magnetostrictive actuators were chosen to move the mirror over piezoelectrics due to their larger driving force and larger range of motion. The actuator and mirror mounts are currently under development.

Taylor, Jaime R.↗

A Determination of the Precession and Rotation-Rate of Mars from Tracking and Altimetry Data

The doppler tracking of the MGS spacecraft and the laser altimetry acquired by the MOLA instrument on MGS have been used to measure the rotation of Mars. By observing the repeating pattern of gravity anomalies and the topography the direction of the rotation axis in both inertial space and body fixed coordinates, and the rate of rotation of the planet on its axis have been derived. The tracking of MGS in orbit about Mars by stations on Earth permit the orbit of MGS to be obtained in an inertial coordinate system enabling the periodic perturbation of the MGS orbit by Mars' gravity field to be used as a sensor to detect the vector rotation of Mars under the orbit of MGS. The altimetry data acquired by MOLA on MGS permit the topography to help define the orbit of MGS about Mars in a body fixed coordinate system thus enabling the pole of rotation of the planet to be determined in a coordinate system defined by the altimetry. Together they enable the precession of the rotation axis to be observed from which the moment of inertia can be derived, changes in length of day to be measured, and a method to check for possible Mars polar motion. Tracking data on the Mariner 9, Viking 1 and 2, from the 1970's have been analyzed with MGS data obtained in the mapping orbit in 1999 to provide a new estimate of the moment of inertia.

Smith, David E.↗

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

compressor↗

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

unsteady↗

An Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. II: Approach Based on the Dyson and Bethe-Salpeter Equations

In this paper, the vector radiative transfer equation is derived by means of the vector integral Foldy equations describing the electromagnetic scattering by a group of particles. By assuming that in a discrete random medium the positions of the particles are statistically independent and by applying the Twersky approximation to the order-of-scattering expansion of the total field, we derive the Dyson equation for the coherent field and the ladder approximated Bethe–Salpeter equation for the dyadic correlation function. Then, under the far-field assumption for sparsely distributed particles, the Dyson equation is reduced to the Foldy integral equation for the coherent field, while the iterated solution of the Bethe–Salpeter equation ultimately yields the vector radiative transfer equation.

Electromagnetic scattering↗

Voyager planetary radio astronomy studies

Analysis of nonthermal radio emission data obtained by the Planetary Radio Astronomy (PRA) spectrometers on the Voyager 1 and 2 spacecraft was performed. This PRA data provided unique insights into the radio emission characteristics of the outer planets because of PRA's unique spectral response below the terrestrial ionospheric plasma frequency and its unprecedented proximity to the source. Of those results which were documented or published, this final report surveys only the highlights and cites references for more complete discussions. Unpublished results for Uranus, Neptune, and theoretical Ionian current distributions are presented at greater length. The most important conclusion to be drawn from these observations is that banded spectral emission is common to the radio emission below 1-2 MHz observed from all four Jovian planets. In every case multiple spectral features evolve on time scales of seconds to minutes. To the extent these features drift in frequency, they appear never to cross one another. The Neptunian spectral features appear to drift little or not at all, their evolution consisting principally of waxing and waning. Since other evidence strongly suggests that most or all of this radio emission is occurring near the local magnetospheric electron cyclotron frequency, this implies that this emission preferentially occurs at certain continually changing planetary radii. It remains unknown why certain radii might be favored, unless radial electric field components or other means serve to differentiate radially the magnetospheric plasma density, particle energy vectors, or particle coherence. Calculation of the spatial distribution and intensity of the Io-generated magnetospheric currents are also presented; these currents may be limited principally by wave impedance and local field strengths.

Staelin, David H.↗

Study of waves in the earth's bow shock.

Results of an examination of the perturbation vectors of waves upstream and downstream from the region of maximum compression in the bow shock on Ogo 5 under particularly steady solar-wind conditions. The polarization of the upstream waves was right-hand circular, and that of the downstream waves left-hand elliptical in the spacecraft frame. By observing that the polarization of the waves remained unchanged as the shock motion swept the wave structure back and forth across the satellite three times in eight minutes, it was found that the waves were not stationary in the shock frame. A study of the methods of determining the shock normal indicates that the normal estimated from a shock model should be superior to the normal based on magnetic coplanarity. The propagation vectors of the waves examined did not coincide with the shock-model normal, the average magnetic field, or the plasma-flow velocity. However, the major axis of the polarization ellipse of the downstream wave was nearly parallel to the upstream propagation vector.

Holzer, R. E.↗

The Geopotential Research Mission - Mapping the near earth gravity and magnetic fields

The Geopotential Research Mission (GRM), NASA's low-level satellite system designed to measure the gravity and magnetic fields of the earth, and its objectives are described. The GRM will consist of two, Shuttle launched, satellite systems (300 km apart) that will operate simultaneously at a 160 km circular-polar orbit for six months. Current mission goals include mapping the global geoid to 10 cm, measuring gravity-field anomalies to 2 mgal with a spatial resolution of 100 km, detecting crustal magnetic anomalies of 100 km wavelength with 1 nT accuracy, measuring the vectors components to + or - 5 arc sec and 5 nT, and computing the main dipole or core field to 5 nT with a 2 nT/year secular variation detection. Resource analysis and exploration geology are additional applications considered.

Taylor, P. T.↗

Expanded Equations for Torque and Force on a Cylindrical Permanent Magnet Core in a Large-Gap Magnetic Suspension System

The expanded equations for torque and force on a cylindrical permanent magnet core in a large-gap magnetic suspension system are presented. The core is assumed to be uniformly magnetized, and equations are developed for two orientations of the magnetization vector. One orientation is parallel to the axis of symmetry, and the other is perpendicular to this axis. Fields and gradients produced by suspension system electromagnets are assumed to be calculated at a point in inertial space which coincides with the origin of the core axis system in its initial alignment. Fields at a given point in the core are defined by expanding the fields produced at the origin as a Taylor series. The assumption is made that the fields can be adequately defined by expansion up to second-order terms. Examination of the expanded equations for the case where the magnetization vector is perpendicular to the axis of symmetry reveals that some of the second-order gradient terms provide a method of generating torque about the axis of magnetization and therefore provide the ability to produce six-degree-of-freedom control.

Groom, Nelson J.↗

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗

Static performance investigation of a skewed-throat multiaxis thrust-vectoring nozzle concept

The static performance of a jet exhaust nozzle which achieves multiaxis thrust vectoring by physically skewing the geometric throat has been characterized in the static test facility of the 16-Foot Transonic Tunnel at NASA Langley Research Center. The nozzle has an asymmetric internal geometry defined by four surfaces: a convergent-divergent upper surface with its ridge perpendicular to the nozzle centerline, a convergent-divergent lower surface with its ridge skewed relative to the nozzle centerline, an outwardly deflected sidewall, and a straight sidewall. The primary goal of the concept is to provide efficient yaw thrust vectoring by forcing the sonic plane (nozzle throat) to form at a yaw angle defined by the skewed ridge of the lower surface contour. A secondary goal is to provide multiaxis thrust vectoring by combining the skewed-throat yaw-vectoring concept with upper and lower pitch flap deflections. The geometric parameters varied in this investigation included lower surface ridge skew angle, nozzle expansion ratio (divergence angle), aspect ratio, pitch flap deflection angle, and sidewall deflection angle. Nozzle pressure ratio was varied from 2 to a high of 11.5 for some configurations. The results of the investigation indicate that efficient, substantial multiaxis thrust vectoring was achieved by the skewed-throat nozzle concept. However, certain control surface deflections destabilized the internal flow field, which resulted in substantial shifts in the position and orientation of the sonic plane and had an adverse effect on thrust-vectoring and weight flow characteristics. By increasing the expansion ratio, the location of the sonic plane was stabilized. The asymmetric design resulted in interdependent pitch and yaw thrust vectoring as well as nonzero thrust-vector angles with undeflected control surfaces. By skewing the ridges of both the upper and lower surface contours, the interdependency between pitch and yaw thrust vectoring may be eliminated and the location of the sonic plane may be further stabilized.

Wing, David J.↗

Wave Normal and Poynting Vector Calculations using the Cassini Radio and Plasma Wave Instrument

Wave normal and Poynting vector measurements from the Cassini radio and plasma wave instrument (RPWS) are used to examine the propagation characteristics of various plasma waves during the Earth flyby on August 18, 1999. Using the five-channel waveform receiver (WFR), the wave normal vector is determined using the Means method for a lightning-induced whistler, equatorial chorus, and a series of low-frequency emissions observed while Cassini was in the magnetosheath. The Poynting vector for these emissions is also calculated from the five components measured by the WFR. The propagation characteristics of the lightning-induced whistler were found to be consistent with the whistler wave mode of propagation, with propagation antiparallel to the magnetic field (southward) at Cassini. The sferic associated with this whistler was observed by both Cassini and the Stanford VLF group at the Palmer Station in Antarctica. Analysis of the arrival direction of the sferic at the Palmer Station suggests that the lightning stroke is in the same sector as Cassini. Chorus was observed very close (within a few degrees) to the magnetic equator during the flyby. The chorus was found to propagate primarily away from the magnetic equator and was observed to change direction as Cassini crossed the magnetic equator. This suggests that the source region of the chorus is very near the magnetic equator. The low-frequency emission in the magnetosheath has many of the characteristics of lion roars. The average value of the angle between the wave normal vector and the local magnetic field was found to be 16 degrees, and the emissions ranged in frequency from 0. 19 to 0.75 f(sub ce), where f(sub ce) is the electron cyclotron frequency. The wave normal vectors of these waves were primarily in one direction for each individual burst (either parallel or antiparallel to the local field) but varied in direction throughout the magnetosheath. This suggests that the sources of the emissions are far from the spacecraft and that there are multiple source regions.

Hospodarsky, G. B.↗

Correlation of the CME Productivity of Solar Active Regions with Measures of their Global Nonpotentiality from Vector Magnetograms: Baseline Results

From conventional magnetograms and chromospheric and coronal images, it is known qualitatively that the fastest coronal mass ejections (CMEs) are magnetic explosions from sunspot active regions in which the magnetic field is globally strongly sheared and twisted from its minimum-energy potential configuration. In this paper, we present measurements from active-region vector magnetograms that begin to quantify the dependence of the CME productivity of an active region on the global nonpotentiality of its magnetic field. From each of 17 magnetograms of 12 bipolar active regions, we obtain a measure of the size of the active region (the magnetic flux content, phi) and three different measures of the global nonpotentiality (L(sub SS), the length of strong-shear, strong-field main neutral line; I(sub N), the net electric current arching from one polarity to the other; and alpha = muI(subN/phi), a flux-normalized measure of the field twist).

Falconer, David A.↗

Analytic Solution to the Problem of Aircraft Electric Field Mill Calibration

It is by no means a simple task to retrieve storm electric fields from an aircraft instrumented with electric field mill sensors. The presence of the aircraft distorts the ambient field in a complicated way. Before retrievals of the storm field can be made, the field mill measurement system must be "calibrated". In other words, a relationship between impressed (i.e., ambient) electric field and mill output must be established. If this relationship can be determined, it is mathematically inverted so that ambient field can be inferred from the mill outputs. Previous studies have primarily focused on linear theories where the relationship between ambient field and mill output is described by a "calibration matrix" M. Each element of the matrix describes how a particular component of the ambient field is enhanced by the aircraft. For example the product M(sub ix), E(sub x), is the contribution of the E(sub x) field to the i(th) mill output. Similarly, net aircraft charge (described by a "charge field component" E(sub q)) contributes an amount M(sub iq)E(sub q) to the output of the i(th) sensor. The central difficulty in obtaining M stems from the fact that the impressed field (E(sub x), E(sub y), E(sub z), E(sub q) is not known but is instead estimated. Typically, the aircraft is flown through a series of roll and pitch maneuvers in fair weather, and the values of the fair weather field and aircraft charge are estimated at each point along the aircraft trajectory. These initial estimates are often highly inadequate, but several investigators have improved the estimates by implementing various (ad hoc) iterative methods. Unfortunately, none of the iterative methods guarantee absolute convergence to correct values (i.e., absolute convergence to correct values has not been rigorously proven). In this work, the mathematical problem is solved directly by analytic means. For m mills installed on an arbitrary aircraft, it is shown that it is possible to solve for a single 2m-vector that provides all other needed variables (i.e., the unknown fair weather field, the unknown aircraft charge, and the unknown matrix M). Numerical tests of the solution, effects of measurement errors, and studies of solution non-uniqueness are ongoing as of this writing.

Koshak, William↗

A plane-polar approach for far-field construction from near-field measurements

The planar configuration with a probe scanning a polar geometry is discussed with reference to its usefulness in the determination of a far field from near-field measurements. The accuracy of the method is verified numerically, using the concept of probe compensation as a vector deconvolution. Advantages of the Jacobi-Bessel series over the fast Fourier transforms for the plane-polar geometry are demonstrated. Finally, the far-field pattern of the Viking high gain antenna is constructed from the plane-polar near-field measured data and compared with the previously measured far-field pattern.

Rahmat-Samii, Y.↗

A Prospective Method for Predicting Coronal Mass Ejections From Vector Magnetograms

Two quantitative measures of the global nonpotentiality of predominantly bipolar active regions are explored for their usefulness as indicators of whether an active region is likely to produce coronal mass ejections (CMEs). The two measures, length of strong-field strong-shear main neutral line (Lss) and global net current (IN), are obtained directly from vector magnetograms of the regions. From measurements of Lss and IN in four active regions for which the global nonpotentiality and CME productivity are assessed independently of Lss and IN, it appears that Lss and IN are each a useful indicator of an active region's likely CME productivity. The results also suggest that Lss and IN may be more reliable predictors of the CME productivity of active regions than is the presence or absence of sigmoidal structure in coronal X-ray images of the regions and that there might be threshold values of Lss and IN above which CME production is likely and below which CME production is unlikely. While these results need to be certified by an expansion of this pilot study to more active regions, they show the promise of Lss and IN as quantitative measures of global nonpotentiality for CME prediction.

Falconer, D. A.↗

Cutaneous microvascular flow in the foot during simulated variable gravities

Our objective was to understand how weight bearing with varying gravitational fields affects blood perfusion in the sole of the foot. Human subjects underwent whole body tilting at four angles: upright [1 gravitational vector from head to foot (Gz)], 22 degrees (0.38 Gz), 10 degrees (0.17 Gz), and supine (0 Gz), simulating the gravitational fields of Earth, Mars, Moon, and microgravity, respectively. Cutaneous capillary blood flow was monitored on the plantar surface of the heel by laser Doppler flowmetry while weight-bearing load was measured. At each tilt angle, subjects increased weight bearing on one foot in graded load increments of 1 kg beginning with zero. The weight bearing at which null flow first occurred was determined as the closing load. Subsequently, the weight bearing was reduced in reverse steps until blood flow returned (opening load). Mean closing loads for simulated Earth gravity, Mars gravity, Moon gravity, and microgravity were 9.1, 4.6, 4.4, and 3.6 kg, respectively. Mean opening loads were 7.9, 4.1, 3.5, and 3.1 kg, respectively. Mean arterial pressures in the foot (MAP(foot)) calculated for each simulated gravitational field were 192, 127, 106, and 87 mmHg, respectively. Closing load and opening load were significantly correlated with MAP(foot) (r =0.70, 0.72, respectively) and were significantly different (P < 0.001) from each other. The data suggest that decreased local arterial pressure in the foot lowers tolerance to external compression. Consequently, the human foot sole may be more prone to cutaneous ischemia during load bearing in microgravity than on Earth.

NASA Discipline Cardiopulmonary↗

Doppler lidar sampling strategies and accuracies: Regional scale

It has been proposed that a Doppler lidar be placed in a polar orbit and scanned to provide estimates of lower tropospheric winds twice per day and with a spatial resolution of 300 km. Initial feasibility studies conducted primarily by NOAA and NASA presented an optimistic outlook for a space based lidar. The technology appeared within reach and initial computer simulations suggested that acceptable accuracies could be obtained. Those early studies exposed, however, several potential problem areas which included: (1) the algorithms for computing the wind vectors did not perform well when there were coherent gradients in the wind fields; and (2) the lifetime and power requirements of the lidar put severe restrictions on the pulse repetition frequency (PRF). These two basic problems are currently being addressed by a Doppler lidar simulation study focussed upon three primary objectives: (1) to develop optimum scan parameters and shot patterns for a satellite-based Doppler lidar; (2) to develop robust algorithms for computing wind vectors from lidar returns; and (3) to evaluate the impact of coherent mesoscale structures (wind gradients, clouds, aerosols) on up-scale wind estimates. An overview is provided of the simulation efforts with particular emphasis upon rationale and methodology. Since this research is currently underway, any results shown are meant only as evidence of progress.

Emmitt, G. D.↗