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

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↗

The fragmentation of 510 MeV/nucleon iron-56 in polyethylene. II. Comparisons between data and a model

The results of a Monte Carlo model for calculating fragment fluences and LET spectra are compared to data taken with 600 MeV/nucleon iron ions incident on an accelerator beamline configured for irradiation of biological samples, with no target and with 2, 5 and 8 cm of polyethylene. The model uses a multi-generation nuclear fragmentation code, coupled with a formulation of ionization energy loss based on the Bethe-Bloch equation. In the region where the data are reliable and the experimental acceptance is well understood, many of the features of the experimental spectra are well replicated by the model. To obtain good agreement with the experimental data, the model must allow for at least two generations of fragment production in the target.

NASA Discipline Number 99-99↗

Study of the Radiative Properties of Inhomogeneous Stratocumulus Clouds

Clouds play an important role in the radiation budget of the atmosphere. A good understanding of how clouds interact with solar radiation is necessary when considering their effects in both general circulation models and climate models. This study examined the radiative properties of clouds in both an inhomogeneous cloud system, and a simplified cloud system through the use of a Monte Carlo model. The purpose was to become more familiar with the radiative properties of clouds, especially absorption, and to investigate the excess absorption of solar radiation from observations over that calculated from theory. The first cloud system indicated that the absorptance actually decreased as the cloud's inhomogeneity increased, and that cloud forcing does not indicate any changes. The simplified cloud system looked at two different cases of absorption of solar radiation in the cloud. The absorptances calculated from the Monte Carlo is compared to a correction method for calculating absorptances and found that the method can over or underestimate absorptances at cloud edges. Also the cloud edge effects due to solar radiation points to a possibility of overestimating the retrieved optical depth at the edge, and indicates a possible way to correct for it. The effective cloud fraction (Ne) for a long time has been calculated from a cloud's reflectance. From the reflectance it has been observed that the N, for most cloud geometries is greater than the actual cloud fraction (Nc) making a cloud appear wider than it is optically. Recent studies we have performed used a Monte Carlo model to calculate the N, of a cloud using not only the reflectance but also the absorptance. The derived Ne's from the absorptance in some of the Monte Carlo runs did not give the same results as derived from the reflectance. This study also examined the inhomogeneity of clouds to find a relationship between larger and smaller scales, or wavelengths, of the cloud. Both Fourier transforms and wavelet transforms were used to analyze the liquid water content of marine stratocumulus clouds taken during the ASTEX project. From the analysis it was found that the energy in the cloud is not uniformly distributed but is greater at the larger scales than at the smaller scales. This was determined by examining the slope of the power spectrum, and by comparing the variability at two scales from a wavelet analysis.

Batey, Michael↗

Monte Carlo transport model comparison with 1A GeV accelerated iron experiment: heavy-ion shielding evaluation of NASA space flight-crew foodstuff

Deep-space manned flight as a reality depends on a viable solution to the radiation problem. Both acute and chronic radiation health threats are known to exist, with solar particle events as an example of the former and galactic cosmic rays (GCR) of the latter. In this experiment Iron ions of 1A GeV are used to simulate GCR and to determine the secondary radiation field created as the GCR-like particles interact with a thick target. A NASA prepared food pantry locker was subjected to the iron beam and the secondary fluence recorded. A modified version of the Monte Carlo heavy ion transport code developed by Zeitlin at LBNL is compared with experimental fluence. The foodstuff is modeled as mixed nuts as defined by the 71st edition of the Chemical Rubber Company (CRC) Handbook of Physics and Chemistry. The results indicate a good agreement between the experimental data and the model. The agreement between model and experiment is determined using a linear fit to ordered pairs of data. The intercept is forced to zero. The slope fit is 0.825 and the R2 value is 0.429 over the resolved fluence region. The removal of an outlier, Z=14, gives values of 0.888 and 0.705 for slope and R2 respectively. c2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

