Engineering PapersSearch

SEARCH · Engineering Papers

Results for “centroid estimation”

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 91 records · Page 5

Aerodynamic Performance and Static Stability at Mach Number 3.3 of an Aircraft Configuration Employing Three Triangular Wing Panels and a Body Equal Length

An aircraft configuration, previously conceived as a means to achieve favorable aerodynamic stability characteristics., high lift-drag ratio, and low heating rates at high supersonic speeds., was modified in an attempt to increase further the lift-drag ratio without adversely affecting the other desirable characteristics. The original configuration consisted of three identical triangular wing panels symmetrically disposed about an ogive-cylinder body equal in length to the root chord of the panels. This configuration was modified by altering the angular disposition of the wing panels, by reducing the area of the panel forming the vertical fin, and by reshaping the body to produce interference lift. Six-component force and moment tests of the modified configuration at combined angles of attack and sideslip were made at a Mach number of 3.3 and a Reynolds number of 5.46 million. A maximum lift-drag ratio of 6.65 (excluding base drag) was measured at a lift coefficient of 0.100 and an angle of attack of 3.60. The lift-drag ratio remained greater than 3 up to lift coefficient of 0.35. Performance estimates, which predicted a maximum lift-drag ratio for the modified configuration 27 percent greater than that of the original configuration, agreed well with experiment. The modified configuration exhibited favorable static stability characteristics within the test range. Longitudinal and directional centers of pressure were slightly aft of the respective centroids of projected plan-form and side area.

James, Carlton S.

Detection of Multi-Layer and Vertically-Extended Clouds Using A-Train Sensors

The detection of mUltiple cloud layers using satellite observations is important for retrieval algorithms as well as climate applications. In this paper, we describe a relatively simple algorithm to detect multiple cloud layers and distinguish them from vertically-extended clouds. The algorithm can be applied to coincident passive sensors that derive both cloud-top pressure from the thermal infrared observations and an estimate of solar photon pathlength from UV, visible, or near-IR measurements. Here, we use data from the A-train afternoon constellation of satellites: cloud-top pressure, cloud optical thickness, the multi-layer flag from the Aqua MODerate-resolution Imaging Spectroradiometer (MODIS) and the optical centroid cloud pressure from the Aura Ozone Monitoring Instrument (OMI). For the first time, we use data from the CloudSat radar to evaluate the results of a multi-layer cloud detection scheme. The cloud classification algorithms applied with different passive sensor configurations compare well with each other as well as with data from CloudSat. We compute monthly mean fractions of pixels containing multi-layer and vertically-extended clouds for January and July 2007 at the OMI spatial resolution (l2kmx24km at nadir) and at the 5kmx5km MODIS resolution used for infrared cloud retrievals. There are seasonal variations in the spatial distribution of the different cloud types. The fraction of cloudy pixels containing distinct multi-layer cloud is a strong function of the pixel size. Globally averaged, these fractions are approximately 20% and 10% for OMI and MODIS, respectively. These fractions may be significantly higher or lower depending upon location. There is a much smaller resolution dependence for fractions of pixels containing vertically-extended clouds (approx.20% for OMI and slightly less for MODIS globally), suggesting larger spatial scales for these clouds. We also find higher fractions of vertically-extended clouds over land as compared with ocean, particularly in the tropics and summer hemisphere.

Joiner, J.

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference

X-ray locations and sizes of clusters of galaxies observed with the HEAO 1 scanning modulation collimator

Observations with the HEAO 1 scanning modulation collimator (MC) have yielded improved locations for the centroids of X-ray sources in clusters of galaxies dominated by a giant central galaxy. For A85, A478, A1795, A2029, A2199, A2319 A, and SC 0316-444, these positions are consistent with those of the clusters' dominant cD galaxies. By comparing the HEAO 1 MC flux with that obtained with the HEAO 1 A-2 experiment, estimates of the X-ray sizes for these clusters have been obtained. The correlation of these sizes with other cluster properties is discussed.

