Search NASA⌕ Search

SEARCH · Search NASA

Results for “invariant domain”

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 163 records · Page 9

Spatial and Spatial Frequency Spreads Of Masking: Measurements and a Contrast-Gain-Control Model

Masked pattern detection is affected by the mask's proximity to the target in both the spatial and the spatial frequency domains. We measured contrast thresholds for detection of a Gabor patch centered in a sample of static band-pass noise with a central noise-free aperture. To examine spatial spread, target and mask frequency were identical, at either 2, 4, or 8 cycles/degree (cpd), and the radius of the aperture was 0, 1/2 or 1 cycle of the noise band's center frequency. The results suggest that the spatial spread of masking is scale invariant, with the largest radius producing little masking. To examine spatial frequency spread, a radius of 0 (no aperture) was used with all pairings of 2,4, and 8 cpd for target and mask. The results suggest that masking is asymmetrical over log frequency: 8 cpd noise masks a 2 cpd target, but 2 cpd noise does not mask an 8 cpd target. We have used these results to calibrate a model of contrast-gain-control.

Solomon, J. A.↗

Efficient Computation of Closed-loop Frequency Response for Large Order Flexible Systems

An efficient and robust computational scheme is given for the calculation of the frequency response function of a large order, flexible system implemented with a linear, time invariant control system. Advantage is taken of the highly structured sparsity of the system matrix of the plant based on a model of the structure using normal mode coordinates. The computational time per frequency point of the new computational scheme is a linear function of system size, a significant improvement over traditional, full-matrix techniques whose computational times per frequency point range from quadratic to cubic functions of system size. This permits the practical frequency domain analysis of systems of much larger order than by traditional, full-matrix techniques. Formulations are given for both open and closed loop loop systems. Numerical examples are presented showing the advantages of the present formulation over traditional approaches, both in speed and in accuracy. Using a model with 703 structural modes, a speed-up of almost two orders of magnitude was observed while accuracy improved by up to 5 decimal places.

Maghami, Peiman G.↗

Coarsening dynamics of Ising-nematic order in a frustrated Heisenberg antiferromagnet

We study the phase ordering dynamics of the classical antiferromagnetic 𝐽 1 −𝐽 2 (nearest-neighbor and next-nearest-neighbor couplings) Heisenberg model on the square lattice in the strong frustration regime (𝐽 2 /𝐽 1 > 1/2). While thermal fluctuations preclude any long-range magnetic order at finite temperatures, the system exhibits a long-range spin-driven nematic phase at low temperatures. The transition into the nematic phase is further shown to belong to the two-dimensional Ising universality class based on the critical exponents near the phase transition. Our large-scale stochastic Landau-Lifshitz-Gilbert simulations find a two-stage phase ordering when the system is quenched from a high-temperature paramagnetic state into the nematic phase. In the early stage, collinear alignments of spins lead to a locally saturated Ising-nematic order. Once domains of well-defined Ising order are developed, the late-stage relaxation is dominated by curvature-driven domain coarsening, as described by the Allen-Cahn equation. The characteristic size of Ising-nematic domains scales as the square root of time, similar to the kinetic Ising model described by the time-dependent Ginzburg-Landau theory. Our results confirm that the late-stage ordering kinetics of the spin-driven nematic, which is a vestigial order of the frustrated Heisenberg model, belongs to the dynamical universality class of a nonconserved Ising order. Interestingly, the system shows no violation of the superuniversality hypothesis under weak bond disorder. The dynamic scaling invariance is preserved in the presence of weak bond disorder. Here, we also discuss possible applications of our results to materials for which vestigial Ising-nematic order is realized.

Antiferromagnets↗

Imaging Analysis of the Hard X-Ray Telescope ProtoEXIST2 and New Techniques for High-Resolution Coded-Aperture Telescopes

