Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “positive definite functions”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Orbital Manuvering System Design and Performance For the Magnetosperic Multiscale Constellation

The Magnetospheric Multiscale (MMS) mission, launched on March 13, 2015, is the fourth mission of NASA's Solar Terrestrial Probe program. The MMS mission consists of four identically instrumented observatories that function as a constellation to provide the first definitive study of magnetic reconnection in space. Since it is frequently desirable to isolate electric and magnetic field sensors from stray effects caused by the spacecraft's core-body, the suite of instruments on MMS includes six radial and two axial instrument-booms with deployed lengths ranging from 5-60 meters (see Figure 1). The observatory is spin-stabilized about its positive z-axis with a nominal rate slightly above 3 rev/min (RPM). The spin is also used to maintain tension in the four radial wire-booms. Each observatory's Attitude Control System (ACS) consists of digital sun sensors, star cameras, accelerometers, and mono-propellant hydrazine thrusters-responsible for orbital adjustments, attitude control, and spin adjustments. The sections that follow describe performance requirements, the hardware and algorithms used for 6-DOF estimation, and then similarly for 6-DOF control. The paper concludes with maneuver performance based on both simulated and on-orbit telem.

Determination and Control↗

Application of multidimensional structural characterization to detect and describe moderate forest disturbance

The study of vegetation community and structural change has been central to ecology for over a century, yet the ways in which disturbances reshape the physical structure of forest canopies remain relatively unknown. Moderate severity disturbances affect different canopy strata and plant species, resulting in variable structural outcomes and ecological consequences. Terrestrial lidar (light detection and ranging) offers an unprecedented view of the interior arrangement and distribution of canopy elements, permitting the derivation of multidimensional measures of canopy structure that describe several canopy structural traits (CSTs) with known links to ecosystem function. We used lidar-derived CSTs within a machine learning framework to detect and describe the structural changes that result from various disturbance agents, including moderate severity fire, ice storm damage, age-related senescence, hemlock woolly adelgid, beech bark disease, and chronic acidification. We found that fire and ice storms primarily affected the amount and position of vegetation within canopies, while acidification, senescence, pathogen, and insect infestation altered canopy arrangement and complexity. Only two of the six disturbance agents significantly reduced leaf area, counter to common assumptions regarding many moderate severity disturbances. While findings are limited in their generalizability due to lack of replication among disturbances, they do suggest that the current limitations of standard disturbance detection methods—such as optical-based remote sensing platforms, which are often above-canopy perspectives—limit our ability to understand the full ecological and structural impacts of disturbance, and to evaluate the consistency of structural patterns within and among disturbance agents. A more broadly inclusive definition of ecological disturbance that incorporates multiple aspects of canopy structural change may potentially improve the modeling, detection, and prediction of functional implications of moderate severity disturbance as well as broaden our understanding of the ecological impacts of disturbance.

54 ENVIRONMENTAL SCIENCES↗

Deconditioning and Reconditioning: Humans in Stressful Environments

Deconditioning is an integrated physiological response of the body to a reduction in metabolic rate; that is, to a reduction in energy use or in exercise level. While it may involve assumption of a horizontal body position, it certainly perturbs bodily homeostasis - at least temporarily. The reduction in physical activity that causes deconditioning is often associated with an increase in the time spent, for whatever reason, in a sitting or horizontal position. As a result, orthostatic factors may also contribute to the deconditioning mechanism. The word decondition may be defined as "1: to cause extinction of (a conditioned response) 2: to cause to lose physical fitness". This definition implies that psychological/emotional factors may accompany physical deconditioning, and it is this interpretation of the word that is used throughout this volume. It is apparent that deconditioning plays a major role in the mechanism of the general adaptive (homeostatic) response that is initiated by exposure to prolonged bed rest (BR). And the total homeostatic response to BR involves more than deconditioning per se. For example, it has been shown that the restoration of plasma volume and maximal work capacity after 4 weeks of BR deconditioning left other bodily functions (submaximal exercise oxygen uptake and cardiac output, leg proprioception and posterior leg muscle thickness and volume, head-up tilt tolerance, and sleep quality) functioning at decreased levels. The precise effect of deconditioning on BR homeostasis is difficult to determine, because the fundamental interactive neuro-endocrine-immune control networks that facilitate conditioning and deconditioning also act to maintain basic wholebody homeostasis. For example, is the mechanism of BR-induced deconditioning independent of the mechanism that provokes concomitant orthostatic intolerance; that is, fainting? Assumption of the recumbent body position for prolonged periods of time, results in a new adaptive-homeostatic state. This state occurs in response to the mutually interactive effects of the change in bodily position (hydrostatic pressure), to the virtual elimination of longitudinal pressure on the bones, to the increased confinement with possible reduction in total daily energy (exercise) expenditure, to the reorientation of stimuli within the vestibular organs, and (often) to altered socio-psychological conditions. The exercise-training (reconditioning) syndrome affects total body homeostasis by facilitating increases in work capacity and endurance, whereas deconditioning decreases physical performance. There are many interrelated factors that influence the control parameters that seek to maintain the adaptive conditioning-deconditioning syndrome. These control parameters can be better elucidated by subjecting otherwise healthy ambulatory people to various stresses, such as exercise training and prolonged spaceflight, bed rest, water immersion, hyperbaria, and isolation and confinement. Changes in control parameters will be manifested in muscle function, orthostatic tolerance, cardiorespiratory responses, musculo-skeletal systems, free-radical processes, and body thermoregulation with overarching effects on the subjects' psycho-sociological states. A discussion of these factors and the control parameters constitutes the substance of this volume. Special emphasis is placed on delineating practical applications of the findings that will be of special interest to physicians, nurses, and other health-care workers.