Schwarz, J.

The interpretation of crustal dynamics data in terms of plate motions and regional deformation near plate boundaries

During our participation in the NASA Crustal Dynamics Project under NASA contract NAS-27339 and grant NAG5-814 for the period 1982-1991, we published or submitted for publication 30 research papers and 52 abstracts of presentations at scientific meetings. In addition, five M.I.T. Ph.D. students (Eric Bergman, Steven Bratt, Dan Davis, Jeanne Sauber, Anne Sheehan) were supported wholly or in part by this project during their thesis research. Highlights of our research progress during this period include the following: application of geodetic data to determine rates of strain in the Mojave block and in central California and to clarify the relation of such strain to the San Andreas fault and Pacific-North American plate motions; application of geodetic data to infer post seismic deformation associated with large earthquakes in the Imperial Valley, Hebgen Lake, Argentina, and Chile; determination of the state of stress in oceanic lithosphere from a systematic study of the centroid depths and source mechanisms of oceanic intraplate earthquakes; development of models for the state of stress in young oceanic regions arising from the differential cooling of the lithosphere; determination of the depth extent and rupture characteristics of oceanic transform earthquakes; improved determination of earthquake slip vectors in the Gulf of California, an important data set for the estimation of Pacific-North American plate motions; development of models for the state of stress and mechanics of fold-and-thrust belts and accretionary wedges; development of procedures to invert geoid height, residual bathymetry, and differential body wave travel time residuals for lateral variations in the characteristic temperature and bulk composition of the oceanic upper mantle; and initial GPS measurements of crustal deformation associated with the Imperial-Cerro Prieto fault system in southern California and northern Mexico. Full descriptions of the research conducted on these topics may be found in the Semi-Annual status Reports submitted regularly to NASA over the course of this project and in the publications listed.

Solomon, Sean C.

Assessment of Error in Synoptic-Scale Diagnostics Derived from Wind Profiler and Radiosonde Network Data

A topic of current practical interest is the accurate characterization of the synoptic-scale atmospheric state from wind profiler and radiosonde network observations. We have examined several related and commonly applied objective analysis techniques for performing this characterization and considered their associated level of uncertainty both from a theoretical and a practical standpoint. A case study is presented where two wind profiler triangles with nearly identical centroids and no common vertices produced strikingly different results during a 43-h period. We conclude that the uncertainty in objectively analyzed quantities can easily be as large as the expected synoptic-scale signal. In order to quantify the statistical precision of the algorithms, we conducted a realistic observing system simulation experiment using output from a mesoscale model. A simple parameterization for estimating the uncertainty in horizontal gradient quantities in terms of known errors in the objectively analyzed wind components and temperature is developed from these results.

Mace, Gerald G.

Geoscience Laser Altimeter System (GLAS) for the ICESat Mission