Wide-field (greater than or approximately equal to 100 degrees squared) hard X-ray coded-aperture telescopes with high angular resolution (greater than or approximately equal to 2 minutes) will enable a wide range of time domain astrophysics. For instance, transient sources such as gamma-ray bursts can be precisely localized without the assistance of secondary focusing X-ray telescopes to enable rapid followup studies. On the other hand, high angular resolution in coded-aperture imaging introduces a new challenge in handling the systematic uncertainty: the average photon count per pixel is often too small to establish a proper background pattern or model the systematic uncertainty in a timescale where the model remains invariant. We introduce two new techniques to improve detection sensitivity, which are designed for, but not limited to, a high-resolution coded-aperture system: a self-background modeling scheme which utilizes continuous scan or dithering operations, and a Poisson-statistics based probabilistic approach to evaluate the significance of source detection without subtraction in handling the background. We illustrate these new imaging analysis techniques in high resolution coded-aperture telescope using the data acquired by the wide-field hard X-ray telescope ProtoEXIST2 during a high-altitude balloon flight in fall 2012. We review the imaging sensitivity of ProtoEXIST2 during the flight, and demonstrate the performance of the new techniques using our balloon flight data in comparison with a simulated ideal Poisson background.

balloons - instrumentation: detectors - techniques↗

Design of non-interacting multivariable discrete time systems

The theory of decoupling a multivariable system with the help of state variable feedback is applied to discrete time systems. The system differential equations in the continuous time domain are converted to a discrete time representation. A discrete time controller is designed to control the overall input to the system such that a specified desired decoupled output response is obtained. The decoupled system behaves as a set of single-input, single-output systems and state variable techniques can be easily applied to control each output individually. Control is achieved by periodically sampling the inputs and all the states and feeding them back through a controller. For linear time invariant systems the control is time invariant. The basic theory of decoupling involves feeding back all the states twice, first to decouple the system and then to achieve a desired response. The states of the system are computed with the help of a state estimator which calculates all the present states of the system so that on-line control can be applied to real time systems. The theory is applied to a model of a nuclear rocket engine.

Agrawal, B. L.↗

System identification of the Large-Angle Magnetic Suspension Test Facility (LAMSTF)

The Large-Angle Magnetic Suspension Test Facility (LAMSTF), a laboratory-scale research project to demonstrate the magnetic suspension of objects over wide ranges of attitudes, has been developed. This system represents a scaled model of a planned Large-Gap Magnetic Suspension System (LGMSS). The LAMSTF system consists of a planar array of five copper electromagnets which actively suspend a small cylindrical permanent magnet. The cylinder is a rigid body and can be controlled to move in five independent degrees of freedom. Five position variables are sensed indirectly by using infra-red light-emitting diodes and light-receiving phototransistors. The motion of the suspended cylinder is in general nonlinear and hence only the linear, time-invariant perturbed motion about an equilibrium state is considered. One of the main challenges in this project is the control of the suspended element over a wide range of orientations. An accurate dynamic model plans an essential role in controller design. The analytical model of the LAMSTF system includes highly unstable real poles (about 10 Hz) and low-frequency flexible modes (about 0.16 Hz). Projection filters are proposed to identify the state space model from closed-loop test data in time domain. A canonical transformation matrix is also derived to transform the identified state space model into the physical coordinate. The LAMSTF system is stabilized by using a linear quadratic regulator (LQR) feedback controller. The rate information is obtained by calculating the back difference of the sensed position signals. The reference inputs contain five uncorrelated random signals. This control input and the system reponse are recorded as input/output data to identify the system directly from the projection filters. The sampling time is 4 ms and the model is fairly accurate in predicting the step responses for different controllers while the analytical model has a deficiency in the pitch axis.

Huang, Jen-Kuang↗

Synthesizing Certified Code

