Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ellipse fits”

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 37 records · Page 2

A Comparative Evaluation of Snowflake Particle Shape Estimation Techniques Used By the Precipitation Imaging Package (PIP), Multi-Angle Snowflake Camera (MASC), and Two-Dimensional Video Disdrometer (2DVD)

Measurements of snowflake particle shape are important for studying snow microphysics. While a number of instruments exist that are capable of measuring particle shape, this study focuses on the measurement techniques of three digital video disdrometers: the Precipitation Imaging Package (PIP), the Multi-Angle Snowflake Camera (MASC), and the Two-Dimensional Video Disdrometer (2DVD). To gain a better understanding of the relative strengths and weaknesses of these instruments and to provide a foundation upon which comparisons can be made between studies using data from different instruments, we perform a comparative analysis of the shape measurement algorithms employed by each of the three instruments by applying the algorithms to snowflake images captured by PIP during the ICE-POP 2018 field campaign. Our analysis primarily focuses on the measurement of the aspect ratio of either the particle itself, in the case of PIP and MASC, or of the particle bounding box, in the case of PIP and 2DVD. Both PIP and MASC use shape-fitting algorithms to measure aspect ratio. While our analysis of the MASC aspect ratio suggests that the measurements are reliable, our findings indicate that both the ellipse and rectangle aspect ratios produced by PIP underperformed considerably due to the shortcomings of the PIP shape-fitting techniques. We also demonstrate that reliable measurements of aspect ratio can be retrieved from PIP by reprocessing the raw PIP images using either the MASC ellipse-fitting algorithm or a tensor-based ellipse-fitting algorithm. Because of differences in instrument design, 2DVD produces measurements of particle horizontal and vertical extent rather than length and width. Furthermore, the 2DVD measurements of particle horizontal extent can be contaminated by horizontal particle motion. Our findings indicate that, although the correction technique used to remove the horizontal motion contamination performs remarkably well with snowflakes despite being designed for use with raindrops, the 2DVD measurements of particle horizontal extent are less reliable than those measured by PIP.

Charles N Helms↗

Acoustically induced oscillation and rotation of a large drop in space

A 2.5 cm diameter water drop was successfully deployed and manipulated in a triaxial acoustic resonance chamber during a 240 sec low-gravity SPAR rocket flight. Oscillation and rotation were induced by modulating and phase shifting the signals to the speakers. Portions of the film record were digitized and analyzed. Spectral analysis brought out the n = 2, 3, 4 free oscillation modes of the drop, its very low-frequency center-of-mass motion in the acoustic potential well, and the forced oscillation frequency. The drop boundaries were least-square fitted to general ellipses, providing eccentricities of the distorted drop. The normalized equatorial area of the rotating drop was plotted vs a rotational parameter, and was in excellent agreement with values derived from the theory of equilibrium shapes of rotating liquid drops.

Jacobi, N.↗

Near-Infrared Interferometric Images of the Solar System Sized Disk Surrounding the Herbig Ae/Be Star MWC 349A

We present images of the Herbig Ae/Be star MWC 349A at 1.65 and 2.27, and 3.08 micrometers, reconstructed from complex visibility data obtained with an aperture masking interferometric technique on the Keck I telescope. These images have an approximately elliptical shape, and are consistent with the expected shape of a nearly edge-on Keplerian disk. Visibility data were fitted with uniform ellipses with major axes 36 +/- 2, 47 +/- 2, and 62 +/- 1 mas, respectively. The axial ratio of the ellipses is approximately 0.5 +/- 0.1, and the major axis is at a position angle of 100 +/- 3 degrees, consistent with the position angle of the dark lane observed previously in the Very Large Array (VLA) radio continuum maps at 8 and 22 GHz, perpendicular to the symmetry axis of the bipolar lobes of H66(alpha) recombination line emission, and consistent with positions of the recombination line maser spots at 1.3 mm. At an assumed distance of 1.2 kpc, the linear sizes of the disk are 44 and 57 AU at 1.65 and 2.2 micrometers, respectively. The disk is the presumed source of ionized material in the bipolar outflow and ultracompact HII region around the star.