The Geoscience Laser Altimeter System (GLAS) is a new generation lidar and is the primary science payload for NASA's ICESat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical distribution of clouds and aerosols on a global scale. GLAS will be integrated onto a small spacecraft built by Ball Aerospace, and will be launched into a polar orbit with a 590-630 km altitude at an inclination of 94 degrees. ICESat is is currently planned to launch in winter 2002/03 and GLAS is designed to operate continuously in space for a minimum of 3 years. GLAS will measure the vertical distance from orbit to the Earth's surface with pulses from a ND:YAG laser at a 40 Hz rate. Each 6 nsec wide 1064 nm laser pulse is used to produce a single range measurement. On the surface, the laser footprints have 66 m diameter and approx. 170 m center-center spacings. The GLAS receiver uses a I m diameter telescope to detect laser backscatter and a Si APD to detect the 1064 nm signals. The detector's output is sampled by a digital ranging receiver, which records each transmitted pulse and surface echo waveform with 1 nsec (15 cm) resolution. Each echo pulse is digitized and is reported to ground with a record length of from 200 to 544 samples, depending on the spacecraft's location . The GLAS location and epoch times are measured by a precision GPS receiver carried on the ICESat spacecraft. Initial processing of the echo waveforms within GLAS permits discrimination between cloud and surface echoes for selecting appropriate waveform samples. This selection is guided by an on-board DEM which is used to set the boundaries for the echo pulse search algorithm. Subsequent ground-based echo pulse analysis, along with GPS-based clock frequency estimates, permit final determination of the range to the surface, degree of pulse spreading, and vertical distribution of any vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when measuring over tilted surfaces, such as near the boundaries of ice sheets. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid angle to better than 10 urad is needed. GLAS uses a stellar reference system (SRS) to measure the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). GLAS will also measure the vertical distributions of clouds and aerosols by recording the vertical profiles of laser pulse backscatter at both 1064 and 532 nm. The 1064 rim measurements use the Si APD detector and will be used to measure the height and echo pulse shape from thicker clouds. The lidar receiver at 532 nm uses a narrow bandwidth etalon filter and highly sensitive photon counting detectors. The 532 nm backscatter profiles will be used to measure the vertical extent of thinner clouds and the atmospheric boundary layer. The GLAS instrument component development is complete and the instrument is undergoing final testing and qualification at NASA-Goddard. The GLAS "as-built" characteristics and its expected measurement performance will be discussed.

Abshire, James B.

1919+479: Big WAT in a poor cluster

New x-ray, optical, and redshift data are presented for the cluster of galaxies associated with the giant, 1 Mpc diameter, wide-angle tailed (WAT) radio galaxy 1919+479. The ROSAT Position-Sensitive Proportional Counter (PSPC) pointed observation shows an x-ray peak on the WAT and elongated diffuse emission tracing the galaxy distribution. In addition, an asymmetric extension of emission exists between the tails of the WAT. The fitting of a Raymond-Smith thermal model to the x-ray spectra suggests an approximately = 2 keV temperature intracluster medium (ICM). The cooling time and irregular morphology rule out a cluster-wide cooling flow. The x-ray luminosity and temperature are consistent with the velocity dispersion, 480 km/s, estimated from 31 galaxy velocities. However, this velocity distribution is significantly non-Gaussian, which along with the x-ray morphology, suggests incomplete virialization in the cluster. Substructure analysis does not reveal significant clumping in the velocities/positions; but, the spatial distribution of galaxies is very elongated. Also, the cD galaxy producing the WAT does not have a significant radial peculiar velocity with respect to the cluster centroid. These characteristics are consistent with a merger scenario in which a subcluster has crossed the cluster core in the plane of the sky and has dispersed. We compare this cluster with the post-merger cluster Abell 2634 (Pinkney et al., 1993), containing the prototype WAT, and with a recent N-Body/Hydro simulation of merging clusters. The similarities indicate that the cluster 1919+479 may be the poor extreme of a class of clusters in which the bulk motion in the ICM, caused by a subcluster merger, in shaping the central WAT.

Pinkney, Jason

Compact Optoelectronic Compass

