Search NASA⌕ Search

SEARCH · Search NASA

Results for “satisfiability problems”

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 271 records · Page 15

Telemetry and command standards

The first phase of the international Consultative Committee for Space Data Systems (CCSDS) efforts toward the definition of standards for space telemetry, spacecraft tracking, and command functions has established a set of standard space communications techniques capable of satisfying almost the entire spectrum of space mission user requirements. This was achieved by focusing on the distinctive problems associated with the space/ground data link, and developing the infrastructural system designated the 'Open Systems Interconnection'. The intrinsically international coordination by CCSDS of development efforts ensures highly flexible mutual support activities by the various national space agencies.

Hooke, Adrian J.↗

Application of Computational Fluid Dynamics to the Study of Vortex Flow Control for the Management of Inlet Distortion

The present study demonstrates that the Reduced Navier-Stokes code RNS3D can be used very effectively to develop a vortex generator installation for the purpose of minimizing the engine face circumferential distortion by controlling the development of secondary flow. The computing times required are small enough that studies such as this are feasible within an analysis-design environment with all its constraints of time and costs. This research study also established the nature of the performance improvements that can be realized with vortex flow control, and suggests a set of aerodynamic properties (called observations) that can be used to arrive at a successful vortex generator installation design. The ultimate aim of this research is to manage inlet distortion by controlling secondary flow through an arrangements of vortex generators configurations tailored to the specific aerodynamic characteristics of the inlet duct. This study also indicated that scaling between flight and typical wind tunnel test conditions is possible only within a very narrow range of generator configurations close to an optimum installation. This paper also suggests a possible law that can be used to scale generator blade height for experimental testing, but further research in this area is needed before it can be effectively applied to practical problems. Lastly, this study indicated that vortex generator installation design for inlet ducts is more complex than simply satisfying the requirement of attached flow, it must satisfy the requirement of minimum engine face distortion.

Anderson, Bernhard H.↗

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

Monitoring the Intrinsic Absorption Complex Toward RX J1230.8+115

There is a very prominent intrinsic absorption complex in the spectrum of the quasar RX J1230.8+115. It was diagnosed as intrinsic due to partial covering of the quasar by the absorber, which can only happen for material very near to the quasar. Measuring a variability timescale for the absorption would allow us to place a more stringent constraint on the distance of the absorber from the central engine. Furthermore, understanding the variability and covering factors in different transitions can help to identify if it is due to changes in the ionizing source or transverse motion of the absorber. The OVI doublet is covered in the FUSE wavelength range. Unfortunately, our initial target, RX J1230.8+0115 could not be observed due to the failure of the first reaction wheel. We determined that another target, H 1821+643, could be observed and would satisfy similar science goals, understanding intrinsic narrow absorption line systems. Our specific goal for that target was to search for Ne VI11 absorption from the associated absorber. This observation was attempted, but the counts in the Sic channel, where the absorption was expected to lie, were too low to satisfy our science goal. This was, at least in part, due to a misalignment of the SiC mirrors. Because of this problem, more time was awarded. However, even the extra time was too little. It appears that the quasar may be too faint to properly align all the mirrors. So with disappointment, I report that we were not able to accomplish our science goals for this program. The funds were used to support my work with a graduate student, Jie Ding, to evaluate possible replacement targets once we learned that Rx J1230.8+115 would not be observable. We considered the intrinsic NAL candidates that we could find through a partial coverage analysis of doublets in the STIS E140M and E230M quasar spectra taken from the HST archive. We then checked the flux in the FUSE wavelength range, and in most cases of interest it was found to be too low to satisfy our S/N requirement in the exposure time we had been awarded. We decided that H1821+643 was the best choice for this project, even though the associated absorber we would be studying might in fact originate in the cluster with the quasar and not in intrinsic gas. Part of the appeal of the target, was in fact the opportunity to consider this issue. Please do not hesitate to contact me if more detail is needed for this report. I found it difficult to write a satisfying report under the circumstance and am uncertain if I've provided the proper level of information.

Charlton, Jane↗

Hybrid Model Based Approaches for Systems Health Management and Prognostics