Greenleaf, John E.↗

Comparing hard and soft prior bounds in geophysical inverse problems

In linear inversion of a finite-dimensional data vector y to estimate a finite-dimensional prediction vector z, prior information about X sub E is essential if y is to supply useful limits for z. The one exception occurs when all the prediction functionals are linear combinations of the data functionals. Two forms of prior information are compared: a soft bound on X sub E is a probability distribution p sub x on X which describeds the observer's opinion about where X sub E is likely to be in X; a hard bound on X sub E is an inequality Q sub x(X sub E, X sub E) is equal to or less than 1, where Q sub x is a positive definite quadratic form on X. A hard bound Q sub x can be softened to many different probability distributions p sub x, but all these p sub x's carry much new information about X sub E which is absent from Q sub x, and some information which contradicts Q sub x. Both stochastic inversion (SI) and Bayesian inference (BI) estimate z from y and a soft prior bound p sub x. If that probability distribution was obtained by softening a hard prior bound Q sub x, rather than by objective statistical inference independent of y, then p sub x contains so much unsupported new information absent from Q sub x that conclusions about z obtained with SI or BI would seen to be suspect.

Backus, George E.↗

Comparing hard and soft prior bounds in geophysical inverse problems

In linear inversion of a finite-dimensional data vector y to estimate a finite-dimensional prediction vector z, prior information about X sub E is essential if y is to supply useful limits for z. The one exception occurs when all the prediction functionals are linear combinations of the data functionals. Two forms of prior information are compared: a soft bound on X sub E is a probability distribution p sub x on X which describes the observer's opinion about where X sub E is likely to be in X; a hard bound on X sub E is an inequality Q sub x(X sub E, X sub E) is equal to or less than 1, where Q sub x is a positive definite quadratic form on X. A hard bound Q sub x can be softened to many different probability distributions p sub x, but all these p sub x's carry much new information about X sub E which is absent from Q sub x, and some information which contradicts Q sub x. Both stochastic inversion (SI) and Bayesian inference (BI) estimate z from y and a soft prior bound p sub x. If that probability distribution was obtained by softening a hard prior bound Q sub x, rather than by objective statistical inference independent of y, then p sub x contains so much unsupported new information absent from Q sub x that conclusions about z obtained with SI or BI would seen to be suspect.

Backus, George E.↗

Fully autonomous navigation for the NASA cargo transfer vehicle

A great deal of attention has been paid to navigation during the close approach (less than or equal to 1 km) phase of spacecraft rendezvous. However, most spacecraft also require a navigation system which provides the necessary accuracy for placing both satellites within the range of the docking sensors. The Microcosm Autonomous Navigation System (MANS) is an on-board system which uses Earth-referenced attitude sensing hardware to provide precision orbit and attitude determination. The system is capable of functioning from LEO to GEO and beyond. Performance depends on the number of available sensors as well as mission geometry; however, extensive simulations have shown that MANS will provide 100 m to 400 m (3(sigma)) position accuracy and 0.03 to 0.07 deg (3(sigma)) attitude accuracy in low Earth orbit. The system is independent of any external source, including GPS. MANS is expected to have a significant impact on ground operations costs, mission definition and design, survivability, and the potential development of very low-cost, fully autonomous spacecraft.

Wertz, James R.↗

Three-Dimensional Field Solutions for Multi-Pole Cylindrical Halbach Arrays in an Axial Orientation

