Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Saturation problem”

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 109 records · Page 6

Integration of airborne Thematic Mapper Simulator (TMS) data and digitized aerial photography via an ISH transformation

A simple method for enhancing the spatial and spectral resolution of disparate data sets is presented. Two data sets, digitized aerial photography at a nominal spatial resolution 3,7 meters and TMS digital data at 24.6 meters, were coregistered through a bilinear interpolation to solve the problem of blocky pixel groups resulting from rectification expansion. The two data sets were then subjected to intensity-saturation-hue (ISH) transformations in order to 'blend' the high-spatial-resolution (3.7 m) digitized RC-10 photography with the high spectral (12-bands) and lower spatial (24.6 m) resolution TMS digital data. The resultant merged products make it possible to perform large-scale mapping, ease photointerpretation, and can be derived for any of the 12 available TMS spectral bands.

Ambrosia, Vincent G.↗

Removing Spectral Diagnostics of Galactic and Stellar X-Ray Emission from Charged Exchange Recombination

Our research uses the electron beam ion trap (EBIT) at the Lawrence Livermore National Laboratory to study X-ray emission from the charge exchange (CX) of highly charged ions with neutral gases. The resulting data help to fill a void in existing experimental and theoretical understanding of this atomic physics process, and are needed to explain all or part of the observed X-ray emission from the soft X-ray background, stellar winds, the Galactic Center and Galactic Ridge, supernova ejecta, and photoionized nebulae. Appreciation of the astrophysical relevance of our work continues to grow with the publication of roughly a dozen papers in the past four years describing Chandra and XMM observations of geocoronal and heliospheric CX emission, the temporal variation of such emission and correlation with X-ray emission enhancements observed by ROSAT, the theoretical spatial distribution of that emission, and CX emission around other stars. A similar number of papers were also published during that time describing CX emission from planets and comets. We expect that the launch of ASTRSE2, with its second-generation XRS microcalo- (with 6-eV resolution), will reveal even more clearly the contributions of CX to astrophysical emission. In our EBIT work we collected CX spectra from such ions as H-like and He-like Ne, Ar, and Fe. Our early measurements were made with a high-purity Ge detector, but during the second year we began operation of the first-generation XRS microcalorimeter (a twin of the XRS on ASTRO-E) and greatly improved the resolution of our measurements from roughly 150 eV (FWHM) with the Ge detectors to 10 eV with the XRS. We found that saturation of the XRS counting apparatus, which we described in our proposal as a potential concern, is not a problem for studying CX. During the course of our research, we expanded the number of injection gases permitted by the LLNL safety team, purchased and eventually operated an atomic H source, and clearly demonstrated the feasibility of our longer-range plan. For example, we successfully injected He into EBIT (not a small feat because of the difficulty of maintaining a good vacuum with He and avoiding electrical breakdown) to collect a H-like oxygen CX spectrum. The highest energy CX spectrum recorded with the XRS to date is that of the Ar K-shell emission. These measurements provided the first observation of the relative intensity ratios of resolved He-like singlet and triplet n=2->1 lines. We also carried out measurements of He-like Ne as a function of collision energy (i.e., ion temperature). Significant differences in the resulting x-ray spectra were noted. In all cases, the intensity of high-n H-like Lyman lines is significantly higher than current theoretical CX models predict.

Wargelin, Brad↗

Lunar Impact Basins: Stratigraphy, Sequence and Ages from Superposed Impact Crater Populations Measured from Lunar Orbiter Laser Altimeter (LOLA) Data

