Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase errors”

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 1,297 records · Page 72

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗

Synthesis and Processing by Design of High-Nickel Cathode Materials

The high demand of lightweight, high energy density batteries for energy storage promotes new materials discovery and development. Despite the large number of battery materials being discovered, very few of them have been commercially deployed, mostly bottlenecked by synthesis and processing – namely, making certain phases with the desired structure and properties to meet the multifaceted performance requirements. Alternative to the traditional trial and error, we present here an in situ study aided synthesis- and processing-by-design approach. With specific examples, we illustrate how to use the approach to identify reaction pathways in synthesis and processing of high-Nickel (Ni) cathode materials for next-generation lithium-ion batteries, thereby ensuring precise control of their structure, morphology, and surface properties. Furthermore, perspectives are provided on the wide applicability of the approach to solving critical issues inherent to high-Ni cathodes and the new directions and opportunities in the area.

36 MATERIALS SCIENCE↗

Towards Quantum Computing Phase Diagrams of Gauge Theories with Thermal Pure Quantum States

The phase diagram of strong interactions in nature at finite temperature and chemical potential remains largely theoretically unexplored due to inadequacy of Monte-Carlo–based computational techniques in overcoming a sign problem. Quantum computing offers a sign-problem-free approach, but evaluating thermal expectation values is generally resource intensive on quantum computers. To facilitate thermodynamic studies of gauge theories, we propose a generalization of the thermal-pure-quantum-state formulation of statistical mechanics applied to constrained gauge-theory dynamics, and numerically demonstrate that the phase diagram of a simple low-dimensional gauge theory is robustly determined using this approach, including mapping a chiral phase transition in the model at finite temperature and chemical potential. Quantum algorithms, resource requirements, and algorithmic and hardware error analysis are further discussed to motivate future implementations. Thermal pure quantum states, therefore, may present a suitable candidate for efficient thermal simulations of gauge theories in the era of quantum computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A0535+26 - Refined position measurement and new pulse period data

The hard, pulsing, transient X-ray source A0535+26 has been observed with SAS 3 on three occasions during 1977-1978. These observations have yielded a precise position measurement (20 arcsec error radius), which renders the identification of A0535+26 with the Be star HDE 245770 virtually certain. The pulse phase was tracked for about 9 days in April 1978 and clearly showed both first and second derivatives in the pulse period. An analysis of these timing data, combined with data from previous observations, leads to the following conclusions: (1) a significant fraction of the observed changes in pulse period is probably intrinsic to the compact X-ray star (e.g., accretion torques on a neutron star), and (2) conservative limits on binary orbital parameters tend to further confirm a long orbital period (at least about 20 days).

Li, F.↗

Terminal area automatic navigation, guidance, and control research using the Microwave Landing System (MLS). Part 2: RNAV/MLS transition problems for aircraft

The problems in navigation and guidance encountered by aircraft in the initial transition period in changing from distance measuring equipment, VORTAC, and barometric instruments to the more precise microwave landing system data type navaids in the terminal area are investigated. The effects of the resulting discontinuities on the estimates of position and velocity for both optimal (Kalman type navigation schemes) and fixed gain (complementary type) navigation filters, and the effects of the errors in cross track, track angle, and altitude on the guidance equation and control commands during the critical landing phase are discussed. A method is presented to remove the discontinuities from the navigation loop and to reconstruct an RNAV path designed to land the aircraft with minimal turns and altitude changes.

Pines, S.↗

Design and Performance of Sampled Data Loops for Subcarrier and Carrier Tracking

Design parameters and resulting performance are presented for the sampled data analogies of continuous time phase locked loops of second and third order containing perfect integrators. Expressions for noise equivalent bandwidth and steady state errors are given. Stability and gain margin are investigated using z plane root loci. Finally, an application is presented for Voyager subcarrier and carrier tracking under the dynamics of the encounters with Uranus and Neptune. For carrier tracking, loop bandwidth narrow enough for satisfactory loop signal to noise ratios can be achieved using third order loops without rate aiding, whereas second order loops would require aiding. For subcarrier tracking, third order loops can be used when the sampling rate is limited to approximately once per second, as in the Baseband Assembly, whereas second order loops sufficiently wide to track the dynamics have stability problems at that sampling rate.

Aguirre, S.↗

Method of resolving radio phase ambiguity in satellite orbit determination

