Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vector fitting”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Binning is Sinning (Supernova Version): The Impact of Self-calibration in Cosmological Analyses with Type Ia Supernovae

Recent cosmological analyses (e.g., JLA, Pantheon) of Type Ia supernovae (SNe Ia) have propagated systematic uncertainties into a covariance matrix and either binned or smoothed the systematic uncertainty vectors in redshift space. We demonstrate that systematic error budgets of these analyses can be improved by a factor of ~ 1.5 × with the use of unbinned and unsmoothed covariance matrices. To understand this, we employ a separate approach that simultaneously fits for cosmological parameters and additional self-calibrating scale parameters that constrain the size of each systematic. We show that the covariance-matrix approach and scale-parameter approach indeed yield equivalent results, implying that in both cases the data can self-calibrate certain systematic uncertainties, but that this ability is hindered when information is binned or smoothed in redshift space. We review the top systematic uncertainties in current analyses and find that the reduction of systematic uncertainties in the unbinned case depends on whether a systematic is solely degenerate with the cosmological model in redshift space or whether it can be described by additional correlations between supernova properties and luminosity. Furthermore, we show that the power of self-calibration increases with the size of the data set, which presents a tremendous opportunity for upcoming analyses of photometrically classified samples, like those of Legacy Survey of Space and Time (LSST) and the Nancy Grace Roman Telescope (NGRST). However, to take advantage of self-calibration in large, photometrically classified samples, we must first address the issue that binning is required in currently used photometric analysis methods.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2008

Visual Data Analysis for Satellites

The Visual Data Analysis Package is a collection of programs and scripts that facilitate visual analysis of data available from NASA and NOAA satellites, as well as dropsonde, buoy, and conventional in-situ observations. The package features utilities for data extraction, data quality control, statistical analysis, and data visualization. The Hierarchical Data Format (HDF) satellite data extraction routines from NASA's Jet Propulsion Laboratory were customized for specific spatial coverage and file input/output. Statistical analysis includes the calculation of the relative error, the absolute error, and the root mean square error. Other capabilities include curve fitting through the data points to fill in missing data points between satellite passes or where clouds obscure satellite data. For data visualization, the software provides customizable Generic Mapping Tool (GMT) scripts to generate difference maps, scatter plots, line plots, vector plots, histograms, timeseries, and color fill images.

Lau, Yee↗
DOE OSTI2025

Constraining Absolute Higgs Couplings using Vector Boson Fusion H to WW with FCC-hh

The absolute normalisation scale of Higgs couplings is difficult to determine in hadron colliders. An experimentally and theoretically clean way to establish the coupling scale is to measure the ratio of resonant Vector Boson Fusion $H \rightarrow WW^{*}$ production to non-resonant $qq \rightarrow qqWW$ because both the theoretical and detector response systematics all cancel to first order in the ratio. To set the scale of the sensitivity, a simple selection and fit is presented with gives a 1.8\% (1.9\%) statistical uncertainty on the ratio for a 100 (84) TeV FCC-hh at 30\,ab$^{-1}$, which would give a corresponding uncertainty on the $\kappa$ scale of 0.5-1\% depending on systematic effects and potential impacts of more sophisticated analysis.

Lipeles, Elliot [University of Pennsylvania] (ORCI↗
NASA NTRS2012

Comparison of PARASOL Observations with Polarized Reflectances Simulated Using Different Ice Habit Mixtures

Insufficient knowledge of the habit distribution and the degree of surface roughness of ice crystals within ice clouds is a source of uncertainty in the forward light scattering and radiative transfer simulations required in downstream applications involving these clouds. The widely used MODerate Resolution Imaging Spectroradiometer (MODIS) Collection 5 ice microphysical model assumes a mixture of various ice crystal shapes with smooth-facets except aggregates of columns for which a moderately rough condition is assumed. When compared with PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) polarized reflection data, simulations of polarized reflectance using smooth particles show a poor fit to the measurements, whereas very rough-faceted particles provide an improved fit to the polarized reflectance. In this study a new microphysical model based on a mixture of 9 different ice crystal habits with severely roughened facets is developed. Simulated polarized reflectance using the new ice habit distribution is calculated using a vector adding-doubling radiative transfer model, and the simulations closely agree with the polarized reflectance observed by PARASOL. The new general habit mixture is also tested using a spherical albedo differences analysis, and surface roughening is found to improve the consistency of multi-angular observations. It is suggested that an ice model incorporating an ensemble of different habits with severely roughened surfaces would potentially be an adequate choice for global ice cloud retrievals.

