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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Simulation of DNA Damage in Human Cells from Space Radiation Using a Physical Model of Stochastic Particle Tracks and Chromosomes

The formation of double-strand breaks (DSBs) and chromosomal aberrations (CAs) is of great importance in radiation research and, specifically, in space applications. We are presenting a recently developed model, in which chromosomes simulated by NASARTI (NASA Radiation Tracks Image) is combined with nanoscopic dose calculations performed with the Monte-Carlo simulation by RITRACKS (Relativistic Ion Tracks) in a voxelized space. The model produces the number of DSBs, as a function of dose for high-energy iron, oxygen, and carbon ions, and He ions. The combined model calculates yields of radiation-induced CAs and unrejoined chromosome breaks in normal and repair deficient cells. The merged computational model is calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. The model considers fractionated deposition of energy to approximate dose rates of the space flight environment. The merged model also predicts of the yields and sizes of translocations, dicentrics, rings, and more complex-type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation.

Ponomarev, Artem↗
NASA NTRS2012

Thermal Phase Variations of WASP-12b: Defying Predictions

We report Warm Spitzer full-orbit phase observations of WASP-12b at 3.6 and 4.5 micrometers. This extremely inflated hot Jupiter is thought to be overflowing its Roche lobe, undergoing mass loss and accretion onto its host star, and has been claimed to have a C/O ratio in excess of unity. We are able to measure the transit depths, eclipse depths, thermal and ellipsoidal phase variations at both wavelengths. The large-amplitude phase variations, combined with the planet's previously measured dayside spectral energy distribution, are indicative of non-zero Bond albedo and very poor day-night heat redistribution. The transit depths in the mid-infrared-(R(sub p)/R(sub ∗))(sup 2) = 0.0123(3) and 0.0111(3) at 3.6 and 4.5 micrometers, respectively-indicate that the atmospheric opacity is greater at 3.6 than at 4.5 micrometers, in disagreement with model predictions, irrespective of C/O ratio. The secondary eclipse depths are consistent with previous studies: F(sub day)/F(sub ∗) = 0.0038(4) and 0.0039(3) at 3.6 and 4.5 micrometers, respectively. We do not detect ellipsoidal variations at 3.6 micrometers, but our parameter uncertainties-estimated via prayer-bead Monte Carlo-keep this non-detection consistent with model predictions. At 4.5 micrometers, on the other hand, we detect ellipsoidal variations that are much stronger than predicted. If interpreted as a geometric effect due to the planet's elongated shape, these variations imply a 3:2 ratio for the planet's longest:shortest axes and a relatively bright day-night terminator. If we instead presume that the 4.5 micrometer ellipsoidal variations are due to uncorrected systematic noise and we fix the amplitude of the variations to zero, the best-fit 4.5 micrometer transit depth becomes commensurate with the 3.6 micrometer depth, within the uncertainties. The relative transit depths are then consistent with a solar composition and short scale height at the terminator. Assuming zero ellipsoidal variations also yields a much deeper 4.5 micrometer eclipse depth, consistent with a solar composition and modest temperature inversion. We suggest future observations that could distinguish between these two scenarios.

planetary↗
NASA NTRS2017

A Fast Hybrid (3-D/1-D) Model for Thermal Radiative Transfer in Cirrus via Successive Orders of Scattering