This article presents three-dimensional B field solutions for the cylindrical Halbach array in an axial orientation. This arrangement has applications in the design of axial motors and passive axial magnetic bearings and couplers. The analytical model described here assumes ideal magnets with fixed and uniform magnetization. The field component functions are expressed as sums of 2-D definite integrals that are easily computed by a number of mathematical analysis software packages. The analysis is verified with sample calculations and the results are compared to equivalent results from traditional finite-element analysis (FEA). The field solutions are then approximated for use in flux linkage and induced EMF calculations in nearby stator windings by expressing the field variance with angular displacement as pure sinusoidal function whose amplitude depends on radial and axial position. The primary advantage of numerical implementation of the analytical approach presented in the article is that it lends itself more readily to parametric analysis and design tradeoffs than traditional FEA models.

Thompson, William K.↗

Feedback on Forty-year Long Clean-up Operations of a Contaminated Soil for Environmental Purpose - 20026

One of the CEA facilities in France is the place of miscellaneous mapping and clean-up operations since the 1980's. The final purpose of the conducted work deals with the environmental remediation of the building. In France, the absence of a regulatory framework for the management of sites polluted by radioactive substances, especially as the lack of release thresholds, has led to develop different approaches of soil radiological characterization. The geostatistical approach is part of them. Preliminary investigations are an unavoidable step prior to clearance and remediation. First, historical analysis enables to define site perimeters to investigate as well as to precisely locate expected contaminated parts. For radioactive pollutants, the media (water, soil, air) and transfer routes participate in the definition of the investigation perimeter and thus also need to be identified through a geological study. At the same time, functional analysis leads to choose the best measurement means which could provide useful information when combined with visual inspection. An equipment well adapted to easy-to-measure radionuclides as gamma emitters, such as gamma probe or on-site gamma spectrometry, is usually selected as a non-destructive assay. Coupled with spatial positions, processed on-site radiological data can be mapped. Data processing consists in mapping and kriging before a geostatistical analysis in order to identify the zones of interest to be targeted. Selected zones are not only linked to high level count rates but also to low level ones so that a pollution-free reference is known in the area. This data processing also enables to categorize waste types according to their origin and contamination levels. An expected volume of different waste then follows as well as the waste characterization equipments. In-depth investigations, e.g. core drillings and samplings, naturally ensue to aim at remediation optimization through the assessment of environmental impact. After clean-up operations, non-destructive assays as well as destructive samplings enable to check the effective clearance. Following this methodology in the present remediation work, results of the first non-destructive assay campaign highlighted a long-lived transuranic radionuclides contamination in the concrete flagstone. Successive non-destructive assays were then performed in the defined zones of interest. Count rates measured with a surface probe were first mapped. The best localization of destructive assays, i.e. core drillings at a 1 meter depth, came from this surface mapping followed by a geostatistical analysis. The concrete flagstone was totally removed. Then, the depth of soil to remove under the flagstone was optimized from the results of radiological activities measured in core drillings samples. A new surface mapping was then drawn after this partial soil excavation. The sand and demolition rubble samples, homogeneously constituted from the flagstone fragments, were analyzed by gamma spectrometry. Meanwhile, waste drums and other containers were also measured by gamma spectrometry and passive neutron measurement devices. All these results have brought miscellaneous pieces of information. In this context, major issues appear in terms of physical constraints such as the premises tininess as well as of soil sampling techniques and in terms of measurement performances assessment and the representativeness of homogeneous samples. The choice of radiological soil and waste characterization devices is described in the light of enhanced performances by keeping in mind radiation protection requirements and a willing to always optimize waste categories (very low level waste) and volume. The improvement of the characterization methods appears through this forty-year long work while sharing the feedback of encountered difficulties. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Definition study of land/sea civil user navigational location monitoring systems for NAVSTAR GPS: User requirements and systems concepts

A low-cost GPS civil-user mobile terminal whose purchase cost is substantially an order of magnitude less than estimates for the military counterpart is considered with focus on ground station requirements for position monitoring of civil users requiring this capability and the civil user navigation and location-monitoring requirements. Existing survey literature was examined to ascertain the potential users of a low-cost NAVSTAR receiver and to estimate their number, function, and accuracy requirements. System concepts are defined for low cost user equipments for in-situ navigation and the retransmission of low data rate positioning data via a geostationary satellite to a central computing facility.

Devito, D. M.↗

Journal bearing impedance descriptions for rotordynamic applications

