Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte-Carlo methods”

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 19 records

Use of the Wigner representation in scattering problems

The basic equations of quantum scattering were translated into the Wigner representation, putting quantum mechanics in the form of a stochastic process in phase space, with real valued probability distributions and source functions. The interpretative picture associated with this representation is developed and stressed and results used in applications published elsewhere are derived. The form of the integral equation for scattering as well as its multiple scattering expansion in this representation are derived. Quantum corrections to classical propagators are briefly discussed. The basic approximation used in the Monte-Carlo method is derived in a fashion which allows for future refinement and which includes bound state production. Finally, as a simple illustration of some of the formalism, scattering is treated by a bound two body problem. Simple expressions for single and double scattering contributions to total and differential cross-sections as well as for all necessary shadow corrections are obtained.

Bemler, E. A.↗

Remote sensing of chlorophyll concentration: State-of-the-art, 1975

Remote measurement of chlorophyll concentration of the world's oceans from satellite observations could potentially be extremely useful for assessments of productivity in large areas for which measurements by other means would be impractical. The basis of these measurements rests with the physics of the interaction of light with material dissolved and suspended in the water. It is theoretically possible to predict the nature of light upwelled from the ocean surface from a solution to the radiative transfer equation. Practically, however, this is difficult. Monte-Carlo methods presently are thought to be the most viable method to treat the general theoretical problem. With restrictive assumptions of the nature of scattering, it is possible to construct simpler models. Algorithms developed to relate chlorophyll concentration (or some other parameter, i.e., seechi depth) to the upwelled light spectrum are discussed.

Atwell, B. H.↗

Sensitivity analysis as an aid in modelling and control of (poorly-defined) ecological systems

A literature review of the use of sensitivity analyses in modelling nonlinear, ill-defined systems, such as ecological interactions is presented. Discussions of previous work, and a proposed scheme for generalized sensitivity analysis applicable to ill-defined systems are included. This scheme considers classes of mathematical models, problem-defining behavior, analysis procedures (especially the use of Monte-Carlo methods), sensitivity ranking of parameters, and extension to control system design.

Hornberger, G. M.↗

New results in the theory of dust grain alignment

Two complementary approaches are used in an attempt to propose an appropriate formulation of the solution to the problem of magnetic alignment of grains in the diffuse and/or the more denser clouds, whatever the mechanism of rotational excitation can be. The interest of such a unified formulation is mainly that the same theoretical expression for polarization can be used everywhere, allowing for easier comparisons between regions where the physical conditions are highly different. The first consists in applying a Monte-Carlo method to a limited number of representative cases, for which all the torques acting on the grain are taken into account: impulsive random torques due to direct collisions with gas atoms, to evaporation of atoms from the surface, and to exo-energetic recombinations forming hydrogen molecules, followed by violent ejections from peculiar sites; magnetic torques. Three characteristic times are associated with these torques: the collisional damping time, the time necessary to change completely the actual sites configuration narrowly bound to the correlation time of the suprathermal torque; and the magnetic damping time. The second approach starts from a heuristic point of view. It consists in a generalization of results (Cugnon, 1983; see also Purcell and Spitzer, 1971; Greenberg, 1978) obtained for thermal alignment to the suprathermal case. It appears indeed that in two extreme cases, the thermal formulation may be used after redefinition of involved times and temperatures.

Cugnon, Pierre↗

Typical motions in multiple systems

In very old times, people counted - one, two, many. The author wants to show that they were right. Consider the motions of isolated bodies: (1) N = 1 - simple motion; (2) N = 2 - Keplerian orbits; and (3) N = 3 - this is the difficult problem. In general, this problem can be studied only by computer simulations. The author studied this problem over many years (see, e.g., Agekian and Anosova, 1967; Anosova, 1986, 1989 a,b). The principal result is that two basic types of dynamics take place in triple systems. The first special type is the stable hierarchical systems with two almost Keplerian orbits. The second general type is the unstable triple systems with complicated motions of the bodies. By random choice of the initial conditions, by the Monte-Carlo method, the stable systems comprised about approx. 10% of the examined cases; the unstable systems comprised the other approx. 90% of cases under consideration. In N greater than 3, the studies of dynamics of such systems by computer simulations show that we have in general also the motions roughly as at the cases 1 - 3 with the relative negative or positive energies of the bodies. In the author's picture, the typical trajectories of the bodies in unstable triple systems of the general type of dynamics are seen. Such systems are disrupted always after close triple approaches of the bodies. These approaches play a role like the gravitational slingshot. Often, the velocities of escapers are very large. On the other hand, the movie also shows the dynamical processes of a formation, dynamical evolution and disruption of the temporary wide binaries in triples and a formation of final hard massive binaries in the final evolution of triples.