We investigate the impact of cirrus cloud heterogeneity on the direct emission by cloud or surface and on the scattering by ice particles in the thermal infrared (TIR). Realistic 3-D cirri are modeled with the 3DCLOUD code, and top-of-atmosphere radiances are simulated by the 3-D Monte Carlo radiative transfer (RT) algorithm 3DMCPOL for two (8.65 micrometers and 12.05 micrometers) channels of the Imaging Infrared Radiometer on CALIPSO. At nadir, comparisons of 1-D and 3-D RT show that 3-D radiances are larger than their 1-D counterparts for direct emission but smaller for scattered radiation. For our cirrus cases, 99% of the 3-D total radiance is computed by the third scattering order, which corresponds to 90% of the total computational effort, but larger optical thicknesses need more scattering orders. To radically accelerate the 3-D RT computations (using only few percent of 3-D RT time with a Monte Carlo code), even in the presence of large optical depths, we develop a hybrid model based on exact 3-D direct emission, the first scattering order from 1-D in each homogenized column, and an empirical adjustment linearly dependent on the optical thickness to account for higher scattering orders. Good agreement is found between the hybrid model and the exact 3-D radiances for two very different cirrus models without changing the empirical parameters. We anticipate that a future deterministic implementation of the hybrid model will be fast enough to process multiangle thermal imagery in a practical tomographic reconstruction of 3-D cirrus fields.

radiative transfer↗
NASA NTRS1991

Composite load spectra for select space propulsion structural components

The work performed to develop composite load spectra (CLS) for the Space Shuttle Main Engine (SSME) using probabilistic methods. The three methods were implemented to be the engine system influence model. RASCAL was chosen to be the principal method as most component load models were implemented with the method. Validation of RASCAL was performed. High accuracy comparable to the Monte Carlo method can be obtained if a large enough bin size is used. Generic probabilistic models were developed and implemented for load calculations using the probabilistic methods discussed above. Each engine mission, either a real fighter or a test, has three mission phases: the engine start transient phase, the steady state phase, and the engine cut off transient phase. Power level and engine operating inlet conditions change during a mission. The load calculation module provides the steady-state and quasi-steady state calculation procedures with duty-cycle-data option. The quasi-steady state procedure is for engine transient phase calculations. In addition, a few generic probabilistic load models were also developed for specific conditions. These include the fixed transient spike model, the poison arrival transient spike model, and the rare event model. These generic probabilistic load models provide sufficient latitude for simulating loads with specific conditions. For SSME components, turbine blades, transfer ducts, LOX post, and the high pressure oxidizer turbopump (HPOTP) discharge duct were selected for application of the CLS program. They include static pressure loads and dynamic pressure loads for all four components, centrifugal force for the turbine blade, temperatures of thermal loads for all four components, and structural vibration loads for the ducts and LOX posts.

Newell, J. F.↗
NASA NTRS1996

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗
NASA NTRS1995

Can standard cosmological models explain the observed Abell cluster bulk flow?

Lauer and Postman (LP) observed that all Abell clusters with redshifts less than 15,000 km/s appear to be participating in a bulk flow of 689 km/s with respect to the cosmic microwave background. We find this result difficult to reconcile with all popular models for large-scale structure formation that assume Gaussian initial conditions. This conclusion is based on Monte Carlo realizations of the LP data, drawn from large particle-mesh N-body simulations for six different models of the initial power spectrum (standard, tilted, and Omega(sub 0) = 0.3 cold dark matter, and two variants of the primordial baryon isocurvature model). We have taken special care to treat properly the longest-wavelength components of the power spectra. The simulations are sampled, 'observed,' and analyzed as identically as possible to the LP cluster sample. Large-scale bulk flows as measured from clusters in the simulations are in excellent agreement with those measured from the grid: the clusters do not exhibit any strong velocity bias on large scales. Bulk flows with amplitude as large as that reported by LP are not uncommon in the Monte Carlo data stes; the distribution of measured bulk flows before error bias subtraction is rougly Maxwellian, with a peak around 400 km/s. However the chi squared of the observed bulk flow, taking into account the anisotropy of the error ellipsoid, is much more difficult to match in the simulations. The models examined are ruled out at confidence levels between 94% and 98%.

Strauss, Michael A.↗
NASA NTRS1975

The structure and evolution of Jupiter - The fluid contraction stage