A compact optoelectronic sensor unit measures the apparent motion of the Sun across the sky. The data acquired by this chip are processed in an external processor to estimate the relative orientation of the axis of rotation of the Earth. Hence, the combination of this chip and the external processor finds the direction of true North relative to the chip: in other words, the combination acts as a solar compass. If the compass is further combined with a clock, then the combination can be used to establish a threeaxis inertial coordinate system. If, in addition, an auxiliary sensor measures the local vertical direction, then the resulting system can determine the geographic position. This chip and the software used in the processor are based mostly on the same design and operation as those of the unit described in Micro Sun Sensor for Spacecraft (NPO-30867) elsewhere in this issue of NASA Tech Briefs. Like the unit described in that article, this unit includes a small multiple-pinhole camera comprising a micromachined mask containing a rectangular array of microscopic pinholes mounted a short distance in front of an image detector of the active-pixel sensor (APS) type (see figure). Further as in the other unit, the digitized output of the APS in this chip is processed to compute the centroids of the pinhole Sun images on the APS. Then the direction to the Sun, relative to the compass chip, is computed from the positions of the centroids (just like a sundial). In the operation of this chip, one is interested not only in the instantaneous direction to the Sun but also in the apparent path traced out by the direction to the Sun as a result of rotation of the Earth during an observation interval (during which the Sun sensor must remain stationary with respect to the Earth). The apparent path of the Sun across the sky is projected on a sphere. The axis of rotation of the Earth lies at the center of the projected circle on the sphere surface. Hence, true North (not magnetic North), relative to the chip, can be estimated from paths of the Sun images across the APS. In a test, this solar compass has been found to yield a coarse estimate of the North (within tens of degrees) in an observation time of about ten minutes. As expected, the accuracy was found to increase with observation time: after a few hours, the estimated direction of the rotation axis becomes accurate to within a small fraction of a degree.

Christian, Carl

GHRS Observations of Cool, Low-Gravity Stars: The Outer Atmosphere and Wind of the Nearby K Supergiant Lambda Velorum - 5

UV spectra of lambda Velorum taken with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope are used to probe the structure of the outer atmospheric layers and wind and to estimate the mass-loss rate from this K5 lb-II supergiant. VLA radio observations at lambda = 3.6 cm are used to obtain an independent check on the wind velocity and mass-loss rate inferred from the UV observations, Parameters of the chromospheric structure are estimated from measurements of UV line widths, positions, and fluxes and from the UV continuum flux distribution. The ratios of optically thin C II] emission lines indicate a mean chromospheric electron density of log N(sub e) approximately equal 8.9 +/- 0.2 /cc. The profiles of these lines indicate a chromospheric turbulence (v(sub 0) approximately equal 25-36 km/s), which greatly exceeds that seen in either the photosphere or wind. The centroids of optically thin emission lines of Fe II and of the emission wings of self-reversed Fe II lines indicate that they are formed in plasma approximately at rest with respect to the photosphere of the star. This suggests that the acceleration of the wind occurs above the chromospheric regions in which these emission line photons are created. The UV continuum detected by the GHRS clearly traces the mean flux-formation temperature as it increases with height in the chromosphere from a well-defined temperature minimum of 3200 K up to about 4600 K. Emission seen in lines of C III] and Si III] provides evidence of material at higher than chromospheric temperatures in the outer atmosphere of this noncoronal star. The photon-scattering wind produces self-reversals in the strong chromospheric emission lines, which allow us to probe the velocity field of the wind. The velocities to which these self-absorptions extend increase with intrinsic line strength, and thus height in the wind, and therefore directly map the wind acceleration. The width and shape of these self-absorptions reflect a wind turbulence of approximately equal 9-21 km/s. We further characterize the wind by comparing the observations with synthetic profiles generated with the Lamers et al. Sobolev with Exact Integration (SEI) radiative transfer code, assuming simple models of the outer atmospheric structure. These comparisons indicate that the wind in 1994 can be described by a model with a wind acceleration parameter beta approximately 0.9, a terminal velocity of 29-33 km/s, and a mass-loss rate approximately 3 x 10(exp -9) solar M/yr. Modeling of the 3.6 cm radio flux observed in 1997 suggests a more slowly accelerating wind (higher beta) and/or a higher mass-loss rate than inferred from the UV line profiles. These differences may be due to temporal variations in the wind or from limitations in one or both of the models. The discrepancy is currently under investigation.

Carpenter, Kenneth G.

MTUQ: a framework for estimating moment tensors, point forces, and their uncertainties

