Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “least square approximations”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Theory of wide-angle photometry from standard stars

Wide angle celestial structures, such as bright comet tails and nearby galaxies and clusters of galaxies, rely on photographic methods for quantified morphology and photometry, primarily because electronic devices with comparable resolution and sky coverage are beyond current technological capability. The problem of the photometry of extended structures and of how this problem may be overcome through calibration by photometric standard stars is examined. The perfect properties of the ideal field of view are stated in the guise of a radiometric paraxial approximation, in the hope that fields of view of actual telescopes will conform. Fundamental radiometric concepts are worked through before the issue of atmospheric attenuation is addressed. The independence of observed atmospheric extinction and surface brightness leads off the quest for formal solutions to the problem of surface photometry. Methods and problems of solution are discussed. The spectre is confronted in the spirit of standard stars and shown to be chimerical in that light, provided certain rituals are adopted. After a brief discussion of Baker-Sampson polynomials and the vexing issue of saturation, a pursuit is made of actual numbers to be expected in real cases. While the numbers crunched are gathered ex nihilo, they demonstrate the feasibility of Newton's method in the solution of this overdetermined, nonlinear, least square, multiparametric, photometric problem.

Usher, Peter D.↗
DOE OSTI2025

Uncertainty Quantification Enabled by Automatic Differentiation for Hydrodynamic Simulation of Shock‐to‐Detonation Transition in High Explosives

Quantifying the effects of uncertainty in a reactive burn model on the run-to-detonation time in high explosives (HEs) provides a robust methodology for assessing the probability of an HE failing the IHE qualification standard. Moreover, uncertainty quantification helps evaluate whether the model calibration accurately represents data outside the calibration set. This study uses a specialized hydrodynamic simulation code for modeling detonation to determine the run-to-detonation time of the HE PBX 9502 for various impact velocities. To quickly approximate uncertainties in the model, a surrogate was constructed using a Taylor series expansion centered at the mean of the input parameters. To obtain the sensitivities required for constructing the Taylor series, HYP-percomplex Automatic Differentiation (HYPAD) was implemented. HYPAD is a methodology for infusing existing codes with automatic differentiation capabilities by augmenting variables with one or more imaginary units to compute step-size independent partial derivatives. These derivatives are accurate to machine precision with respect to the implemented numerical algorithm, meaning their accuracy reflects that of the underlying method (e.g., integration or discretization schemes). Using reduced order modeling techniques, the mean and standard deviation of the run-to-detonation time of a shock within PBX 9502 were computed for a number of initial impact velocities. A weighted least squares regression was then performed to obtain a best fit curve and prediction interval for the computed statistics. Historical data points from explosively driven wedge tests were utilized to validate the prediction interval, ensuring its reliability in predicting future outcomes. With this prediction interval and a known safety constraint curve, the most probable point of failure and the probability of failure for the HE PBX 9502 were determined.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗
NASA NTRS1994

Compact radio sources in the starburst galaxy M82 and the Sigma-D relation for supernova remnants

We have obtained an 8.4 GHz Very Large Array (VLA) A-array map of the starburst galaxy M82 with a resolution Full Width at Half Maximum (FWHM) approximately 0.182 sec. About 50 compact radio sources in the central region of M82 were detected with a peak surface brightness approximately greater than 10(exp -17) W/Hz/sq m/sr. Comparison with previous observations shows that most sources are declining in flux. Three previously visible sources have faded into the background of our map (approximately less than 0.2 mJy/beam), while a few sources, including the second and third brightest radio sources in M82, may have increased slightly in flux over the last decade. No new radio supernova was found. The birth rate of the compact radio sources is estimated to be 0.11 + or - 0.05/yr. We attribute the population of such bright, small supernova remnants (SNRs) in M82 to the high pressure in the central region that can truncate the mass loss during a red supergiant phase or allow dense ionized clouds to be present. The compact radio sources obey a Sigma(radio surface brightness) - D(diameter) relation which is remarkably similar to that followed by supernova remnants in the Galaxy and the Magellanic Clouds and by two of the strongest known extragalactic radio supernovae: SN 1986J and SN 1979C. A least-squares fit to the SNR data gives: Sigma(sub 8.4 GHz) (W/Hz/sq m/sr) = 4.4 x 10(exp -16) D(sub pc)(exp -3.5 +/- 0.1) covering seven orders of magnitude in Sigma. Possible selection effects are discussed and a theoretical discussion of the correlation is presented.