Code certification is a lightweight approach to demonstrate software quality on a formal level. Its basic idea is to require producers to provide formal proofs that their code satisfies certain quality properties. These proofs serve as certificates which can be checked independently. Since code certification uses the same underlying technology as program verification, it also requires many detailed annotations (e.g., loop invariants) to make the proofs possible. However, manually adding theses annotations to the code is time-consuming and error-prone. We address this problem by combining code certification with automatic program synthesis. We propose an approach to generate simultaneously, from a high-level specification, code and all annotations required to certify generated code. Here, we describe a certification extension of AUTOBAYES, a synthesis tool which automatically generates complex data analysis programs from compact specifications. AUTOBAYES contains sufficient high-level domain knowledge to generate detailed annotations. This allows us to use a general-purpose verification condition generator to produce a set of proof obligations in first-order logic. The obligations are then discharged using the automated theorem E-SETHEO. We demonstrate our approach by certifying operator safety for a generated iterative data classification program without manual annotation of the code.

Whalen, Michael↗

Invariant two-dimensional structures of vorticity field in nonviscous fluid

It is shown that a fluid flow which are invariant under translations along their vorticity field directions can be described in terms of a formally 2D system. Its nonviscous behavior of the scalar vorticity is governed by Hamiltonian dynamics, similar to the 2D line vortex case. Applying a statistical mechanical method, the equilibrium distribution of the vorticity field for such a general 2D system is shown to obey a nonlinear partial differential equation that is a generalization of the sinh-Poisson equation derived by Montgomery and Joyce (1973, 1974). Thus, if the domain containing the fluid is finite, similar negative temperature states may occur and the distribution of the vorticity intensity may become localized in space. The results are believed to have applications in studying large-scale partial and temporal properties of flows having small but nonzero viscosity. The analysis can probably be extended to other fluid systems with curved magnetic field lines.

Chen, Hudong↗

Automatic Échelle Spectrograph Wavelength Calibration

Time-domain astronomy and the increasing number of exoplanet candidates call for reliable, robust, and automatic wavelength calibration. We present an algorithm for wavelength-calibrating échelle spectrographs that uses order-by-order extracted spectra and a list of laboratory wavelengths. Our approach is fully automatic and does not need the pixel locations of certain spectral features with which to anchor the wavelength solution, nor the true order number of each diffraction order. We use spectral features that are duplicated in adjacent orders to establish the scale-invariant component of the wavelength solution. We then match the central wavelengths of spectral features to laboratory wavelengths to establish the scale and higher-order components of the wavelength solution. We demonstrate our method on the four spectrographs of Las Cumbres Observatory’s Network of Robotic Échelle Spectrographs (NRES), on the High Accuracy Radial Velocity Planet Searcher (HARPS) spectrograph, and on synthetic data. We obtain a velocity-equivalent precision of ∼10 m s{sup −1} on NRES. We achieve ∼1 m s{sup −1} on HARPS, which agrees with the precision reported by the HARPS team. On synthetic data, we achieve the velocity precision set by Gaussian centroiding errors. Our algorithm likely holds for a wide range of spectrographs beyond the five presented here. We provide an open-source Python package, xwavecal, which outputs wavelength-calibrated spectra as well as the wavelengths of spectral features.

79 ASTRONOMY AND ASTROPHYSICS↗

System identification of the Large-Angle Magnetic Suspension Test Fixture (LAMSTF)