SUMMARY We introduce MTUQ, an open-source Python package for seismic source estimation and uncertainty quantification, emphasizing flexibility and operational scalability. MTUQ provides MPI-parallelized grid search and global optimization capabilities, compatibility with 1-D and 3-D Green’s function database formats, customizable data processing, C-accelerated waveform and first-motion polarity misfit functions, and utilities for plotting seismic waveforms and visualizing misfit and likelihood surfaces. Applicability to a range of full- and constrained-moment tensor, point force, and centroid inversion problems is possible via a documented application programming interface, accompanied by example scripts and integration tests. We demonstrate the software using three different types of seismic events: (1) a 2009 intraslab earthquake near Anchorage, Alaska; (2) an episode of the 2021 Barry Arm landslide in Alaska; and (3) the 2017 Democratic People’s Republic of Korea underground nuclear test. With these events, we illustrate the well-known complementary character of body waves, surface waves, and polarities for constraining source parameters. We also convey the distinct misfit patterns that arise from each individual data type, the importance of uncertainty quantification for detecting multimodal or otherwise poorly constrained solutions, and the software’s flexible, modular design.

58 GEOSCIENCES

Optical Navigation Attitude Estimation and Calibration Performance Improvement using Outlier Rejection

Spacecraft optical navigation (OpNav) systems process a sequence of images of celestial bodies against a starfield background to estimate the position and velocity of the vehicle. While attitude is sometimes available from an onboard star tracker, it is often desirable to recognize the background stars in the OpNav images to better align the image. While many image processing algorithms exist for finding stars, efficiency and reliability remain key issues in the presence of extended bodies(e.g. the Moon, Earth), especially when attempting to solve the full lost-in-space problem. Some star outliers(stars identified with high residuals)could appear in the camera field of view, however using them in the attitude estimation or camera calibration would lead to less accurate results. Therefore, we require new and robust approaches to remove these outliers before any further processing. The emphasis of the work is on developing a simple and robust iterative technique to detect and reject the outliers which could be found in any frame during the lost in space attitude determination or during the camera calibration. These outliers are determined based on the residuals of the centroids of the detected stars and the corresponding location using the star catalog. If the residuals exceed a predetermined threshold value, the object will be detected as an outlier and will be removed before another attitude determination and calibration iteration is performed. The performance for both attitude determination and on-orbit camera calibration are improved by an almost two-fold increase in accuracy when applying this outlier rejection technique.

OpNav

Weather, Climate, and Society: New Demands on Science and Services

A new algorithm has been constructed to estimate the path length of lightning channels for the purpose of improving the model predictions of lightning NOx in both regional air quality and global chemistry/climate models. This algorithm was tested and applied to VHF signals detected by the North Alabama Lightning Mapping Array (NALMA). The accuracy of the algorithm was characterized by comparing algorithm output to the plots of individual discharges whose lengths were computed by hand. Several thousands of lightning flashes within 120 km of the NALMA network centroid were gathered from all four seasons, and were analyzed by the algorithm. The mean, standard deviation, and median statistics were obtained for all the flashes, the ground flashes, and the cloud flashes. Channel length distributions were also obtained for the different seasons.

Source record

Fast Simulators for Satellite Cloud Optical Centroid Pressure Retrievals, 1. Evaluation of OMI Cloud Retrievals

The cloud Optical Centroid Pressure (OCP), also known as the effective cloud pressure, is a satellite-derived parameter that is commonly used in trace-gas retrievals to account for the effects of clouds on near-infrared through ultraviolet radiance measurements. Fast simulators are desirable to further expand the use of cloud OCP retrievals into the operational and climate communities for applications such as data assimilation and evaluation of cloud vertical structure in general circulation models. In this paper, we develop and validate fast simulators that provide estimates of the cloud OCP given a vertical profile of optical extinction. We use a pressure-weighting scheme where the weights depend upon optical parameters of clouds and/or aerosol. A cloud weighting function is easily extracted using this formulation. We then use fast simulators to compare two different satellite cloud OCP retrievals from the Ozone Monitoring Instrument (OMI) with estimates based on collocated cloud extinction profiles from a combination of CloudS at radar and MODIS visible radiance data. These comparisons are made over a wide range of conditions to provide a comprehensive validation of the OMI cloud OCP retrievals. We find generally good agreement between OMI cloud OCPs and those predicted by CloudSat. However, the OMI cloud OCPs from the two independent algorithms agree better with each other than either does with the estimates from CloudSat/MODIS. Differences between OMI cloud OCPs and those based on CloudSat/MODIS may result from undetected snow/ice at the surface, cloud 3-D effects, low altitude clouds missed by CloudSat, and the fact that CloudSat only observes a relatively small fraction of an OMI field-of-view.

