Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Transport Simulations”

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 289 records · Page 16

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Molecular Simulations Probing the Adsorption and Diffusion of Ammonia, Nitrogen, Hydrogen, and Their Mixtures in Bulk MFI Zeolite and MFI Nanosheets at High Temperature and Pressure

Recent advances in the synthesis of MFI zeolite nanosheets have led to highly selective membranes that are promising candidates for small-scale ammonia separation from ammonia/nitrogen/hydrogen mixtures in distributed green ammonia production plants. Using force-field-based molecular simulations, we evaluate the performance of bulk all-silica MFI zeolite and a 3 nm thick all-silica MFI nanosheet for ammonia/nitrogen/hydrogen separations at moderate temperature and pressure conditions (T = 373 K, p = 5 bar) as well as conditions relevant for a membrane-based reactor–separator process (T = 523 or 623 K, p = 80 bar). Isobaric–isothermal Gibbs ensemble Monte Carlo (GEMC) simulations were carried out to understand the selective adsorption behavior for these mixtures. Molecular dynamics simulations in the canonical ensemble were performed to examine the transport behavior of ammonia in a mixture with nitrogen or hydrogen at loadings obtained from the GEMC simulations. Our results show that both bulk MFI and the nanosheet with explicit surface silanols are highly selective toward ammonia adsorption, but the adsorption selectivity decreases by factors of 4 and 10 for ammonia/nitrogen and ammonia/hydrogen mixtures as the temperature is increased from 373 to 623 K. Conversely, the diffusion selectivities toward ammonia are more favorable at process-relevant temperatures and for the MFI nanosheet. Here, at T = 523 K and p = 80 bar, our simulations predict overall separation factors of 3.8 ± 0.2 and 4.5 ± 0.3 for ammonia/nitrogen and ammonia/hydrogen mixtures, respectively, in the bulk MFI zeolite, and separation factors of 8.6 ± 0.4 and 12.5 ± 0.8, respectively, for the MFI nanosheet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Grain Boundary-Limited Thermal Transport in Suspended Thin Graphite across an Unexplored Thickness Regime

Here, we present systematic thermal conductivity (κ) measurements of suspended thin graphite ribbons, 234–527 nm thick, using a four-probe 3ω method. Unlike recent reports of phonon hydrodynamics and exceptionally high κ in micrometer-thick graphite ( Science, 2020), we observe significantly lower κ and no signatures of collective phonon flow in this intermediate thickness regime. Instead, our measured κ lies between few-layer graphene and bulk graphite. These results agree with a first-principles-informed Peierls–Boltzmann transport model with spatially resolved Monte Carlo sampling. Additionally, the temperature for the peak κ shifts lower with increasing thickness, due to the interplay of phonon-boundary and phonon-isotope scattering. Incorporating grain boundary scattering into simulations is necessary to replicate the experimental trends. These findings delineate the boundary between ballistic, hydrodynamic, and diffusive transport regimes in graphite and underscore the dominant role of disorder and geometry in phonon transport in quasi-two-dimensional materials, offering insights for nanoscale thermal management.

Boltzmann transport model↗

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

22 Ne Phase Separation as a Solution to the Ultramassive White Dwarf Cooling Anomaly

The precise astrometric measurements of the Gaia Data Release 2 have opened the door to detailed tests of the predictions of white dwarf cooling models. Significant discrepancies between theory and observations have been identified, the most striking affecting ultramassive white dwarfs. Cheng et al. found that a small fraction of white dwarfs on the so-called Q branch must experience an extra cooling delay of ~8 Gyr not predicted by current models. 22 Ne phase separation in a crystallizing C/O white dwarf can lead to a distillation process that efficiently transports 22 Ne toward its center, thereby releasing a considerable amount of gravitational energy. Here, using state-of-the-art Monte Carlo simulations, we show that this mechanism can largely resolve the ultramassive cooling anomaly if the delayed population consists of white dwarfs with moderately above-average 22 Ne abundances. We also argue that 22 Ne phase separation can account for the smaller cooling delay currently missing for models of white dwarfs with more standard compositions.

79 ASTRONOMY AND ASTROPHYSICS↗

Orbit IMU alignment: Error analysis