For satellite orbit determination, the most accurate observable available today is microwave radio phase, which can be differenced between observing stations and between satellites to cancel both transmitter- and receiver-related errors. For maximum accuracy, the integer cycle ambiguities of the doubly differenced observations must be resolved. To perform this ambiguity resolution, a bootstrapping strategy is proposed. This strategy requires the tracking stations to have a wide ranging progression of spacings. By conventional 'integrated Doppler' processing of the observations from the most widely spaced stations, the orbits are determined well enough to permit resolution of the ambiguities for the most closely spaced stations. The resolution of these ambiguities reduces the uncertainty of the orbit determination enough to enable ambiguity resolution for more widely spaced stations, which further reduces the orbital uncertainty. In a test of this strategy with six tracking stations, both the formal and the true errors of determining Global Positioning System satellite orbits were reduced by a factor of 2.

Councelman, Charles C., III↗

A radio optical reference frame. I - Precise radio source positions determined by Mark III VLBI - Observations from 1979 to 1988 and a tie to the FK5

Observations from 600 Mark III VLBI experiments from 1979 to 1988, resulting in 237,681 acceptable pairs of group delay and phase delay rate observations, have been used to derive positions of 182 extragalactic radio sources with typical formal standard errors less than 1 mas. The sources are distributed fairly evenly above delta = -30 deg, and 70 sources have delta greater than 0 deg. Analysis with different troposphere models, as well as internal and external comparisons, indicates that a coordinate frame defined by this set of radio sources should be reliable at the 1 mas level. The right ascension zero point of this reference frame has been aligned with the FK5 by using the optical positions of 28 extragalactic radio sources whose positions are on the FK5 system. Because of known defects in the knowledge of astronomical constants, daily nutation offsets in longitude and obliquity were determined relative to an arbitrary reference day in the set of experiments.

Ma, C.↗

Doppler tracking of planetary spacecraft

This article concerns the measurement of Doppler shift on microwave links that connect planetary spacecraft with the Deep Space Network. Such measurements are made by tracking the Doppler effect with phase-locked loop receivers. A description of equipment and techniques as well as a summary of the appropriate mathematical models are given. The two-way Doppler shift is measured by transmitting a highly-stable microwave (uplink) carrier from a ground station, having the spacecraft coherently transpond this carrier, and using a phase-locked loop receiver at the ground station to track the returned (downlink) carrier. The largest sources of measurement error are usually plasma noise and thermal noise. The plasma noise, which may originate in the ionosphere or the solar corona, is discussed; and a technique to partially calibrate its effect, involving the use of two simultaneous downlink carriers that are coherently related, is described. Range measurements employing Doppler rate-aiding are also described.

Kinman, Peter W.↗

Measuring rapid ocean tidal earth orientation variations with very long baseline interferometry

Ocean tidal effects on universal time and polar motion (UTPM) are investigated at four nearly diurnal (K(sub 1), P(sub 1), O(sub 1), and Q(sub 1)) and four nearly semidiurnal (K(sub 2), S(sub 2), M(sub 2), and N(sub 2)) frequencies by analyzing very long baseline interferometry (VLBI) data extending from 1978 to 1992. We discuss limitations of comparisons between experiment and theory for the retograde nearly diurnal polar motion components due to their degeneracy with prograde components of the nutation model. Estimating amplitudes of contributions to the modeled VLBI observables at these eight frequencies produces a statistically highly significant improvement of 7 mm to the residuals of a fit to the observed delays. Use of such an improved UTPM model also reduces the 14-30 mm scatter of baseline lengths about a time-linear model of tectonic motion by 3-14 mm, also withhigh significance levels. A total of 28 UTPM ocean tidal amplitudes can be unambiguously estimated from the data, with resulting UTI and PM magnitudes as large as 21 micro secs and 270 microarc seconds and formal uncertainties of the order of 0.3 micro secs and 5 microarc secs for UTI and PM, respectively. Empirically determined UTPM amplitudes and phases are com1pared to values calculated theoretically by Gross from Seiler's global ocean tide model. The discrepancy between theory and experiment is larger by a factor of 3 for UTI amplitudes (9 micro secs) than for prograde PM amplitudes (42 microarc secs). The 14-year VLBI data span strongly attenuates the influence of mismodeled effects on estimated UTPM amplitudes and phases that are not coherent with the eight frequencies of interest. Magnitudes of coherent and quasi-coherent systematic errors are quantified by means of internal consistency tests. We conclude that coherent systematic effects are many times larger than the formal uncertainties and can be as large as 4 micro secs for UTI and 60 microarc secs for polar motion. On the basis of such ealistic error estimates, 22 of the 31 fitted UTPM ocean tidal amplitudes differ from zero by more than 2 sigma.