Joiner, J.

Secondary Fragmentation of Comet Shoemaker-Levy 9 and the Ramifications for the Progenitor's Breakup in July 1992

Comprehensive analysis of discrete events of secondary fragmentation leads to a conceptually new understanding of the process of disintegration of comet Shoemaker-Levy 9. We submit that the jovian tidal forces inflicted extensive cracks throughout the interior of the original nucleus but did not split it apart. The initial disruption was apparently accomplished by stresses exerted on the cracked object by its fast rotation during the early post-perijove period of time. We argue that this disruption was in fact a rapid sequence of episodes during July 1992 that gave birth to the 12 on-train, or primary, fragments: A, C, D, E, G, H, K, L, Q (later Q(sub 1)), R, S, and W. The discrete events of secondary fragmentation, which gave birth to the off-train fragments, are understood in this scenario as stochastic manifestations of the continuing process of progressive disintegration. Of the 13 off-train fragments considered, nine were secondary--B, F, G(sub 2), M, N, P (later P(sub 2) or P(sub 2a)), Q(sub 2), U, and V--and four tertiary (J, P(sub 1), P(sub 2b), and T). The separation parameters of 11 off-train fragments were determined. The vectorial distribution of separation velocities of these fragments shows a strong concentration toward a great circle, unquestionably an effect of the approximately conserved angular momentum of the progenitor comet since the time of its initial disruption. Also apparent is their clumping (except for P(sub 1)) to a segment along the great circle, implying that the fragments were consistently released from one side of their parents, thus explaining for the first time why the off-train fragments preferentially appeared on one side of the nuclear train. In order to obtain a consistent solution, our model requires that the points of separation be on the antisolar side of the parent fragments, where thermal stresses are likely to enhance the effect of rotation. The episodes of secondary fragmentation are found to have nine months after the close encounter with Jupiter in early July 1992, and the separation velocities ranged between 0.36 and 1.7 m/s. The spin-axis position is determined to have been nearly in the jovicentric orbit plane, which rules out the Asphaug-Benz-Solem strengthless aggregate model as a plausible breakup hypothesis. Since the separation velocities are rotational in nature, they cannot substantially exceed the critical limit for centrifugal breakup and offer an estimate for the original nuclear dimensions. The comet's nucleus is found to have been approximately 10 km in diameter and spinning rapidly. With the exception of P(sub 1) and apparently also P(sub 2) and F, no nongravitational deceleration was detected in the motions of the off-train fragments. Serious doubts are cast on continuing appreciable activity of any of these fragments. Indeed, when it was necessary to introduce a deceleration into the equations of motion, the effect appears to have been due to the action of solar radiation pressure on the centroid of centimeter-sized particulates in the disintegrating condensations.

Sekanina, Z.

Q-Cluster: Quantum Error Mitigation Through Noise-Aware Unsupervised Learning