Cole, Benjamin H.↗
NASA NTRS2016

Modeling in the State Flow Environment to Support Launch Vehicle Verification Testing for Mission and Fault Management Algorithms in the NASA Space Launch System

Analysis methods and testing processes are essential activities in the engineering development and verification of the National Aeronautics and Space Administration's (NASA) new Space Launch System (SLS). Central to mission success is reliable verification of the Mission and Fault Management (M&FM) algorithms for the SLS launch vehicle (LV) flight software. This is particularly difficult because M&FM algorithms integrate and operate LV subsystems, which consist of diverse forms of hardware and software themselves, with equally diverse integration from the engineering disciplines of LV subsystems. M&FM operation of SLS requires a changing mix of LV automation. During pre-launch the LV is primarily operated by the Kennedy Space Center (KSC) Ground Systems Development and Operations (GSDO) organization with some LV automation of time-critical functions, and much more autonomous LV operations during ascent that have crucial interactions with the Orion crew capsule, its astronauts, and with mission controllers at the Johnson Space Center. M&FM algorithms must perform all nominal mission commanding via the flight computer to control LV states from pre-launch through disposal and also address failure conditions by initiating autonomous or commanded aborts (crew capsule escape from the failing LV), redundancy management of failing subsystems and components, and safing actions to reduce or prevent threats to ground systems and crew. To address the criticality of the verification testing of these algorithms, the NASA M&FM team has utilized the State Flow environment6 (SFE) with its existing Vehicle Management End-to-End Testbed (VMET) platform which also hosts vendor-supplied physics-based LV subsystem models. The human-derived M&FM algorithms are designed and vetted in Integrated Development Teams composed of design and development disciplines such as Systems Engineering, Flight Software (FSW), Safety and Mission Assurance (S&MA) and major subsystems and vehicle elements such as Main Propulsion Systems (MPS), boosters, avionics, Guidance, Navigation, and Control (GN&C), Thrust Vector Control (TVC), liquid engines, and the astronaut crew office. Since the algorithms are realized using model-based engineering (MBE) methods from a hybrid of the Unified Modeling Language (UML) and Systems Modeling Language (SysML), SFE methods are a natural fit to provide an in depth analysis of the interactive behavior of these algorithms with the SLS LV subsystem models. For this, the M&FM algorithms and the SLS LV subsystem models are modeled using constructs provided by Matlab which also enables modeling of the accompanying interfaces providing greater flexibility for integrated testing and analysis, which helps forecast expected behavior in forward VMET integrated testing activities. In VMET, the M&FM algorithms are prototyped and implemented using the same C++ programming language and similar state machine architectural concepts used by the FSW group. Due to the interactive complexity of the algorithms, VMET testing thus far has verified all the individual M&FM subsystem algorithms with select subsystem vendor models but is steadily progressing to assessing the interactive behavior of these algorithms with LV subsystems, as represented by subsystem models. The novel SFE applications has proven to be useful for quick look analysis into early integrated system behavior and assessment of the M&FM algorithms with the modeled LV subsystems. This early MBE analysis generates vital insight into the integrated system behaviors, algorithm sensitivities, design issues, and has aided in the debugging of the M&FM algorithms well before full testing can begin in more expensive, higher fidelity but more arduous environments such as VMET, FSW testing, and the Systems Integration Lab7 (SIL). SFE has exhibited both expected and unexpected behaviors in nominal and off nominal test cases prior to full VMET testing. In many findings, these behavioral characteristics were used to correct the M&FM algorithms, enable better test coverage, and develop more effective test cases for each of the LV subsystems. This has improved the fidelity of testing and planning for the next generation of M&FM algorithms as the SLS program evolves from non-crewed to crewed flight, impacting subsystem configurations and the M&FM algorithms that control them. SFE analysis has improved robustness and reliability of the M&FM algorithms by revealing implementation errors and documentation inconsistencies. It is also improving planning efficiency for future VMET testing of the M&FM algorithms hosted in the LV flight computers, further reducing risk for the SLS launch infrastructure, the SLS LV, and most importantly the crew.