Sovers, O. J.↗

Heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(s)

The heterogeneous reactions of HNO3(g) + NaCl(s) yields HCl(g) + NaNO3(s) (eq 1) and N2O5(g) + NaCl(s) yields ClNO2(g) + NaNO3(S) (eq 2) were investigated over the temperature range 223-296 K in a flow-tube reactor coupled to a quadrupole mass spectrometer. Either a chemical ionization mass spectrometer (CIMS) or an electron-impact ionization mass spectrometer (EIMS) was used to provide suitable detection sensitivity and selectivity. In order to mimic atmospheric conditions, partial pressures of HNO3 and N2O5 in the range 6 x 10(exp -8) - 2 x 10(exp -6) Torr were used. Granule sizes and surface roughness of the solid NaCl substrates were determined by using a scanning electron microscope. For dry NaCl substrates, decay rates of HNO3 were used to obtain gamma(1) = 0.013 +/- 0.004 (1sigma) at 296 K and > 0.008 at 223 K, respectively. The error quoted is the statistical error. After all corrections were made, the overall error, including systematic error, was estimated to be about a factor of 2. HCl was found to be the sole gas-phase product of reaction 1. The mechanism changed from heterogeneous reaction to predominantly physical adsorption when the reactor was cooled from 296 to 223 K. For reaction 2 using dry salts, gamma(2) was found to be less than 1.0 x 10(exp -4) at both 223 and 296 K. The gas-phase reaction product was identified as ClNO2 in previous studies using an infrared spectrometer. An enhancement in reaction probability was observed if water was not completely removed from salt surfaces, probably due to the reaction of N2O5(g) + H2O(s) yields 2HNO3(g). Our results are compared with previous literature values obtained using different experimental techniques and conditions. The implications of the present results for the enhancement of the hydrogen chloride column density in the lower stratosphere after the El Chichon volcanic eruption and for the chemistry of HCl and HNO3 in the marine troposphere are discussed.

Leu, Ming-Taun↗

Task Analytic Models to Guide Analysis and Design: Use of the Operator Function Model to Represent Pilot-Autoflight System Mode Problems

Task-analytic models structure essential information about operator interaction with complex systems, in this case pilot interaction with the autoflight system. Such models serve two purposes: (1) they allow researchers and practitioners to understand pilots' actions; and (2) they provide a compact, computational representation needed to design 'intelligent' aids, e.g., displays, assistants, and training systems. This paper demonstrates the use of the operator function model to trace the process of mode engagements while a pilot is controlling an aircraft via the, autoflight system. The operator function model is a normative and nondeterministic model of how a well-trained, well-motivated operator manages multiple concurrent activities for effective real-time control. For each function, the model links the pilot's actions with the required information. Using the operator function model, this paper describes several mode engagement scenarios. These scenarios were observed and documented during a field study that focused on mode engagements and mode transitions during normal line operations. Data including time, ATC clearances, altitude, system states, and active modes and sub-modes, engagement of modes, were recorded during sixty-six flights. Using these data, seven prototypical mode engagement scenarios were extracted. One scenario details the decision of the crew to disengage a fully automatic mode in favor of a semi-automatic mode, and the consequences of this action. Another describes a mode error involving updating aircraft speed following the engagement of a speed submode. Other scenarios detail mode confusion at various phases of the flight. This analysis uses the operator function model to identify three aspects of mode engagement: (1) the progress of pilot-aircraft-autoflight system interaction; (2) control/display information required to perform mode management activities; and (3) the potential cause(s) of mode confusion. The goal of this paper is twofold: (1) to demonstrate the use of the operator functio model methodology to describe pilot-system interaction while engaging modes And monitoring the system, and (2) to initiate a discussion of how task-analytic models might inform design processes. While the operator function model is only one type of task-analytic representation, the hypothesis of this paper is that some type of task analytic structure is a prerequisite for the design of effective human-automation interaction.

Degani, Asaf↗

Top-of-Atmosphere Albedo Estimation from Angular Distribution Models using Scene Identification from Satellite Cloud Property Retrievals