The Large-Angle Magnetic Suspension Test Fixture (LAMSTF), a laboratory-scale research project to demonstrate the magnetic suspension of objects over wide ranges of attitudes, has been developed. This system represents a scaled model of a planned Large-Gap Magnetic Suspension System (LGMSS). The LAMSTF consists of a small cylindrical permanent magnet suspended element which is levitated above a planar array of five electromagnets mounted in a circular configuration. The cylinder is a rigid body and can be controlled to move in five independent degrees of freedom. Six position variables are sensed indirectly by using infrared light-emitting diodes and light-receiving phototransistors. The motion of the suspended cylinder is in general nonlinear and hence only the linear, time-invariant perturbed motion about an equilibrium state is considered. One of the main challenges in this project is the control of the suspended element over a wide range of orientations. An accurate dynamic model plays an essential role in controller design. The analytical model is first derived and open-loop characteristics discussed. The system is shown to be highly unstable and requires feedback control for system identification. Projection filters are first proposed to identify the state space model from closed-loop input/output test data in the time domain. This method is then extended to identify linear systems from the frequency test data. A canonical transformation matrix is also derived to transform the identified state space model into the physical coordinate. The LAMSTF system is stabilized by using a linear quadratic regulator (LQR) feedback controller for closed-loop identification. The rate information is obtained by calculating the back difference of the sensed position signals. Only the closed-loop random input/output data are recorded. Preliminary results from numerical simulations demonstrate that the identified system model is fairly accurate from either time domain or frequency-domain data. Experiments will be performed to validate the proposed closed-loop identification algorithms.

Huang, Jen-Kuang↗

Assimilation of SMOS Brightness Temperatures or Soil Moisture Retrievals into a Land Surface Model

Three different data products from the Soil Moisture Ocean Salinity (SMOS) mission are assimilated separately into the Goddard Earth Observing System Model, version 5 (GEOS-5) to improve estimates of surface and root-zone soil moisture. The first product consists of multi-angle, dual-polarization brightness temperature (Tb) observations at the bottom of the atmosphere extracted from Level 1 data. The second product is a derived SMOS Tb product that mimics the data at a 40 degree incidence angle from the Soil Moisture Active Passive (SMAP) mission. The third product is the operational SMOS Level 2 surface soil moisture (SM) retrieval product. The assimilation system uses a spatially distributed ensemble Kalman filter (EnKF) with seasonally varying climatological bias mitigation for Tb assimilation, whereas a time-invariant cumulative density function matching is used for SM retrieval assimilation. All assimilation experiments improve the soil moisture estimates compared to model-only simulations in terms of unbiased root-mean-square differences and anomaly correlations during the period from 1 July 2010 to 1 May 2015 and for 187 sites across the US. Especially in areas where the satellite data are most sensitive to surface soil moisture, large skill improvements (e.g., an increase in the anomaly correlation by 0.1) are found in the surface soil moisture. The domain-average surface and root-zone skill metrics are similar among the various assimilation experiments, but large differences in skill are found locally. The observation-minus-forecast residuals and analysis increments reveal large differences in how the observations add value in the Tb and SM retrieval assimilation systems. The distinct patterns of these diagnostics in the two systems reflect observation and model errors patterns that are not well captured in the assigned EnKF error parameters. Consequently, a localized optimization of the EnKF error parameters is needed to further improve Tb or SM retrieval assimilation.

SMOS↗

Implementing distributed Ada for real-time applications

The discussion of applying a (distributed) High Order Language (HOL) to applications which require real-time performance invariably invokes the subject of excessive overhead. In a related paper, some of the basic language-specific issues involved in distributing a High Order Language, with special attention paid to the Ada language, are discussed. In a traditional implementation, several of these issues imply considerable, if not prohibitive overhead. An implementation strategy is introduced which promises to deal with these issues in a manner that will provide significant performance improvements. These improvements should in fact be sufficient to make use of distributed Ada feasible even in highly-constrained application domains. Additionally, the general approach should be applicable to nondistributed implementations as well.

Rogers, Patrick↗

The Characterization of Deep Convective Clouds as an Invariant Calibration Target and as a Visible Calibration Technique