Trevino, Luis↗
NASA NTRS1994

Root-sum-square structural strength verification approach

Utilizing a proposed fixture design or some variation thereof, this report presents a verification approach to strength test space flight payload components, electronics boxes, mechanisms, lines, fittings, etc., which traditionally do not lend themselves to classical static loading. The fixture, through use of ordered Euler rotation angles derived herein, can be mounted on existing vibration shakers and can provide an innovative method of applying single axis flight load vectors. The versatile fixture effectively loads protoflight or prototype components in all three axes simultaneously by use of a sinusoidal burst of desired magnitude at less than one-third the first resonant frequency. Cost savings along with improved hardware confidence are shown. The end product is an efficient way to verify experiment hardware for both random vibration and strength.

Lee, Henry M.↗
NASA NTRS1997

Kinematic principles of primate rotational vestibulo-ocular reflex. II. Gravity-dependent modulation of primary eye position

The kinematic constraints of three-dimensional eye positions were investigated in rhesus monkeys during passive head and body rotations relative to gravity. We studied fast and slow phase components of the vestibulo-ocular reflex (VOR) elicited by constant-velocity yaw rotations and sinusoidal oscillations about an earth-horizontal axis. We found that the spatial orientation of both fast and slow phase eye positions could be described locally by a planar surface with torsional variation of <2.0 +/- 0.4 degrees (displacement planes) that systematically rotated and/or shifted relative to Listing's plane. In supine/prone positions, displacement planes pitched forward/backward; in left/right ear-down positions, displacement planes were parallel shifted along the positive/negative torsional axis. Dynamically changing primary eye positions were computed from displacement planes. Torsional and vertical components of primary eye position modulated as a sinusoidal function of head orientation in space. The torsional component was maximal in ear-down positions and approximately zero in supine/prone orientations. The opposite was observed for the vertical component. Modulation of the horizontal component of primary eye position exhibited a more complex dependence. In contrast to the torsional component, which was relatively independent of rotational speed, modulation of the vertical and horizontal components of primary position depended strongly on the speed of head rotation (i.e., on the frequency of oscillation of the gravity vector component): the faster the head rotated relative to gravity, the larger was the modulation. Corresponding results were obtained when a model based on a sinusoidal dependence of instantaneous displacement planes (and primary eye position) on head orientation relative to gravity was fitted to VOR fast phase positions. When VOR fast phase positions were expressed relative to primary eye position estimated from the model fits, they were confined approximately to a single plane with a small torsional standard deviation ( approximately 1.4-2.6 degrees). This reduced torsional variation was in contrast to the large torsional spread (well >10-15 degrees ) of fast phase positions when expressed relative to Listing's plane. We conclude that primary eye position depends dynamically on head orientation relative to space rather than being fixed to the head. It defines a gravity-dependent coordinate system relative to which the torsional variability of eye positions is minimized even when the head is moved passively and vestibulo-ocular reflexes are evoked. In this general sense, Listing's law is preserved with respect to an otolith-controlled reference system that is defined dynamically by gravity.

NASA Discipline Neuroscience↗
DOE OSTI2023

How to identify different new neutrino oscillation physics scenarios at DUNE

Next generation neutrino oscillation experiments are expected to measure the remaining oscillation parameters with very good precision. They will have unprecedented capabilities to search for new physics that modify oscillations. DUNE, with its broad band beam, good particle identification, and relatively high energies will provide an excellent environment to search for new physics. If deviations from the standard three-flavor oscillation picture are seen however, it is crucial to know which new physics scenario is found so that it can be verified elsewhere and theoretically understood. We investigate several benchmark new physics scenarios by looking at existing long-baseline accelerator neutrino data from NOvA and T2K and determine at what sensitivity DUNE can differentiate among them. We consider sterile neutrinos and both vector and scalar non-standard neutrino interactions, all with new complex phases, the latter of which could conceivably provide absolute neutrino mass scale information. We find that, in many interesting cases, DUNE will have good model discrimination. We also perform a new fit to NOvA and T2K data with scalar NSI.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Noninteracting electrons in a prototypical one-dimensional sinusoidal potential