The next generation of Earth radiation budget satellite instruments will routinely merge estimates of global top-of-atmosphere radiative fluxes with cloud properties. This information will offer many new opportunities for validating radiative transfer models and cloud parameterizations in climate models. In this study, five months of POLarization and Directionality of the Earth's Reflectances (POLDER) 670 nm radiance measurements are considered in order to examine how satellite cloud property retrievals can be used to define empirical Angular Distribution Models (ADMs) for estimating top-of-atmosphere (TOA) albedo. ADMs are defined for 19 scene types defined by satellite retrievals of cloud fraction and cloud optical depth. Two approaches are used to define the ADM scene types: The first assumes there are no biases in the retrieved cloud properties and defines ADMs for fixed discrete intervals of cloud fraction and cloud optical depth (fixed-tau approach). The second approach involves the same cloud fraction intervals, but uses percentile intervals of cloud optical depth instead (percentile-tau approach). Albedos generated using these methods are compared with albedos inferred directly from the mean observed reflectance field. Albedos based on ADMs that assume cloud properties are unbiased (fixed-tau approach) show a strong systematic dependence on viewing geometry. This dependence becomes more pronounced with increasing solar zenith angle, reaching approximately equals 12% (relative) between near-nadir and oblique viewing zenith angles for solar zenith angles between 60 deg and 70 deg. The cause for this bias is shown to be due to biases in the cloud optical depth retrievals. In contrast, albedos based on ADMs built using percentile intervals of cloud optical depth (percentile-tau approach) show very little viewing zenith angle dependence and are in good agreement with albedos obtained by direct integration of the mean observed reflectance field (less than 1% relative error). When the ADMs are applied separately to populations consisting of only liquid water and ice clouds, significant biases in albedo with viewing geometry are observed (particularly at low sun elevations), highlighting the need to account for cloud phase both in cloud optical depth retrievals and in defining ADM scene types. ADM-derived monthly mean albedos determined for all 5 deg x 5 deg latitude/longitude regions over ocean are in good agreement (regional RMS relative errors less than 2%) with those obtained by direct integration when ADM albedos inferred from specific angular bins are averaged together. Albedos inferred from near-nadir and oblique viewing zenith angles are the least accurate, with regional RMS errors reaching approximately 5-10% (relative). Compared to an earlier study involving ERBE ADMs, regional mean albedos based on the 19 scene types considered here show a factor of 4 reduction in bias error and a factor of 3 reduction in RMS error.

Loeb, N. G.↗

Remote Sensing of Cloud Top Heights Using the Research Scanning Polarimeter

Clouds cover roughly two thirds of the globe and act as an important regulator of Earth's radiation budget. Of these, multilayered clouds occur about half of the time and are predominantly two-layered. Changes in cloud top height (CTH) have been predicted by models to have a globally averaged positive feedback, however observational changes in CTH have shown uncertain results. Additional CTH observations are necessary to better and quantify the effect. Improved CTH observations will also allow for improved sub-grid parameterizations in large-scale models and accurate CTH information is important when studying variations in freezing point and cloud microphysics. NASA's airborne Research Scanning Polarimeter (RSP) is able to measure cloud top height using a novel multi-angular contrast approach. RSP scans along the aircraft track and obtains measurements at 152 viewing angles at any aircraft location. The approach presented here aggregates measurements from multiple scans to a single location at cloud altitude using a correlation function designed to identify the location-distinct features in each scan. During NASAs SEAC4RS air campaign, the RSP was mounted on the ER-2 aircraft along with the Cloud Physics Lidar (CPL), which made simultaneous measurements of CTH. The RSPs unique method of determining CTH is presented. The capabilities of using single and combinations of channels within the approach are investigated. A detailed comparison of RSP retrieved CTHs with those of CPL reveal the accuracy of the approach. Results indicate a strong ability for the RSP to accurately identify cloud heights. Interestingly, the analysis reveals an ability for the approach to identify multiple cloud layers in a single scene and estimate the CTH of each layer. Capabilities and limitations of identifying single and multiple cloud layers heights are explored. Special focus is given to sources of error in the method including optically thin clouds, physically thick clouds, multi-layered clouds as well as cloud phase. When determining multi-layered CTHs, limits on the upper clouds opacity are assessed.

Cloud height indicators↗

Two-Phase Convection Heat Transfer Correlations for Liquid Hydrogen Pipe Chilldown

Recently, heat transfer correlations based on liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching data were developed to improve the predictive accuracy of lumped node codes like SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP). After implementing these correlations into both programs, updated model runs showed strong improvement in LN2 pipe chilldown modeling but only modest improvement in LH2 modeling. Due to large differences in thermal and fluid properties between the two fluids, results indicated a need to develop a separate set of LH2-only correlations to improve the accuracy of the simulations. This paper presents a new set of two-phase convection heat transfer correlations based on LH2 pipe quenching data. A correlation to predict the bulk vapor temperature was developed after analysis showed that high amounts of thermal nonequilibrium of the liquid and vapor phases occurred during film boiling of LH2. Implemented in a numerical model, the new correlations achieve a mean absolute error of 19.5 K in the predicted wall temperature when compared to recent LH2 pipe chilldown data, an improvement of 40% over recent GFSSP predictions. This correlation set can be implemented in simulations of the transient LH2 chilldown process. Such simulations are useful for predicting the chilldown time and boil-off mass of LH2 for applications such as the transfer of LH2 from a ground storage tank to the rocket vehicle propellant tank, or through a rocket engine feedline during engine startup.