Danchi, W.C.↗

Robust pupil center detection using a curvature algorithm

Determining the pupil center is fundamental for calculating eye orientation in video-based systems. Existing techniques are error prone and not robust because eyelids, eyelashes, corneal reflections or shadows in many instances occlude the pupil. We have developed a new algorithm which utilizes curvature characteristics of the pupil boundary to eliminate these artifacts. Pupil center is computed based solely on points related to the pupil boundary. For each boundary point, a curvature value is computed. Occlusion of the boundary induces characteristic peaks in the curvature function. Curvature values for normal pupil sizes were determined and a threshold was found which together with heuristics discriminated normal from abnormal curvature. Remaining boundary points were fit with an ellipse using a least squares error criterion. The center of the ellipse is an estimate of the pupil center. This technique is robust and accurately estimates pupil center with less than 40% of the pupil boundary points visible.

Non-NASA Center↗

A model for the characterization of the spatial properties in vestibular neurons

Quantitative study of the static and dynamic response properties of some otolith-sensitive neurons has been difficult in the past partly because their responses to different linear acceleration vectors exhibited no "null" plane and a dependence of phase on stimulus orientation. The theoretical formulation of the response ellipse provides a quantitative way to estimate the spatio-temporal properties of such neurons. Its semi-major axis gives the direction of the polarization vector (i.e., direction of maximal sensitivity) and it estimates the neuronal response for stimulation along that direction. In addition, the semi-minor axis of the ellipse provides an estimate of the neuron's maximal sensitivity in the "null" plane. In this paper, extracellular recordings from otolith-sensitive vestibular nuclei neurons in decerebrate rats were used to demonstrate the practical application of the method. The experimentally observed gain and phase dependence on the orientation angle of the acceleration vector in a head-horizontal plane was described and satisfactorily fit by the response ellipse model. In addition, the model satisfactorily fits neuronal responses in three-dimensions and unequivocally demonstrates that the response ellipse formulation is the general approach to describe quantitatively the spatial properties of vestibular neurons.

Non-NASA Center↗

On the cost of approximating and recognizing a noise perturbed straight line or a quadratic curve segment in the plane

Approximation of noisy data in the plane by straight lines or elliptic or single-branch hyperbolic curve segments arises in pattern recognition, data compaction, and other problems. The efficient search for and approximation of data by such curves were examined. Recursive least-squares linear curve-fitting was used, and ellipses and hyperbolas are parameterized as quadratic functions in x and y. The error minimized by the algorithm is interpreted, and central processing unit (CPU) times for estimating parameters for fitting straight lines and quadratic curves were determined and compared. CPU time for data search was also determined for the case of straight line fitting. Quadratic curve fitting is shown to require about six times as much CPU time as does straight line fitting, and curves relating CPU time and fitting error were determined for straight line fitting. Results are derived on early sequential determination of whether or not the underlying curve is a straight line.

Cooper, D. B.↗

Eccentric features in Saturn's outer C ring

The present search for possible eccentric and inclined features in the outer C ring of Saturn measured all sharp-edged feature radii in Voyager C ring data. The Maxwell ringlet and two other narrow ringlets, 1.470R(s) and 1.495R(s) are found to be eccentric; the latter is best fitted by a model describing a freely precessing Keplerian ellipse, while the former is not conclusively fitted by either a resonant forcing or a free precession model. These two eccentric ringlets are compared with the Titan and Maxwell ringlets.

Porco, Carolyn C.↗

Precession of the epsilon ring of Uranus

It is noted that the outer and inner boundaries of the epsilon ring of Uranus can be fitted by aligned Keplerian ellipses. Four possible mechanisms for maintaining uniform precession in the epsilon ring are considered: the ring's self-gravity, precession due to a satellite, smooth pressure gradients, and shocklike phenomena. It is proposed that apse alignment is maintained by the self-gravity of the ring. In this case, a ring mass of approximately 5 x 10 to the 18th g and a mean surface density at quadrature of about 25 g/sq cm are estimated.

Goldreich, P.↗

Magnetotail views at 33R(sub E): IMP 8 magnetometer observations

