SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS2007

Three-dimensional Model of Tissue and Heavy Ions Effects

A three-dimensional tissue model was incorporated into a new Monte Carlo algorithm that simulates passage of heavy ions in a tissue box . The tissue box was given as a realistic model of tissue based on confocal microscopy images. The action of heavy ions on the cellular matrix for 2- or 3-dimensional cases was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix), and as a multi-layer obtained from confocal microscopy (3-d case). Image segmentation was used to identify cells with precise areas/volumes in an irradiated cell culture monolayer, and slices of tissue with many cell layers. The cells were then inserted into the model box of the simulated physical space pixel by pixel. In the case of modeled tissues (3-d), the tissue box had periodic boundary conditions imposed, which extrapolates the technique to macroscopic volumes of tissue. For the real tissue (3-d), specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated from current experimental data. A spatial correlation function indicating a higher spatial concentration of damaged cells from heavy ions relative to the low-LET radiation cell damage pattern is presented. The spatial correlation effects among necrotic cells can help studying microlesions in organs, and probable effects of directionality of heavy ion radiation on epithelium and endothelium.

Ponomarev, Artem L.↗

NASA NTRS1985

Modeling a high energy gamma-ray telescope

A Monte Carlo program was developed to simulate a high energy gamma-ray telescope using a coded aperture mask (CAMTRAC). The ultimate purpose of the calculation will be to determine the optimum design parameters for such an instrument. The model can also used to ascertain under what conditions CAMTRAC performance is superior to that of conventional telescopes which employ electron-positron pair direction as the sole means of determining gamma ray direction.

Jenkins, T. L.↗

NASA NTRSDate not supplied

A Modern Load Relief Guidance Scheme for Space Launch Vehicles

Launch vehicle load relief algorithms are concerned with realizing a reduction of transient bending moments near maximum dynamic pressure. Traditional approaches to load relief typically use inner-loop acceleration feedback to reduce the wind-induced angle of attack. When implemented in the inner loop, load relief bandwidth is necessarily limited by the achievable stability margins, and when acceleration feedback is employed, by the uncertainty associated with structural modes that couple with the body-mounted accelerometer. The structure of inner loop load relief increases the dimensionality of the flight control gain and filter optimization problem. Most importantly, classical load relief laws do not take advantage of high-rate and high-accuracy GPS-aided inertial velocity data that is readily available from modern strap down IMUs. In this paper, a novel load relief guidance scheme is described that uses direct angle-of-attack feedback in a clever mechanization. The steering commands are determined by examining the wind-perturbed dynamics of a launch vehicle with respect to a gravity turn ascent trajectory. An angle of attack estimate is derived from GPS-aided inertial data and pre-launch range wind measurements, and it is shown that a reduction worst-case rigid-body loads can be realized without requiring air data. The algorithm also includes a high-rate navigation data preprocessing scheme that operates directly on the IMU delta-theta and delta-velocity measurements in order to produce a filtered acceleration estimate at the vehicle center of mass. The outer-loop guidance scheme simplifies the design process for the classical inner-loop autopilot. Algorithm performance is demonstrated using Monte Carlo analysis of a representative liquid booster in a production high fidelity launch vehicle simulation.

NESC↗

NASA NTRS1995

LES, DNS and RANS for the analysis of high-speed turbulent reacting flows

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES), direct numerical simulation (DNS), and Reynolds averaged Navier Stokes (RANS) methods for the computational analysis of high-speed reacting turbulent flows. In the second phase of this work, covering the period 1 Aug. 1994 - 31 Jul. 1995, we have focused our efforts on two programs: (1) developments of explicit algebraic moment closures for statistical descriptions of compressible reacting flows and (2) development of Monte Carlo numerical methods for LES of chemically reacting flows.

Adumitroaie, V.↗

NASA NTRS2007

Spatial Pattern of Cell Damage in Tissue from Heavy Ions

