Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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 433 records · Page 24

Numerical Descriptions of Cosmic-Ray Transport

The behavior of energetic particles in the solar system is described by a well known Fokker-Planck equation. Although analytic methods yield insight into the nature of its solutions, especially in the diffusion regime, calculations that go beyond diffusion are very complicated. The reliability of these calculations is of concern, because numerical methods are notorious for their errors and artifacts. The well known Milne problem of classical transport theory was analyzed with the aid of three different numerical methods. These are: (1) The method of eigenfunctions in which the distribution function is approximated by a sum of eigenfunctions of the scattering operator, (2) Numerical solutions of a finite difference aquation; and (3) Direct simulation of the scattering and streaming of individual particles with the aid of Monte Carlo methods.

Earl, J. A.↗

Effects of Raman scattering on the water-leaving radiance

The contribution of Raman scattering to the water-leaving radiance is examined using Monte Carlo simulations. Exit angle information is retained, allowing a comparison of different satellite viewing directions. Chlorophyll values of 0.0, 0.01, 0.1, and 1.0 mg Chl/cu m are simulated. Little directional variability is found, with the exception of the direct solar backscatter direction. The wavelength variability is greatest for low chlorophyll concentrations and is negligible for 1.0 mg Chl/cu m. At 550 nm the Raman contribution ranges from approximately 18% of the total water-leaving radiance for pure water to 3% for 1.0 mg Chl/cu m. At 440 nm the range is from 6% to 2%, indicating that Raman scattering will impact radiance ratios for ocean color satellite algorithms.

Waters, Kirk J.↗

Covariance Analysis Tool (G-CAT) for Computing Ascent, Descent, and Landing Errors

G-CAT is a covariance analysis tool that enables fast and accurate computation of error ellipses for descent, landing, ascent, and rendezvous scenarios, and quantifies knowledge error contributions needed for error budgeting purposes. Because GCAT supports hardware/system trade studies in spacecraft and mission design, it is useful in both early and late mission/ proposal phases where Monte Carlo simulation capability is not mature, Monte Carlo simulation takes too long to run, and/or there is a need to perform multiple parametric system design trades that would require an unwieldy number of Monte Carlo runs. G-CAT is formulated as a variable-order square-root linearized Kalman filter (LKF), typically using over 120 filter states. An important property of G-CAT is that it is based on a 6-DOF (degrees of freedom) formulation that completely captures the combined effects of both attitude and translation errors on the propagated trajectories. This ensures its accuracy for guidance, navigation, and control (GN&C) analysis. G-CAT provides the desired fast turnaround analysis needed for error budgeting in support of mission concept formulations, design trade studies, and proposal development efforts. The main usefulness of a covariance analysis tool such as G-CAT is its ability to calculate the performance envelope directly from a single run. This is in sharp contrast to running thousands of simulations to obtain similar information using Monte Carlo methods. It does this by propagating the "statistics" of the overall design, rather than simulating individual trajectories. G-CAT supports applications to lunar, planetary, and small body missions. It characterizes onboard knowledge propagation errors associated with inertial measurement unit (IMU) errors (gyro and accelerometer), gravity errors/dispersions (spherical harmonics, masscons), and radar errors (multiple altimeter beams, multiple Doppler velocimeter beams). G-CAT is a standalone MATLAB- based tool intended to run on any engineer's desktop computer.

Boussalis, Dhemetrios↗

Monte Carlo modeling of atomic oxygen attack of polymers with protective coatings on LDEF

Characterization of the behavior of atomic oxygen interaction with materials on the Long Duration Exposure Facility (LDEF) will assist in understanding the mechanisms involved, and will lead to improved reliability in predicting in-space durability of materials based on ground laboratory testing. A computational simulation of atomic oxygen interaction with protected polymers was developed using Monte Carlo techniques. Through the use of assumed mechanistic behavior of atomic oxygen and results of both ground laboratory and LDEF data, a predictive Monte Carlo model was developed which simulates the oxidation processes that occur on polymers with applied protective coatings that have defects. The use of high atomic oxygen fluence-directed ram LDEF results has enabled mechanistic implications to be made by adjusting Monte Carlo modeling assumptions to match observed results based on scanning electron microscopy. Modeling assumptions, implications, and predictions are presented, along with comparison of observed ground laboratory and LDEF results.

Banks, Bruce A.↗

Determination of Barometric Altimeter Errors for the Orion Exploration Flight Test-1 Entry