The paper deals with the development of analytic descriptions for plain circumferentially-symmetric fluid journal bearings, which are suitable for use in rotor dynamic analysis. The bearing impedance vector is introduced, which defines the bearing reaction force components as a function of the bearing motion. Impedances are derived directly for the Ocvirk (short) and Sommerfeld (long) bearings, and the relationships between the impedance vector and the more familiar mobility vector are developed and used to derive analytic impedance for finite-length bearings. The static correctness of the finite-length cavitating impedance is verified. Analytic stiffness and damping coefficient definitions are derived in terms of an impedance vector for small motion around an equilibrium position and demonstrated for the finite-length cavitating impedance. Nonlinear transient rotordynamic simulations are presented for the short pi and 2-pi impedances and the finite-length cavitating impedance. It is shown that finite-length impedance yields more accurate results for substantially less computer time than the short-bearing numerical-pressure-integration approach.

Childs, D.↗

Implementation of a Wavefront-Sensing Algorithm

A computer program has been written as a unique implementation of an image-based wavefront-sensing algorithm reported in "Iterative-Transform Phase Retrieval Using Adaptive Diversity" (GSC-14879-1), NASA Tech Briefs, Vol. 31, No. 4 (April 2007), page 32. This software was originally intended for application to the James Webb Space Telescope, but is also applicable to other segmented-mirror telescopes. The software is capable of determining optical-wavefront information using, as input, a variable number of irradiance measurements collected in defocus planes about the best focal position. The software also uses input of the geometrical definition of the telescope exit pupil (otherwise denoted the pupil mask) to identify the locations of the segments of the primary telescope mirror. From the irradiance data and mask information, the software calculates an estimate of the optical wavefront (a measure of performance) of the telescope generally and across each primary mirror segment specifically. The software is capable of generating irradiance data, wavefront estimates, and basis functions for the full telescope and for each primary-mirror segment. Optionally, each of these pieces of information can be measured or computed outside of the software and incorporated during execution of the software.

Smith, Jeffrey S.↗

G2Aero: A Python package for separable shape tensors

G2Aero is a Python package for the design and deformation of discrete planar curves and tubular surfaces using a geometric data-driven approach. G2Aero utilizes a topology of product manifolds: the Grassmannian, $\mathcal{G}$($\mathcal{n}$, 2) - the set of 2-dimensional subspaces in $\mathbb{R}$ $\mathcal{n}$ - and the symmetric positive-definite (SPD) manifold, $\mathcal{S}^{2}_{++}$ - the set of 2x2 SPD matrices. The package provides a novel framework for representing separable deformations to shapes, which consist of stretching, scaling, rotating, and translating - also known as affine deformations - and a set of complementary deformations - which we refer to as undulation-type deformations. We focus on airfoil and blade design applications to emphasize the utility of the methods in an environment where the separation of affine and undulation-type deformations is critical. Notable functionalities of the framework for blade design include: 1) generating novel 2D (airfoil) shapes informed by a database of physically relevant airfoils, 2) building 3D blades by interpolating sequences of 2D airfoil cross-sections, and 3) generating blades with consistent perturbations along the blade span. We discuss the framework and provide examples in the context of wind energy applications, specifically wind turbine blade design. Figure 1 shows the wire frame obtained by interpolating airfoils defining the IEA 15-MW wind turbine blade and applying affine transformations corresponding to twist, chordal scaling, and bending. This, and all other figures in the paper, can be reproduced following examples and referencing supporting documentation provided in the G2Aero package.

97 MATHEMATICS AND COMPUTING↗

A note on correlation description

The paper suggests some modifications to definitions of correlation functions made by Freeman and Bennett (1974) specifically for comparing measured quantities with those predicted by rotor theories, with reference to rotorcraft flight simulation. Special attention is paid to the quality of correlation, determined by a quality factor relating the measured and computed entities. Correlation criteria (mean load and rotor performance, overall amplitude and phase, harmonic altitude and harmonic phase) are examined, along with phase angle correlation criteria (harmonic phase defined by sine-cosine relation and harmonic phase defined by a first positive peak rotor azimuth position)

Mccloud, J. L., III↗

Concurrent Cholesky factorization of positive definite banded Hermitian matrices

First, the Cholesky factorization is extended to cover uniformly partitioned banded positive definite matrices of rank n which may be real symmetric or Hermitian. Then, two stratagems are given for the use of the algorithm in concurrent machines where the number of processing elements is less than required to factor the matrix in as few serial steps as possible, and where uniformly high efficiency is expected from all processing elements. Expressions are given for the efficiency factor e appearing in the speed-up expression q = eN, and these are specialized for the N node hypercube machine as a function of partition size s, the number N of processing elements of the hypercube machine, and the cost mu of interelement transmission relative to computation. It is shown that the efficiency factor e is inversely proportional to mu/s, and that e is almost independent of N when N is large and mu/s = 0. The task is completed in n/s serial steps with no limit on n. The half bandwidth b of the matrix is 2 Ns.

