Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low earth orbit for scientific, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with 8 Teledyne dosimeters over varying shielding effectiveness. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 –AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat 3 unit (U) models. Using NOVICE model radiation analysis coding, the shielding effectiveness's, based on a generic CubeSat 3U structure with 4 electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cmto4.44 cm (0.550 g/cm2–12.0 g/cm2) thick aluminum. For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID)increases by nearly a magnitude for thin-walled aluminum 0.550g/cm2-0.686 g/cm2(0.204 cm –0.254 cm) typical CubeSat structures. The shielding effectiveness by NOVICE Sigma estimates, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from Earth's magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast, to a thin walled generic 3U CubeSat, Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm2from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the spacecraft’s mission value versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗
NASA NTRSDate not supplied

Shielding Considerations for CubeSat Structures During Solar Maximum

The purpose of this lessons learned paper is to communicate the utility of shielding in small spacecraft planning for the support of mission assurance and reliability. Numerous SmallSats have been flying in polar low Earth orbit for science, communications, technology demonstrations, and imaging with academic, commercial, and government interests. Shielding has been part of mission assurance and reliability from the advent of long duration spacecraft missions. The Shields-1 CubeSat has been operating in polar low Earth orbit since 16 December 2018 with atomic number (Z)-grade radiation shielding and demonstrates shielding effectiveness. Shields-1 has collected a representative example of solar minimum data in 2019 with eight Teledyne µDosimeters over varying shielding effectivenesses. It serves as current experimental data and has been compared with NOVICE Shielding estimates using the AP8 – AE8 trapped radiation model with the Shields-1 CAD and generic CubeSat three unit (3U) models. Using NOVICE model radiation analysis coding, the shielding effectivenesses, based on a generic CubeSat 3U structure with four electronic boards, were estimated for aluminum wall thicknesses ranging from 0.204 cm to 4.44 cm (0.550 g/cm 2 to 12.0 g/cm 2 ). For modeled polar orbiting spacecraft, solar maximum total ionizing dose (TID) increases by nearly a magnitude for thin-walled aluminum 0.550-g/cm 2 - 0.686-g/cm 2 (0.204-cm – 0.254-cm) typical CubeSat structures. The shielding effectiveness estimated by NOVICE Sigma, which is a shielding sphere approximation around a detector, showed a linear relationship with wall thickness, which increased over the wall thickness by a ratio of 1.43 determined by linear regression analysis. Using NOVICE Adjoint Monte-Carlo Modeling of solar minimum and solar maximum with the inclusion of a worst-case solar particle event over a one-year mission without geomagnetic shielding, the TID for minimum and maximum conditions for a generic 3U with a wall thickness of 0.254 cm is 158 RAD and 1540 RAD, respectively. The modeled total solar maximum TID is over estimated, because at low orbital latitudes a spacecraft will have shielding from the Earth’s magnetic field. However, TID will still be significant at high latitudes over the poles, where a spacecraft is exposed in a solar particle event. In contrast to a thin-walled generic 3U CubeSat, the Shields-1 electronics enclosure has a shielding effectiveness of 21.3 g/cm 2 from NOVICE Sigma modeling and is expected to show reduced total ionizing dose increases during the present active Solar Cycle 25 period. Because solar particle events during solar maximum increase TID on electronic parts with thin-walled shielding in short periods of time, it is a mission assurance and reliability consideration on the mission value of the spacecraft versus adding shielding for risk reduction of premature spacecraft or instrument payload loss. Since the volumes of many instruments and system electronics have reduced with small spacecraft, shielding material costs and weight penalties have diminished. A small spacecraft project budget and schedule may limit traditional radiation-hardened part use and radiation testing requirements, where shielding can contribute to mission assurance and reliability with reduced costs.

Shields-1↗
NASA NTRS2021

Recent Improvements in the Dart Model for Atmosphere, Topography, Large Landscape, Chlorophyll Fluorescence, Satellite Image Inversion