Deep convective clouds (DCCs) are ideal visible calibration targets because they are bright nearly isotropic solar reflectors located over the tropics and they can be easily identified using a simple infrared threshold. Because all satellites view DCCs, DCCs provide the opportunity to uniformly monitor the stability of all operational sensors, both historical and present. A collective DCC anisotropically corrected radiance calibration approach is used to construct monthly probability distribution functions (PDFs) to monitor sensor stability. The DCC calibration targets were stable to within 0.5% and 0.3% per decade when the selection criteria were optimized based on Aqua MODerate Resolution Imaging Spectroradiometer 0.65-micrometer-band radiances. The Tropical Western Pacific (TWP), African, and South American regions were identified as the dominant DCC domains. For the 0.65-micrometer band, the PDF mode statistic is preferable, providing 0.3%regional consistency and 1%temporal uncertainty over land regions. It was found that the DCC within the TWP had the lowest radiometric response and DCC over land did not necessarily have the highest radiometric response. For wavelengths greater than 1 micrometer, the mean statistic is preferred, and land regions provided a regional variability of 0.7%with a temporal uncertainty of 1.1% where the DCC land response was higher than the response over ocean. Unlike stratus and cirrus clouds, the DCC spectra were not affected by water vapor absorption.

Doelling, David R.↗

Data-Driven Dynamic Modeling in Power Systems: A Fresh Look on Inverter-Based Resource Modeling

Here, this article introduces ways to identify dynamic system models using measurement data. In power system analysis, a static model represents the time-invariant input-output relationship of a system, while a dynamic model describes the behavior of the system over time. For example, how will a system transit from one steady-state operation point to another?

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stability Assessment and Tuning of an Adaptively Augmented Classical Controller for Launch Vehicle Flight Control