A comprehensive accuracy analysis of orbit inertial measurement unit (IMU) alignments using the shuttle star trackers was completed and the results are presented. Monte Carlo techniques were used in a computer simulation of the IMU alignment hardware and software systems to: (1) determine the expected Space Transportation System 1 Flight (STS-1) manual mode IMU alignment accuracy; (2) investigate the accuracy of alignments in later shuttle flights when the automatic mode of star acquisition may be used; and (3) verify that an analytical model previously used for estimating the alignment error is a valid model. The analysis results do not differ significantly from expectations. The standard deviation in the IMU alignment error for STS-1 alignments was determined to the 68 arc seconds per axis. This corresponds to a 99.7% probability that the magnitude of the total alignment error is less than 258 arc seconds.

Corson, R. W.↗

On the transport of ions released in the magnetotail by the AMPTE-IRM satellite

The Ba and Li ions releasd into the magnetotail in spring 1985 by the AMPTE-IRM satellite were not observed subsequently in the inner magnetosphere with the AMPTE-CCE satellite. These results were studied by using a Monte Carlo code to compute the transport of the ions. For each release several hundred ion guiding-center trajectories were computed under simulated magnetospheric conditions, using the Tsyganenko-Usmanov (1982) magnetic-field model and the Millstone Hill convection-electric-field model (Oliver et al., 1983). The corotation and convection electric fields were mapped to altitudes above the ionosphere, assuming the magnetic-field lines to be equipotentials. The initial conditions of the ions, at the times at which the ions were picked up by the electric field, were estimated by taking into consideration the release conditions and the early-time collective effects. The results indicate that the Ba(+) ions were not observed because the CCE satellite was not along the drift paths of the ions, and the Li(+) ions were not observed because their fluxes at the satellite were too low.

Cladis, J. B.↗

Transport of photons produced by lightning in clouds

The optical effects of the light produced by lightning are of interest to atmospheric scientists for a number of reasons. Two techniques are mentioned which are used to explain the nature of these effects: Monte Carlo simulation; and an equivalent medium approach. In the Monte Carlo approach, paths of individual photons are simulated; a photon is said to be scattered if it escapes the cloud, otherwise it is absorbed. In the equivalent medium approach, the cloud is replaced by a single obstacle whose properties are specified by bulk parameters obtained by methods due to Twersky. Herein, Boltzmann transport theory is used to obtain photon intensities. The photons are treated like a Lorentz gas. Only elastic scattering is considered and gravitational effects are neglected. Water droplets comprising a cuboidal cloud are assumed to be spherical and homogeneous. Furthermore, it is assumed that the distribution of droplets in the cloud is uniform and that scattering by air molecules is neglible. The time dependence and five dimensional nature of this problem make it particularly difficult; neither analytic nor numerical solutions are known.

Solakiewicz, Richard↗

Modeling the Impact of Drizzle and 3D Cloud Structure on Remote Sensing of Effective Radius

Remote sensing of cloud particle size with passive sensors like MODIS is an important tool for cloud microphysical studies. As a measure of the radiatively relevant droplet size, effective radius can be retrieved with different combinations of visible through shortwave infrared channels. MODIS observations sometimes show significantly larger effective radii in marine boundary layer cloud fields derived from the 1.6 and 2.1 pm channel observations than for 3.7 pm retrievals. Possible explanations range from 3D radiative transport effects and sub-pixel cloud inhomogeneity to the impact of drizzle formation on the droplet distribution. To investigate the potential influence of these factors, we use LES boundary layer cloud simulations in combination with 3D Monte Carlo simulations of MODIS observations. LES simulations of warm cloud spectral microphysics for cases of marine stratus and broken stratocumulus, each for two different values of cloud condensation nuclei density, produce cloud structures comprising droplet size distributions with and without drizzle size drops. In this study, synthetic MODIS observations generated from 3D radiative transport simulations that consider the full droplet size distribution will be generated for each scene. The operational MODIS effective radius retrievals will then be applied to the simulated reflectances and the results compared with the LES microphysics.

Platnick, Steven↗

Monte-Carlo Estimation of the Inflight Performance of the GEMS Satellite X-Ray Polarimeter

We report a Monte-Carlo estimation of the in-orbit performance of a cosmic X-ray polarimeter designed to be installed on the focal plane of a small satellite. The simulation uses GEANT for the transport of photons and energetic particles and results from Magboltz for the transport of secondary electrons in the detector gas. We validated the simulation by comparing spectra and modulation curves with actual data taken with radioactive sources and an X-ray generator. We also estimated the in-orbit background induced by cosmic radiation in low Earth orbit.

Monte-Carlo simulation↗

PUFFIn – A user friendly fast interface for calculating and visualizing the dose distribution in materials

