Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “metric 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 325 records · Page 18

High-Level Performance Modeling of SAR Systems

SAUSAGE (Still Another Utility for SAR Analysis that s General and Extensible) is a computer program for modeling (see figure) the performance of synthetic- aperture radar (SAR) or interferometric synthetic-aperture radar (InSAR or IFSAR) systems. The user is assumed to be familiar with the basic principles of SAR imaging and interferometry. Given design parameters (e.g., altitude, power, and bandwidth) that characterize a radar system, the software predicts various performance metrics (e.g., signal-to-noise ratio and resolution). SAUSAGE is intended to be a general software tool for quick, high-level evaluation of radar designs; it is not meant to capture all the subtleties, nuances, and particulars of specific systems. SAUSAGE was written to facilitate the exploration of engineering tradeoffs within the multidimensional space of design parameters. Typically, this space is examined through an iterative process of adjusting the values of the design parameters and examining the effects of the adjustments on the overall performance of the system at each iteration. The software is designed to be modular and extensible to enable consideration of a variety of operating modes and antenna beam patterns, including, for example, strip-map and spotlight SAR acquisitions, polarimetry, burst modes, and squinted geometries.

Chen, Curtis↗

Calabi-Yau CFTs and random matrices

Using numerical methods for finding Ricci-flat metrics, we explore the spectrum of local operators in two-dimensional conformal field theories defined by sigma models on Calabi-Yau targets at large volume. Focusing on the examples of K3 and the quintic, we show that the spectrum, averaged over a region in complex structure moduli space, possesses the same statistical properties as the Gaussian orthogonal ensemble of random matrix theory.

superstring vacua↗

Turbulent Navier-Stokes Flow Analysis of an Advanced Semispan Diamond-Wing Model in Tunnel and Free Air at High-Lift Conditions

Turbulent Navier-Stokes computational results are presented for an advanced diamond wing semispan model at low-speed, high-lift conditions. The numerical results are obtained in support of a wind-tunnel test that was conducted in the National Transonic Facility at the NASA Langley Research Center. The model incorporated a generic fuselage and was mounted on the tunnel sidewall using a constant-width non-metric standoff. The computations were performed at to a nominal approach and landing flow conditions.The computed high-lift flow characteristics for the model in both the tunnel and in free-air environment are presented. The computed wing pressure distributions agreed well with the measured data and they both indicated a small effect due to the tunnel wall interference effects. However, the wall interference effects were found to be relatively more pronounced in the measured and the computed lift, drag and pitching moment. Although the magnitudes of the computed forces and moment were slightly off compared to the measured data, the increments due the wall interference effects were predicted reasonably well. The numerical results are also presented on the combined effects of the tunnel sidewall boundary layer and the standoff geometry on the fuselage forebody pressure distributions and the resulting impact on the configuration longitudinal aerodynamic characteristics.

Ghaffari, Farhad↗

From Fundamental Interfacial Reaction Kinetics to Macroscopic Current–Voltage Characteristics: Case Study of Solid Acid Fuel Cell Limitations and Possibilities

The unique properties of solid acid electrolytes, in particular CsH 2 PO 4 , are in many ways ideal for fuel cell operation. However, the technology is constrained by high cathode overpotentials. Here a simplified cathode geometry is employed to obtain the fundamental electrochemical parameters (exchange current density and charge transfer coefficient) describing the oxygen reduction reaction (ORR) at the CsH 2 PO 4 -Pt-gas interface. The parameters are incorporated into a 1D model of the voltage–current characteristics of realistic SAFC cathodes, which reproduced the measured polarization behavior of such cathodes without recourse to fitting adjustable parameters. Following this validation, the model is utilized to evaluate the impact of changes to cathode properties, microstructure, and operating conditions. Of these, the charge transfer coefficient, measured to have a value of ≈0.6 for ORR on Pt in the SAFC cathode environment, is found to have the greatest impact on power output. Nevertheless, even without material modifications, a combination of microstructural and operational modifications are identified with projected performance metrics meeting Department of Energy targets (0.8 V at 300 mA cm –2 , and peak power density of 1 W cm –2 ), albeit at high Pt loadings. However, the analysis indicates that truly meaningful advances will likely necessitate the discovery of alternative ORR catalysts.