This is a previously approved and published presentation. To accurately predict the future state of any system, it is required to possess knowledge of its current health state and future operational conditions. Present achievements of data-driven algorithms in regression of complex nonlinear functions and classification tasks have generated a growing interest in artificial intelligence for industrial applications. Complex multi-physics models as well as digital twins, once purely built on physics and corresponding simplified lumped parameter iterations, can now benefit from machine learning algorithms to mitigate the lack of understanding of some complex behavior. Given models of the current and future system behavior, a general approach of model-based prognostics can solve the prediction problem and further decision-making. In principle, data-driven approaches can replace expensive experimental test-setups as well as reduce the number of simulations needed to explore, e.g., the parametric space of a multi-parameter model. Nonetheless, the limitations of pure data-driven methods came to light rather quickly, at least for some industries. In many industrial applications, data acquisition is costly, and the volume of data that can be collected does not satisfy the requirements for effective model training and cross-validation. Therefore, some recent works in the area of machine learning is focusing on blending physics with data-driven algorithms, thus mitigating the drawbacks of the two approaches and emphasizing respective advantages. Partial physical knowledge of the problem can aid the learning process by “guiding” the algorithm towards efficient solutions that satisfy the physics driving the system behavior. The result is a hybrid modeling approach combining physical knowledge as well data-driven methods to develop a unified hybrid approach. A hybrid framework for fusing information from physics-based performance models along with deep learning algorithms for prognostics of complex safety-critical systems is presented. In this framework, physics-based performance models infer unobservable model parameters related to the system's components health solving a calibration problem in the deep learning approach.

Prognostics↗

Hybrid Model Based Approaches for Systems Health Management and Prognostics

To facilitate and solve the prediction problem, awareness of the current health state of the system is key, since it is necessary to perform condition-based predictions. To accurately predict the future state of any system, it is required to possess knowledge of its current health state and future operational conditions. Latest achievements of data-driven algorithms in regression of complex nonlinear functions and classification tasks have generated a growing interest in artificial intelligence for industrial applications. Complex multi-physics models as well as digital twins, once purely built on physics and corresponding simplified lumped parameter iterations, can now benefit from machine learning algorithms to mitigate the lack of understanding of some complex behavior. Given models of the current and future system behavior, a general approach of model-based prognostics can solve the prediction problem and further decision making. In principle, data driven approaches can replace expensive experimental test-setups as well as reduce the number of simulations needed to explore, e.g., the parametric space of a multi-parameter model. Nonetheless, the limitations of pure data-driven methods came to light rather quickly, at least for some industries. In many industrial applications, data acquisition is costly, and the volume of data that can be collected does not satisfy the requirements for an effective model training and cross-validation. Therefore, some recent works in the area of machine learning is focusing on blending physics with data-driven algorithms, thus mitigating the drawbacks of the two approaches and emphasizing respective advantages. Partial physical knowledge of the problem can aid the learning process by “guiding” the algorithm towards efficient solutions that satisfy the physics driving the system behavior. The result is a hybrid modeling approach combining physical knowledge as well data driven methods to develop a unified hybrid approach. A hybrid framework for fusing information from physics-based performance models along with deep learning algorithms for prognostics of complex safety critical systems is presented. In this framework, physics-based performance models infer unobservable model parameters related to the system's components health solving a calibration problem in the deep learning approach.

Hybrid Modeling↗

Hybrid Approaches to Systems Health Management and Prognostics