Huang, Z. P.↗
NASA NTRS2005

Experimental and Analytical Determination of the Geometric Far Field for Round Jets

An investigation was conducted at the NASA Glenn Research Center using a set of three round jets operating under unheated subsonic conditions to address the question: "How close is too close?" Although sound sources are distributed at various distances throughout a jet plume downstream of the nozzle exit, at great distances from the nozzle the sound will appear to emanate from a point and the inverse-square law can be properly applied. Examination of normalized sound spectra at different distances from a jet, from experiments and from computational tools, established the required minimum distance for valid far-field measurements of the sound from subsonic round jets. Experimental data were acquired in the Aeroacoustic Propulsion Laboratory at the NASA Glenn Research Center. The WIND computer program solved the Reynolds-Averaged Navier-Stokes equations for aerodynamic computations; the MGBK jet-noise prediction computer code was used to predict the sound pressure levels. Results from both the experiments and the analytical exercises indicated that while the shortest measurement arc (with radius approximately 8 nozzle diameters) was already in the geometric far field for high-frequency sound (Strouhal number >5), low-frequency sound (Strouhal number <0.2) reached the geometric far field at a measurement radius of at least 50 nozzle diameters because of its extended source distribution.

Koch, L. Danielle↗
NASA NTRS1993

Birkeland currents in the plasma sheet

A search was conducted for the signatures of Birkeland currents in the Earth's magnetic tail, using observed values of B(sub x) and B(sub y) from large sets of spacecraft data. The data were binned by x and y for -10 greater than x(sub GSM) greater than -35 and absolute value of y(sub GSM) less than or equal to 20 R(sub E) (less than or equal to 30 R(sub E) for x(sub GSM) less than or equal to -25 R(sub E)) and in each bin their distribution in the (B(sub x), B(sub y)) plane was fitted by least squares to a piecewise linear function. That gave average x-y distributions of the flaring angle between B(sub xy) and the x direction, as well as that angle's variation across the thickness of the plasma sheet. Angles obtained in the central plasma sheet differed from those derived near the lobe boundary. That is the expected signature if earthward or tailward Birkeland current sheets are embedded in the plasma sheet, and from this dfiference we derived the dawn-dusk profiles of the tail Birkeland currents for several x(sub GSM) intervals. It was found that (1) the Birkeland currents have the sense of region 1 currents, when mapped to the ionosphere; (2) both the linear current density (kiloamperes/R(sub E)) and the net magnitude of the field-aligned currents decrease rapidly down the tail; (3) the total Birkeland current at x approximately equals -10 R(sub E) equals approximately equals 500-700 kA, which is approx. 30% of the net region 1 current observed at ionospheric altitudes, in agreement with model mapping results; and (4) the B(sub z) and B(sub y) components of the interplanetary magnetic field influence the distribution of Birkeland currents in the tail.

Tsyganenko, Nikolai A.↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Magnetic flux rope versus the spheromak as models for interplanetary magnetic clouds