36 MATERIALS SCIENCE↗

Assessment of Radiated Fan Noise Prediction Capabilities Using Static Engine Test Data

This paper describes further assessment of the CDUCT-LaRC code via comparison with static engine test data. In an effort to improve confidence in the use of CDUCT-LaRC for liner optimization studies addressing realistic three-dimensional geometries, inlet radiated fan noise predictions were performed at 54% and 87% engine speed settings. Predictions were then compared with far-field measurements to assess the approach and implementation. The particular configurations were chosen to exercise the three-dimensional capability of CDUCT-LaRC and it s applicability to realistic configurations and conditions. At the 54% engine speed setting, the predictions capture the general directivity and acoustic treatment effects quite well. Comparisons of the predicted and measured directivity at the 87% power setting were more problematic. This was likely due in part to the difficulties in source specification and possibly the nonlinear nature of buzz-saw tones at this engine operating condition. Overall, the approach captured the basic trends and provided a conservative estimate of liner effects from which relative performance metrics could be inferred.

Nark, Douglas M.↗

Dataset describing two reference models for full-spectral lighting and daylight simulations together with implementations for two software systems

A dataset of two spectral lighting simulation reference models - one office and one factory hall - is presented. It aims to demonstrate and support full-spectral daylight and electric lighting simulations and facilitate evaluation of non-visual effects of light. The dataset includes Rhino CAD geometry, comprehensive spectral material and light source data and window system BSDF data. Example implementations in the two software tools, Radiance and OWL, enable reproducible workflows and support adoption in other software. The dataset is openly available on Zenodo. The office model reproduces Room 518 at the University of Innsbruck, including a west-facing façade and interior furnishings. The factory hall model follows the proposed geometry in the European standard 15193 for building energy performance. Interior reflectances in the office were measured in-situ using a handheld spectrometer. Exterior spectra and factory hall materials matching specified reflectances were obtained from an online spectral materials database. Glazing transmittance was derived from IGDB data using LBNL Optics/WINDOW. BSDFs for venetian blinds at various tilt angles, and for a diffusing pane adapted from the Complex Glazing Database, were generated in WINDOW. Luminaires in both models are specified with photometric files (Eulumdat/IES) and lamp spectra (Fluorescent 840, 4000 K LED). The provided example implementations (Radiance, OWL) include prepared input data and scripts to run first spectral simulations; example results are also included. The dataset is prepared to support reuse by researchers, designers and software developers for method validation, software engineering and comparison, and development of spectral metrics and controls.

Geisler-Moroder, David↗

Post-puff SOL broadening on MAST-U under high-recycling conditions: evidence consistent with cross-field transport changes

Transient broadening of the scrape-off layer (SOL) density profile can modify main-chamber first-wall particle fluxes and divertor loading, yet its control parameters remain debated between divertor-regime transitions, neutral dynamics and changes in cross-field transport. We investigate fueling-driven SOL density-profile evolution and post-fueling relaxation on MAST-Upgrade (MAST-U) in ohmic L-mode, using two otherwise similar double-null Conventional Divertor discharges (I p = 450 kA, B T = 0.33 T) with identical 50 ms low-field side gas puffs; the only intentional difference is the puff start time. Upstream Thomson scattering shows that both discharges develop a transient far-SOL density profile modification, expressed as an increased SOL-width metric λ n e and a far-SOL enhancement consistent with a shoulder-like signature. In the earlier-puff case, the SOL-width metric remains elevated after puff termination and the post-puff decay time scales are systematically longer across the analysed radii. Outer-target Langmuir probes indicate high-recycling conditions during the analysed window (few-eV T e with radially peaked j sat,∥ and q ∥ profiles, without signatures of deep detachment). The outer-divertor collisionality proxy Λ div is elevated in both cases and does not discriminate between the different post-puff persistence. The target-integrated ion flux proxy ∫ J sat dA evolves nearly identically in both discharges when aligned to the puff start. Taken together, these observations suggest that the late/post-puff upstream SOL evolution is not set by parallel exhaust to the outer targetalone and is more consistent with upstream cross-field redistribution, with a possible role for plasma–neutral coupling and fueling geometry.