This article describes a new graphical user interface that uses the PENELOPE Monte Carlo code to calculate dose distributions in materials. It is named PUFFIn, for the Penelope User-Friendly Fast Interface and was developed as an educational and scoping analysis tool for non-experts in radiation modeling. PUFFIn enables the user to visualize and compare the dose distributions in objects (e.g., sterilized healthcare products) irradiated with cobalt-60 gamma-rays, electron beam (E-beam) or X-rays. From such comparisons, the user can determine the most efficient product and/or packaging designs for any given radiation field – whether for conceptual or existing products and packaging. PUFFIn is distributed in a complete, self-contained package of software, including the PENELOPE radiation transport code, a standard graphics package, and a set of simple exercises. The package is available to any user at no cost and requires a minimal amount of training compared to other similar software. PUFFIn's capabilities are described, as well as validation measurements performed at Texas A&M University E-beam facility and the Aerial E-beam Facility in France.

97 MATHEMATICS AND COMPUTING↗

Light Transfers Through a Koch Shape Cloud

Abstract Modeling radiative transfer in a 3D cloudy atmosphere is critical to climate projections. A recently developed fast 3D radiation parameterization scheme gains some success in quantifying horizontal radiative transfer through cloud sides using cloud area fraction. Based on 3D Monte Carlo simulations of radiative transfer through an idealized single‐layer cloud with Koch‐shaped fractal geometry edges, here we show that radiative energy transport through cloud sides correlates more significantly with cloud area fraction than with cloud perimeter length. The results exemplify the importance of accounting for the horizontal radiative energy exchanges between cloud‐free and cloudy regions with cloud area fraction. Results from additional sensitivity simulations show that increased cloud vertical extent often enhances cloud‐side sunlight leak more significantly than cloud‐side sunlight interception. At low sun elevations, cloud‐side sunlight interception is enhanced more than cloud‐side sunlight leak does with the increase of cloud mass.

58 GEOSCIENCES↗

Modeling of the Advanced Test Reactor Using OpenMC, Cubit and Griffin

In the pursuit of the ability to perform multiphysics simulations of the Advanced Test Reactor, it is crucial to have a fast and highly accurate deterministic model. To achieve this, a contemporary two-step method is employed. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor multiphysics application based on the Multiphysics Object-Oriented Simulation Environment. To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material IDs are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Preliminary comparisons indicate good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 50 pcm in the two-dimensional geometry configuration. However, in three-dimensional calculations, an unacceptably large error is found in the Griffin solution. Future work is planned to resolve this discrepancy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Nucleosynthesis in Outflows from Black Hole–Neutron Star Merger Disks with Full GR(ν)RMHD

Abstract Along with binary neutron star mergers, the inspiral and merger of a black hole and a neutron star is a predicted site of r -process nucleosynthesis and associated kilonovae. For the right mass ratio, very large amounts of neutron-rich material (relative to the dynamical ejecta) may become unbound from the post-merger accretion disk. We simulate a suite of four post-merger disks with three-dimensional general-relativistic magnetohydrodynamics with time-dependent Monte Carlo neutrino transport. We find that within 10 4 GM BH / c 3 (∼200–500 ms), the outflows from these disks are very close to the threshold conditions for robust r -process nucleosynthesis. For these conditions, the detailed properties of the outflow determine whether a full r -process can or cannot occur, implying that a wide range of observable phenomena is possible. We show that on average the disk outflow lanthanide fraction is suppressed relative to the solar isotopic pattern. In combination with the dynamical ejecta, these outflows imply a kilonova with both blue and red components.

79 ASTRONOMY AND ASTROPHYSICS↗

Unraveling small-scale defects in irradiated ThO 2 using kinetic Monte Carlo simulations

Point defects and their clusters generated through irradiation can have significant impact on the physical and mechanical properties of materials. However, direct experimental visualization of these small-scale defects using high-resolution scanning transmission electron microscopy remains a challenging task. In this study, using thorium dioxide (ThO 2 ) with the fluorite structure as a model system, we demonstrate the use of ab initio basin-hopping simulations in synergy with object kinetic Monte Carlo simulations as a powerful tool for identifying small defect complexes in irradiated materials. In addition to providing quantitative insights into defect evolution in ThO 2 under irradiation, our study reveals an unexpected role of bound anti-Schottky defect clusters in mediating defect transport. Remarkably, despite their poor thermal stability against dissociation at high temperatures, the transient formation of bound anti-Schottky defects under irradiation and their subsequent migration provide the dominant mechanism for the growth of large interstitial loops that have been experimentally observed in ThO 2 .

36 MATERIALS SCIENCE↗

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