NASA Discipline Radiation Health↗

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3, documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation is to describe the application of the hydraulic property fields and recharge parameters documented in ECF-HANFORD-20-0027, Null Space Monte Carlo Evaluation of the Plateau to River Model, to the CA flow and fate and transport simulation results to quantify the uncertainty in the simulated results due to input parameter selection. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, models lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. The result of NSCM analysis is a set of F&T simulations that provide an estimate of the range of possible outcomes that are used to quantify the uncertainty in simulated concentrations produced using the base case simulations. The simulated concentrations from all simulations will support calculation of the uncertainty of the total dose calculated in a separate calculation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Charge order in the kagome lattice Holstein model: a hybrid Monte Carlo study

The Holstein model is a paradigmatic description of the electron-phonon interaction, in which electrons couple to local dispersionless phonon modes, independent of momentum. The model has been shown to host a variety of ordered ground states such as charge density wave (CDW) order and superconductivity on several geometries, including the square, honeycomb, and Lieb lattices. In this work, we study CDW formation in the Holstein model on the kagome lattice, using a recently developed hybrid Monte Carlo simulation method. We present evidence for √3 × √3 CDW order at an average electron filling of < n > = 2/3 per site, with an ordering wavevector at the K-points of the Brillouin zone. We estimate a phase transition occurring at T c ≈ t/18, where t is the nearest-neighbor hopping parameter. Our simulations find no signature of CDW order at other electron fillings or ordering momenta for temperatures T ≥ t/20.

36 MATERIALS SCIENCE↗

Reaction Rate Distribution Measurements in ZPPR-15

Small metal foils were irradiated in ZPPR-15 to measure distributions for 239 Pu(n,f), 235 U(n,f), 238 U(n,f), 238 U(n,γ), 115 In(n,n’) and 115 In(n,γ) reaction rates in ZPPR-15. The measured data were reported as basic reaction rate distributions for individual foils, cell-averaged reaction rate distributions and spectral indices or reaction rate ratios. All of the measured data were analyzed and tabulated in such a manner that the data can be readily compared with computed values. The published uncertainties for ZPPR-15 the reaction rate measurements include the statistical uncertainties in the measurements, reproducibility uncertainty and small uncertainties related to corrections for items such as foil compositions and post-irradiation decay between reactor shutdown and foil measurement. There are additional uncertainty components related to detector calibration that are not included in the published uncertainties. A full uncertainty analysis was performed for each foil measurement and for quantities derived from the basic foil measurements. The significant uncertainties were quantified, and a total uncertainty was derived for each measured value. The ZPPR-15 experimental records were used to create detailed as-built Monte Carlo models for the foil reaction rate measurements in ZPPR-15 loadings 123, 134 and 203 and for the fission chamber measurements of 237 Np and 242 Pu fission rates in ZPPR-15 loading 141. Because there were many foils distributed in the ZPPR-15 assembly and because an individual foil occupies a negligible volume in the ZPPR-15 assembly, direct Monte Carlo calculations of the measured foil reaction rates would require an impractical computational effort. The as-built Monte Carlo models and related files were prepared for conversion to deterministic models that can be calculated with realistic computational effort.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SiO(X) coatings for atomic oxygen protection of polyimide Kapton in low earth orbit

Sputter-deposited SiO(X) (where X between 1.9 and 2.0) thin film coatings have been found to be durable to atomic oxygen. Such coatings will be used to protect polyimide Kapton photovoltaic array blankets from atomic oxygen attack in low earth orbit (LEO) on the Space Station Freedom (SSF). Monte Carlo modeling of atomic oxygen attack at defect sites in protected Kapton exposed in laboratory RF plasma ashers and on solar tracking photovoltaic arrays in space has been conducted to enable understanding of degradation processes relevant to the durability of the SSF solar array blanket. SiO(X) protective coating performance data from RF plasma asher tests will be presented, along with Monte Carlo modeling considerations, to enable the projection of in-space durability of the SSF solar array blankets based on ground laboratory test results.

Banks, Bruce A.↗

An Efficient Monte Carlo Method for Modeling Radiative Transfer in Protoplanetary Disks

