Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo modeling”

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 703 records · Page 39

STAT7 v2.0 User Guide

The STAT7 software was developed to perform steady-state, single-phase thermal hydraulics analysis of plate-fueled reactors based on statistical propagation of uncertainties. Application of the software includes non-power research and test reactor analysis, and it has been used for the conversion to low- enriched uranium fuel of U.S. High Performance Research Reactors such as the Massachusetts Institute of Technology Research Reactor. Since it can be necessary to repeat reactor safety analysis, such as during fuel reloading, STAT7 accommodates flexibility in analyzing many practical aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model uncertainty in common fuel fabrication parameters and other key reactor operating parameters required for reactor thermal hydraulics analysis. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux, or onset of flow instability. STAT7 supports water properties based on the IAPWS-IF97 functions (The International Association for the Properties of Water and Steam Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam) in addition to fitted functions. STAT7 predicts axial profiles of fuel, cladding, and coolant temperature along a lateral stripe that runs the full length of the fuel plate from the bottom to the top. STAT7 can simultaneously analyze every axial node in each lateral stripe of all fuel plates and coolant channels in every fuel element of an entire reactor core. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. In a single execution, STAT7 can be used to perform full core analysis by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each fuel element in the core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulation Study of High-Precision Characterization of MeV Electron Interactions for Advanced Nano-Imaging of Thick Biological Samples and Microchips

The resolution of a mega-electron-volt scanning transmission electron microscope (MeV-STEM) is primarily governed by the properties of the incident electron beam and angular broadening effects that occur within thick biological samples and microchips. A precise understanding and mitigation of these constraints require detailed knowledge of beam emittance, aberrations in the STEM column optics, and energy-dependent elastic and inelastic critical angles of the materials being examined. This simulation study proposes a standardized experimental framework for comprehensively assessing beam intensity, divergence, and size at the sample exit. This framework aims to characterize electron-sample interactions, reconcile discrepancies among analytical models, and validate Monte Carlo (MC) simulations for enhanced predictive accuracy. Our numerical findings demonstrate that precise measurements of these parameters, especially angular broadening, are not only feasible but also essential for optimizing imaging resolution in thick biological samples and microchips. By utilizing an electron source with minimal emittance and tailored beam characteristics, along with amorphous ice and silicon samples as biological proxies and microchip materials, this research seeks to optimize electron beam energy by focusing on parameters to improve the resolution in MeV-STEM/TEM. This optimization is particularly crucial for in situ imaging of thick biological samples and for examining microchip defects with nanometer resolutions. Our ultimate goal is to develop a comprehensive mapping of the minimum electron energy required to achieve a nanoscale resolution, taking into account variations in sample thickness, composition, and imaging mode.

36 MATERIALS SCIENCE↗

Martian doublet craters

Monte Carlo cratering simulation model to show nonrandomness of formation of Mars tangential meteorite craters

Aoyagi, M.↗

Digital simulation of hybrid loop operation in RFI backgrounds.

A digital computer model for Monte-Carlo simulation of an imperfect second-order hybrid phase-locked loop (PLL) operating in radio-frequency interference (RFI) and Gaussian noise backgrounds has been developed. Characterization of hybrid loop performance in terms of cycle slipping statistics and phase error variance, through computer simulation, indicates that the hybrid loop has desirable performance characteristics in RFI backgrounds over the conventional PLL or the costas loop.

Ziemer, R. E.↗

Martian doublet craters.

A large number of Mars craters are nearly tangential to other craters. They occur in clusters or as isolated crater doublets. Results of probability calculations and a Monte Carlo cratering simulation model show conclusively that many of the Mars craters could not have resulted from random single-body impact. The possibility that these craters are calderas is considered possible only if calderas on Mars form by mechanisms different from those on earth. However, clusters and doublets could be caused by meteoroid breakup resulting from stresses induced in the meteoroid by the gravitational field of Mars. It is concluded that, under certain conditions, doublets should be produced on Mars as a direct result of breakup of an impacting meteoroid. The impact process can yield nonrandom crater distributions that should be observed in different degrees of development on different planetary surfaces.

Oberbeck, V. R.↗

Computed laser backscattering from turbid liquids - Comparison with laboratory results

A recently developed semianalytic Monte Carlo radiative transfer model applicable to oceanographic lidar systems (SALMON) has been used to simulate a series of laboratory experiments studying the backscatter of laser light from dispersions of Teflon spheres. Results obtained with SALMON, using best estimates of pertinent optical parameters, compared quite well with experimental data in both a qualitative and quantitative sense, with the largest relative discrepancies being approximately 30%. The results firmly establish the validity of SALMON in studying the performance of real oceanographic lidar systems.

Poole, L. R.↗

Segmentation and pulse shape discrimination techniques for rejecting background in germanium detectors

The possibility of rejecting the internal beta-decay background in coaxial germanium detectors by distinguishing between the multi-site energy losses characteristic of photons and the single-site energy losses of electrons in the range 0.2 - 2 MeV is examined. The photon transport was modeled with a Monte Carlo routine. Background rejection by both multiple segmentation and pulse shape discrimination techniques is investigated. The efficiency of a six 1 cm-thick segment coaxial detector operating in coincidence mode alone is compared to that of a two-segment (1 cm and 5 cm) detector employing both front-rear coincidence and PSD in the rear segment to isolate photon events. Both techniques can provide at least 95 percent rejection of single-site events while accepting at least 80 percent of the multi-site events above 500 keV.

Roth, J.↗

Vector computer memory bank contention