critical heat flux↗

Covariance Analysis Tool (G-CAT) for Computing Ascent, Descent, and Landing Errors

G-CAT is a covariance analysis tool that enables fast and accurate computation of error ellipses for descent, landing, ascent, and rendezvous scenarios, and quantifies knowledge error contributions needed for error budgeting purposes. Because GCAT supports hardware/system trade studies in spacecraft and mission design, it is useful in both early and late mission/ proposal phases where Monte Carlo simulation capability is not mature, Monte Carlo simulation takes too long to run, and/or there is a need to perform multiple parametric system design trades that would require an unwieldy number of Monte Carlo runs. G-CAT is formulated as a variable-order square-root linearized Kalman filter (LKF), typically using over 120 filter states. An important property of G-CAT is that it is based on a 6-DOF (degrees of freedom) formulation that completely captures the combined effects of both attitude and translation errors on the propagated trajectories. This ensures its accuracy for guidance, navigation, and control (GN&C) analysis. G-CAT provides the desired fast turnaround analysis needed for error budgeting in support of mission concept formulations, design trade studies, and proposal development efforts. The main usefulness of a covariance analysis tool such as G-CAT is its ability to calculate the performance envelope directly from a single run. This is in sharp contrast to running thousands of simulations to obtain similar information using Monte Carlo methods. It does this by propagating the "statistics" of the overall design, rather than simulating individual trajectories. G-CAT supports applications to lunar, planetary, and small body missions. It characterizes onboard knowledge propagation errors associated with inertial measurement unit (IMU) errors (gyro and accelerometer), gravity errors/dispersions (spherical harmonics, masscons), and radar errors (multiple altimeter beams, multiple Doppler velocimeter beams). G-CAT is a standalone MATLAB- based tool intended to run on any engineer's desktop computer.

Boussalis, Dhemetrios↗

Quantum eigenvalue estimation via time series analysis

We present an efficient method for estimating the eigenvalues of a Hamiltonian H from the expectation values of the evolution operator for various times. For a given quantum state ρ, our method outputs a list of eigenvalue estimates and approximate probabilities. Each probability depends on the support of ρ in those eigenstates of H associated with eigenvalues within an arbitrarily small range. The complexity of our method is polynomial in the inverse of a given precision parameter ε, which is the gap between eigenvalue estimates. Unlike the well-known quantum phase estimation algorithm that uses the quantum Fourier transform, our method does not require large ancillary systems, large sequences of controlled operations, or preserving coherence between experiments, and is therefore more attractive for near-term applications. The output of our method can be used to estimate spectral properties of H and other expectation values efficiently, within additive error proportional to ε.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Possible Detection of an Emission Cyclotron Resonance Scattering Feature from the Accretion-Powered Pulsar 4U 1626-67

We present analysis of 4U 1626-67, a 7.7 s pulsar in a low-mass X-ray binary system, observed with the hard X-ray detector of the Japanese X-ray satellite Suzaku in 2006 March for a net exposure of ~88 ks. The source was detected at an average 10-60 keY flux of approx 4 x 10-10 erg / sq cm/ s. The phase-averaged spectrum is reproduced well by combining a negative and positive power-law times exponential cutoff (NPEX) model modified at approx 37 keY by a cyclotron resonance scattering feature (CRSF). The phase-resolved analysis shows that the spectra at the bright phases are well fit by the NPEX with CRSF model. On the other hand. the spectrum in the dim phase lacks the NPEX high-energy cutoff component, and the CRSF can be reproduced by either an emission or an absorption profile. When fitting the dim phase spectrum with the NPEX plus Gaussian model. we find that the feature is better described in terms of an emission rather than an absorption profile. The statistical significance of this result, evaluated by means of an F test, is between 2.91 x 10(exp -3) and 1.53 x 10(exp -5), taking into account the systematic errors in the background evaluation of HXD-PIN. We find that the emission profile is more feasible than the absorption one for comparing the physical parameters in other phases. Therefore, we have possibly detected an emission line at the cyclotron resonance energy in the dim phase.

Iwakiri, W. B.↗