A prototypical model of a one-dimensional metallic monatomic solid containing noninteracting electrons is studied, where the argument of the cosine potential energy, periodic with the lattice, contains the first reciprocal lattice vector G 1 =2π/a, where a is the lattice constant. The time-independent Schrödinger equation can be written in reduced variables as a Mathieu equation for which numerically exact solutions for the band structure and wave functions are obtained. The band structure has band gaps that increase with increasing amplitude q of the cosine potential. In the extended-zone scheme, the energy gaps decrease with increasing index n of the Brillouin-zone boundary ka=nπ, where k is the crystal momentum of the electron. The wave functions of the band electron are derived for various combinations of k and q as complex combinations of the real Mathieu functions with even and odd parity, and the normalization factor is discussed. The wave functions at the bottoms and tops of the bands are found to be real or imaginary, respectively, corresponding to standing waves at these energies. Irrespective of the wave vector k within the first Brillouin zone, the electron probability density is found to be periodic with the lattice. The Fourier components of the wave functions are derived versus q, which reveal multiple reciprocal-lattice-vector components with variable amplitudes in the wave functions unless q = 0. The magnitudes of the Fourier components are found to decrease exponentially as a power of n for n~3 to 45 for ka = π/2 and q = 2, and a precise fit is obtained to the data. The probability densities and probability currents obtained from the wave functions are also discussed. The probability currents are found to be zero for crystal momenta at the tops and bottoms of the energy bands, because the wave functions for these crystal momenta are standing waves. To conlcude, the band structure is calculated from the central equation and compared to the numerically exact band structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1973

Coordinated observations of chemical releases from the ground and from aircraft at high latitudes

The ground observations of the Na-Li trail released from a Nike-Apache rocket obtained by the Geophysical Institute are discussed. By using the nominal trajectory for a 60 pound payload and the particular rocket, a best fit trajectory was determined based on the Ester Dome photographic data, launch time and earth-sun geometrical shadow height. From these calculations, the height of obvious features along the trail were determined and their velocity estimated. A clockwise rotation of the velocity vector with increasing height was observed. Velocities deduced at various altitudes were then compared to meter radar data also obtained during this period. The comparisons of these two neutral wind measurements techniques are satisfactory.

Romick, G. J.↗
NASA NTRS1996

The alpha Centauri Line of Sight: D/H Ratio, Physical Properties of Local Interstellar Gas, and Measurement of Heated Hydrogen (The 'Hydrogen Wall') Near the Heliopause