Anosova, Joanna P.↗

Maximum likelihood estimation for life distributions with competing failure modes

Systems which are placed on test at time zero, function for a period and die at some random time were studied. Failure may be due to one of several causes or modes. The parameters of the life distribution may depend upon the levels of various stress variables the item is subject to. Maximum likelihood estimation methods are discussed. Specific methods are reported for the smallest extreme-value distributions of life. Monte-Carlo results indicate the methods to be promising. Under appropriate conditions, the location parameters are nearly unbiased, the scale parameter is slight biased, and the asymptotic covariances are rapidly approached.

Sidik, S. M.↗

Planetary atmosphere representation for Monte-Carlo simulation of aerospace vehicle trajectories

Methods of simulating large scale random atmospheric variations are examined for use in Monte Carlo simulations of hypervelocity vehicles. Random density, temperature, and pressure variations have a significant impact on the performance of such vehicles. Methods for generating random atmosphere profiles are described and compared in the context of a Monte Carlo simulation of the Aeroassist Flight Experiment (AFE). It is seen that these techniques result in statistically significant differences in key simulation point statistics.

Dutton, K.↗

Use of Polarization Lidar for Investigation of Meteorological Formations

This paper presents the results of theoretical and experimental, investigations of depolarization characteristics of different meteorological formations. Experimental investigations are carried out with a monostatic lidar. The ruby laser radiation is polarized in a vertical plane. The radiation reflected is accepted by a lens system of 150 rom in diameter and a viewing angle of 12' and further it is divided by Wollaston prism into the components polarized orthogonally. In this case the principal plane of the prism is exposed parallel with the laser polarization plane. Investigations show the degree of radiation polarization, reflected from water clouds, to be changed within 1/0.7 (seldom up to 0.6) depending on their density. In most cases a signal reflected from the cloud leading edge is polarized completely. The time shift is observed between polarized and crosspolarized components of a Fignal, reflected from a cloud, depending on the density of a meteorological object. While penetrating into the cloud depth a degree of polarization decreases up to 0.8-0.7, and the character of this decrease is different for various types of clouds. For crystal clouds the shift between the components of the reflected signal is not observed and the magnitude of polarization degree amounts to 0.1/0.3 in comparison with water clouds. The polarization degree of radiation reflected by fog is not less than 0.6, and that in the rains of average intensity (about 5 mm/h) is always about 1. The authors have suggested an algorithm of numerical solution of nonstationary transfer equation in the vector form to forecast the influence of multiple scattering effects on polarization characteristics of the lidar light signal. The method of statistical simulation (Monte-Carlo technique) forms the basis of the algorithm. Numerical estimates obtained for a model of stratocumulus at lambda = 0.6943 microns under boundary conditions close to the conditions of natural experiment being discussed proved to be in a good agreement with the results of observation. Specifically, Fig. 1 shows the profiles of polarization (p) versus depth (L) of the following drop formations: fog (curve 1) with horizontal meteorological visibility of 4 km two stratocumulus at a height of 1100 m with the attenuation factors delta = 0.01m(exp -1) (curve 2) and delta = 0.05m(exp -1) (curve 3). Curve 3 shows the results of numerical estimates and the value of their statistical error.

Balin, Yu. S.↗

Fault-tolerant system optimization

The paper describes the decisions to be made in the design of fault-tolerant systems and provides details of a comprehensive model developed to cost optimize such systems. Economical use of replication is making fault-tolerant systems possible and more applications for safety crucial systems such as active flight controls can be expected. In turn, the use of massive redundancy, fault-tolerance, and reconfigurable systems in stimulating the development of new analytical tools for establishing the cost and effectiveness of the safety and cost effectiveness of the levels of replication will increase. Closed-form analytical solutions for the reliability and maintainability analysis of fault-tolerant systems are complex, and Monte-Carlo simulation appears to be a more desirable method of establishing the reliability and maintainability of such systems.

Rose, J.↗

Evaluation of concepts and methods for extraction of flow parameters from laser transit anemometry data

A Monte-Carlo type simulation system was developed to test concepts and associated algorithms proposed by Mayo and Smart (1984) for the extraction of two-dimensional flow parameters from laser transit anemometry data sets. The Monte-Carlo LTA simulation program models a flow field as well as the sample volume geometry of the LTA system, thereby permitting simulation of two-dimensional data acquisition in a three-dimensional flow field. The results of these tests on the proposed processing concepts suggest that the instrument is capable of measuring mean velocities to less than +0.15 percent error and flow angles to less then +0.3 percent for turbulence intensities of up to 15.0 percent.