The complete evolution of a contracting star of Jovian mass consisting of a convective adiabatic homogeneous fluid is determined using stellar structure methods, improved model atmosphere calculations, and substantially improved thermodynamic properties for hydrogen and hydrogen-helium fluids. The model atmospheres are calculated in the form of time-averaged vertical temperature structures, including all relevant sources of opacity and a solar energy deposition component, and the thermodynamic properties are modified to obtain better agreement with Monte Carlo results for metallic fluids. The resultant gravitationally contracting evolutionary models are found to have two phases: an early stellar phase similar to a typical low-mass pre-main-sequence body and a later phase constituting an approach to a degenerate-dwarf cooling curve. The first phase is shown to have high luminosities and internal temperatures, while the second gives excellent agreement with the observed radius and luminosity of Jupiter. Analysis indicates that the equation of state and superadiabaticity have the strongest influence on evolution over planetary time scales.

Graboske, H. C., Jr.↗
NASA NTRSDate not supplied

Algorithmic Detection of Elemental Biosignatures

Machine learning models that classify a sample as indicative or non-indicative of life could play an important role in life-detection missions. Their predictions result from agnostic algorithms and thereby add redundancy to judgements resulting from human expertise. Additionally, their important features can reveal the most informative measurements within the operational constraints of a life-detection mission. The Ladder of Life Detection (Neveu 2018) identifies the need for an understanding of how combinations of multiple biosignatures affect overall confidence. The present work provides a starting point to answer this need, and future work will expand the data types to obtain even more predictive combinations of features. Elemental abundance was chosen as a starting set of features due to its availability in diverse sample types, which are needed to train a generalizable model. A standardized dataset was collected, including 35 non-indicative, e.g., lunar rock, basalt; 19 indicative mixed, e.g., seawater, agricultural soil; 46 indicative non-alive, e.g., coal, chalk; and 10 indicative alive, e.g., biofilm, bacteria. This dataset could be valuable for complementary biosignature research. The samples were standardized to the same limit of detection of a simulated mission scenario. Four classification models were used: k-nearest neighbors (KNN), logistic regression (LR), linear support vector machines (SVM), and Gaussian naïve Bayes (GNB). To obtain feature importances, KNN was run on three principal components of the training data and LR and SVM were run with L1 and L2 regularization. The performances and feature importances of the six model variants on 40:60 train to validation ratios were assessed with Monte Carlo simulations. ROC AUC and mean accuracy scores ranged between 82% - 94%, with sensitivity greater than specificity. For indicative of life predictors, all models had C and Ca as strong and Cl as medium; a majority of models had N, K, and P as medium. For non-indicative of life predictors, all models had Si as strong, and a majority of models had Mg, Al, and Ti as medium. Varied elements were Fe (slightly non-indicative), H (slightly indicative), O (widely varied), Na, Mn, and S. These results serve as a proof of concept and suggest important elemental signals beyond merely the CHNOPS of Earth-based life.

Algorithmic↗
NASA NTRS2009

Proton Upset Monte Carlo Simulation

The Proton Upset Monte Carlo Simulation (PROPSET) program calculates the frequency of on-orbit upsets in computer chips (for given orbits such as Low Earth Orbit, Lunar Orbit, and the like) from proton bombardment based on the results of heavy ion testing alone. The software simulates the bombardment of modern microelectronic components (computer chips) with high-energy (.200 MeV) protons. The nuclear interaction of the proton with the silicon of the chip is modeled and nuclear fragments from this interaction are tracked using Monte Carlo techniques to produce statistically accurate predictions.

O'Neill, Patrick M.↗
NASA NTRS1971

The earth and the moon /Harold Jeffreys Lecture/.

The internal structures of the earth and the moon are compared in the light of the latest extensive data on the earth structure, mobility of the earth outer layers, and the properties of lunar crust. The Monte Carlo method is applied to develop an earth model by a stepwise process beginning with a random distribution of two elastic velocities and the density as a function of de pth. Lunar seismic, magnetic, and rock analysis data are used to infer the properties of the moon. The marked planetological contrast between the earth and the moon is shown to consist in that the earth is highly differentiated and still undergoes a large-scale differentiation, while the moon has lost its volatiles in its early history and has a cold dynamically inactive shell which has been without basic changes for three billion years.