A new Monte Carlo algorithm was developed that can model passage of heavy ions in a tissue, and their action on the cellular matrix for 2- or 3-dimensional cases. The build-up of secondaries such as projectile fragments, target fragments, other light fragments, and delta-rays was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix). A simple model of tissue was given as abstract spheres with close approximation to real cell geometries (3-d cell matrix), as well as a realistic model of tissue was proposed based on microscopy images. Image segmentation was used to identify cells in an irradiated cell culture monolayer, or slices of tissue. The cells were then inserted into the model box pixel by pixel. In the case of cell monolayers (2-d), the image size may exceed the modeled box size. Such image was is moved with respect to the box in order to sample as many cells as possible. In the case of the simple tissue (3-d), the tissue box is modeled with periodic boundary conditions, which extrapolate the technique to macroscopic volumes of tissue. For real tissue, specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated based on BNL data, and other experimental data.

Ponomarev, Artem L.↗

NASA NTRS2018

Heasim and Skyback Simulation Tools and Their Application to the Hitomi Mission

We present an introduction to the heasim multimission observation and skyback background, high-energy pseudo Monte Carlo astrophysical simulation tools. Heasim may be used to accurately and efficiently construct flexible image transport system (FITS) event files for simple or composite sources with a wide range of standard and user-defined spatial, spectral, and temporal characteristics. Skyback is designed to enable users to assess the impact of background discrete and diffuse emission on prospective observations, and skyback output may be directly input into heasim. We present a brief overview of heasim and skyback input, algorithms, usage, and output. We also introduce the sxsbranch tool that computes Hitomi soft X-ray spectrometer resolution grade branching ratios, emphasizing its application to simulations. We include several examples of particular relevance to the Hitomi mission.

Loewenstein, Michael↗

NASA NTRS2010

Monte-Carlo Simulation of Heavy Ion Track Structure Calculation of Local Dose and 3D Time Evolution of Radiolytic Species

Heavy ions have gained considerable importance in radiotherapy due to their advantageous dose distribution profile and high Relative Biological Effectiveness (RBE). Heavy ions are difficult to produce on Earth, but they are present in space and it is impossible at this moment to completely shield astronauts from them. The risk of these radiations is poorly understood, which is a concern for a 3-years Mars mission. The effects of radiation are mainly due to DNA damage such as DNA double-strand breaks (DSBs), although non-targeted effects are also very important. DNA can be damaged by the direct interaction of radiation and by reactions with chemical species produced by the radiolysis of water. The energy deposition is of crucial importance to understand biological effects of radiation. Therefore, much effort has been done recently to improve models of radiation tracks.

Plante, Ianik↗

NASA NTRS1993

Parallelizing Timed Petri Net simulations