The Exploration Flight Test 1 (EFT-1) mission is the unmanned flight test for the upcoming Multi-Purpose Crew Vehicle (MPCV). During entry, the EFT-1 vehicle will trigger several Landing and Recovery System (LRS) events, such as parachute deployment, based on on-board altitude information. The primary altitude source is the filtered navigation solution updated with GPS measurement data. The vehicle also has three barometric altimeters that will be used to measure atmospheric pressure during entry. In the event that GPS data is not available during entry, the altitude derived from the barometric altimeter pressure will be used to trigger chute deployment for the drogues and main parachutes. Therefore it is important to understand the impact of error sources on the pressure measured by the barometric altimeters and on the altitude derived from that pressure. The error sources for the barometric altimeters are not independent, and many error sources result in bias in a specific direction. Therefore conventional error budget methods could not be applied. Instead, high fidelity Monte-Carlo simulation was performed and error bounds were determined based on the results of this analysis. Aerodynamic errors were the largest single contributor to the error budget for the barometric altimeters. The large errors drove a change to the altitude trigger setpoint for FBC jettison deploy.

Brown, Denise L.↗

Reflectance of broken cloud fields - Simulation and parameterization

The transfer of solar irradiance in plane parallel and broken cloud fields is simulated using a Monte Carlo method. The angular distribution pattern of radiances exiting the cloud layer is studied with varying cloud geometries, optical thicknesses, cloudiness, and solar zenith angles. A rather large anisotropy of the reflected flux is found, usually increasing with solar zenith angle and with patterns that strongly depend on cloud geometry. The main features are (1) a local maximum of reflected intensity in the forward direction for all cases, (2) a limb darkening for the plane parallel case, and (3) a limb brightening and a local maximum of reflected intensity in the backward direction for broken clouds. A parameterization for the azimuth-average reflectance function is developed. It reproduces the Monte Carlo simulation with a reasonable accuracy and illustrates that, when azimuthally averaged, the reflectance function is dominated by side viewing and intercloud shadowing effects.

Breon, Francois-Marie↗