MAST-U↗

Search for Effects of an Electrostatic Potential on Clocks in the Frame of Reference of a Charged Particle

Results of experiments to confirm a theory that links classical electromagnetism with the geometry of spacetime are described. The theory, based on the introduction of a Torsion tensor into Einstein s equations and following the approach of Schroedinger, predicts effects on clocks attached to charged particles, subject to intense electric fields, analogous to the effects on clocks in a gravitational field. We show that in order to interpret this theory, one must re-interpret all clock changes, both gravitational and electromagnetic, as arising from changes in potential energy and not merely potential. The clock is provided naturally by proton spins in hydrogen atoms subject to Nuclear Magnetic Resonance trials. No frequency change of clocks was observed to a resolution of 6310(exp -9). A new "Clock Principle" was postulated to explain the null result. There are two possible implications of the experiments: (a) The Clock Principle is invalid and, in fact, no metric theory incorporating electromagnetism is possible; (b) The Clock Principle is valid and it follows that a negative rest mass cannot exist.

Ringermacher, Harry I.↗

An Efficient, Multi-Layered Crown Delineation Algorithm for Mapping Individual Tree Structure Across Multiple Ecosystems

Deriving individual tree information from discrete return, small footprint LiDAR data may improve forest above ground biomass estimates, and provide tree-level information that is important in many ecological studies. Several crown delineation algorithms have been developed to extract individual tree information from LiDAR point clouds or rasterized canopy height models (CHM), but many of these algorithms have difficulty discriminating between overlapping crowns, and also may fail to detect understory trees. Our approach uses a watershed based delineation of a CHM, which is subsequently refined using the LiDAR point cloud. Individual tree detection was validated with stem mapped field data from the Smithsonian Environmental Research Center (SERC), Maryland, and on a plot and stand level through comparisons of stem density and basal area to delineated metrics at both SERC and a study area in the Sierra Nevada, California. For individual tree detection, the algorithm correctly identified 70% of dominant trees, 58% of co-dominant trees, 35% of intermediate trees and 21% of suppressed trees at SERC. The algorithm had difficulty distinguishing between crowns of small, dense understory trees of approximately the same height. Delineated crown volume alone explained 53% and 84% of the variability in basal area at the SERC and Sierra Nevada sites, respectively. The algorithm produced crown area distributions comparable to diameter at breast height (DBH) size class distributions observed in the field in both study sites. The algorithm detected understory crowns better in the conifer-dominated Sierra Nevada site than in the closed-canopy deciduous site in Maryland. The ability for the algorithm to reproduce both accurate tree size distributions and individual crown geometries in two dissimilar and complex forests suggests great promise for applicability to a wide range of forest systems.

LiDAR↗

Extreme Environment Hot Fire Durability of Post-Processed Additively Manufactured GRCop-Alloy Combustion Chambers in LOX/Hydrogen and LOX/Methane

Extreme environment survivability of metal additive manufactured (AM) GRCop-alloy thrust chambers has been demonstrated in different bi-propellants at near stoichiometric and even oxygen rich combustion. GRCop-alloy chambers tested at NASA Marshall Space Flight Center (MSFC) have accumulated over 26,000 seconds of hot fire duration and over 500 starts. These chambers are produced using an AM process called laser powder bed fusion (L-PBF). A major feature of this process is high wall surface roughness which can be customized in post processing to leverage various performance advantages. Post-processing can include heat treatment, final machining, polishing, and welding and is key to hardware survivability in extreme environments. Surface finish enhancement techniques were applied to the hot wall and coolant channels to reduce the overall total heat load to the chamber walls. Performance optimization of various thrust class TCA’s is a strategic technology goal of NASA MSFC. Three different chamber geometries using cryogenic methane and de-ionized water as coolants were hot fire tested to obtain their life cycle, pressure drop, and heat load performances. Several 1.2K lbf LOX/H2 chambers, 1K lbf LOX/CH4 chambers, and 7K lbf LOX/CH4 chambers were tested. All post-processed configurations performed remarkably well when subjected to extreme hot fire test conditions. Streaking and blanching due to localized oxygen rich conditions were observed on some test articles in their as-built surface finish state. However, this is actually a function of the injector mixing and only serves to further establish the durability of L-PBF GRCop-alloy thrust chambers. Several different polishing techniques were applied to the hot wall and integrated coolant channels prior to hot fire testing and their performances assessed. Overall, AM produced GRCop chambers are extremely reliable, durable, and customizable to the desired performance metrics.