Magnetic clouds form a subset of interplanetary ejecta with well-defined magnetic and thermodynamic properties. Observationally, it is well established that magnetic clouds expand as they propagate antisunward. The aim of this paper is to compare and contrast two models which have been proposed for the global magnetic field line topology of magnetic clouds: a magnetic flux tube geometry, on the one hand, and a spheromak geometry (including possible higher multiples), on the other. Traditionally, the magnetic structure of magnetic clouds has been modeled by force-free configurations. In a first step, we therefore analyze the ability of static force-free models to account for the asymmetries observed in the magnetic field profiles of magnetic clouds. For a cylindrical flux tube the magnetic field remains symmetric about closest approach to the magnetic axis on all spacecraft orbits intersecting it, whereas in a spheromak geometry one can have asymmetries in the magnetic field signatures along some spacecraft trajectories. The duration of typical magnetic cloud encounters at 1 AU (1 to 2 days) is comparable to their travel time from the Sun to 1 AU and thus magnetic clouds should be treated as strongly nonstationary objects. In a second step, therefore, we abandon the static approach and model magnetic clouds as self-similarly evolving MHD configurations. In our theory, the interaction of the expanding magnetic cloud with the ambient plasma is taken into account by a drag force proportional to the density and the velocity of expansion. Solving rigorously the full set of MHD equations, we demonstrate that the asymmetry in the magnetic signature may arise solely as a result of expansion. Using asymptotic solutions of the MHD equations, we least squares fit both theoretical models to interplanetary data. We find that while the central part of the magnetic cloud is adequately described by both models, the 'edges' of the cloud data are modeled better by the magnetic flux tube. Further comparisons of the two models necessarily involve thermodynamic properties, since real magnetic configurations are never exactly force-free and gas pressure plays an essential role. We consider a polytropic gas. Our theoretical analysis shows that the self-similar expansion of a magnetic flux tube requires the polytropic index gamma to be less than unity. For the spheromak, however, self-similar, radially expanding solutions are known only for gamma equal to 4/3. This difference, therefore, yields a good way of distinguishing between the two geometries. It has been shown recently that the polytropic relationship is applicable to magnetic clouds and that the corresponding polytropic index is approximately 0.5. This observational result is consistent with the self-similar model of the magnetic flux rope but is in conflict with the self-similar spheromak model.

Farrugia, C. J.↗
NASA NTRS2003

High Resolution Spectroscopy to Support Atmospheric Measurements

Spectroscopic parameters (such as line position, intensity, broadening and shifting coefficients and their temperature dependences, line mixing coefficients etc.) for various molecular species of atmospheric interest are determined. In order to achieve these results, infrared spectra of several molecular bands are obtained using high-resolution recording instruments such as tunable diode laser spectrometer and Fourier transform spectrometers. Using sophisticated analysis routines (Multispectrum nonlinear least squares technique) these high-resolution infrared spectra are processed to determine the various spectral line parameters that are cited above. Spectra were taken using the McMath-Pierce Fourier transform spectrometer (FTS) at the National Solar Observatory on Kitt Peak, Arizona as well as the Bruker FTS at the Pacific Northwest National Laboratory (PNNL) at Richland, Washington. Most of the spectra are acquired not only at room temperature, but also at several different cold temperatures. This procedure is necessary to study the variation of the spectral line parameters as a function of temperature in order to simulate the Earth's and other planetary atmospheric environments. Depending upon the strength or weakness of the various bands recorded and analyzed, the length(s) of the absorption cells in which the gas samples under study are kept varied from a few centimeters up to several meters and the sample temperatures varied from approximately +30 C to -63 C. Research on several infrared bands of various molecular species and their isotopomers are undertaken. Those studies are briefly described.

Venkataraman, Malathy Devi↗
NASA NTRS1999

Interplanetary Fast Shocks and Associated Drivers Observed through the Twenty-Third Solar Minimum by WIND Over its First 2.5 Years