Monte Carlo methods have been shown to be effective and versatile in modeling radiative transfer processes to calculate model temperature profiles for protoplanetary disks. Temperatures profiles are important for connecting physical structure to observation and for understanding the conditions for planet formation and migration. However, certain areas of the disk such as the optically thick disk interior are under-sampled, or are of particular interest such as the snow line (where water vapor condenses into ice) and the area surrounding a protoplanet. To improve the sampling, photon packets can be preferentially scattered and reemitted toward the preferred locations at the cost of weighting packet energies to conserve the average energy flux. Here I report on the weighting schemes developed, how they can be applied to various models, and how they affect simulation mechanics and results. We find that improvements in sampling do not always imply similar improvements in temperature accuracies and calculation speeds.

radiative transfer↗

Atomic oxygen interaction with solar array blankets at protective coating defect sites

Atomic oxygen in the low-Earth-orbital environment oxidizes SiOx protected polyimide Kapton solar array blankets at sites which are not protected such as pin windows or scratches in the protective coatings. The magnitude and shape of the atomic oxygen undercutting which occurs at these sites is dependent upon the exposure environment details such as arrival direction and reaction probability. The geometry of atomic oxygen undercutting at defect sites exposed to atomic oxygen in plasma asher was used to develop a Monte Carlo model to simulate atomic oxygen erosion processes at defect sites in protected Kapton. Comparisons of Monte Carlo predictions and experimental results are presented for plasma asher atomic oxygen exposures for large and small defects as well as for protective coatings on one or both sides of Kapton. The model is used to predict in-space exposure results at defect sites for both directed and sweeping atomic oxygen exposure. A comparison of surface textures predicted by the Monte Carlo model and those experimentally observed from both directed space ram and laboratory plasma asher atomic oxygen exposure indicate substantial agreement.

Banks, Bruce A.↗

Imaging and spectroscopy of Comet P/Halley

The goals of this investigation are the analysis of a large set of high-resolution echelle/reticon spectra, and the reduction and analysis of a set of narrow-band-filtered charge coupled device (CCD) images of Comet Halley taken during the preperihelion period at Oak Ridge Observatory by Dr. R. E. McCrosky. The scientific objectives associated with these goals are the determination of the spatial distributions of several important radicals, atoms and ions in the coma. These include C2, CN, C3, H2O(+) and CO(+) from the image data. The analysis of the neutral species distribution with Monte Carlo models will aid in the understanding of their production and decay mechanisms as well as serve as an important indicator of the physical conditions in the inner coma. The spatial distributions of the ions will serve as a guide to constrain the complex model necessary for understanding the interaction of the solar wind and the cometary ions. Work during this past year has been devoted largely to the reduction of the standard star photometry for the CCD image data set, as well as the re-flat-fielding of a number of the comet images. We are pleased to report that despite a number of setbacks and the small effort devoted to this work (2 1/2 months for the PI and a generous share of completely unsupported time by Dr. McCrosky) that this portion of the work has been successfully completed. The goals for the upcoming final year of this project (under a new project number) are to complete the calibration of the CCD image data for inclusion in the IHW archive, to analyze a select portion of the neutral radical images with our Monte Carlo models, and to present the results of the 6300/region spectra as a guide to low-resolution spectral observers in order to yield the unambiguous separation of the contributions of cometary O(1D), airglow O(1D), and the numerous NH2 lines in that region of the spectrum.

Combi, Michael R.↗

Differential Measurements of Jet Substructure and Partonic Energy Loss in Au+Au Collisions at √ sNN = 200 GeV