Impact basin formation is a fundamental process in the evolution of the Moon and records the history of impactors in the early solar system. In order to assess the stratigraphy, sequence, and ages of impact basins and the impactor population as a function of time, we have used topography from the Lunar Orbiter Laser Altimeter (LOLA) on the Lunar Reconnaissance Orbiter (LRO) to measure the superposed impact crater size-frequency distributions for 30 lunar basins (D ≥ 300 km). These data generally support the widely used Wilhelms sequence of lunar basins, although we find significantly higher densities of superposed craters on many lunar basins than derived by Wilhelms (50% higher densities). Our data also provide new insight into the timing of the transition between distinct crater populations characteristic of ancient and young lunar terrains. The transition from a lunar impact flux dominated by Population 1 to Population 2 occurred before the mid-Nectarian. This is before the end of the period of rapid cratering, and potentially before the end of the hypothesized Late Heavy Bombardment. LOLA-derived crater densities also suggest that many Pre-Nectarian basins, such as South Pole-Aitken, have been cratered to saturation equilibrium. Finally, both crater counts and stratigraphic observations based on LOLA data are applicable to specific basin stratigraphic problems of interest; for example, using these data, we suggest that Serenitatis is older than Nectaris, and Humboldtianum is younger than Crisium. Sample return missions to specific basins can anchor these measurements to a Pre-Imbrian absolute chronology.

Fassett, C. I.↗

Design challenges encountered in the F-15 PCA flight test program

The NASA Dryden Flight Research Center conducted flight tests of a propulsion-controlled aircraft system on an F-15 airplane. This system was designed to explore the feasibility of providing safe emergency landing capability using only the engines to provide flight control in the event of a catastrophic loss of conventional flight controls. Control laws were designed to control the flight path and bank angle using only commands to the throttles. While the program was highly successful, this paper concentrates on the challenges encountered using engine thrust as the only control effector. Compared to conventional flight control surfaces, the engines are slow, nonlinear, and have limited control effectiveness. This increases the vulnerability of the system to outside disturbances and changes in aerodynamic conditions. As a result, the PCA system had problems with gust rejection. Cross coupling of the longitudinal and lateral axis also occured, primarily as a result of control saturation. The normally negligible effects of inlet airframe interactions became significant with the engines as the control effector. Flight and simulation data are used to illustrate these difficulties.

Maine, Trindel A.↗

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Method for Exploiting Redundancy to Accommodate Actuator Limits in Multivariable Systems

This paper introduces a new method for accommodating actuator saturation in a multivariable system with actuator redundancy. Actuator saturation can cause significant deterioration in control system performance because unmet demand may result in sluggish transients and oscillations in response to setpoint changes. To help compensate for this problem, a technique has been developed which takes advantage of redundancy in multivariable systems to redistribute the unmet control demand over the remaining useful effectors. This method is not a redesign procedure, rather it modifies commands to the unlimited effectors to compensate for those which are limited, thereby exploiting the built-in redundancy. The original commands are modified by the increments due to unmet demand, but when a saturated effector comes off its limit, the incremental commands disappear and the original unmodified controller remains intact. This scheme provides a smooth transition between saturated and unsaturated modes as it divides up the unmet requirement over any available actuators. This way, if there is sufficiently redundant control authority, performance can be maintained.

Litt, Jonathan↗

Circuit complexity and functionality: A statistical thermodynamics perspective

Circuit complexity, defined as the minimum circuit size required for implementing a particular Boolean computation, is a foundational concept in computer science. Determining circuit complexity is believed to be a hard computational problem. Recently, in the context of black holes, circuit complexity has been promoted to a physical property, wherein the growth of complexity is reflected in the time evolution of the Einstein-Rosen bridge (“wormhole”) connecting the two sides of an anti-de Sitter “eternal” black hole. Here, we are motivated by an independent set of considerations and explore links between complexity and thermodynamics for functionally equivalent circuits, making the physics-inspired approach relevant to real computational problems, for which functionality is the key element of interest. In particular, our thermodynamic framework provides an alternative perspective on the obfuscation of programs of arbitrary length—an important problem in cryptography—as thermalization through recursive mixing of neighboring sections of a circuit, which can be viewed as the mixing of two containers with “gases of gates.” This recursive process equilibrates the average complexity and leads to the saturation of the circuit entropy, while preserving functionality of the overall circuit. The thermodynamic arguments hinge on ergodicity in the space of circuits which we conjecture is limited to disconnected ergodic sectors due to fragmentation. The notion of fragmentation has important implications for the problem of circuit obfuscation as it implies that there are circuits of same size and functionality that cannot be connected via a polynomial number of local moves. Furthermore, we argue that fragmentation is unavoidable unless the complexity classes NP and coNP coincide, a statement that implies the collapse of the polynomial hierarchy of computational complexity theory to its first level.