We analyze high-resolution spectra of the nearby (1.34 pc) stars alpha Cen A (G2 V) and alpha Cen B (K1 V), which were obtained with the Goddard High Resolution Spectrograph on the Hubble Space Telescope. The observations consist of echelle spectra of the Mg II 2800 A and Fe II 2599 A resonance lines and the Lyman-alpha lines of hydrogen and deuterium. The interstellar gas has a velocity (v = - 18.0 +/- 0.2 km/s) consistent with the local flow vector proposed for this line of sight by Lailement & Berlin (1992). The temperature and nonthermal velocity inferred from the Fe II, Mg II, and D I line profiles are T = 5400 +/- 500 K and xi = 1.20 +/- 0.25 km/s, respectively. However, single-component fits to the H I Lyman-alpha lines yield a Doppler parameter (b(sub HI) = 11.80 km/s) that implies a significantly warmer temperature of 8350 K, and the velocity of the H I absorption (v = - 15.8 +/- 0.2 km/s) is redshifted by about 2.2 km/s with respect to the Fe II, Mg II, and D I lines. The one-component model of the interstellar gas suggests natural logarithm N base HI = 18.03 +/- 0.01 and D/H = (5.7 +/- 0.2) x 10(exp -6) . These parameters lead to a good fit to the observed spectra, but this model does not explain the higher temperature and redshift of H I relative to the other interstellar lines. The most sensible way to resolve the discrepancy between H(I) and the other lines is to add a second absorption component to the H(I) lines. This second component is hotter (T approx. equals 30,000 K), is redshifted relative to the primary component by 2-4 km/s, and has a column density too low to be detected in the Fe(II), Mg(II), and D(I) lines. We propose that the gas responsible for this component is located near the heliopause, consisting of the heated H I gas from the interstellar medium that is compressed by the solar wind. This so-called 'hydrogen wall' is predicted by recent multifluid gasdynamical models of the interstellar gas and solar wind interaction. Our data provide the first measurements of the temperature and column density of H(I) in the hydrogen wall. After considering the effects that a corresponding hydrogen wall around alpha Cen would have on our analysis, our best estimates for the parameters of the solar hydrogen wall are natural log N(sup (2))(H(I)) = 14.74 +/- 0.24, b(sup (2))(H(I)) = 21.9 +/- 1.7 km/s (corresponding to T = 29,000 +/- 5000 K), and v(sup (2))(H(I)) greater than -16km/s. Unfortunately, the existence of this heated H(I) reduces our ability to compute the H(I) column density of the interstellar medium accurately because, with slight alterations to our assumed stellar Lyman-alpha profiles, we discovered that acceptable two-component fits also exist with natural log N(H(I))approx. 17.6. We, therefore, quote large error bars for the H I column density along the alpha Cen line of sight, natural log N(H(I)) = 17.80 +/- 0.30. For this range in N(H(I)), n(H(I)) = 0.15 /cu.cm (+/- a factor of 2) and D/H = (0.5-1.9) x 10(exp -5). This is the first direct measurement of the H(I) density in a local cloud and allows us to predict the distance from the Sun to the edge of the local cloud along various lines of sight. This range in D/H is consistent with the value D/H = 1.6 x 10(exp -5) previously derived for the Capella and Procyon lines of sight. We cannot tell whether D/H ratio varies or is constant in the local interstellar medium, but we do find that the D(I)/Mg(II) ratio for the alpha Cen line of sight is about 4 times smaller than for the Capella and Procyon lines of sight. Therefore, either D/H or the Mg depletion varies significantly over distance scales of only a few parsecs.

Linsky, Jeffrey L.↗
NASA NTRS1999

Novel, Miniature Multi-Hole Probes and High-Accuracy Calibration Algorithms for their use in Compressible Flowfields

Two new calibration algorithms were developed for the calibration of non-nulling multi-hole probes in compressible, subsonic flowfields. The reduction algorithms are robust and able to reduce data from any multi-hole probe inserted into any subsonic flowfield to generate very accurate predictions of the velocity vector, flow direction, total pressure and static pressure. One of the algorithms PROBENET is based on the theory of neural networks, while the other is of a more conventional nature (polynomial approximation technique) and introduces a novel idea of local least-squares fits. Both algorithms have been developed to complete, user-friendly software packages. New technology was developed for the fabrication of miniature multi-hole probes, with probe tip diameters all the way down to 0.035". Several miniature 5- and 7-hole probes, with different probe tip geometries (hemispherical, conical, faceted) and different overall shapes (straight, cobra, elbow probes) were fabricated, calibrated and tested. Emphasis was placed on the development of four stainless-steel conical 7-hole probes, 1/16" in diameter calibrated at NASA Langley for the entire subsonic regime. The developed calibration algorithms were extensively tested with these probes demonstrating excellent prediction capabilities. The probes were used in the "trap wing" wind tunnel tests in the 14'x22' wind tunnel at NASA Langley, providing valuable information on the flowfield over the wing. This report is organized in the following fashion. It consists of a "Technical Achievements" section that summarizes the major achievements, followed by an assembly of journal articles that were produced from this project and ends with two manuals for the two probe calibration algorithms developed.

Rediniotis, Othon K.↗
DOE OSTI2025

Small-𝑥 gluon GPD constrained from deeply virtual 𝐽/𝜓 production and gluon PDF through universal-moment parametrization