Utku, S.↗

Delaying Trains of Short Light Pulses in WGM Resonators

Suitably configured whispering-gallery-mode (WGM) optical resonators have been proposed as delay lines for trains of short light pulses. Until now, it has been common practice to implement an optical delay line as a coiled long optical fiber, which is bulky and tends to be noisy. An alternative has been to implement an optical delay line as a coupled-resonator optical waveguide (a chain of coupled optical resonators), which is compact but limits the width of the pulse spectrum to the width of an optical resonance and thereby places a lower limit on the duration of a pulse. In contrast, a delay line according to the proposal could be implemented as a single WGM resonator, and the pulses delayed by the resonator could be so short that their spectral widths could greatly exceed the spectral width of any single resonance. The proposal emerged from theoretical and experimental studies of the propagation of a pulse train in a WGM resonator. An important element of the theoretical study was recognition that the traditional definition of group velocity in effect, the velocity of a single pulse comprising a packet of waves propagating in a medium, the responsivity of which is a monotonous function of frequency does not necessarily apply in the case of a WGM resonator or other medium having a spectrum consisting of discrete resonance peaks at different frequencies. A new definition of group velocity, applicable to a train of pulses propagating in such a medium, was introduced and found to lead to the discovery of previously unknown features of propagation. Notably, it was found that in a micro - sphere optical resonator that supports a suitable combination of WGM modes, the group velocity for a train of light pulses could be positive, zero, or negative. A positive group velocity could be so small that the delay could be much longer than the ring-down time of the resonator; a delay of such great length is impossible for a single pulse interacting with either a linearly responding lossless resonator or a coupled- resonator optical waveguide.

Matsko, Andrey↗

Selective Enhancement of Spectroscopic Features by Quantum Optimal Control

Tailored light can be used to steer atomic motions into selected quantum pathways. In optimal control theory (OCT), the target is usually expressed in terms of the molecular wave function, a quantity that is not directly observable in experiment. In this study, we present simulations using OCT that optimize the spectroscopic signal itself. By shaping the optical pump, the x-ray stimulated Raman signal, which occurs solely during the passage through conical intersections, is temporally controlled and amplified by up to 2 orders of magnitude. This enhancement can be crucial in order to bring small coherence-based signatures above the detectable threshold. Our approach is applicable to any signal that depends on the expectation value of a positive definite operator.

74 ATOMIC AND MOLECULAR PHYSICS↗

Note on the definitions of branching ratios of overlapping resonances

Branching ratios for the decay of hadrons with large width or near thresholds depend on their definition. We test different definitions and show that rather different branching ratios can be obtained. For wide resonances and for sequential decays with wide intermediate resonances, integration over the spectral functions is mandatory. The tests are performed exploiting the latest solution of the Bonn-Gatchina multi-channel analysis and published values for residues of light scalar mesons. For a resonance overlapping with a threshold, in case its pole lies in a non-adjacent sheet, we show how the total width, needed for the branching ratios, does not correspond to the imaginary part of the pole position. We use the Madrid-Krakow dispersive parameterizations to illustrate this situation with the f 0 (980).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Formally Verified Floating-Point Implementation of the Compact Position Reporting Algorithm

The Automatic Dependent Surveillance-Broadcast (ADS-B) system allows aircraft to communicate their current state, including position and velocity information, to other aircraft in their vicinity and to ground stations. The Compact Position Reporting (CPR) algorithm is the ADS-B module responsible for the encoding and decoding of aircraft positions. CPR is highly sensitive to computer arithmetic since it heavily relies on functions that are intrinsically unstable such as floor and modulo. In this paper, a formally-verified double-precision floating-point implementation of the CPR algorithm is presented. The verification proceeds in three steps. First, an alternative version of CPR, which reduces the floating-point rounding error is proposed. Then, the Prototype Verification System (PVS) is used to formally prove that the ideal real-number counterpart of the improved algorithm is mathematically equivalent to the standard CPR definition. Finally, the static analyzer Frama-C is used to verify that the double-precision implementation of the improved algorithm is correct with respect to its operational requirement. The alternative algorithm is currently being considered for inclusion in the revised version of the ADS-B standards document as the reference implementation of the CPR algorithm.

Laura Titolo↗