To facilitate and solve the prediction problem, awareness of the current health state of the system is key, since it is necessary to perform condition-based predictions. To accurately predict the future state of any system, it is required to possess knowledge of its current health state and future operational conditions. Latest achievements of data-driven algorithms in regression of complex nonlinear functions and classification tasks have generated a growing interest in artificial intelligence for industrial applications. Complex multi-physics models as well as digital twins, once purely built on physics and corresponding simplified lumped parameter iterations, can now benefit from machine learning algorithms to mitigate the lack of understanding of some complex behavior. Given models of the current and future system behavior, a general approach of model-based prognostics can solve the prediction problem and further decision making. In principle, data driven approaches can replace expensive experimental test-setups as well as reduce the number of simulations needed to explore, e.g., the parametric space of a multi-parameter model. Nonetheless, the limitations of pure data-driven methods came to light rather quickly, at least for some industries. In many industrial applications, data acquisition is costly, and the volume of data that can be collected does not satisfy the requirements for an effective model training and cross-validation. Therefore, some recent works in the area of machine learning is focusing on blending physics with data-driven algorithms, thus mitigating the drawbacks of the two approaches and emphasizing respective advantages. Partial physical knowledge of the problem can aid the learning process by “guiding” the algorithm towards efficient solutions that satisfy the physics driving the system behavior. The result is a hybrid modeling approach combining physical knowledge as well data driven methods to develop a unified hybrid approach. A hybrid framework for fusing information from physics-based performance models along with deep learning algorithms for prognostics of complex safety critical systems is presented. In this framework, physics-based performance models infer unobservable model parameters related to the system's components health solving a calibration problem in the deep learning approach.

Systems Health Management↗

An algorithm for the weighting matrices in the sampled-data optimal linear regulator problem

The sampled-data optimal linear regulator problem provides a means whereby a control designer can use an understanding of continuous optimal regulator design to produce a digital state variable feedback control law which satisfies continuous system performance specifications. A basic difficulty in applying the sampled-data regulator theory is the requirement that certain digital performance index weighting matrices, expressed as complicated functions of system matrices, be computed. Infinite series representations are presented for the weighting matrices of the time-invariant version of the optimal linear sampled-data regulator problem. Error bounds are given for estimating the effect of truncating the series expressions after a finite number of terms, and a method is described for their computer implementation. A numerical example is given to illustrate the results.

Armstrong, E. S.↗

A new solution to the accretion problem.

The Bondi-Hoyle theory of line-accretion is briefly reviewed, and it is shown that there are an infinity of steady-state solutions satisfying the requirements originally imposed as boundary conditions. Questions are raised as to the applicability of these solutions to physical reality. The existence of an entirely different type of steady-state solution that implies slow velocity in the accretion stream beyond the neutral point is demonstrated. Accurate numerical values of this solution are obtained for a number of cases. The braking action given by this solution, for any selected cut-off distance, is far stronger than in the Bondi-Hoyle solution. Consideration of the physical state in the accretion stream shows that for this slow solution the appropriate cutoff range would be only of the order of 100 AU rather than of interstellar-distance order.

Lyttleton, R. A.↗

Calculation of free fall trajectories based on numerical optimization techniques

The development of a means of computing free-fall (nonthrusting) trajectories from one specified point in the solar system to another specified point in the solar system in a given amount of time was studied. The problem is that of solving a two-point boundary value problem for which the initial slope is unknown. Two standard methods of attack exist for solving two-point boundary value problems. The first method is known as the initial value or shooting method. The second method of attack for two-point boundary value problems is to approximate the nonlinear differential equations by an appropriate linearized set. Parts of both boundary value problem solution techniques described above are used. A complete velocity history is guessed such that the corresponding position history satisfies the given boundary conditions at the appropriate times. An iterative procedure is then followed until the last guessed velocity history and the velocity history obtained from integrating the acceleration history agree to some specified tolerance everywhere along the trajectory.

Source record↗

Interplanetary plasma concentration measurements by use of dual frequencies from separated sites

In the interplanetary plasma studies considered, signals at two frequencies are sent simultaneously from a single site on the earth to a spacecraft, or vice versa. Changes in group and phase velocity are proportional to the number density of the plasma and inversely proportional to the square of the signal frequency. Plasma-induced radio timing changes can be measured through comparisons of the changes at widely different frequencies. Such a system has, however, certain weaknesses, which can be overcome by making use of an approach which uses separated sites for sources of the low and high frequencies. The practical advantages gained by site separation are balanced somewhat by several technical problems, but each problem can be solved. The main barrier appears to be a very exacting timekeeping requirement. These requirements can be satisfied with the aid of the cesium clock.

Croft, T. A.↗

Geologic utility of improved orbital measurements