Press, F.↗
NASA NTRS1972

The effect of random fluctuations in atmospheric density on significant space shuttle reentry parameters.

The effects of realistic atmospheric density variations on reentry vehicles are investigated with the aid of a random model of atmospheric density and a series of Monte Carlo reentry trajectories generated for a typical delta-wing orbiter design. The random behavior of five major performance indicators is studied, namely, maximum dynamic pressure, maximum g load, maximum stagnation heating rate, down-range distance, and cross-range distance.

Campbell, J. W.↗
NASA NTRS1972

Degree and direction of polarization of multiple scattered light. II - Earth's atmosphere with aerosols.

The radiance, polarization, and direction of polarization of the radiation reflected and transmitted through the atmosphere are calculated by a Monte Carlo method on the basis of a realistic model of the atmosphere. The calculated polarization is shown to depend on the aerosol amount and to have normally a value intermediate between that for pure Rayleigh and pure aerosol scattering. The polarization for multiple scattered photons is usually less than the value calculated for single scattering, but may be larger in low-polarization regions near other regions of high polarization. The direction of polarization shows little variation with either the aerosol amount or the value of the surface albedo.

Plass, G. N.↗
NASA NTRS1974

Atmospheric spreading of protons in auroral arcs

A model is developed to calculate the effect of atmospheric spreading on the flux and angular distribution of protons in homogeneous auroral arcs. An expression is derived that indicates the angular distribution in the atmosphere as a function of distance from arc center, neutral scale height, arc width, and initial angular distribution. The results of the model agree favorably with those based on Monte-Carlo calculations. From these results the enhancement factors needed to compute the original proton current above the atmosphere are obtained. A technique is indicated for determining the incident angular distribution from rocket-based measurements of the arc width and angular distribution.

Iglesias, G. E.↗
NASA NTRS1988

Nonequilibrium thermal radiation for an aeroassist flight experiment vehicle

The direct-simulation Monte Carlo method incorporating a dissociating and ionizing gas model for air with thermal radiation is used to characterize the hypersonic flow about an axisymmetric representation of an aeroassist flight experiment (AFE) vehicle, whose freestream conditions correspond to selected points along the entry, aerobraking, and exit phases of the trajectory. Calculations for two trajectory conditions indicate that the radiative heating of the AFE forebody is lower than the convective heating, but becomes significant as the maximum convective heating rate condition is approached.

Moss, James N.↗
NASA NTRS1990

Efficient Monte Carlo simulation of rarefied flow in a small nozzle

The flow of carbon dioxide through a small conical nozzle is computed using two different numerical approaches. The first of these is the direct simulation Monte Carlo method in which the physical phenomena are modeled at the molecular level. The second technique solves the Navier-Stokes equations of fluid flow. The purpose of the paper is to compare the solutions generated by these two methods for the conditions prevailing in a small, low density nozzle. The comparison between the two sets of solutions reveals only small differences. However, the flow properties in the thick viscous layer close to the nozzle at the exit plane show important deviations which have implications for the subsequent determination of the plume backflow.

Boyd, Iain D.↗
NASA NTRS1995

Orientation and velocity dependence of the nonequilibrium partition coefficient

Monte Carlo simulations based on a Spin-1 Ising Model for binary alloys have been used to investigate the non-equilibrium partition coefficient (k(sub neq)) as a function of solid-liquid interface velocity and orientation. In simulations of Si with a second component k(sub neq) is greater in the [111] direction than the [100] direction in agreement with experimental results reported by Azlz et al. The simulated partition coefficient scales with the square of the step velocity divided by the diffusion coefficient of the secondary component in the liquid.

Beatty, K. M.↗