This paper presents magnetic field vector (B) maps, electric current vector (curl B) maps, magnetic force (JxB) contour maps, and total field contour maps covering the full tail cross section in the yz plane. The maps are based on 16 years of 5-min averages of Interplanetary Monitoring Platform (IMP) 8 magnetic field data. During this time, IMP 8 traversed the tail between -25R(sub E) and -40R(sub E) in the x direction. Its average x distance was -33R(sub E). For this average distance we show separate maps for low and high dipole tilts, corresponding to equinox and northern hemisphere summer seasons. The low-tilt (equinox) maps show symmetrical field and current patterns; the high-tilt (solstice) maps show the cross-tail current sheet arcing above the equatorial diagonal in the center and dipping below it on the flanks. The shape of warped current sheet fits Fairfield's (1980) displaced ellipse model fairly well. The distance at which the current sheet is hinged to the magnetic equator is found to be 9.88R(sub E) and is independent of Kp. The z profile of current density shows a central peak, 3R(sub E) full width at half maximum, and smaller, flanking shoulders. A Harris sheet profile with a 7R(sub E) thickness fits the B(sub x) profile. Though these are magnetic field data, the JxB maps clearly outline the plasma sheet. This approach also gives 7R(sub E) thickness. Many of the average field and current features inferred and demonstrated in earlier studies are confirmed here; some of them are seen for the first time in full cross-section view. Among new features revealed are a large current vortex in the winter hemisphere lobe, a dawn-dusk asymmetry in the JxB force in the plasma sheet (it is stronger on the duskside), and a separation of the cross-tail current sheet into core and wing parts.

Kaymaz, Zerefsan↗

A Methodology to Assess the Capability of Engine Designs to Meet Closed-loop Performance and Operability Requirements

Designing a closed-loop controller for an engine requires balancing trade-offs between performance and operability of the system. One such trade-off is the relationship between the 95% response time and minimum high-pressure compressor (HPC) surge margin (SM) attained during acceleration from idle to takeoff power. Assuming a controller has been designed to meet some specification on response time and minimum HPC SM for a mid-life (nominal) engine, there is no guarantee that these limits will not be violated as the engine ages, particularly as it reaches the end of its life. A characterization for the uncertainty in this closed-loop system due to aging is proposed that defines elliptical boundaries to estimate worst-case performance levels for a given control design point. The results of this characterization can be used to identify limiting design points that bound the possible con- troller designs yielding transient results that do not exceed specified limits in response time or minimum HPC SM. This characterization involves performing Monte Carlo simulation of the closed-loop system with controller constructed for a set of trial design points and developing curve fits to describe the size and orientation of each ellipse; a binary search procedure is then employed that uses these fits to identify the limiting design point. The method is demonstrated through application to a generic turbofan engine model in closed- loop with a simplified controller; it is found that the limit for which each controller was designed was exceeded by less than 4.76%. Extension of the characterization to another trade-off, that between the maximum high-pressure turbine (HPT) entrance temperature and minimum HPC SM, showed even better results: the maximum HPT temperature was estimated within 0.76%. Because of the accuracy in this estimation, this suggests another limit that may be taken into consideration during design and analysis. It also demonstrates the extension of the characterization to other attributes that contribute to the performance or operability of the engine. Metrics are proposed that, together, provide information on the shape of the trade-off between response time and minimum HPC SM, and how much each varies throughout the life cycle, at the limiting design points. These metrics also facilitate comparison of the expected transient behavior for multiple engine models.

engine control↗

A Methodology to Assess the Capability of Engine Designs to Meet Closed-Loop Performance and Operability Requirements