Science & Technology - Other Topics↗

Theory of wide-angle photometry from standard stars

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

Usher, Peter D.↗

The R-matrix bootstrap for the 2d O( N ) bosonic model with a boundary

The S-matrix bootstrap is extended to a 1+1d theory with O(N) symmetry and a boundary in what we call the R-matrix bootstrap since the quantity of interest is the reflection matrix (R-matrix). Given a bulk S-matrix, the space of allowed R-matrices is an infinite dimensional convex space from which we plot a two dimensional section given by a convex domain on a 2d plane. In certain cases, at the boundary of the domain, we find vertices corresponding to integrable R-matrices with no free parameters. In other cases, when there is a one-parameter family of integrable R-matrices, the whole boundary represents integrable theories. We also consider R-matrices which are analytic in an extended region beyond the physical cuts, thus forbidding poles (resonances) in that region. In certain models, this drastically reduces the allowed space of R-matrices leading to new vertices that again correspond to integrable theories. We also work out the dual problem, in particular in the case of extended analyticity, the dual function has cuts on the physical line whenever unitarity is saturated. For the periodic Yang-Baxter solution that has zero transmission, we computed the R-matrix initially using the bootstrap and then derived its previously unknown analytic form.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Assessing the difficulty of capturing the distribution function of neutrinos in neutron star merger simulations

The collision of two neutron stars is a rich source of information about nuclear physics. In particular, the kilonova signal following a merger can help us elucidate the role of neutron stars in nucleosynthesis, and informs us about the properties of matter above nuclear saturation. Approximate modeling of neutrinos remains an important limitation to our ability to make predictions for these observables. Part of the problem is the fermionic nature of neutrinos. By the exclusion principle, the expected value 𝑓 𝜈 for the number of neutrinos in a quantum state is at most 1. Any process producing neutrinos is suppressed by a blocking factor (1 −𝑓 𝜈 ). Recent simulations focused on neutrino physics mostly use a gray two-moment scheme to evolve neutrinos. This evolves integrals of 𝑓 𝜈 over momentum space, preventing direct calculations of blocking factors. Monte Carlo methods may be an attractive alternative, providing access to the full distribution of neutrinos. Their current implementation is, however, inadequate to estimate 𝑓 𝜈 : in our most recent simulations, a single Monte Carlo packet causes, in the worst cases, estimates of 𝑓 𝜈 to jump from 𝑓 𝜈 =0 to 𝑓 𝜈 ∼10 5 . While this is concerning, this brazen violation of the fermionic nature of neutrinos has been largely inconsequential, as the interactions used in simulations avoid direct calculations of 𝑓 𝜈 . We are, however, reaching a level of modeling at which this problem can no longer be ignored. Here, we discuss the relatively simple origin of this issue. We then show that very rough estimates of 𝑓 𝜈 can in theory be obtained in merger simulations, but that they will require a combination of unintuitive weighting schemes for Monte Carlo packets and smoothing of the neutrino distribution at coarser resolution than what the merger simulation uses.

79 ASTRONOMY AND ASTROPHYSICS↗

Lower Bounds on Quantum Annealing Times