From a geologists point of view, the most important use of orbital cameras is probably to provide spatial or context information in remote areas. This information can be extracted even from broad-band, single channel images provided that the spatial resolution is adequate for the problem at hand. Topographic information, usually provided by stereoscopic images, is next in importance. Any satellite system designed to satisfy general geologic requirements must produce stereo images, preferably with large base/height ratios of 0.5 or more. Satellite coverage at different illumination geometries could provide low Sun angles for regions of low relief where topography must be accentuated, and high Sun angles for areas of high relief or when multispectral data are being collected. Thus a different polar orbit than that of LANDSAT would be advantageous. The extended spectral sensitivity of scanner permits interferences of the composition of rocks in the scene to be drawn directly, rather than from contextual clues.

Gillespie, A. R.↗

Discrete approximation methods for parameter identification in delay systems

Approximation schemes for parameter identification problems in which the governing state equation is a linear functional differential equation of retarded type are constructed. The basis of the schemes is the replacement of the parameter identification problem having an infinite dimensional state equation by a sequence of approximating parameter identification problems in which the states are given by finite dimensional discrete difference equations. The difference equations are constructed using linear semigroup theory and rational function approximations to the exponential. Sufficient conditions are given for the convergence of solutions to the approximating problems, which can be obtained using conventional methods, to solutions to the original parameter identification problem. Finite difference and spline based schemes using Paderational function approximations to the exponential are constructed, and shown to satisfy the sufficient conditions for convergence. A discussion and analysis of numerical results obtained through the application of the schemes to several examples is included.

Rosen, I. G.↗

Dynamic ground effects on a two-dimensional flat plate

The effect of time-variant vortex shedding is simulated by a sequence of discrete vortices convecting downstream in the wake of a two-dimensional flat plate whose lifting condition is modeled by means of the quasi-vortex lattice method. The boundary condition of this problem is specified in such a way that the tangency condition on the surface of the flat plate is satisfied; the boundary condition also takes into account the effect of airfoil motion relative to the ground. Significant lift changes are shown to occur, due to the dynamic ground effect, that are crucial in aircraft takeoff and landing transitions.

Chen, Y.-S.↗

The effects of strain heating in lithospheric stretching models

The deformation by stretching of a continental type lithosphere has been formulated so that the problem can be solved by a continuum mechanical approach. The deformation, stress state, and temperature distribution are constrained to satisfy the physical laws of conservation of mass, energy, momentum, and an experimentally defined rheological response. The conservation of energy equation including a term of strain energy dissipation is given. The continental lithosphere is assumed to have the rheology of an isotropic, incompressible, nonlinear viscous, two layered solid.

Stanton, M.↗

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗

Analysis of spectral projectors in one-dimensional domains

A class of projection operators with values in a subspace of polynomials is analyzed. These projection operators are related to the Hilbert spaces involved in the numerical analysis of spectral methods. They are, in the first part of the paper, the standard Sobolev spaces and, in the second part, some weighted Sobolev spaces, the weight of which is related to the orthoronality relation satisfied by the Chebyshev polynomials. These results are used to study the approximation of a model fourth-order problem.

Maday, Y.↗

A Variational Formalism for the Radiative Transfer Equation: Prelude to Model 3

The MODEL III variational data assimilation model is the third of four general assimilation models designed to blend weather data measured from space based platforms in the meteorological data mainstream in a way that maximizes the information content of the satellite data. Because there are many different observation locations and there are many instruments with different measurement error characteristics, it is also necessary to require that the blending be done to maximize the information content of the data and simultaneously to retain a dynamically consistent and reasonably accurate description of the state of the atmosphere. This is ideally a variational problem for which the data receive relative weights that are inversely proportional to measurement error and are adjusted to satisfy a set of dynamical equations that govern atmospheric processes. The advantage of MODEL III over the previous two models is that radiance, the atmospheric variable measured by satellite, becomes a dependent variable. In the previous versions, mean layer temperatures that had been retrieved from the radiances by some method, were included in the assimilation by substituting them in place of the rawinsonde temperatures. Now both rawinsonde temperatures and satellite radiances are included independently in the assimilation.

Achtemeier, Gary L.↗