The possibility of using parallel processing to accelerate the simulation of Timed Petri Nets (TPN's) was studied. It was recognized that complex system development tools often transform system descriptions into TPN's or TPN-like models, which are then simulated to obtain information about system behavior. Viewed this way, it was important that the parallelization of TPN's be as automatic as possible, to admit the possibility of the parallelization being embedded in the system design tool. Later years of the grant were devoted to examining the problem of joint performance and reliability analysis, to explore whether both types of analysis could be accomplished within a single framework. In this final report, the results of our studies are summarized. We believe that the problem of parallelizing TPN's automatically for MIMD architectures has been almost completely solved for a large and important class of problems. Our initial investigations into joint performance/reliability analysis are two-fold; it was shown that Monte Carlo simulation, with importance sampling, offers promise of joint analysis in the context of a single tool, and methods for the parallel simulation of general Continuous Time Markov Chains, a model framework within which joint performance/reliability models can be cast, were developed. However, very much more work is needed to determine the scope and generality of these approaches. The results obtained in our two studies, future directions for this type of work, and a list of publications are included.

Nicol, David M.↗

NASA NTRS1994

Predictions of a population of cataclysmic variables in globular clusters

We have studied the number of cataclysmic variables (CVs) that should be active in globular clusters during the present epoch as a result of binary formation via two-body tidal capture. We predict the orbital period and luminosity distributions of CVs in globular clusters. The results arebased on Monte Carlo simulations combined with evolution calculations appropriate to each system formed during the lifetime of two specific globular clusters, omega Cen and 47 Tuc. From our study of these two clusters, which represent the range of core densities and states of mass segregation that are likely to be interesting, we extrapolate our results to the Galactic globlular cluster system. Although there is at present little direct observational evidence of CVs in globular clusters, we find that there should be a large number of active systems. We predict that there should be more than approximately 100 CVs in both 47 Tuc and omega Cen and several thousand in the Galactic globular cluster system. These numbers are based on two-body processes alone and represent a lower bound on the number of systems that may have been formed as a result of stellar interaction within globular clusters. The relation between these calculations and the paucity of optically detected CVs in globular clusters is discussed. Should future observations fail to find convincing evidence of a substantial population of cluster CVs, then the two-body tidal capture scenario is likely to be seriously constrained. Of the CVs we espect in 47 Tuc and omega Cen, approximately 45 and 20, respectively, should have accretion luminosities above 10(exp 33) ergs/s. If one utilizes a relation for converting accretion luminosity to hard X-ray luminosity that is based on observations of Galactic plane CVs, even these sources will not exhibit X-ray luminosities above 10(exp 33) ergs/s. While we cannot account directly for the most luminous subset of the low-luminosity globular cluster X-ray sources without assuming an evolutionary pattern that is different from that of the majority of CVs in the disk, we are able to account for all of the observed lower luminosity subset of these sources, many of which have been recently discovered through ROSAT observations. In order for our predicted integrated cluster X-ray luminosities to be consistent with observational upper limits, the relation between accretion and X-ray luminosities should be something like that inferred from the Galactic plane population of CVs. Our calculations predict a large number of systems with L(sub acc) is less than 10(exp 32) ergs/s. Although our calculations imply that globular clusters should have an enhancement of CVs relative to the number thought to be present in the Galactic disk, this enhancement is at most roughly an order of magnitude, not comparable to the factor of approximately 100 for low-mass X-ray binaries (LMXBs).

Di Stefano, R.↗

NASA NTRS2007

A Neutron Detector for the Electron Calorimeter (ECAL) Long Duration Balloon Experiment

The highest energy measurements of cosmic ray electrons extend just beyond 1 TeV. High energy electrons are of particular interest because energy losses during interstellar propagation insure that they arrive primarily from nearby sources. This may produce observable structure in their spectrum. Further, it is predicted that electrons and positrons result from the annihilation of many exotic particles deposited as dark matter candidates. These electrons may appear as excesses in the cosmic ray electron spectrum from 200 GeV to 1000 GeV. A new long duration balloon experiment, ECAL, is being planned to provide direct cosmic ray electron measurements from approx.50 GeV to >1 TeV. To make these measurements ECAL must discriminate strongly against showers from protons and heavier ions. One of the techniques used to make this discrimination may be based on measuring the secondary neutrons produced by events in the instrument. The neutron detector configuration and technique will be discussed along with its expected performance based on Monte Carlo simulations.

Adams, J. H., Jr.↗

NASA NTRS2004

Ionization Cross Sections and Dissociation Channels of DNA Bases by Electron Collisions

Free secondary electrons are the most abundant secondary species in ionizing radiation. Their role in DNA damage, both direct and indirect, is an active area of research. While indirect damage by free radicals, particularly by the hydroxyl radical generated by electron collision with water. is relatively well studied, damage by direct electron collision with DNA is less well understood. Only recently Boudaiffa et al. demonstrated that electrons at energies well below ionization thresholds can induce substantial yields of single- and double-strand breaks in DNA by a resonant, dissociative attachment process. This study attracted renewed interest in electron collisions with DNA, especially in the low energy region. At higher energies ionization becomes important. While Monte Carlo track simulations of radiation damage always include ionization, the probability of dissociative ionization, i.e., simultaneous ionization and dissociation, is ignored. Just like dissociative attachment, dissociative ionization may be an important contributor to double-strand breaks since the radicals and ions produced by dissociative ionization, located in the vicinity of the DNA coil, can readily interact with other parts of the DNA. Using the improved binary-encounter dipole (iBED) formulation, we calculated the ionization cross sections of the four DNA bases, adenine, cytosine, guanine, and thymine, by electrons at energies from threshold to 1 KeV. The present calculation gives cross sections approximately 20% lower than the results by Bemhardt and Paretzke using the Deutsch-Mark and Binary-Encounter-Bethe (BEB) formalisms. The difference is most likely due to the lack of a shielding term in the dipole potential used in the Deutsch-Mark and BEB formalisms. The dissociation channels of ionization for the bases are currently being studied.

Huo, Winifred M.↗

NASA NTRS2018

A Component-Level Model of Automatic Dependent Surveillance - Broadcast (ADS-B)

Automatic Dependent Surveillance – Broadcast (ADS-B) is being employed in numerous peer-to-peer initiatives attempting to expand the capacity of the National Airspace System (NAS) or enable mixed operations of manned and unmanned vehicles. Safety assessments of these initiatives rely, in part, on modeling the accuracy of ADS-B in reporting the position and direction of an ownship and surrounding traffic. Frequently, these initiatives utilize a position uncertainty model that applies a reported ADS-B estimation position uncertainty (EPU) value to a Rayleigh distribution and uses a Gauss-Markov random walk to add error to the ADS-B output of a vehicle. This model of ADS-B state error is easy to implement and apply to numerous problems. However, it has a couple of draw-backs. First, the ADS-B state errors are equally probable in all directions. This is a good assumption in situations where aircraft maneuvering is not constrained. However, in situations where the aircraft maneuvering is constrained such as landing, the error distribution is likely to exhibit directionality and the non-directional model may skew results especially when assessing very low probabilities (e.g., 10(exp -9)) of catastrophic encounters. Second, the model does not account for processing latency in the receiving aircraft. NASA Langley Research Center (LaRC) recently examined the feasibility of decreasing the spacing of aircraft on parallel approaches to runways separated by as little as 700 feet. For Monte-Carlo analysis using a high-fidelity simulation of a large transport, LaRC started with a Gauss-Markov model of ADS-B error but then developed a component level model of ADS-B error to increase the fidelity of results.

Madden, Michael M.↗

NASA NTRSDate not supplied

Performance Risk Model (PRisM) Proof-of-Concept: An Operational Decision Support Tool to Predict Crew Performance in Space from Available Performance Tests

The Crew Health and Performance-Probabilistic Risk Assessment (CHP-PRA) team at NASA Glenn Research Center has developed a range of tools to evaluate astronaut health during spaceflight and to optimize the medical set required for missions. Among these, the Performance Risk Model (PRisM) represents a novel advancement, extending CHP-PRA’s focus beyond medical systems into the domain of human performance. Such tool could be pivotal in optimizing astronaut capabilities during space travel, thereby enhancing overall mission success. PRisM leverages data from well-established performance assessments conducted during previous crewed space missions and Earth analogs to strategically predict outcomes for planned tasks, even when direct performance testing has not been conducted for those specific tasks. To evaluate performance, PRisM references the various metrics outlined in NASA-STD-3001 and employs a methodology to integrate different performance scales. This framework analyzes the contribution of various human system task categories (HSTCs) to task execution and compares these contributions to the HSTCs’ involvement in other known tests. The model further employs a Monte Carlo simulation to sample performance scores from their distribution in operationally relevant tests such as those in Mulavara et al. (2018) and, by leveraging similarities in HSTC involvement, transfers this knowledge to actual mission tasks, such as those outlined in the "Human Exploration of Mars: Preliminary List of Crew Tasks”. The current PRisM proof-of-concept includes analyses of the impact of exercise and specific medications on astronauts’ performance, with provisions to incorporate additional countermeasures as data becomes available. Furthermore, the tool is customizable to include any system necessary to fully encompass the domain of human systems and can be adapted to evaluate performance for any spaceflight activities as requested by operational stakeholders. PRisM has the potential to assist the Human Research Program in exploring the capabilities trade space for optimized crew performance.

performance modeling↗

NASA NTRS2003

Validation of an In-Water, Tower-Shading Correction Scheme

Large offshore structures used for the deployment of optical instruments can significantly perturb the intensity of the light field surrounding the optical measurement point, where different portions of the visible spectrum are subject to different shadowing effects. These effects degrade the quality of the acquired optical data and can reduce the accuracy of several derived quantities, such as those obtained by applying bio-optical algorithms directly to the shadow-perturbed data. As a result, optical remote sensing calibration and validation studies can be impaired if shadowing artifacts are not fully accounted for. In this work, the general in-water shadowing problem is examined for a particular case study. Backward Monte Carlo (MC) radiative transfer computations- performed in a vertically stratified, horizontally inhomogeneous, and realistic ocean-atmosphere system are shown to accurately simulate the shadow-induced relative percent errors affecting the radiance and irradiance data profiles acquired close to an oceanographic tower. Multiparameter optical data processing has provided adequate representation of experimental uncertainties allowing consistent comparison with simulations. The more detailed simulations at the subsurface depth appear to be essentially equivalent to those obtained assuming a simplified ocean-atmosphere system, except in highly stratified waters. MC computations performed in the simplified system can be assumed, therefore, to accurately simulate the optical measurements conducted under more complex sampling conditions (i.e., within waters presenting moderate stratification at most). A previously reported correction scheme, based on the simplified MC simulations, and developed for subsurface shadow-removal processing of in-water optical data taken close to the investigated oceanographic tower, is then validated adequately under most experimental conditions. It appears feasible to generalize the present tower-specific approach to solve other optical sensor shadowing problems pertaining to differently shaped deployment platforms, and also including surrounding structures and instrument casings.

Hooker, Stanford B.↗

NASA NTRS2020

Development of a Numeric Predictor-Corrector Aerocapture Guidance for Direct Force Control

Direct force control, where the angle of attack and sideslip angle are modulated, has been proposed as an alternative to bank angle control for aero-assist maneuvers. This paper reimplements the current state-of-the-art aerocapture guidance for bank angle control, Fully Numeric Predictor-corrector Aerocapture Guidance (FNPAG), for direct force control. The optimal control theory underlying the structure of FNPAG is shown to not be applicable to the direct force control approach. Several solution structures for the longitudinal channel are compared by simulating dispersed three-degree-of-freedom trajectories for a reference mission consisting of a low lift-to-drag vehicle and a highly elliptical, 1-sol target orbit around Mars. The equations of motion for the lateral channel are derived, and a controller is designed to target a specified orbital plane. Finally, a Monte Carlo is used to demonstrate the performance of the new guidance.

Matz, Daniel A.↗

NASA NTRS1978

Radiative transfer model for remote sensing of suspended sediments in water

A Monte Carlo simulation model of radiative transfer in turbid water is discussed. The model can be used to calculate characteristics of the backscattered signal from an illuminated body of water as a function of the turbidity level and spectral properties of the suspended particulates. The dependence of remote sensing applications on the concentration and spectral properties of sediments in the environmental waters is considered in terms of the model. Attention is directed to the effects of various inputs for the volume-scattering function on backscattered radiance from natural waters. The wavelength dependence of single scattering albedo is investigated.

Ghovanlou, A. H.↗

NASA NTRS1989

Interference effects on the hypersonic, rarefied flow about a flat plate

The Direct Simulatiaon Monte Carlo method is used to study the hypersonic, rarefied flow interference effects on a flat plate caused by nearby surfaces. Calculations focus on shock-boundary-layer and shock-lip interactions in hypersonic inlets. Results are presented for geometries consisting of a flat plate with different leading-edge shapes over a flat lower wall and a blunt-edge flat plate over a 5-degree wedge. The problems simulated correspond to a typical entry flight condition of 7.5 km/s at altitudes of 75 to 90 km. The results show increases in predicted local heating rates for shock-boundary-layer and shock-lip interactions that are quantitatively similar to those observed experimentally at much higher densities.

Wilmoth, Richard G.↗

NASA NTRS1988

Power spectral density analysis of wind-shear turbulence for related flight simulations

Meteorological phenomena known as microbursts can produce abrupt changes in wind direction and/or speed over a very short distance in the atmosphere. These changes in flow characteristics have been labelled wind shear. Because of its adverse effects on aerodynamic lift, wind shear poses its most immediate threat to flight operations at low altitudes. The number of recent commercial aircraft accidents attributed to wind shear has necessitated a better understanding of how energy is transferred to an aircraft from wind-shear turbulence. Isotropic turbulence here serves as the basis of comparison for the anisotropic turbulence which exists in the low-altitude wind shear. The related question of how isotropic turbulence scales in a wind shear is addressed from the perspective of power spectral density (psd). The role of the psd in related Monte Carlo simulations is also considered.

Laituri, Tony R.↗