Designing a closed-loop controller for an engine requires balancing trade-offs between performance and operability of the system. One such trade-off is the relationship between the 95 percent response time and minimum high-pressure compressor (HPC) surge margin (SM) attained during acceleration from idle to takeoff power. Assuming a controller has been designed to meet some specification on response time and minimum HPC SM for a mid-life (nominal) engine, there is no guarantee that these limits will not be violated as the engine ages, particularly as it reaches the end of its life. A characterization for the uncertainty in this closed-loop system due to aging is proposed that defines elliptical boundaries to estimate worst-case performance levels for a given control design point. The results of this characterization can be used to identify limiting design points that bound the possible controller designs yielding transient results that do not exceed specified limits in response time or minimum HPC SM. This characterization involves performing Monte Carlo simulation of the closed-loop system with controller constructed for a set of trial design points and developing curve fits to describe the size and orientation of each ellipse; a binary search procedure is then employed that uses these fits to identify the limiting design point. The method is demonstrated through application to a generic turbofan engine model in closed-loop with a simplified controller; it is found that the limit for which each controller was designed was exceeded by less than 4.76 percent. Extension of the characterization to another trade-off, that between the maximum high-pressure turbine (HPT) entrance temperature and minimum HPC SM, showed even better results: the maximum HPT temperature was estimated within 0.76 percent. Because of the accuracy in this estimation, this suggests another limit that may be taken into consideration during design and analysis. It also demonstrates the extension of the characterization to other attributes that contribute to the performance or operability of the engine. Metrics are proposed that, together, provide information on the shape of the trade-off between response time and minimum HPC SM, and how much each varies throughout the life cycle, at the limiting design points. These metrics also facilitate comparison of the expected transient behavior for multiple engine models.

systems analysis↗

Scalable Freeze-Tape-Casting Fabrication and Pore Structure Analysis of 3D LLZO Solid-State Electrolytes

Nonflammable solid-state electrolytes can potentially address the reliability and energy density limitations of lithium-ion batteries. Garnet-structured oxides such as Li 7 La 3 Zr 2 O 12 (LLZO) are some of the most promising candidates for solid-state devices. Here, three-dimensional (3D) solid-state LLZO frameworks with low tortuosity pore channels are proposed as scaffolds, into which active materials and other components can be infiltrated to make composite electrodes for solid-state batteries. To make the scaffolds, we employed aqueous freeze tape casting (FTC), a scalable and environmentally friendly method to produce porous LLZO structures. Using synchrotron radiation hard X-ray microcomputed tomography, we confirmed that LLZO films with porosities of up to 75% were successfully fabricated from slurries with a relatively wide concentration range. The acicular pore size and shape at different depths of scaffolds were quantified by fitting the pore shapes with ellipses, determining the long and short axes and their ratios, and investigating the equivalent diameter distribution. The results show that relatively homogeneous pore sizes and shapes were sustained over a long range along the thickness of the scaffold. Additionally, these pores had low tortuosity and the wall thickness distributions were found to be highly homogeneous. These are desirable characteristics for 3D solid electrolytes for composite electrodes, in terms of both the ease of active material infiltration and also minimization of Li diffusion distances in electrodes. The advantages of the FTC scaffolds are demonstrated by the improved conductivity of LLZO scaffolds infiltrated with poly(ethylene oxide)/lithium bis(trifluoromethanesulfonyl)imide (PEO/LITFSI) compared to those of PEO/LiTFSI films alone or composites containing LLZO particles.

36 MATERIALS SCIENCE↗

A model for the motion of the Philippine Sea plate consistent with NUVEL-1 and geological data

We investigate angular velocity vectors of the Philippine Sea (PH) plate relative to the adjacent major plates, Eurasia (EU) and Pacific (PA), and the smaller Caroline (CR) plate. Earthquake slip vector data along the Philippine Sea plate are inverted, subject to the constraint that EU-PA motion equals that predicted by the global relative plate model NUVEL-1. The resulting solution fails to satisfy geological constraints along the Caroline-Pacific boundary: convergence along the Mussau Trench and divergence along the Sorol Trough. We then seek solutions satisfying both the CR-PA boundary conditions and the Philippine Sea slip vector data, by adjusting the PA-PH and EU-PH best fitting poles within their error ellipses. We also consider northern Honshu to be part of the North American plate and impose the constraint that the Philippine Sea plate subducts beneath northern Honshu along the Sagmi Trough in a NNW-NW direction. Of the solutions satisfying these conditions, we select the best EU-PH as 48.2 deg N, 157.0 deg E, 1.09 deg/my, corresponding to a pole far from Japan and south of Kamchatka, and PA-PH, 1.2 deg N, 134.2 deg E, 1.00 deg/my. Predicted NA-PH and EU-PH convergence rates in central Honshu are consistent with estimated seismic slip rates. Previous estimates of the EU-PH pole close to central Honshu are inconsistent with extension within the Bonin backarc implied by earthquake slip vectors and NNW-NW convergence of the Bonin forearc at the Sagami Trough.