The adiabatic theorem provides sufficient conditions for the time needed to prepare a target ground state. While it is possible to prepare a target state much faster with more general quantum annealing protocols, rigorous results beyond the adiabatic regime are rare. Here, we provide such a result, deriving lower bounds on the time needed to successfully perform quantum annealing. The bounds are asymptotically saturated by three toy models where fast annealing schedules are known: the Roland and Cerf unstructured search model, the Hamming spike problem, and the ferromagnetic p-spin model. Our bounds demonstrate that these schedules have optimal scaling. Herein, our results also show that rapid annealing requires coherent superpositions of energy eigenstates, singling out quantum coherence as a computational resource.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Integrated sample-handling and mounting system for fixed-target serial synchrotron crystallography

Serial synchrotron crystallography (SSX) is enabling the efficient use of small crystals for structure–function studies of biomolecules and for drug discovery. An integrated SSX system has been developed comprising ultralow background-scatter sample holders suitable for room and cryogenic temperature crystallographic data collection, a sample-loading station and a humid `gloveless' glovebox. The sample holders incorporate thin-film supports with a variety of designs optimized for different crystal-loading challenges. These holders facilitate the dispersion of crystals and the removal of excess liquid, can be cooled at extremely high rates, generate little background scatter, allow data collection over >90° of oscillation without obstruction or the risk of generating saturating Bragg peaks, are compatible with existing infrastructure for high-throughput cryocrystallography and are reusable. The sample-loading station allows sample preparation and loading onto the support film, the application of time-varying suction for optimal removal of excess liquid, crystal repositioning and cryoprotection, and the application of sealing films for room-temperature data collection, all in a controlled-humidity environment. The humid glovebox allows microscope observation of the sample-loading station and crystallization trays while maintaining near-saturating humidities that further minimize the risks of sample dehydration and damage, and maximize working times. This integrated system addresses common problems in obtaining properly dispersed, properly hydrated and isomorphous microcrystals for fixed-orientation and oscillation data collection. Its ease of use, flexibility and optimized performance make it attractive not just for SSX but also for single-crystal and few-crystal data collection. Fundamental concepts that are important in achieving desired crystal distributions on a sample holder via time-varying suction-induced liquid flows are also discussed.

47 OTHER INSTRUMENTATION↗

Multiobjective Dynamic Optimization for Optimal Load-Following of Natural Gas Combined Cycle Power Plants under Stress Constraints

With the increasing penetration of intermittent renewable energy sources into the electric grid, there is an associated need for conventional thermal power plants that are designed to operate at base-load conditions to cycle their load rapidly and frequently and operate under low-load conditions. Fast load-following operations lead to plant efficiency loss and the reduction of equipment life. A high-fidelity dynamic model of a natural gas combined cycle power plant, with rigorous equipment level submodels, is developed to capture the plant transient performance and the thermomechanical stress evolution at the high-pressure drum is computed to assess the drum life consumption. Stress evolution is modeled for the location of the edge at the drum/downcomer junction that experiences high circumferential stress amplitude during the fast load-following operation, leading to higher fatigue damage than locations considered under existing design standards. A dynamic optimization problem is solved for optimal load-following operation. Depending on the value of the stress constraint and the desired ramp rate, satisfying the desired stress constraint may be infeasible without relaxing the ramp rate. A multiobjective dynamic optimization problem is solved using a lexicographic approach that minimizes ramp rate relaxation and maximizes efficiency while satisfying stress constraints, avoiding spraying to saturation, and maintaining main steam and reheat steam temperature within bounds. It was observed that the optimal ramp rate can be highly nonlinear as opposed to the industry-standard linear ramp rate. Nonlinear optimal ramp-rate profiles not only help to avoid stress constraints that may be unavoidable by using the linear profile but also result in higher efficiency than the linear profile. Furthermore, the study shows that there is a strong tradeoff between the relaxation in the plant ramp rate and the time-average thermal efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GSoFa: Scalable Sparse Symbolic LU Factorization on GPUs

Decomposing a matrix $\mathbf {A}$ into a lower matrix $\mathbf {L}$ and an upper matrix $\mathbf {U}$, which is also known as LU decomposition, is an essential operation in numerical linear algebra. For a sparse matrix, LU decomposition often introduces more nonzero entries in the $\mathbf {L}$ and $\mathbf {U}$ factors than in the original matrix. A symbolic factorization step is needed to identify the nonzero structures of $\mathbf {L}$ and $\mathbf {U}$ matrices. Attracted by the enormous potentials of the Graphics Processing Units (GPUs), an array of efforts have surged to deploy various LU factorization steps except for the symbolic factorization, to the best of our knowledge, on GPUs. This article introduces gSoFa, the first GPU-based symbolic factorization design with the following three optimizations to enable scalable LU symbolic factorization for nonsymmetric pattern sparse matrices on GPUs. First, here we introduce a novel fine-grained parallel symbolic factorization algorithm that is well suited for the Single Instruction Multiple Thread (SIMT) architecture of GPUs. Second, we tailor supernode detection into a SIMT friendly process and strive to balance the workload, minimize the communication and saturate the GPU computing resources during supernode detection. Third, we introduce a three-pronged optimization to reduce the excessive space consumption problem faced by multi-source concurrent symbolic factorization. Taken together, gSoFa achieves up to 31× speedup from 1 to 44 Summit nodes (6 to 264 GPUs) and outperforms the state-of-the-art CPU project, on average, by 5×. Notably, gSoFa also achieves up to 47 percent of the peak memory throughput of a V100 GPU in the Summit Supercomputer.

97 MATHEMATICS AND COMPUTING↗

Factors affecting snow assessment from LANDSAT data

Problems associated with using LANDSAT as a snow monitoring satellite are studied. Data cover problems of slop hinderance in thematic mapping of snow and detection of snow in forested areas. It was concluded that if detector saturation threshold is raised and the upper spectral limit of the MSS is extended, snow extent mapping will be improved along with the likelihood of monitoring some aspects of the physical condition of snow pack.

Mcginnis, D. F., Jr.↗

Output statistics of laser anemometers in sparsely seeded flows

It is noted that until very recently, research on this topic concentrated on the particle arrival statistics and the influence of the optical parameters on them. Little attention has been paid to the influence of subsequent processing on the measurement statistics. There is also controversy over whether the effects of the particle statistics can be measured. It is shown here that some of the confusion derives from a lack of understanding of the experimental parameters that are to be controlled or known. A rigorous framework is presented for examining the measurement statistics of such systems. To provide examples, two problems are then addressed. The first has to do with a sample and hold processor, the second with what is called a saturable processor. The sample and hold processor converts the output to a continuous signal by holding the last reading until a new one is obtained. The saturable system is one where the maximum processable rate is arrived at by the dead time of some unit in the system. At high particle rates, the processed rate is determined through the dead time.

Edwards, R. V.↗

Maturing Electron Multiplying Charge Coupled Device Photon-Counting with Variable Multiplication Gain Imaging for a Coronagraphic Instrument

This is an introduction to a US Government Program that conducted high-contrast imaging experiments with an Electron Multiplying Charge Coupled Device (EMCCD) in an interferometric coronagraph. This report will introduce the concepts of “Charge Blooming” and “Starlight Saturation” in the context of high-contrast astronomical imaging. These phenomena adversely affect the performance of high-contrast photon-counting instruments that do not use a mask to physically block starlight in the science channel of the coronagraph. The problems will be presented with the help of images taken with a commercial EMCCD camera in the Visible Nulling Coronagraph at the Goddard Space Flight Center (GSFC). A new clocking scheme for EMCCDs – Variable Multiplication Gain Clocking – will be proposed as a means for suppressing horizontal blooming and starlight saturation in an astronomical camera. This opening report from the program will conclude with an introduction to a new controller for high-contrast imaging with EMCCDs in coronagraphs. This controller is being designed to allow a single frame from an EMCCD to be scanned in multiple modes – photon counting and digitization – to enable direct imaging of an exo-planet and wavefront control of a coronagraph, simultaneously.

Udayan Mallik↗