The STAR collaboration presents jet substructure measurements related to both the momentum fraction and the opening angle within jets in p+p and Au+Au collisions at √sNN =200GeV. The substructure observables include SoftDrop groomed momentum fraction (zg), groomed jet radius (Rg), and subjet momentum fraction (zSJ) and opening angle (θSJ). The latter observable is introduced for the first time. Fully corrected subjet measurements are presented for p+p collisions and are compared to leading-order Monte Carlo models. The subjet θSJ distributions reflect the jets leading opening angle and are utilized as a proxy for the resolution scale of the medium in Au+Au collisions. We compare data from Au+Au collisions to those from p+p which are embedded in minimum-bias Au+Au events in order to include the effects of detector smearing and the heavy-ion collision underlying event. The subjet observables are shown to be more robust to the background than zg and Rg. We observe no significant modifications of the subjet observables within the two highest-energy, back-to-back jets, resulting in a distribution of opening angles and the splittings that are vacuumlike. We also report measurements of the differential dijet momentum imbalance (AJ) for jets of varying θSJ. We find no qualitative differences in energy loss signatures for varying angular scales in the range 0.1< θSJ<0.3, leading to the possible interpretation that energy loss in this population of high-momentum dijet pairs, is due to soft medium-induced gluon radiation from a single color charge as it traverses the medium.

Abdallah, M. S.↗

Low-energy enhancement in the magnetic dipole γ-ray strength functions of heavy nuclei

A low-energy enhancement (LEE), which was observed experimentally in the gamma-ray strength function (γSF) describing the decay of compound nuclei, would have profound effects on r-process nucleosynthesis if it persists in heavy neutron-rich nuclei. The LEE was shown to be a feature of the magnetic dipole (M1) strength function in configuration-interaction shell-model calculations in medium-mass nuclei. However, its existence in heavy open-shell nuclei remains an open question. Here, using a combination of many-body methods, we identify a LEE in the M1 γSFs of heavy samarium nuclei. In particular, we use the static-path plus random-phase approximation (SPA+RPA), which includes static and small-amplitude quantal fluctuations beyond the mean field. Using the SPA+RPA strength as a prior, we apply the maximum-entropy method (MEM) to obtain finite-temperature M1 γSFs from exact imaginary-time response functions calculated with the shell model Monte Carlo (SMMC) method. We find that the slope of the LEE in samarium isotopes is roughly independent of the average initial energy over a wide range below the neutron separation energy. As the neutron number increases, strength transfers to a low-energy excitation, which we interpret as the scissors mode built on top of excited states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Low-energy enhancement of the magnetic dipole radiation in odd-mass lanthanides

We compute the magnetic dipole (M1) $\gamma$-ray strength functions ($\gamma$SF) for the odd-mass lanthanides $^{\textrm{143-151}}$Nd and $^{\textrm{147-153}}$Sm using the shell-model Monte Carlo method in combination with the static-path approximation and the maximum-entropy method. In particular, we quantify the statistical uncertainties in the calculated M1 $\gamma$SFs and show that they are under control for the excitation energies relevant to the experiments despite a Monte Carlo sign problem that originates in the projection onto an odd number of neutrons. We identify a low-energy enhancement (LEE) in the M1 $\gamma$SFs of these odd-mass lanthanides, which was recently observed experimentally in some of them. We also find a scissors mode resonance (SR) in the strongly deformed isotopes. We observe that the decrease in the LEE strength with neutron number along an isotopic chain is compensated for by an increase in the SR strength in the deformed nuclei. Furthermore, we compare our results with recent experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau-to-River (P2R) Model is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by Central Plateau Cleanup Company (CPCCo) at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the composite analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Monte Carlo lattice models for adsorbed polymer conformation

The adhesion between a polymer film and a metal surface is of great technological interest. However, the prediction of adhesion and wear properties of polymer coated metals is quite difficult because a fundamental understanding of the polymer surface interaction does not yet exist. A computer model for the conformation of a polymer molecule adsorbed on a surface is discussed. The chain conformation is assumed to be described by a partially directed random walk on a three dimensional simple cubic lattice. An attractive surface potential is incorporated into the model through the use of a random walk step probability distribution that is anisotropic in the direction normal to the attractive surface. The effects of variations in potential characteristics are qualitatively included by varying both the degree of anisotropy of the step distribution and the range of the anisotropy. Polymer conformation is characterized by the average end to end distance, average radius of gyration, and average number of chain segments adsorbed on the surface.

Good, B. S.↗