grcop↗

Correlation Calculations for the Russian Pu Metal Fast Experiments

Nuclear criticality experiments are often conducted in campaigns with multiple variations. These experiments reuse the same basic components, like the fuel, moderator, or positioning machines. The components have uncertainties in their geometry and composition that propagate to models of the experiments. Shared components create shared uncertainty between the $k_{eff}$ of benchmarks. The shared uncertainty is commonly quantified with a covariance, or correlation coefficient. These covariances can impact criticality safety and nuclear data validation applications. While benchmark evaluations tabulate an experiment’s uncertainty, they often lack a detailed calculation of correlations between experiments. Even some very commonly used benchmarks, like the Russian Pu Metal Fast (PMF) experiments, have missing correlations. This paper presents our approach to calculate the correlations for five of the Russian PMF experiments. The experiments share hemispherical Pu shells that induce a correlation between modeled $k_{eff}$ values. We estimated the correlations with simplified and detailed models of the experiments through linear-perturbation theory. The correlations between the experiments vary significantly between the detailed vs. simplified models. We also investigate how the correlations affect validation metrics of the experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Visual Object Recognition and Tracking of Tools

A method has been created to automatically build an algorithm off-line, using computer-aided design (CAD) models, and to apply this at runtime. The object type is discriminated, and the position and orientation are identified. This system can work with a single image and can provide improved performance using multiple images provided from videos. The spatial processing unit uses three stages: (1) segmentation; (2) initial type, pose, and geometry (ITPG) estimation; and (3) refined type, pose, and geometry (RTPG) calculation. The image segmentation module files all the tools in an image and isolates them from the background. For this, the system uses edge-detection and thresholding to find the pixels that are part of a tool. After the pixels are identified, nearby pixels are grouped into blobs. These blobs represent the potential tools in the image and are the product of the segmentation algorithm. The second module uses matched filtering (or template matching). This approach is used for condensing synthetic images using an image subspace that captures key information. Three degrees of orientation, three degrees of position, and any number of degrees of freedom in geometry change are included. To do this, a template-matching framework is applied. This framework uses an off-line system for calculating template images, measurement images, and the measurements of the template images. These results are used online to match segmented tools against the templates. The final module is the RTPG processor. Its role is to find the exact states of the tools given initial conditions provided by the ITPG module. The requirement that the initial conditions exist allows this module to make use of a local search (whereas the ITPG module had global scope). To perform the local search, 3D model matching is used, where a synthetic image of the object is created and compared to the sensed data. The availability of low-cost PC graphics hardware allows rapid creation of synthetic images. In this approach, a function of orientation, distance, and articulation is defined as a metric on the difference between the captured image and a synthetic image with an object in the given orientation, distance, and articulation. The synthetic image is created using a model that is looked up in an object-model database. A composable software architecture is used for implementation. Video is first preprocessed to remove sensor anomalies (like dead pixels), and then is processed sequentially by a prioritized list of tracker-identifiers.

English, James↗

Assessment of Noise Reduction Concepts for Leading-Edge Slat Noise

The leading-edge slat of a high-lift airfoil can be a significant noise contributor during aircraft landing. This paper summarizes the effects of several passive noise reduction devices on the 30P30N high-lift airfoil. Experiments are conducted on a two-dimensional multi-element high-lift airfoil with leading-edge slat extensions, gap filler, and cove filler in an anechoic wind tunnel to evaluate the effect of passive flow control on the acoustics generated by the unsteady flow field. Slat geometry modifications associated with the treatments alter the flow field in the region that dominates the generation of the acoustic field. Three angles of attack (a(k) = 8°, 10°, and 15.5°) and three different Reynolds numbers (Re(c) = 1.2e6, 1.5e6, and 1.71e6) are selected as the test conditions. Steady surface pressure measurements are conducted to assess the effect of the treatments on the lift and drag. Unsteady surface pressure measurements along with the far-field acoustic array measurements are performed to evaluate the changes in near- and far-field pressure fluctuations, respectively. Delay and Sum (DAS) beamforming method is applied to locate the noise sources on the model and provide integrated spectra. Implementation difficulties with the gap filler led to structural integration deficiencies that prevented a fair assessment of this technology. Among the other passive devices, the cove filler s the most effective noise reduction, along with a negligible change in the aerodynamic metrics.