A number of vector supercomputers feature very large memories. Unfortunately the large capacity memory chips that are used in these computers are much slower than the fast central processing unit (CPU) circuitry. As a result, memory bank reservation times (in CPU ticks) are much longer than on previous generations of computers. A consequence of these long reservation times is that memory bank contention is sharply increased, resulting in significantly lowered performance rates. The phenomenon of memory bank contention in vector computers is analyzed using both a Markov chain model and a Monte Carlo simulation program. The results of this analysis indicate that future generations of supercomputers must either employ much faster memory chips or else feature very large numbers of independent memory banks.

Bailey, D. H.↗

The gamma-ray telescope Gamma-1

French and Soviet specialists have designed and built the gamma-ray telescope GAMMA-1 to detect cosmic gamma rays above 50 MeV. The sensitive area of the detector is 1400 sq cm, energy resolution is 30% at 300 MeV, and angular resolution 1.2 deg at 300 MeV (and less than 20' arc when a coded aperture mask is used). Results on calibration of the qualification model and Monte-Carlo calculations are presented.

Akimov, V. V.↗

Cloud shadow effects on remote sensing

On relatively clear days, the shadow cast by an isolated, optically dense cloud produces an effect which is similar in extent and magnitude to the adjacency effects associated with high-contrast ground albedo boundaries. A Monte Carlo radiative transfer model has been used to examine the intensity and its components, comparing the cloud shadow case with that for a Lambertian albedo boundary. Although the behavior of the total intensity is nearly the same for both cases, the behaviors of the intensity components differ significantly.

Pearce, W. A.↗

A laboratory investigation of the variability of cloud reflected radiance fields

A method to determine the radiative properties of complex cloud fields was developed. A Cloud field optical simulator (CFOS) was constructed to simulate the interaction of cloud fields with visible radiation. The CFOS was verified by comparing experimental results from it with calculations performed with a Monte Carlo radiative transfer model. A software library was developed to process, reduce, and display CFOS data. The CFSOS was utilized to study the reflected radiane patterns from simulated cloud fields.

Mckee, T. B.↗

Effects of multiple scattering on time- and depth-resolved signals in airborne lidar systems

A semianalytic Monte Carlo radiative transfer model (SALMON) is employed to probe the effects of multiple-scattering events on the time- and depth-resolved lidar signals from homogeneous aqueous media. The effective total attenuation coefficients in the single-scattering approximation are determined as functions of dimensionless parameters characterizing the lidar system and the medium. Results show that single-scattering events dominate when these parameters are close to their lower bounds and that when their values exceed unity multiple-scattering events dominate.

Punjabi, A.↗

Vector computer memory bank contention

A number of vector supercomputers feature very large memories. Unfortunately the large capacity memory chips that are used in these computers are much slower than the fast central processing unit (CPU) circuitry. As a result, memory bank reservation times (in CPU ticks) are much longer than on previous generations of computers. A consequence of these long reservation times is that memory bank contention is sharply increased, resulting in significantly lowered performance rates. The phenomenon of memory bank contention in vector computers is analyzed using both a Markov chain model and a Monte Carlo simulation program. The results of this analysis indicate that future generations of supercomputers must either employ much faster memory chips or else feature very large numbers of independent memory banks.

Bailey, David H.↗

The role of hydrogen thermalization in shaping the Lyman-alpha coma

Two interrelated aspects of the thermalization of cometary hydrogen were addressed using a three-dimensional Monte Carlo particle trajectory model: photochemical heating of the coma and the velocity distribution of cometary hydrogen leaving the inner coma that produces the extended Lyman-alpha coma. Observed cometary conditions were examined, ranging from productive comets like Kohoutek at very small heliocentric distances to intrinsically small comets like Giacobini-Zinner at 1 AU. The collisional decoupling of hydrogen yields bulk flow speeds and gas temperatures which are smaller than generally predicted in hydrodynamic models. By properly including these conditions the hydrogen velocity distribution required to produce the observed Lyman-alpha coma can be derived naturally. The two-dimensional morphology of the H coma is diagnostic of the physical conditions in the inner coma.

Combi, Michael R.↗

A satellite-based radar wind sensor

The objective is to investigate the application of Doppler radar systems for global wind measurement. A model of the satellite-based radar wind sounder (RAWS) is discussed, and many critical problems in the designing process, such as the antenna scan pattern, tracking the Doppler shift caused by satellite motion, and backscattering of radar signals from different types of clouds, are discussed along with their computer simulations. In addition, algorithms for measuring mean frequency of radar echoes, such as the Fast Fourier Transform (FFT) estimator, the covariance estimator, and the estimators based on autoregressive models, are discussed. Monte Carlo computer simulations were used to compare the performance of these algorithms. Anti-alias methods are discussed for the FFT and the autoregressive methods. Several algorithms for reducing radar ambiguity were studied, such as random phase coding methods and staggered pulse repitition frequncy (PRF) methods. Computer simulations showed that these methods are not applicable to the RAWS because of the broad spectral widths of the radar echoes from clouds. A waveform modulation method using the concept of spread spectrum and correlation detection was developed to solve the radar ambiguity. Radar ambiguity functions were used to analyze the effective signal-to-noise ratios for the waveform modulation method. The results showed that, with suitable bandwidth product and modulation of the waveform, this method can achieve the desired maximum range and maximum frequency of the radar system.

Xin, Weizhuang↗

Bi-directional reflectance of finite and infinite clouds

The transfer of solar irradiance in plane parallel and broken cloud fields is modeled 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: a local maximum of reflected intensity in the forward direction for all cases; a limb darkening for the plane parallel case; and a limb brightening and a local maximum of reflected intensity in the backward direction for broken clouds. The method is used to study the cloud reflectance sensitivity to various parameters. A more precise description of cloud field internal and external structure is needed in order to obtain accurate bi directional reflectance diagrams.

Breon, Francois-Marie↗