Recently, a robust and practical adaptive control scheme for launch vehicles [ [1] has been introduced. It augments a classical controller with a real-time loop-gain adaptation, and it is therefore called Adaptive Augmentation Control (AAC). The loop-gain will be increased from the nominal design when the tracking error between the (filtered) output and the (filtered) command trajectory is large; whereas it will be decreased when excitation of flex or sloshing modes are detected. There is a need to determine the range and rate of the loop-gain adaptation in order to retain (exponential) stability, which is critical in vehicle operation, and to develop some theoretically based heuristic tuning methods for the adaptive law gain parameters. The classical launch vehicle flight controller design technics are based on gain-scheduling, whereby the launch vehicle dynamics model is linearized at selected operating points along the nominal tracking command trajectory, and Linear Time-Invariant (LTI) controller design techniques are employed to ensure asymptotic stability of the tracking error dynamics, typically by meeting some prescribed Gain Margin (GM) and Phase Margin (PM) specifications. The controller gains at the design points are then scheduled, tuned and sometimes interpolated to achieve good performance and stability robustness under external disturbances (e.g. winds) and structural perturbations (e.g. vehicle modeling errors). While the GM does give a bound for loop-gain variation without losing stability, it is for constant dispersions of the loop-gain because the GM is based on frequency-domain analysis, which is applicable only for LTI systems. The real-time adaptive loop-gain variation of the AAC effectively renders the closed-loop system a time-varying system, for which it is well-known that the LTI system stability criterion is neither necessary nor sufficient when applying to a Linear Time-Varying (LTV) system in a frozen-time fashion. Therefore, a generalized stability metric for time-varying loop=gain perturbations is needed for the AAC.

VanZwieten, Tannen↗

Recoil polarization in 𝐾 + ⁢𝑌 electroproduction in the nucleon resonance region with CLAS12

Hyperon recoil polarization measurements for the exclusive electroproduction of 𝐾 + ⁢Λ and 𝐾 + ⁢Σ 0 final states from an unpolarized proton target have been carried out using the CLAS12 spectrometer at Jefferson Laboratory. The measurements at beam energies of 6.535 and 7.546 GeV span the range of four-momentum transfer 𝑄 2 from 0.3 to 4.5GeV 2 and invariant mass 𝑊 from 1.6 to 2.4 GeV, while covering the full center-of-mass angular range of the 𝐾 + . These new Λ polarization observables extend the existing data in a similar kinematic range but from a significantly larger dataset. However, they represent the first electroproduction measurements of this observable for the Σ 0 . These data will allow for better exploration of the reaction mechanism in strangeness production, for further understanding of the spectrum and structure of excited nucleon states that couple to 𝐾⁢𝑌, and for improved insight into the strong interaction in the nonperturbative domain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Applications of Fractal Analytical Techniques in the Estimation of Operational Scale

The observational scale and the resolution of remotely sensed imagery are essential considerations in the interpretation process. Many atmospheric, hydrologic, and other natural and human-influenced spatial phenomena are inherently scale dependent and are governed by different physical processes at different spatial domains. This spatial and operational heterogeneity constrains the ability to compare interpretations of phenomena and processes observed in higher spatial resolution imagery to similar interpretations obtained from lower resolution imagery. This is a particularly acute problem, since longterm global change investigations will require high spatial resolution Earth Observing System (EOS), Landsat 7, or commercial satellite data to be combined with lower resolution imagery from older sensors such as Landsat TM and MSS. Fractal analysis is a useful technique for identifying the effects of scale changes on remotely sensed imagery. The fractal dimension of an image is a non-integer value between two and three which indicates the degree of complexity in the texture and shapes depicted in the image. A true fractal surface exhibits self-similarity, a property of curves or surfaces where each part is indistinguishable from the whole, or where the form of the curve or surface is invariant with respect to scale. Theoretically, if the digital numbers of a remotely sensed image resemble an ideal fractal surface, then due to the self-similarity property, the fractal dimension of the image will not vary with scale and resolution, and the slope of the fractal dimension-resolution relationship would be zero. Most geographical phenomena, however, are not self-similar at all scales, but they can be modeled by a stochastic fractal in which the scaling properties of the image exhibit patterns that can be described by statistics such as area-perimeter ratios and autocovariances. Stochastic fractal sets relax the self-similarity assumption and measure many scales and resolutions to represent the varying form of a phenomenon as the pixel size is increased in a convolution process. We have observed that for images of homogeneous land covers, the fractal dimension varies linearly with changes in resolution or pixel size over the range of past, current, and planned space-borne sensors. This relationship differs significantly in images of agricultural, urban, and forest land covers, with urban areas retaining the same level of complexity, forested areas growing smoother, and agricultural areas growing more complex as small pixels are aggregated into larger, mixed pixels. Images of scenes having a mixture of land covers have fractal dimensions that exhibit a non-linear, complex relationship to pixel size. Measuring the fractal dimension of a difference image derived from two images of the same area obtained on different dates showed that the fractal dimension increased steadily, then exhibited a sharp decrease at increasing levels of pixel aggregation. This breakpoint of the fractal dimension/resolution plot is related to the spatial domain or operational scale of the phenomenon exhibiting the predominant visible difference between the two images (in this case, mountain snow cover). The degree to which an image departs from a theoretical ideal fractal surface provides clues as to how much information is altered or lost in the processes of rescaling and rectification. The measured fractal dimension of complex, composite land covers such as urban areas also provides a useful textural index that can assist image classification of complex scenes.

Emerson, Charles W.↗

Low-latency Jet Tagging for HL-LHC Using Transformer Architectures

Transformers are the state-of-the-art model architectures and widely used in application areas of machine learning. However the performance of such architectures is less well explored in the ultra-low latency domains where deployment on FPGAs or ASICs is required. Such domains include the trigger and data acquisition systems of the LHC experiments. We present a transformer-based algorithm for jet tagging built with the HGQ2 framework, which is able to produce a model with heterogeneous bitwidths for fast inference on FPGAs, as required in the trigger systems at the LHC experiments. The bitwidths are acquired during training by minimizing the total bit operations as an additional parameter. By allowing a bitwidth of zero, the model is pruned in-situ during training. Using this quantization-aware approach, our algorithm achieves state-of-the-art performance while also retaining permutation invariance which is a key property for particle physics applications. Due to the strength of transformers in representation learning, our work also serves as a stepping stone for the development of a larger foundation model for trigger applications.

Laatu, Lauri [Imperial Coll., London]↗