A list of the interplanetary (IP) shocks observed by WIND from its launch (in November 1994) to May 1997 is presented. Forty two shocks were identified. The magnetohydrodynamic nature of the shocks is investigated, and the associated shock parameters and their uncertainties are accurately computed using a practical scheme which combines two techniques. These techniques are a combination of the "pre-averaged" magnetic-coplanarity, velocity-coplanarity, and the Abraham-Schrauner-mixed methods, on the one hand, and the Vinas and Scudder [1986] technique for solving the non-linear least-squares Rankine-Hugoniot shock equations, on the other. Within acceptable limits these two techniques generally gave the same results, with some exceptions. The reasons for the exceptions are discussed. It is found that the mean strength and rate of occurrence of the shocks appears to correlated with the solar cycle. Both showed a decrease in 1996 coincident with the time of the lowest ultraviolet solar radiance, indicative of solar minimum and start of solar cycle 23, which began around June 1996. Eighteen shocks appeared to be associated with corotating interaction regions (CIRs). The distribution of their shock normals showed a mean direction peaking in the ecliptic plane and with a longitude (phi(sub n)) in that plane between perpendicular to the Parker spiral and radial from the Sun. When grouped according to the sense of the direction of propagation of the shocks the mean azimuthal (longitude) angle in GSE coordinates was approximately 194 deg for the fast-forward and approximately 20 deg for the fast-reverse shocks. Another 16 shocks were determined to be driven by solar transients, including magnetic clouds. These shocks had a broader distribution of normal directions than those of the CIR cases with a mean direction close to the Sun-Earth line. Eight shocks of unknown origin had normal orientation well off the ecliptic plane. No shock propagated with longitude phi(sub n) >= 220 +/- 10 deg, this would suggest strong hindrance to the propagation of shocks contra a rather tightly winding Parker spiral. Examination of the obliquity angle theta(sub Bn) (that between the shock normal and the upstream interplanetary magnetic field) for the full set of shocks revealed that about 58% was quasi-perpendicular, and some were very nearly perpendicular. About 32% of the shocks were oblique, and the rest (only 10%) were quasi-parallel, with one on Dec. 9, 1996 that showed field pulsations. Small uncertainty in the estimated angle theta(sub Bn) was obtained for about 10 shocks with magnetosonic Mach numbers between 1 and 2, hopefully significantly contributing to studies researching particle acceleration mechanisms at IP shocks, and to investigations where accurate values of theta(sub Bn) are crucial.

Mariani, F.↗
DOE OSTI2025

Chemometrics and visible diffuse reflectance spectroscopy to classify plutonium dioxide

Diffuse reflectance (DR) spectra in the Vis-NIR (∼380–1050 nm) region were acquired for a series of PuO 2 samples with a spot size of about 10 × 10 μm. Two batches of six PuO 2 samples, synthesized approximately 7.5 months apart, were prepared using both Pu(III) and Pu(IV) oxalate precursors at three distinct calcination temperatures (450, 650, and 950 °C). This yielded a total of 12 PuO 2 samples and 433 DR spectra. The DR spectrum of PuO 2 contained numerous peaks in the visible region, and characteristic features were identified with respect to calcination temperature and chemistry. A distinct peak multiplet near 615 nm was observed for samples prepared at low calcination temperatures, and a peak near 660 nm was observed for higher calcination temperatures. A multivariate classification strategy based on principal component analysis (PCA) was developed to distinguish PuO 2 calcination temperatures of 450, 650, and 950 °C with 100 % accuracy. Classification results also indicate the potential to distinguish chemical processing history (i.e., Pu(III) or Pu(IV)) based on the spectra with 72 % accuracy based on k-nearest neighbors applied to the PCA scores. Partial least squares discriminant analysis was used to identify variation among batches with 88 % accuracy and found that peaks near 669, 681, 811, and 970 nm were the most useful for predicting the batch identity. Here, this work demonstrates how micro-diffuse reflectance spectroscopy and chemometrics can be used to classify PuO 2 processing history based on Vis-NIR spectral features. Combining the chemometric approach with mapping sequences could provide a rapid, nondestructive approach to classify Pu oxide materials for environmental, forensics, and nonproliferation applications.

Actinide↗
NASA NTRS1984

Center of mass estimation in closed vortices - A verification in principle and practice

The problem of tracking closed mesoscale vortices using center of mass estimation techniques is studied. Three estimators are evaluated using data from a warm core Gulf Stream ring. The comparisons show that a method based on the intersection of perpendicular bisectors and one using a least-squares fit of a conic section perform comparably. The perpendicular bisector algorithm is used in conjunction with a Gaussian ring model and a star-shaped survey pattern to produce an expected error curve as a function of vortex translation, survey speed and vortex size. For typical ring parameters, center estimation is usually possible to within + or - 5 km. The feasibility of using differing data sets to construct a history of ring motion based on a coordinate system moving with the ring is also investigated. In this way, the validity of using satellite-derived data and drifter trajectories to estimate the center of mass of a mesoscale feature is assessed. The results of the analysis demonstrate that the location of the deeper structure of the ring and the surface expression are sufficiently well correlated to permit dynamically relevant calculations based on surface measurements. It is shown that satellite-derived data can be used to approximate the center of mass trajectory to within the error in the individual center estimates for the period analyzed. The Lagrangian-drifter-derived centers are offset from the center of mass trajectory in a manner consistent with kinematic arguments.