Zhang, Yang↗

Nuclear Data Assessment for Advanced Reactors

Advanced reactor concepts being developed throughout the industry are significantly different from light-water reactor (LWR) designs with respect to geometry, materials, and operating conditions, and consequently, with respect to their reactor physics behavior. Given the limited operating experience with non-LWRs, the accurate simulation of reactor physics and the quantification of associated uncertainties are important for ensuring that the nuclear design for advanced reactor concepts include appropriate margins. Nuclear data are a major source of input uncertainties in reactor physics analysis. As part of a project sponsored by the US Nuclear Regulatory Commission at Oak Ridge National Laboratory (ORNL), key nuclear data relevant to reactor safety analysis in selected advanced reactor technologies were identified, and their impacts on important key figures of merit were assessed based on (1) a review of available advanced reactor specifications, (2) analysis of previous studies performed at ORNL and other research institutions, and (3) sensitivity and uncertainty analyses performed for six selected benchmarks—three experimental and three computational—to quantify the impacts of the identified key nuclear data on several key metrics. This report summarizes the key nuclear data—nominal data and nuclear data uncertainties— considering the most important nuclear reactions in the fuel and in various materials for the moderator, coolant, and structure of the considered advanced reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Probabilistic Nanomagnetic Memories for Uncertain and Robust Machine Learning

This project evaluated the use of emerging spintronic memory devices for robust and efficient variational inference schemes. Variational inference (VI) schemes, which constrain the distribution for each weight to be a Gaussian distribution with a mean and standard deviation, are a tractable method for calculating posterior distributions of weights in a Bayesian neural network such that this neural network can also be trained using the powerful backpropagation algorithm. Our project focuses on domain-wall magnetic tunnel junctions (DW-MTJs), a powerful multi-functional spintronic synapse design that can achieve low power switching while also opening the pathway towards repeatable, analog operation using fabricated notches. Our initial efforts to employ DW-MTJs as an all-in-one stochastic synapse with both a mean and standard deviation didn’t end up meeting the quality metrics for hardware-friendly VI. In the future, new device stacks and methods for expressive anisotropy modification may make this idea still possible. However, as a fall back that immediately satisfies our requirements, we invented and detailed how the combination of a DW-MTJ synapse encoding the mean and a probabilistic Bayes-MTJ device, programmed via a ferroelectric or ionically modifiable layer, can robustly and expressively implement VI. This design includes a physics-informed small circuit model, that was scaled up to perform and demonstrate rigorous uncertainty quantification applications, up to and including small convolutional networks on a grayscale image classification task, and larger (Residual) networks implementing multi-channel image classification. Lastly, as these results and ideas all depend upon the idea of an inference application where weights (spintronic memory states) remain non-volatile, the retention of these synapses for the notched case was further interrogated. These investigations revealed and emphasized the importance of both notch geometry and anisotropy modification in order to further enhance the endurance of written spintronic states. In the near future, these results will be mapped to effective predictions for room temperature and elevated operation DW-MTJ memory retention, and experimentally verified when devices become available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Are Solar Energetic Particle Events and Type II Bursts Associated with Fast and Narrow Coronal Mass Ejections?

Gradual solar energetic ( E>10 MeV ) particle (SEP) events and metric through kilometric wavelength type II radio bursts are usually associated with shocks driven by fast ( V > 900 kms-1 ) and wide ( W>60deg ) coronal mass ejections (FW CMEs). This criterion was established empirically by several studies from solar cycle 23. The characteristic Alfven speed in the corona, which ranges over 500-1500 kms-1 at heights > 2 Ro, provides the minimum V requirement for a CME to drive a shock, but the general absence of SEP events or type II bursts with fast and narrow ( W<60deg ) CMEs has not been explained. We review and confirm the earlier studies with a more comprehensive comparison of SEP events and type II bursts with fast and narrow (FN) CMEs. We offer an explanation for the lack of SEP event and type II burst associations with FN CMEs in terms of recent heuristic arguments and modeling that show that the response of a magnetized plasma to the propagation of a CME depends on the CME geometry as well as on its speed. A clear distinction is made between a projectile that propagates through the medium to produce a bow shock, and a 3D piston that everywhere accumulates material to produce a broad shock and sheath. The bow shock is unfavorable for producing SEP events and type II bursts, but the 60 deg cut-off is not explained.