Seno, Tetsuzo↗

Dynamic Differential Image Circle Diameter Measurement Precision Assessment: Application to Burning Droplets

Dynamic measurement precision assessment has been achieved for a differential circle measurement application. Differential circle diameter measurement, in image analysis, typically requires fitting a circle model that optimizes for image distortions, defects or occlusions. The differential task occurs when precise measurements of diameter change are required given object size variation with time. An automated system was designed to provide diameter measurements and associated measurement precision of images of a fuel droplet undergoing combustion in zero gravity for the FLEX-2 dataset. An image gradient-based, least-squares boundary point fitting method to a circle or ellipse model is used for diameter measurement. The presence of soot aggregates poses significant challenges for diameter measurements when it occludes part of the droplet boundary. The precision of the diameter measurements depends upon the image quality. Using synthetic image simulations that model the soot behavior, we developed a model based on image quality measures that assesses the measurement precision for each individual diameter measurement. Thus, diameter measurements with precision assessments were made available for follow-up scientific analysis. As a result, the algorithm's success rate for measurable runs was 98%. In cases of limited occlusion, a measurement precision of ±0.2 pixels for the FLEX-2 dataset was achieved.

42 ENGINEERING↗

Laser transit anemometer software development program

Algorithms were developed for the extraction of two components of mean velocity, standard deviation, and the associated correlation coefficient from laser transit anemometry (LTA) data ensembles. The solution method is based on an assumed two-dimensional Gaussian probability density function (PDF) model of the flow field under investigation. The procedure consists of transforming the data ensembles from the data acquisition domain (consisting of time and angle information) to the velocity space domain (consisting of velocity component information). The mean velocity results are obtained from the data ensemble centroid. Through a least squares fitting of the transformed data to an ellipse representing the intersection of a plane with the PDF, the standard deviations and correlation coefficient are obtained. A data set simulation method is presented to test the data reduction process. Results of using the simulation system with a limited test matrix of input values is also given.

Abbiss, John B.↗

The rings of Uranus

Stellar occultation data on the Uranian rings are summarized. Spatial resolution is 0.1 km, for positions of ring segments, 4 km for structural details. The system consists of 9 narrow rings with average geometric albedo of 0.02. The data are well fitted by a kinematic model that describes the rings as ellipses, inclined to the equatorial plane of Uranus and precess due to harmonic terms in the Uranian potential fields. The main characteristics of the Uranian ring system that need a confirmed theoretical explanation are: narrow rings; sharp edges; uniform orbital precession; origin of the eccentricities and inclinations; the structure of the epsilon ring; the structure of the eta ring; and the origin of the ring system.

Elliot, J. L.↗

The rings of Uranus - Results of the 10 April 1978 occultation

Observations of the April 10, 1978, stellar occultation by the rings of Uranus are presented. Nine rings were observed, and their radii and widths are calculated. Rings eta, gamma, and delta are found to be most likely circular and coplanar, in agreement with previous analyses; the remaining rings are either noncircular or slightly inclined. The width of the epsilon ring is a linear function of its radius from the center of Uranus, projected onto the satellites' orbital plane; this suggests that it forms one continuous noncircular ring. The optical-depth profile of the epsilon ring has not changed significantly since March 1977. A model of this ring which fits all available observations adequately is that of a uniformly precessing Keplerian ellipse coplanar with the satellites' orbits. This model permits predictions of the radius and width of the epsilon ring for future occultations. The precession rate is used to determine J2 for Uranus, on the assumption that precession is caused solely by the planetary oblateness and not by satellite-ring interactions.

Nicholson, P. D.↗