Hooker, S. B.↗
NASA NTRS2013

The NuSTAR Extragalactic Survey: A First Sensitive Look at the High-Energy Cosmic X-Ray Background Population

We report on the first 10 identifications of sources serendipitously detected by the Nuclear Spectroscopic Telescope Array (NuSTAR) to provide the first sensitive census of the cosmic X-ray background source population at approximately greater than 10 keV. We find that these NuSTAR-detected sources are approximately 100 times fainter than those previously detected at approximately greater than 10 keV and have a broad range in redshift and luminosity (z = 0.020-2.923 and L(sub 10-40 keV) approximately equals 4 × 10(exp 41) - 5 × 10(exp 45) erg per second; the median redshift and luminosity are z approximately equal to 0.7 and L(sub 10-40 keV) approximately equal to 3 × 10(exp 44) erg per second, respectively. We characterize these sources on the basis of broad-band approximately equal to 0.5 - 32 keV spectroscopy, optical spectroscopy, and broad-band ultraviolet-to-mid-infrared spectral energy distribution analyses. We find that the dominant source population is quasars with L(sub 10-40 keV) greater than 10(exp 44) erg per second, of which approximately 50% are obscured with N(sub H) approximately greater than 10(exp 22) per square centimeters. However, none of the 10 NuSTAR sources are Compton thick (N(sub H) approximately greater than 10(exp 24) per square centimeters) and we place a 90% confidence upper limit on the fraction of Compton-thick quasars (L(sub 10-40 keV) greater than 10(exp 44) erg per second) selected at approximately greater than 10 keV of approximately less than 33% over the redshift range z = 0.5 - 1.1. We jointly fitted the rest-frame approximately equal to 10-40 keV data for all of the non-beamed sources with L(sub 10-40 keV) greater than 10(exp 43) erg per second to constrain the average strength of reflection; we find R less than 1.4 for gamma = 1.8, broadly consistent with that found for local active galactic nuclei (AGNs) observed at approximately greater than 10 keV. We also constrain the host-galaxy masses and find a median stellar mass of approximately 10(exp 11) solar mass, a factor approximately 5 times higher than the median stellar mass of nearby high-energy selected AGNs, which may be at least partially driven by the order of magnitude higher X-ray luminosities of the NuSTAR sources. Within the low source-statistic limitations of our study, our results suggest that the overall properties of the NuSTAR sources are broadly similar to those of nearby high-energy selected AGNs but scaled up in luminosity and mass.

Spectroscopic↗
NASA NTRS2013

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗
NASA NTRS2014

Grail Refinements to Lunar Seismic Structure

To probe a planet's interior, seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL (Gravity Recovery and Interior Laboratory) mission's high-resolution measurements of the lunar gravity field provide constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten). These data complement seismic investigations, and joint interpretation permits improved constraints on the Moon's internal structure. Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗
NASA NTRS1991

The control network of Mars: April 1991

The modern geodetic control network of Mars was first established based on Mariner 9 images with 1-2 km/pixel resolutions and covered almost the entire Martian surface. The introduction of higher resolution (10-200 meter/pixel) Viking Orbiter images greatly improved the accuracy and density of points in the control network. Analysis of the Viking Lander radio tracking data led to more accurate measurements of Mars' rotation period, spin axis direction, and the lander coordinates relative to the inertial reference frame. The prime meridian on Mars was defined by the Geodesy/Cartography Group of the Mariner 9 Television Team as the crater Airy-0, located about 5 degrees south of the equator. The Viking 1 Lander site was identified on a high resolution Viking frame. The control point measurements form the basis of a least squares solution determined by analytical triangulation after the pixel measurements are corrected for geometric distortions and converted to millimeter coordinates in the camera focal plane. Photogrammetric strips encircling Mars at the equator and at 60 degree north south were used to strengthen the overall net and improve the accuracy of the coordinates of points. In addition, photogrammetric strips along 0, 90, 180, and 270 degrees longitude to the Viking 1 Lander site have all significantly strengthened the control network. Most recently, photogrammetric strips were added to the net along 30 degrees north latitude between 0 and 180 degrees, and along 30 degrees between 180 and 360 degrees. The Viking 1 Lander site and Airy-0 are linked through photogrammetric strips occurring along the 0 degree meridian from Airy-0 to 65 degrees north, from that point through the Viking 1 Lander site to the equator, and along the equator to 180 degrees longitude. The Viking 1 lander site is thus a well calibrated area with coordinates of points accurate to approximately 200 meters relative to the J2000 inertial coordinate system. This will be a useful calibration region for upcoming missions. The current status of the control network calculations is presented.