Kahler, S.W.↗

Assessing Utility of Different Orbits for a Lunar PNT Constellation

"As the Artemis program establishes a long-term presence on the lunar surface, lunar surface users will need an accurate source of Position, Navigation, and Timing (PNT). Under the LunaNet framework, NASA has proposed the Lunar Navigation System (LNS), a constellation of satellites in lunar orbit which will broadcast LunaNet’s Augmented Forward Signal (AFS) [1] to the lunar surface. LNS will have the objective of providing user position accuracies of 10 m and timing accuracies of 10 microseconds on the lunar surface and up to 200 km above the surface. Several classes of orbits are possibilities for the LNS constellation, including circular orbits, eccentric orbits, and various families of 3-body orbits. Each of these orbit classes have both advantages and disadvantages in their utility as PNT constellations. In this study we examine a wide variety of potential LNS constellations and evaluate their performance against several metrics, including: a. Dilution of Precision (DOP) performance over the lunar service volume and the South Pole region, during both the constellation buildup and after the constellation is completed. b. Number of satellites required to achieve required DOP performance. c. Relative dynamics of LNS satellites with respect to surface users. d. Signal in space error (SISE) performance, assuming orbit determination that includes the Global Positioning System (GPS), crosslinks, and monitor sites on the surface of the moon. e. Constellation insertion, station keeping, and disposal costs. f. Specific Concept of Operations (CONOPS) and engineering challenges associated with establishing and operating the constellations. (E.g., frequent station keeping, variable signal strength, divergent CONOPS for different planes of the constellation.) The system’s performance is considered both independently and in combination with other proposed systems, such as NASA’s Lunar Communication Relay and Navigation System (LCRNS) and ESA’s Moonlight constellation. The study concludes that a 3-plane Walker constellation at 6 lunar radii (slightly higher than GPS’s 4.2 earth radii from a geometry perspective, but significantly closer in an absolute sense) is a good candidate for the LNS constellation, providing the user with good SISE and DOP performance, while not requiring excessive station keeping or complex CONOPS. However, other viable concepts exist and could be chosen if early South Pole coverage is a driving factor in the system design."

navigation↗

Assessment of ICESat-2 Sea Ice Surface Classification with Sentinel-2 Imagery: Implications for Freeboard and New Estimates of Lead and Floe Geometry

NASA's Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) mission launched in September 2018 and is now providing high-resolution surface elevation profiling across the entire globe, including the sea ice cover of the Arctic and Southern Oceans. For sea ice applications, successfully discriminating returns between sea ice and open water is key for accurately determining freeboard (the extension of sea ice above local sea level) and new information regarding the geometry of sea ice floes and leads. We take advantage of near-coincident optical imagery obtained from the European Space Agency (ESA) Sentinel-2 (S-2) satellite over the Western Weddell Sea of the Southern Ocean in March 2019 and the Lincoln Sea of the Arctic Ocean in May 2019 to evaluate the surface classification scheme in the ICESat-2 ATL07 and ATL10 sea ice products. We find a high level of agreement between the ATL07 (specular) lead classification and visible leads in the S-2 imagery in these two coincident images across all six ICESat-2 beams, increasing our confidence in the freeboard products and deriving new estimates of the sea ice state. The S-2 overlays provide additional, albeit limited, evidence of the misclassification of dark leads due to clouds. Dark leads are no longer used to derive sea surface and thus freeboard as of the third release (r003) of the ICESat-2 sea ice products. We show estimates of lead fraction and more preliminary estimates of chord length (a proxy for floe size) using two metrics for classifying sea surface (lead) segments across both the Arctic and Southern Ocean for the first winter season of data collection.

A. A. Petty↗