We phenomenologically constrain the small-𝑥 and small-𝜉 gluon generalized parton distributions (GPDs) with the deeply virtual 𝐽/𝜓 production (DV⁢𝐽/𝜓⁢P) in the framework of GPDs through universal moment parameterization (GUMP). We use a hybrid cross-section formula combining collinear factorization to the next-to-leading order (NLO) accuracy of the strong coupling 𝛼 𝑠 , with corrections from nonrelativistic QCD to account for the power corrections due to the heavy 𝐽/𝜓 mass. We reach reasonable fit to the measured differential cross sections of DV⁢𝐽/𝜓⁢P by H1 at the Hadron-Electron Ring Accelerator as well as forward gluon parton distribution functions (PDFs) from JAM22 global analysis. We find that both NLO and nonrelativistic corrections are significant for heavy vector meson productions. Of course, the gluon GPD we obtain still contains considerable freedom in need of inputs from other constraints, particularly in the distribution-amplitude-like region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Search for Higgs bosons produced in association with a high-energy photon via vector-boson fusion and decaying to a pair of b -quarks in the ATLAS detector

A search for Standard Model Higgs bosons produced in association with a high-energy photon and decaying to $b\bar{b}$ is performed using 133 fb −1 of $\sqrt{s}$ = 13 TeV pp collision data collected with the ATLAS detector at the Large Hadron Collider at CERN. The photon requirement reduces the multijet background, and the H → $b\bar{b}$ decay is the dominant decay mode. Event selection requirements target events produced by vector-boson fusion, the dominant production mode in this channel. Several improvements enhance the search sensitivity compared to previous measurements. These improvements include better background modelling and characterization, the use of a neural-network classifier, and an updated signal extraction strategy adopting a direct binned-likelihood fit to the classifier output. With these improvements, the Higgs boson signal strength is measured to be 0.2 ± 0.7 relative to the Standard Model prediction. This corresponds to an observed significance of 0.3 standard deviations, compared to an expectation of 1.5 standard deviations assuming the Standard Model.

ATLAS experiment↗
DOE OSTI2025

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↗
NASA NTRS1978

Numerical solution of the Navier-Stokes equations for arbitrary blunt bodies in supersonic flows

A time-dependent, two-dimensional Navier-Stokes code employing the body-fitted coordinate technique has been developed for supersonic flows past blunt bodies of arbitrary shape. The computer program is based on the finite-difference approximation of the compressible Navier-Stokes equations transformed to nonorthogonal curvilinear coordinates with the contravariant components of the velocity vector as dependent variables. The bow shock ahead of the body is obtained as part of the solution, by 'shock capturing'. Numerical solutions of the complete equations are presented in detail for free-stream Mach number 4.6, Reynolds number 10,000, and an isothermal wall temperature of 556 K for a circular cylinder with the free-stream outer boundaries forming a hyperbola in the front and a circular arc in the back.

Warsi, Z. U. A.↗
NASA NTRS1980

Design of a recursive vector processor using polynomial splines

The problem of obtaining smoothed estimates of function values, particularly their derivatives, from a finite set of inaccurate measurements is considered. A recursive two-dimensional vector processor is introduced as an approximation to the nonrecursive constrained least-squares estimation. Here, piecewise bicubic Hermite polynomials are extensively used as approximating functions, and the smoothing integral is converted to a discrete quadratic form. This makes it possible to convert the problem of fitting an approximating function to one of estimating the function values and derivatives at the nodes.

Kim, C. S.↗
NASA NTRS1983

Intensity dependence of the pulse profile and polarization of the Vela pulsar

Polarization measurements of the Vela pulsar PSR 0833-45 at 2295 MHz revealed that the radio pulse profile is composed of four components. The four component configuration was found to fit the observed dependence of the pulse profile on the intensity of the pulsar, and the dependence of the modulation index and power spectrum of the intensity fluctuations on the pulse-longitude. Deviations of the observed polarization position angles from those predicted by a rotating vector model were explained by the emission regions for each component being located at different distances from the star. The components are affected by retardation, aberration, and magnetic field curvature effects, all of which contribute to arrival of the components at different longitudes even though the components actually carry the position angle signature predicted by the rotating vector model.

Krishnamohan, S.↗