Davies, Merton E.↗
NASA NTRS1951

Wind-tunnel investigation of six shielded total-pressure tubes at high angles of attack subsonic speeds

The effect of inclination of the air stream on the measured pressures of six shielded total-pressure tubes (all designed for end-mounting on a horizontal boom) has been determined for an angle-of-attack range of approximately -30° to 65°. The tests were conducted at a Mach number of 0.26 in the Langley stability tunnel and at Mach numbers of 0.50, 0.70, 0.90, and 0.95 in the Langley high-speed 7- by 10-foot tunnel. The results of the tests showed that curved venturi entries are less sensitive to inclination of the air stream than the conical entry of the standard Kiel design. Of the three curved entries tested, the tube having the bluntest entry proved the least sensitive to inclination. The critical angle (that is, the sngle.at which the total-pressure error equals 1 percent of the impact pressure) for this tube was ±63.2° at a Mach number of 0.26. The tests also showed that the critical angle of a shielded tube could be extended to higher positive angles of attack by means of a slant profile. The critical angle of a tube having a conical entry with a 10° slant profile, for example, was found to average 9° greater than that of a similar tube having a square profile. Although the critical angle of this tube at angles of yaw was less than that at angles of attack, the critical angle in yaw was equal to that of the tube with the square profile. Tests of the effect of varying the position of the total-pressure probe in the shield showed that, in general, no advantage was to be gained by changing from the position used in the standard Kiel design. Limited tests of the effect of varying the internal diameter of the shield showed that an increase in the diameter of the throat corresponding an increase in the internal area of 50 percent resulted in a relatively small increase (2.5°) in the critical angle. The effect of Mach number on the sensitivity of the tubes to inclination was shown to be appreciable. An increase in Mach number from 0.26 to 0.90, for example, resulted in a decrease in the critical angle of most of the tubes of about 4°.

Walter R Russell↗
NASA NTRS1942

Radiator Design and Installation - II, Special Report

A mathematical analysis of radiator design has been made. The volume of the radiator using least total power has been expressed in a single formula which shows that the optimum radiator volume is independent of the shape of the radiator and which makes possible the construction of design tables that give the optimum radiator volume per 100-horsepower heat dissipation as a function of the speed, of the altitude, and of one parameter involving characteristics of the airplane. Although, for a given set of conditions, the radiator volume using the least total power is fixed, the frontal area, or the length of the radiator needs to be separately specified in order to satisfy certain other requirement such as the ability to cool with the pressure drop available while the airplane is climbing. In order to simplify the specification for the shape of the radiator and in order to reduce the labor involved in calculating the detailed performance of radiators, generalized design curves have been developed for determining the pressure drop, the mass flow of air, and the power expended in overcoming the cooling drag of a radiator from the physical dimensions of the radiator. In addition, a table is derived from these curves, which directly gives the square root of the pressure drop required for ground cooling as a function of the radiator dimensions, of the heat dissipation and of the available temperature difference. Typical calculations using the tables of optimum radiator volume and the design curves are given. The jet power that can be derived from the heated air is proportional to the heat dissipation and is approximately proportional to the square of the airplane speed and to the reciprocal of the absolute temperature of the atmosphere. A table of jet power, per 100 horsepower of heat dissipation at various airplane speeds and altitudes is presented.

Tifford, Arthur N.↗