Quantum error mitigation (QEM) is critical in reducing the impact of noise in the pre-fault-tolerant era, and is expected to complement error correction in fault-tolerant quantum computing (FTQC). In this work, we propose a novel QEM approach, Q-Cluster, that uses unsupervised learning (clustering) to reshape the measured bit-string distribution. Our approach starts with a simplified bit-flip noise model. It first performs clustering on noisy measurement results, i.e., bit-strings, based on the Hamming distance. The centroid of each cluster is calculated using a qubit-wise majority vote. Next, the noisy distribution is adjusted with the clustering outcomes and the bitflip error rates using Bayesian inference. Our simulation results show that Q-Cluster can mitigate high noise rates (up to 40% per qubit) with the simple bit-flip noise model. However, real quantum computers do not fit such a simple noise model. To address the problem, we (a) apply Pauli twirling to tailor the complex noise channels to Pauli errors, and (b) employ a machine learning model, ExtraTrees regressor, to estimate an effective bit-flip error rate using a feature vector consisting of machine calibration data (gate & measurement error rates), circuit features (number of qubits, numbers of different types of gates, etc.) and the shape of the noisy distribution (entropy). Our experimental results show that our proposed Q-Cluster scheme improves the fidelity by a factor of 1.46x, on average, compared to the unmitigated output distribution, for a set of low-entropy benchmarks on five different IBM quantum machines. Our approach outperforms the state-of-art QEM approaches RZNE [28], M3 [24], Hammer [35], and QBEEP [33] by 1.26x,1.29x,1.47x, and 2.65 x, respectively.

42 ENGINEERING

Experimental Investigation of a Hypersonic Glider Configuration at a Mach Number of 6 and at Full-Scale Reynolds Numbers

The aerodynamic characteristics of a hypersonic glider configuration, consisting of a slender ogive cylinder with three highly swept wings, spaced 120 apart, with the wing chord equal to the body length, were investigated experimentally at a Mach number of 6 and at Reynolds numbers from 6 to 16 million. The objectives were to evaluate the theoretical procedures which had been used to estimate the performance of the glider, and also to evaluate the characteristics of the glider itself. A principal question concerned the viscous drag at full-scale Reynolds number, there being a large difference between the total drags for laminar and turbulent boundary layers. It was found that the procedures which had been applied for estimating minimum drag, drag due to lift, lift curve slope, and center of pressure were generally accurate within 10 percent. An important exception was the non-linear contribution to the lift coefficient which had been represented by a Newtonian term. Experimentally, the lift curve was nearly linear within the angle-of-attack range up to 10 deg. This error affected the estimated lift-drag ratio. The minimum drag measurements indicated that substantial amounts of turbulent boundary layer were present on all models tested, over a range of surface roughness from 5 microinches maximum to 200 microinches maximum. In fact, the minimum drag coefficients were nearly independent of the surface smoothness and fell between the estimated values for turbulent and laminar boundary layers, but closer to the turbulent value. At the highest test Reynolds numbers and at large angles of attack, there was some indication that the skin friction of the rough models was being increased by the surface roughness. At full-scale Reynolds number, the maximum lift-drag ratio with a leading edge of practical diameter (from the standpoint of leading-edge heating) was 4.0. The configuration was statically and dynamically stable in pitch and yaw, and the center of pressure was less than 2-percent length ahead of the centroid of plan-form area.

Seiff, Alvin

An Algorithm for Obtaining the Distribution of 1-Meter Lightning Channel Segment Altitudes for Application in Lightning NOx Production Estimation

An algorithm has been developed to estimate the altitude distribution of one-meter lightning channel segments. The algorithm is required as part of a broader objective that involves improving the lightning NOx emission inventories of both regional air quality and global chemistry/climate models. The algorithm was tested and applied to VHF signals detected by the North Alabama Lightning Mapping Array (NALMA). The accuracy of the algorithm was characterized by comparing algorithm output to the plots of individual discharges whose lengths were computed by hand; VHF source amplitude thresholding and smoothing were applied to optimize results. Several thousands of lightning flashes within 120 km of the NALMA network centroid were gathered from all four seasons, and were analyzed by the algorithm. The mean, standard deviation, and median statistics were obtained for all the flashes, the ground flashes, and the cloud flashes. One-meter channel segment altitude distributions were also obtained for the different seasons.

Peterson, Harold