Physical models simulating the radiative budget (RB) and remote sensing (RS) observation of three-dimensional (3D) landscapes are critical to better understand human and natural components of the Earth system and further develop RS technology. DART is one of the most comprehensive 3D models of Earth-atmosphere optical radiative transfer (RT), from ultraviolet (UV) to thermal infrared (TIR). It simulates the optical signal of proximal, aerial and satellite imaging spectrometers and laser scanners, the 3D RB and solar induced chlorophyll fluorescence (SIF) signal, for any urban or natural landscape and any experimental or instrument configuration. It is freely available for research and teaching activities (https://dart.omp.eu). Here, five recent advances are presented. 1) Atmosphere RT. 2) RT in non repetitive topography. 3) Monte Carlo modelling for fast RS image simulation of large landscapes. 4) SIF modelling for vegetation simulated as facets and turbid cells. 5) RS image inversion for mapping the optical properties of urban material and the urban radiative budget.

DART↗
NASA NTRS1992

Origin and evolution of planetary atmospheres

This report concerns several research tasks related to the origin and evolution of planetary atmospheres and the large-scale distribution of volatile elements in the Solar System. These tasks and their present status are as follows: (1) we have conducted an analysis of the volatility and condensation behavior of compounds of iron, aluminum, and phosphorus in the atmosphere of Venus in response to publish interpretations of the Soviet Venera probe XRF experiment data, to investigate the chemistry of volcanic gases, injection of volatiles by cometary and asteroidal impactors, and reactions in the troposphere; (2) we have completed and are now writing up our research on condensation-accretion modeling of the terrestrial planets; (3) we have laid the groundwork for a detailed study of the effects of water transport in the solar nebula on the bulk composition, oxidation state, and volatile content of preplanetary solids; (4) we have completed an extensive laboratory study of cryovolcanic materials in the outer solar system; (5) we have begun to study the impact erosion and shock alteration of the atmosphere of Mars resulting from cometary and asteroidal bombardment; and (6) we have developed a new Monte Carlo model of the cometary and asteroidal bombardment flux on the terrestrial planets, including all relevant chemical and physical processes associated with atmospheric entry and impact, to assess both the hazards posed by this bombardment to life on Earth and the degree of cross-correlation between the various phenomena (NO(x) production, explosive yield, crater production, iridium signature, etc.) that characterize this bombardment. The purpose of these investigations has been to contribute to the developing understanding of both the dynamics of long-term planetary atmosphere evolution and the short-term stability of planetary surface environments.

Lewis, John S.↗
NASA NTRS2016

Recurrence Rate and Magma Effusion Rate for the Latest Volcanism on Arsia Mons, Mars

Magmatism and volcanism have evolved the Martian lithosphere, surface, and climate throughout the history of Mars. Constraining the rates of magma generation and timing of volcanism on the surface clarifies the ways in which magma and volcanic activity have shaped these Martian systems. The ages of lava flows on other planets are often estimated using impact crater counts, assuming that the number and size-distribution of impact craters per unit area reflect the time the lava flow has been on the surface and exposed to potential impacts. Here we show that impact crater age model uncertainty is reduced by adding stratigraphic information observed at locations where neighboring lavas abut each other, and demonstrate the significance of this reduction in age uncertainty for understanding the history of a volcanic field comprising 29 vents in the 110-kilometer-diameter caldera of Arsia Mons, Mars. Each vent within this caldera produced lava flows several to tens of kilometers in length; these vents are likely among the youngest on Mars, since no impact craters in their lava flows are larger than 1 kilometer in diameter. First, we modeled the age of each vent with impact crater counts performed on their corresponding lava flows and found very large age uncertainties for the ages of individual vents, often spanning the estimated age for the entire volcanic field. The age model derived from impact crater counts alone is broad and unimodal, with estimated peak activity in the field around 130Ma (megaannum, 1 million years). Next we applied our volcano event age model (VEAM), which uses a directed graph of stratigraphic relationships and random sampling of the impact crater age determinations to create alternative age models. Monte Carlo simulation was used to create 10,000 possible vent age sets. The recurrence rate of volcanism is calculated for each possible age set, and these rates are combined to calculate the median recurrence rate of all simulations. Applying this approach to the 29 volcanic vents, volcanism likely began around 200-300Ma then first peaked around 150Ma, with an average production rate of 0.4 vents per Myr (million years). The recurrence rate estimated including stratigraphic data is distinctly bimodal, with a second, lower peak in activity around 100Ma. Volcanism then waned until the final vents were produced 10-90Ma. Based on this model, volume flux is also bimodal, reached a peak rate of 1-8 cubic kilometers per million years by 150Ma and remained above half this rate until about 90Ma, after which the volume flux diminished greatly. The onset of effusive volcanism from 200-150Ma might be due to a transition of volcanic style away from explosive volcanism that emplaced tephra on the western flank of Arsia Mons, while the waning of volcanism after the 150Ma peak might represent a larger-scale diminishing of volcanic activity at Arsia Mons related to the emplacement of flank apron lavas.

Richardson, Jacob A.↗
NASA NTRS1987

On the possibility of observing cosmic ray sources in high energy gamma rays

If cosmic rays are accelerated by strong shocks, then cosmic ray sources should be characterized by spectra, dN/dE alpha E exp -(2.0-2.2), reflecting the strength of those shocks. This is expected from the 'standard leaky box' model of cosmic ray propagation in which the source spectra are harder than the observed spectra because higher energy particles have shorter residence times in the galactic magnetic fields. Furthermore, data on cosmic ray nucleons suggest that these sources might be surrounded by material. If the latter is true, such sources should be observable in gamma rays at energies beyond 1 GeV where the angular resolution of gamma-ray telescopes is optimized and the background is significantly reduced. For identified sources, the source location accuracy can be shown to improve with increasing energy in spite of the decreasing statistics, as long as the gamma-ray spectrum is harder than dN/dE alpha E exp -gamma. A Monte Carlo model is used to predict the photon spectra which would be expected from cosmic ray sources under varying assumptions about the strength of the shocks in the acceleration region.

Ormes, J. F.↗
NASA NTRS1999

External Sources of Water for Mercury's Putative Ice Deposits

Radar images have revealed the possible presence of ice deposits in Mercury's polar regions. Although thermal models indicate that water ice can be stable in permanently shaded regions near Mercury's poles, the ultimate source of the water remains unclear. We use stochastic models and other theoretical methods to investigate the role of external sources in supplying Mercury with the requisite amount of water. By extrapolating the current terrestrial influx of interplanetary dust particles to that at Mercury, we find that continual micrometeoritic bombardment of Mercury over the last 3.5 byr could have resulted in the delivery of (3-60) x 10(exp 16) grams of water ice to the permanently shaded regions at Mercury's poles (equivalent to an average ice thickness of 0.8-20 m). Erosion by micrometeoritic impact on exposed ice deposits could reduce the above value by about a half. For comparison, the current ice deposits on Mercury are believed to be somewhere between approx. 2 and 20 m thick. Using a Monte Carlo model to simulate the impact history of Mercury, we find that asteroids and comets can also deliver an amount of water consistent with the observations. Impacts from Jupiter-family comets over the last 3.5 billion years can supply (0.1-200) x 10(exp 16) grams of water to Mercury's polar regions (corresponding to ice deposits 0.05-60 m thick), Halley-type comets can supply (0.2-20) x 10(exp 16) grams of water to the poles (0.07-7 m of ice), and asteroids can provide (0.4-20) x 10(exp 16) grams of water to the poles (0.1-8 m of ice). Although all these external sources are nominally sufficient to explain the estimated amount of ice currently at Mercury's poles, impacts by a few large comets and/or asteroids seem to provide the best explanation for both the amount and cleanliness of the ice deposits on Mercury. Despite their low population estimates in the inner solar system, Jupiter-family comets are particularly promising candidates for delivering water to Mercury because they have a larger volatile content than asteroids and more favorable orbital and impact characteristics than Halley-type comets.

Moses, Julianne I.↗
NASA NTRS1992

Modeling, design, fabrication, and testing of InP Gunn devices in the D-band (110 GHz - 170 GHz)

The development of fundamental Gunn sources for D-band frequencies requires improvements of doping profiles, processing technology, and circuit design. We have developed a technology for fabricating InP Gunn diodes using an InGaAs etch-stop layer between the InP substrate and the device layers. The epitaxial layers were grown by CBE. During device processing, the substrate is completely removed. Substrateless devices with an n(+) InGaAs cap layer are expected to have reduced contact and series resistances, and skin effect losses. This technology gives better uniformity and control of the device geometry across the processed chip. InP Gunn devices with a 1.7 micron long active region (doping : 9 x 10(exp 15) cm(exp -3)) have been mounted on copper heat sinks. Two tapered leads were then bonded to the diode and to four quartz standoffs. As a preliminary result, an output power of 13 mW at 82 GHz was obtained. Based on these RF measurements, we determine appropriate material parameters to be used in the Ensemble Monte Carlo model. Subsequently, we use this model to design and evaluate the performance of InP Gunn Devices for D-band frequencies. Using the same technology, we are currently processing Gunn devices with a 1 micron long active region for operation at higher frequencies.

Kamoua, R.↗
NASA NTRS1991

Low altitude plume impingement handbook

Plume Impingement modeling is required whenever an object immersed in a rocket exhaust plume must survive or remain undamaged within specified limits, due to thermal and pressure environments induced by the plume. At high altitudes inviscid plume models, Monte Carlo techniques along with the Plume Impingement Program can be used to predict reasonably accurate environments since there are usually no strong flowfield/body interactions or atmospheric effects. However, at low altitudes there is plume-atmospheric mixing and potential large flowfield perturbations due to plume-structure interaction. If the impinged surface is large relative to the flowfield and the flowfield is supersonic, the shock near the surface can stand off the surface several exit radii. This results in an effective total pressure that is higher than that which exists in the free plume at the surface. Additionally, in two phase plumes, there can be strong particle-gas interaction in the flowfield immediately ahead of the surface. To date there have been three levels of sophistication that have been used for low altitude plume induced environment predictions. Level 1 calculations rely on empirical characterizations of the flowfield and relatively simple impingement modeling. An example of this technique is described by Piesik. A Level 2 approach consists of characterizing the viscous plume using the SPF/2 code or RAMP2/LAMP and using the Plume Impingement Program to predict the environments. A Level 3 analysis would consist of using a Navier-Stokes code such as the FDNS code to model the flowfield and structure during a single calculation. To date, Level 1 and Level 2 type analyses have been primarily used to perform environment calculations. The recent advances in CFD modeling and computer resources allow Level 2 type analysis to be used for final design studies. Following some background on low altitude impingement, Level 1, 2, and 3 type analysis will be described.

Smith, Sheldon D.↗
NASA NTRS1971

Measurement of the nuclear electromagnetic cascade development in glass at energies above 200 GeV

The longitudinal development of nuclear-electromagnetic cascades with energies greater than 200 GeV was measured in a low-Z (glass) absorber. This was done in the course of operating an ionization spectrometer at mountain altitude in an experiment to study the properties of gamma rays emitted from individual interactions at energies around 10,000 GeV. The ionization produced by a cascade is sampled by 20 sheets of plastic scintillator spaced uniformly in depth every 2.2 radiation lengths. Adjacent pairs of scintillators are viewed by photomultipliers which measure the mean ionization produced by an individual cascade in 10 layers each 1.1 interaction length (4.4 radiation lengths) thick. The longitudinal development of the cascades was measured for about 250 cascades having energies ranging from 200 GeV to 2500 GeV. The observations are compared with the predictions of calculations made for this specific spectrometer using a three-dimensional Monte Carlo model of the nuclear-electromagnetic cascade.

Gillespie, C. R.↗
NASA NTRS1974

Secondary electron background produced by heavy nuclei in a multiwire proportional counter hodoscope

The secondary electron background produced by heavy nuclei in a multiwire proportional counter hodoscope is calculated using both a simplified and a more complete Monte Carlo model. These results are compared with experimental data from a small multiwire proportional counter hodoscope operated in a 530 MeV/nucleon accelerator beam of nitrogen nuclei. Estimates of the secondary electron background produced by heavy relativistic nuclei are presented along with the detailed results from calculations of energy deposition in the hodoscope counter cells.

Morgan, S. H., Jr.↗
NASA NTRS1973

Effects of changing canopy directional reflectance on feature selection

The use of a Monte Carlo model for generating sample directional reflectance data for two simplified target canopies at two different solar positions is reported. Successive iterations through the model permit the calculation of a mean vector and covariance matrix for canopy reflectance for varied sensor view angles. These data may then be used to calculate the divergence between the target distributions for various wavelength combinations and for these view angles. Results of a feature selection analysis indicate that different sets of wavelengths are optimum for target discrimination depending on sensor view angle and that the targets may be more easily discriminated for some scan angles than others. The time-varying behavior of these results is also pointed out.

Smith, J. A.↗
NASA NTRS1975

Error analysis of Dobson spectrophotometer measurements of the total ozone content

A study of techniques for measuring atmospheric ozone is reported. This study represents the second phase of a program designed to improve techniques for the measurement of atmospheric ozone. This phase of the program studied the sensitivity of Dobson direct sun measurements and the ozone amounts inferred from those measurements to variation in the atmospheric temperature profile. The study used the plane - parallel Monte-Carlo model developed and tested under the initial phase of this program, and a series of standard model atmospheres.

Holland, A. C.↗
NASA NTRS1975

Early scattering by Jupiter and its collision effects in the terrestrial zone

When Jupiter was on the order of three to ten earth masses in size, there undoubtedly was a considerably larger mass of condensed matter in its zone, since Jupiter would have perturbed most of it to other parts of the solar system. Monte Carlo studies indicate a significant portion would have crossed the earth's orbit. If the earth and moon had not yet fully formed, the probability of earth-zone planetesimals being hit by this Jupiter-scattered material was high. Further Monte Carlo models of these collisions and their products indicate a significant portion of matter was heated to melting, even if less than 5% of the relative kinetic energy went into heat. The models include capture probabilities by an embryo earth and a protolunar swarm. Because heat energy is correlated with comminution energy, and because the capture probability of the swarm is mass-dependent while embryo's is not, the protolunar material suffered much higher heating on the average than did the proto-earth material.

Kaula, W. M.↗
NASA NTRS1976

Sample size and sampling errors as the source of dispersion in chemical analyses

The paper describes a Monte Carlo model for simulation of two-dimensional representations of thin sections of some of the more common igneous rock textures. These representations are extrapolated to three dimensions to develop a volume of 'rock'. The model (here applied to a medium-grained high-Ti basalt) can be used to determine a statistically significant sample for a lunar rock or to predict the probable errors in the oxide contents that can occur during the analysis of a sample that is not representative of the parent rock.

Clanton, U. S.↗
NASA NTRS1977

Initial energy and perihelion distributions of Oort-cloud comets

A Monte Carlo model of stellar perturbations of the Oort cloud is used to study the distributions in energy and perihelion of comets entering the planetary region for the first time. The model is run for a variety of initial states and a range of velocity perturbations. In all cases the resulting orbits are uniformly distributed in perihelion distance in the planetary region, q less than 20 AU. Most orbits are confined to a fairly narrow range in 1/a and hyperbolic orbits are rare.

Weissman, P. R.↗