Humphreys, W. M., Jr.↗

Transverse momentum distribution of pion deg in the fragmentation region of super high energy interactions

The lateral distribution of gamma-families observed by emulsion chamber is sensitive to test transverse momentum of high energy interaction. But most part gamma-families are succesive interation's results which makes it necessary to analyze the propagation of gamma-ray in atmosphere. A gamma-ray with energy Er and transverse momentum Pt is produced at the altitude h. In the emulsion chamber experiment, Eob is estimated by decascade method. There are 30 gamma-families observed by Ganbala emulsion Chamber. A group of Monte-Carlo simulation gamma-families are used to compare with experimental data with the same treatment. It is found that both distributions are consistent.

Ren, J. R.↗

Design, Development, and Test of the Advanced Apollo Orbital Assembly System

As spaceflight moves toward commercial solutions for Crew Vehicles and Space Stations, opportunity exists to lower costs with novel designs. Probe and cone docking systems provide a lightweight, low cost, and high-performance docking solution. This work revisits the Apollo probe and cone design and modifies it for the requirements of today’s computer-controlled spacecraft. This new system is called the Advanced Apollo Orbital Assembly (APOA) system, and is intended to support satellite, cargo, and space station assembly missions. A simulation of the APOA was programmed, components were sized, and a design was created by scaling the probe to the size of an EELV Secondary Payload Adapter (ESPA) tunnel. Prototype test hardware was 3D printed using Fused Deposition Modeling (FDM) methods with Polylactic Acid (PLA) material. Testing of the APOA-ESPA was conducted at Marshall Space Flight Center’s (MSFC) Flat Floor, and a test-correlated simulation is used to evaluate a Monte-Carlo of Initial Contact Conditions (ICC’s) to establish baseline performance. The successful development, test, and correlation of the APOA-ESPA proves the design validity and increases Technology Readiness Level (TRL) from 2 to 4. This work opens the door to construction of an APOA-ESPA from flight like materials, and to develop a larger scale prototype APOA. When the full scale APOA is incorporated with the Common Berthing Mechanism (CBM), becomes the Hybrid Berthing System (HBS), which allows for berthing without a robotic arm.

Berthing↗

Translating expert system rules into Ada code with validation and verification

The purpose of this ongoing research and development program is to develop software tools which enable the rapid development, upgrading, and maintenance of embedded real-time artificial intelligence systems. The goals of this phase of the research were to investigate the feasibility of developing software tools which automatically translate expert system rules into Ada code and develop methods for performing validation and verification testing of the resultant expert system. A prototype system was demonstrated which automatically translated rules from an Air Force expert system was demonstrated which detected errors in the execution of the resultant system. The method and prototype tools for converting AI representations into Ada code by converting the rules into Ada code modules and then linking them with an Activation Framework based run-time environment to form an executable load module are discussed. This method is based upon the use of Evidence Flow Graphs which are a data flow representation for intelligent systems. The development of prototype test generation and evaluation software which was used to test the resultant code is discussed. This testing was performed automatically using Monte-Carlo techniques based upon a constraint based description of the required performance for the system.

Becker, Lee↗

Performance evaluation of integrated guidance and navigation systems

The method of linear covariance analysis is applied to the study of an aerospace vehicle integrated guidance-and-navigation system. This requires that the rms trajectory deviations of the vehicle be computed along with the rms errors of the navigation system. The technique is used to examine the performance of the shuttle orbiter vehicle's guidance-and-navigation system during the terminal phase of the landing. Essentially the same rms trajectory deviations were obtained from a single covariance run as from a 30-sample Monte-Carlo run set. Of particular interest in the study was the result that rms trajectory deviations for the integrated-system were generally significantly larger than the rms navigation errors.

Tao, Y. C.↗

Monte-Carlo modeling of polar wind photoelectron distributions with anomalous heat flux

In situ measurements above the polar cap by the DE satellites show asymmetric field-aligned electron velocity distributions in the photoelectron energy range 5-60 eV. The formation of these anisotropic distributions is described by a collisional kinetic equation that includes the effects of Coulomb collisions (via a Fokker-Planck operator), those of a parallel electric field, and the magnetic mirror effects of the geomagnetic field. This equation is solved using a Monte Carlo simulation method. The simulation illustrates the formation of the portion of the observed electron distributions originating in the ionosphere and makes it possible to evaluate the significance of the suprathermal electrons at higher altitude. Both observations and calculations indicate that these electrons possess an energy flux which, by its effect on the polar-wind electric field, may influence the dynamics of the polar-wind outflow.

Yasseen, F.↗