Implementation and Testing of the First Prompt Search for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (OW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient OW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize OW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 172009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify OW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency OW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with approx.50% or better probability with a few pointings of wide-field telescopes. Key words. gravitational waves - methods: observational

Camp, Jordan↗

Implementation and Testing of the First Prompt for Electromagnetic Counterparts to Gravitational Wave Transients

A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW-event candidates and to reconstruct-maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline s ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with 50% or better probability with a few pointings of wide-field telescopes.

Abadie, J.↗

Implementation and Testing of the First Prompt Search for Gravitational Wave Transients with Electromagnetic Counterparts

Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec. 17, 2009 to Jan. 8, 2010 and Sep. 2 to Oct. 20, 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and MilkyWay globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ∼50% or better probability with a few pointings of wide-field telescopes.

testing↗

CFD 2030 Grand Challenge: CFD-in-the-Loop Monte Carlo Flight Simulation for Space Vehicle Design

Flight qualification of space vehicles is markedly different from those typically employed for aircraft. The concept of an extensive flight test campaign for a space vehicle does not exist, and vehicle designers must look to alternative techniques for demonstrating robust and reliable performance of their vehicles prior to operational flight. A space vehicle may undergo only a handful of flight tests in its development cycle, with each flight representing a drastically different flight phase or flight configuration. For instance, NASA’s Space Launch System (SLS) launch vehicle and Orion spacecraft will only see a total of four flight demonstrations before flying a crew on its first operational mission, and each flight demonstrates a unique vehicle configuration and/or set of flight conditions. The SLS will be flown only one time before it becomes operational (Artemis 1). The Orion spacecraft Crew Module (CM) will have been tested twice, once on a Delta IV launch vehicle (Exploration Flight Test 1) and once as a fully integrated system with the SLS launch vehicle (Artemis 1). The Orion Launch abort system will have been tested twice, once in a pad abort scenario (Pad Abort 1) and once in an inflight abort scenario (Ascent Abort 2) on a modified Peacekeeper booster. Both of these latter tests involve only a boiler plate CM, not a functional Orion spacecraft. Thus, unlike aircraft, there is very little opportunity for engineers to assess and evaluate their preflight predictions. Instead, space vehicle designers rely on Monte Carlo flight simulations with detailed dispersions of predicted nominal flight behavior to determine how robust their design is to errors and uncertainties in the flight conditions their vehicle may encounter. These Monte Carlo analyses entail thousands of trajectory simulations to demonstrate that the vehicle can meet design requirements at a specified level of reliability. From an aerodynamics and aerothermodynamics perspective, these trajectory simulations are fueled by an extensive aerodynamic database that covers the complete range of expected flight conditions, vehicle configurations, and flight attitudes expected in a given mission. Today, these databases amount to a table of engineering parameters that can be quickly interrogated by the trajectory simulator. The aerodynamic and aerothermodynamic databases are assembled via a series of ground tests, empirical and analytical analysis, physics-based computational analysis, applicable past flight performance data, and in some cases, engineering judgment. These databases generally take years to assemble for a new space vehicle system and in the case of SLS/Orion, over a decade of test and analysis have been expended to develop the extensive databases required to cover the myriad of configurations and potential flight conditions required for the system. Recently, it has been proposed that Computational Fluid Dynamic (CFD) and computing capability may be reaching a point where it is foreseeable that CFD could be integrated directly into the production trajectory simulation tools used to design NASA’s space vehicles. To demonstrate this, NASA has embarked on two demonstrations of this type of capability, one where six degree of freedom flight trajectory simulation equations are embedded in an existing CFD solver and another where a production CFD solver is loosely coupled with a production trajectory simulation tool. These efforts represent an initial demonstration of a future approach to flight trajectory simulation, but they are a far cry from the capability required to perform a full-up CFD-in-the-loop Monte Carlo trajectory simulation. Therefore, this represents a viable grand challenge for computational methods addressing space vehicle design and development. The final paper/presentation will discuss the many hurdles, beyond simply raw computational power, to realizing this grand challenge and how they map directly to the CFD Vision 2030 ojectives. Among these are the wide range of flight conditions, including accelerating/decelerating flight, encountered by a space vehicle during launch and/or entry. The vehicle can also encounter numerous configuration changes, some of which can be quite drastic, during the course of its flight, so robust, automated geometry modeling, grid generation, and adaptation will play a huge role in reaching this goal. Multiply this by 1000’s of trajectory simulations occurring simultaneously in a given Monte Carlo analysis, and the problem readily scales to absorb virtually any size of supercomputer envisioned today. The concept of CFD-in-the-loop Monte Carlo trajectory simulation poses a formidable challenge for emerging and future computing systems, and it has the potential to shave years off the development cycle for aerodynamic and aerothermodynamic performance predictions as compared to today’s space vehicle design approach.

CFD 2030↗

Uncovering the Geometry of the Hot X-Ray Corona in the Seyfert Galaxy NGC 4151 with IXPE

We present an X-ray spectropolarimetric analysis of the bright Seyfert galaxy NGC 4151. The source has been observed with the Imaging X-ray Polarimetry Explorer (IXPE) for 700 ks, complemented with simultaneous XMM–Newton (50 ks) and NuSTAR (100 ks) pointings. A polarization degree Π = 4.9 ± 1.1 per cent and angle Ψ = 86° ± 7° east of north (68 per cent confidence level) are measured in the 2–8 keV energy range. The spectropolarimetric analysis shows that the polarization could be entirely due to reflection. Given the low reflection flux in the IXPE band, this requires, however, a reflection with a very large (>38 per cent) polarization degree. Assuming more reasonable values, a polarization degree of the hot corona ranging from ∼4 to ∼8 per cent is found. The observed polarization degree excludes a ‘spherical’ lamppost geometry for the corona, suggesting instead a slab-like geometry, possibly a wedge, as determined via Monte Carlo simulations. This is further confirmed by the X-ray polarization angle, which coincides with the direction of the extended radio emission in this source, supposed to match the disc axis. NGC 4151 is the first active galactic nucleus with an X-ray polarization measure for the corona, illustrating the capabilities of X-ray polarimetry and IXPE in unveiling its geometry.

NGC 4151↗

Surface entropy of liquids via a direct Monte Carlo approach - Application to liquid Si

Two methods are presented for a direct Monte Carlo evaluation of the surface entropy S(s) of a liquid interacting by specified, volume-independent potentials. The first method is based on an application of the approach of Ferrenberg and Swendsen (1988, 1989) to Monte Carlo simulations at two different temperatures; it gives much more reliable results for S(s) in liquid Si than previous calculations based on numerical differentiation. The second method expresses the surface entropy directly as a canonical average at fixed temperature.

Wang, Z. Q.↗

Cloud Thickness from Diffusion of Lidar Pulses in Clouds

Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Cahalan, Robert F.↗

Modelling of SAR polarisation phase difference from trees

The data for polarization phase difference Delta Phi between the HH- and VV-polarized backscattered waves from tree-covered fields were obtained with an airborne synthetic aperture radar at 1.225 GHz. The mean values over tree-covered fields were derived from the images of the phase difference and were examined as a function of incident beam angle from 15 to 55 deg. A theoretical model for simulating these data, based on the electromagnetic wave scatterings from the tree trunk and its branches, both of which are assumed as very long dielectric cylinders was developed. The radius and direction of a tree branch are taken as random variables and are chosen by a Monte Carlo method to encounter the incident waves in producing the scattering events. The Monte Carlo simulated results are in good agreement with the observations within experimental uncertainty.

Mo, Tsan↗

CFD2030 Grand Challenge: CFD-in-the-Loop Monte Carlo Simulation for Space Vehicle Design

Space vehicle design and certification differs widely from aircraft design relying more on probabilistic approaches than deterministic. Monte Carlo simulation plays an important role in the probabilistic design of space vehicles to ensure robust and reliable operation. Today, Monte Carlo flight simulation requires 1000’s of trajectory simulations that use databases to provide aerosciences models. These databases can be extremely expensive and time consuming to develop. Replacing these databases with unsteady computational fluid dynamics directly in the simulation loop has potential to significantly reduce the time required to analyze space vehicle concepts, improve simulation accuracy, and reduce the cost of space vehicle development. The CFD Vision 2030 Study outlined gaps and roadblocks to meeting the vision described in the study. The geometric, physical, and computational challenges associated with CFD-in-the-loop Monte Carlo simulation for space vehicle design are substantial and serve as an excellent grand challenge to advance the CFD 2030 vision.

CFD2030↗

Optimizing Muscle Parameters in Musculoskeletal Modeling Using Monte Carlo Simulations

Astronauts assigned to long-duration missions experience bone and muscle atrophy in the lower limbs. The use of musculoskeletal simulation software has become a useful tool for modeling joint and muscle forces during human activity in reduced gravity as access to direct experimentation is limited. Knowledge of muscle and joint loads can better inform the design of exercise protocols and exercise countermeasure equipment. In this study, the LifeModeler(TM) (San Clemente, CA) biomechanics simulation software was used to model a squat exercise. The initial model using default parameters yielded physiologically reasonable hip-joint forces. However, no activation was predicted in some large muscles such as rectus femoris, which have been shown to be active in 1-g performance of the activity. Parametric testing was conducted using Monte Carlo methods and combinatorial reduction to find a muscle parameter set that more closely matched physiologically observed activation patterns during the squat exercise. Peak hip joint force using the default parameters was 2.96 times body weight (BW) and increased to 3.21 BW in an optimized, feature-selected test case. The rectus femoris was predicted to peak at 60.1% activation following muscle recruitment optimization, compared to 19.2% activation with default parameters. These results indicate the critical role that muscle parameters play in joint force estimation and the need for exploration of the solution space to achieve physiologically realistic muscle activation.

Hanson, Andrea↗

Multi-Scale Modeling of Liquid Phase Sintering Affected by Gravity: Preliminary Analysis

A multi-scale simulation concept taking into account impact of gravity on liquid phase sintering is described. The gravity influence can be included at both the micro- and macro-scales. At the micro-scale, the diffusion mass-transport is directionally modified in the framework of kinetic Monte-Carlo simulations to include the impact of gravity. The micro-scale simulations can provide the values of the constitutive parameters for macroscopic sintering simulations. At the macro-scale, we are attempting to embed a continuum model of sintering into a finite-element framework that includes the gravity forces and substrate friction. If successful, the finite elements analysis will enable predictions relevant to space-based processing, including size and shape and property predictions. Model experiments are underway to support the models via extraction of viscosity moduli versus composition, particle size, heating rate, temperature and time.

Olevsky, Eugene↗

A theoretical and numerical study of polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum

In this paper, theoretical and numerical results of the polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum are presented for the remote sensing of ocean and soil surfaces. The polarimetric scattered field for rough dielectric surfaces is derived to the second order by the small perturbation method (SPM). It is found that the second-order scattered field is coherent in nature, and its coefficients for different polarizations present the lowest-order corrections to the Fresnel reflection coefficients of the surfaces. In addition, the cross-polarized (HV and VH) components of the coherent fields are reciprocal and not zero for surfaces with anisotropic directional spectrum when the azimuth angle of the incident direction is not aligned with the symmetry directions of surfaces. In order to verify the energy conservation condition of the theoretical results, which is important if the theory is to be applied to the passive polarimetry of rough surfaces, a Monte Carlo simulation is performed to numerically calculate the polarimetric reflectivities of one-dimensional random rough surfaces which are generated with a prescribed power-law spectrum in the spectral domain and transformed to the spatial domain by the FFT. The surfaces simulated by this approach are periodic with the period corresponding to the low-wavenumber cutoff. To calculate the scattering from periodic dielectric surfaces, the authors present a new numerical technique which applies the Floquet theorem to reduce the problem to one period and does not require the evaluation of one-dimensional periodic Green's function used in the conventional method of moment formulation. Once the scattering coefficients are obtained, the polarimetric Stokes vectors for the emission from the random surfaces are then calculated according to the Kirchhoff's law and are illustrated as functions of relative azimuth observation and row directions. The second-order SPM is also reduced to this one-dimensional case, and theoretical results for all four Stokes parameters are in excellent agreement with the numerical results obtained from the Monte Carlo simulation. Moreover, the second order coherent fields are found to be necessary for the theoretical evaluation of the third and fourth Stokes parameters. Without that, the reflectivities of random rough surfaces are significantly over-estimated, and the signs of the third and the fourth Stokes parameters are found to be incorrect, if calculated by using the first-order incoherent fields alone.

Yueh, S. H.↗