Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo method”

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 685 records · Page 38

Functional nanocomposites prepared by self-assembly and polymerization of diacetylene surfactants and silicic acid

Conjugated polymer/silica nanocomposites with hexagonal, cubic, or lamellar mesoscopic order were synthesized by self-assembly using polymerizable amphiphilic diacetylene molecules as both structure-directing agents and monomers. The self-assembly procedure is rapid and incorporates the organic monomers uniformly within a highly ordered, inorganic environment. By tailoring the size of the oligo(ethylene glycol) headgroup of the diacetylene-containing surfactant, we varied the resulting self-assembled mesophases of the composite material. The nanostructured inorganic host altered the diacetylene polymerization behavior, and the resulting nanocomposites show unique thermo-, mechano-, and solvatochromic properties. Polymerization of the incorporated surfactants resulted in polydiacetylene (PDA)/silica nanocomposites that were optically transparent and mechanically robust. Molecular modeling and quantum calculations and (13)C spin-lattice relaxation times (T(1)) of the PDA/silica nanocomposites indicated that the surfactant monomers can be uniformly organized into precise spatial arrangements prior to polymerization. Nanoindentation and gas transport experiments showed that these nanocomposite films have increased hardness and reduced permeability as compared to pure PDA. Our work demonstrates polymerizable surfactant/silica self-assembly to be an efficient, general approach to the formation of nanostructured conjugated polymers. The nanostructured inorganic framework serves to protect, stabilize, and orient the polymer, mediate its performance, and provide sufficient mechanical and chemical stability to enable integration of conjugated polymers into devices and microsystems.

NASA Discipline Life Support Systems↗

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

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

NASA Discipline Radiation Health↗

Noninvasive identification of the total peripheral resistance baroreflex

We propose two identification algorithms for quantitating the total peripheral resistance (TPR) baroreflex, an important contributor to short-term arterial blood pressure (ABP) regulation. Each algorithm analyzes beat-to-beat fluctuations in ABP and cardiac output, which may both be obtained noninvasively in humans. For a theoretical evaluation, we applied both algorithms to a realistic cardiovascular model. The results contrasted with only one of the algorithms proving to be reliable. This algorithm was able to track changes in the static gains of both the arterial and cardiopulmonary TPR baroreflex. We then applied both algorithms to a preliminary set of human data and obtained contrasting results much like those obtained from the cardiovascular model, thereby making the theoretical evaluation results more meaningful. This study suggests that, with experimental testing, the reliable identification algorithm may provide a powerful, noninvasive means for quantitating the TPR baroreflex. This study also provides an example of the role that models can play in the development and initial evaluation of algorithms aimed at quantitating important physiological mechanisms.

Non-NASA Center↗

Radiation breakage of DNA: a model based on random-walk chromatin structure

Monte Carlo computer software, called DNAbreak, has recently been developed to analyze observed non-random clustering of DNA double strand breaks in chromatin after exposure to densely ionizing radiation. The software models coarse-grained configurations of chromatin and radiation tracks, small-scale details being suppressed in order to obtain statistical results for larger scales, up to the size of a whole chromosome. We here give an analytic counterpart of the numerical model, useful for benchmarks, for elucidating the numerical results, for analyzing the assumptions of a more general but less mechanistic "randomly-located-clusters" formalism, and, potentially, for speeding up the calculations. The equations characterize multi-track DNA fragment-size distributions in terms of one-track action; an important step in extrapolating high-dose laboratory results to the much lower doses of main interest in environmental or occupational risk estimation. The approach can utilize the experimental information on DNA fragment-size distributions to draw inferences about large-scale chromatin geometry during cell-cycle interphase.

Chromatin/radiation effects↗

The response of tissue-equivalent proportional counters to heavy ions

The paper presents a theoretical model for the response of a tissue-equivalent proportional counter (TEPC) irradiated with charged particles. Heavy ions and iron ions in particular constitute a significant part of radiation in space. TEPCs are used for all space shuttle and International Space Station (ISS) missions to estimate the dose and radiation quality (in terms of lineal energy) inside spacecraft. The response of the tissue-equivalent proportional counters shows distortions at the wall/cavity interface. In this paper, we present microdosimetric investigation using Monte Carlo track structure calculations to simulate the response of a TEPC to charged particles of various LET (1 MeV protons, 2.4 MeV alpha particles, 46 MeV/nucleon 20Ne, 55 MeV/nucleon 20Ne, 45 MeV/nucleon 40Ar, and 1.05 GeV/nucleon 56Fe). Data are presented for energy lost and energy absorbed in the counter cavity and wall. The model calculations are in good agreement with the results of Rademacher et al. (Radiat. Res. 149, 387-389, 1998), including the study of the interface between the wall and the sensitive region of the counter. It is shown that the anomalous response observed at large event sizes in the experiment is due to an enhanced entry of secondary electrons from the wall into the gas cavity.

NASA Discipline Radiation Health↗

Radiation dosimetry using three-dimensional optical random access memories

Three-dimensional optical random access memories (3D ORAMs) are a new generation of high-density data storage devices. Binary information is stored and retrieved via a light induced reversible transformation of an ensemble of bistable photochromic molecules embedded in a polymer matrix. This paper describes the application of 3D ORAM materials to radiation dosimetry. It is shown both theoretically and experimentally, that ionizing radiation in the form of heavy charged particles is capable of changing the information originally stored on the ORAM material. The magnitude and spatial distribution of these changes are used as a measure of the absorbed dose, particle type and energy. The effects of exposure on 3D ORAM materials have been investigated for a variety of particle types and energies, including protons, alpha particles and 12C ions. The exposed materials are observed to fluoresce when exposed to laser light. The intensity and the depth of the fluorescence is dependent on the type and energy of the particle to which the materials were exposed. It is shown that these effects can be modeled using Monte Carlo calculations. The model provides a better understanding of the properties of these materials. which should prove useful for developing systems for charged particle and neutron dosimetry/detector applications. c2001 Published by Elsevier Science B.V.

Non-NASA Center↗

Response of Silicon-Based Linear Energy Transfer Spectrometers: Implication for Radiation Risk Assessment in Space Flights

There is considerable interest in developing silicon-based telescopes because of their compactness and low power requirements. Three such telescopes have been flown on board the Space Shuttle to measure the linear energy transfer spectra of trapped, galactic cosmic ray, and solar energetic particles. Dosimeters based on single silicon detectors have also been flown on the Mir orbital station. A comparison of the absorbed dose and radiation quality factors calculated from these telescopes with that estimated from measurements made with a tissue equivalent proportional counter show differences which need to be fully understood if these telescopes are to be used for astronaut radiation risk assessments. Instrument performance is complicated by a variety of factors. A Monte Carlo-based technique was developed to model the behavior of both single element detectors in a proton beam, and the performance of a two-element, wide-angle telescope, in the trapped belt proton field inside the Space Shuttle. The technique is based on: (1) radiation transport intranuclear-evaporation model that takes into account the charge and angular distribution of target fragments, (2) Landau-Vavilov distribution of energy deposition allowing for electron escape, (3) true detector geometry of the telescope, (4) coincidence and discriminator settings, (5) spacecraft shielding geometry, and (6) the external space radiation environment, including albedo protons. The value of such detailed modeling and its implications in astronaut risk assessment is addressed. c2001 Elsevier Science B.V. All rights reserved.

manned↗

Monte Carlo predictions of DNA fragment-size distributions for large sizes after HZE particle irradiation

DSBs (double-strand breaks) produced by densely ionizing space radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. DSB clustering at large scales, from >100 Mbp down to approximately 2 kbp, is modeled using a Monte-Carlo algorithm. A random-walk model of chromatin is combined with a track model, that predicts the radial distribution of energy from an ion, and the RLC (randomly-located-clusters) formalism, in software called DNAbreak. This model generalizes the random-breakage model, whose broken-stick fragment-size distribution is applicable to low-LET radiation. DSB induction due to track interaction with the DNA volume depends on the radiation quality parameter Q. This dose-independent parameter depends only weakly on LET. Multi-track, high-dose effects depend on the cluster intensity parameter lambda, proportional to fluence as defined by the RLC formalism. After lambda is determined by a numerical experiment, the model reduces to one adjustable parameter Q. The best numerical fits to the experimental data, determining Q, are obtained. The knowledge of lambda and Q allows us to give biophysically based extrapolations of high-dose DNA fragment-size data to low doses or to high LETs.

NASA Center JSC↗

Monte Carlo track structure for radiation biology and space applications

Over the past two decades event by event Monte Carlo track structure codes have increasingly been used for biophysical modelling and radiotherapy. Advent of these codes has helped to shed light on many aspects of microdosimetry and mechanism of damage by ionising radiation in the cell. These codes have continuously been modified to include new improved cross sections and computational techniques. This paper provides a summary of input data for ionizations, excitations and elastic scattering cross sections for event by event Monte Carlo track structure simulations for electrons and ions in the form of parametric equations, which makes it easy to reproduce the data. Stopping power and radial distribution of dose are presented for ions and compared with experimental data. A model is described for simulation of full slowing down of proton tracks in water in the range 1 keV to 1 MeV. Modelling and calculations are presented for the response of a TEPC proportional counter irradiated with 5 MeV alpha-particles. Distributions are presented for the wall and wall-less counters. Data shows contribution of indirect effects to the lineal energy distribution for the wall counters responses even at such a low ion energy.

NASA Center JSC↗

Extrapolation of the dna fragment-size distribution after high-dose irradiation to predict effects at low doses

The patterns of DSBs induced in the genome are different for sparsely and densely ionizing radiations: In the former case, the patterns are well described by a random-breakage model; in the latter, a more sophisticated tool is needed. We used a Monte Carlo algorithm with a random-walk geometry of chromatin, and a track structure defined by the radial distribution of energy deposition from an incident ion, to fit the PFGE data for fragment-size distribution after high-dose irradiation. These fits determined the unknown parameters of the model, enabling the extrapolation of data for high-dose irradiation to the low doses that are relevant for NASA space radiation research. The randomly-located-clusters formalism was used to speed the simulations. It was shown that only one adjustable parameter, Q, the track efficiency parameter, was necessary to predict DNA fragment sizes for wide ranges of doses. This parameter was determined for a variety of radiations and LETs and was used to predict the DSB patterns at the HPRT locus of the human X chromosome after low-dose irradiation. It was found that high-LET radiation would be more likely than low-LET radiation to induce additional DSBs within the HPRT gene if this gene already contained one DSB.

NASA Discipline Radiation Health↗

Envelope structure of deeply embedded young stellar objects in the Serpens Molecular Cloud

Aperture-synthesis and single-dish (sub-) millimeter molecular-line and continuum observations reveal in great detail the envelope structure of deeply embedded young stellar objects (SMM 1 = FIRS 1, SMM 2, SMM 3, SMM 4) in the densely star-forming Serpens Molecular Cloud. SMM 1, 3, and 4 show partially resolved (>2" = 800 AU) continuum emission in the beam of the Owens Valley Millimeter Array at lambda = 3.4-1.4 mm. The continuum visibilities accurately constrain the density structure in the envelopes, which can be described by a radial power law with slope -2.0 +/- 0.5 on scales of 300 to 8000 AU. Inferred envelope masses within a radius of 8000 AU are 8.7, 3.0, and 5.3 Msolar for SMM 1, 3, and 4, respectively. A point source with 20%-30% of the total flux at 1.1 mm is required to fit the observations on long baselines, corresponding to warm envelope material within approximately 100 AU or a circumstellar disk. No continuum emission is detected interferometrically toward SMM 2, corresponding to an upper limit of 0.2 Msolar assuming Td = 24 K. The lack of any compact dust emission suggests that the SMM 2 core does not contain a central protostar. Aperture-synthesis observations of the 13CO, C18O, HCO+, H13CO+, HCN, H13CN, N2H+ 1-0, SiO 2-1, and SO 2(2)-1(1) transitions reveal compact emission toward SMM 1, 3, and 4. SMM 2 shows only a number of clumps scattered throughout the primary field of view, supporting the conclusion that this core does not contain a central star. The compact molecular emission around SMM 1, 3, and 4 traces 5"-10" (2000-4000 AU) diameter cores that correspond to the densest regions of the envelopes, as well as material directly associated with the molecular outflow. Especially prominent are the optically thick HCN and HCO+ lines that show up brightly along the walls of the outflow cavities. SO and SiO trace shocked material, where their abundances may be enhanced by 1-2 orders of magnitude over dark-cloud values. A total of 31 molecular transitions have been observed with the James Clerk Maxwell and Caltech Submillimeter telescopes in the 230, 345, 490, and 690 GHz atmospheric windows toward all four sources, containing, among others, lines of CO, HCO+, HCN, H2CO, SiO, SO, and their isotopomers. These lines show 20-30 km s-1 wide line wings, deep and narrow (1-2 km s-1) self-absorption, and 2-3 km s-1 FWHM line cores. The presence of highly excited lines like 12CO 4-3 and 6-5, 13CO 6-5, and several H2CO transitions indicates the presence of material with temperatures > or approximately 100 K. Monte Carlo calculations of the molecular excitation and line transfer show that the envelope model derived from the dust emission can successfully reproduce the observed line intensities. The depletion of CO in the cold gas is modest compared to values inferred in objects like NGC 1333 IRAS 4, suggesting that the phase of large depletions through the entire envelope is short lived and may be influenced by the local star formation density. Emission in high-excitation lines of CO and H2CO requires the presence of a small amount of approximately 100 K material, comprising less than 1% of the total envelope mass and probably associated with the outflow or the innermost region of the envelope. The derived molecular abundances in the warm (Tkin > 20 K) envelope are similar to those found toward other class 0 YSOs like IRAS 16293-2422, though some species appear enhanced toward SMM 1. Taken together, the presented observations and analysis provide the first comprehensive view of the physical and chemical structure of the envelopes of deeply embedded young stellar objects in a clustered environment on scales between 1000 and 10,000 AU.

NASA Discipline Exobiology↗

An analysis of energy deposition in a tissue equivalent proportional counter onboard the space shuttle

An improved prediction for space radiations in the lower earth orbits measured by the shuttle TEPC is obtained when energy loss straggling and chord length distribution of the detector are considered. A generalized analytic model is used to describe the energy deposition of direct ion interaction events in a micron-size detector. The transport calculation accounting for the shuttle configuration is accomplished by using a new version of HZETRN that has been extensively verified with laboratory and flight data. The agreement of predicted and measured lineal energy spectra is within 70% for the region above 2 keV/micrometer but within a factor of 2.3 underpredicted for the region below this value. The inclusion of indirect delta ray events in the model is needed before possible causes for the underprediction below 2 keV/micrometer can be assessed.

Flight Experiment↗

Kinetic dissection of individual steps in the poly(C)-directed oligoguanylate synthesis from guanosine 5'-monophosphate 2-methylimidazolide

A kinetic study of oligoguanylate synthesis on a polycytidylate template, poly(C), as a function of the concentration of the activated monomer, guanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpG, is reported. Reactions were run with 0.005-0.045 M 2-MeImpG in the presence of 0.05 M poly(C) at 23 degrees C. The kinetic results are consistent with a reaction scheme (eq 1) that consists of a series of consecutive steps, each step representing the addition of one molecule of 2-MeImpG to the growing oligomer. This scheme allows the calculation of second-order rate constants for every step by analyzing the time-dependent growth of each oligomer. Computer simulations of the course of reaction based on the determined rate constants and eq 1 are in excellent agreement with the product distributions seen in the HPLC profiles. In accord with an earlier study (Fakhrai, H.; Inoue, T.; Orgel, L. E. Tetrahedron 1984, 40, 39), rate constants, ki, for the formation of the tetramer and longer oligomers up to the 16-mer were found to be independent of length and somewhat higher than k3 (formation of trimer), which in turn is much higher than k2 (formation of dimer). The ki (i > or = 4), k3, and k2 values are not true second-order rate constants but vary with monomer concentration. Mechanistic models for the dimerization (Scheme I) and elongation reactions (Scheme II) are proposed that are consistent with our results. These models take into account that the monomer associates with the template in a cooperative manner. Our kinetic analysis allowed the determination of rate constants for the elementary processes of covalent bond formation between two monomers (dimerization) and between an oligomer and a monomer (elongation) on the template. A major conclusion from our study is that bond formation between two monomer units or between a primer and a monomer is assisted by the presence of additional next-neighbor monomer units. This is consistent with recent findings with hairpin oligonucleotides (Wu, T.; Orgel, L. E. J. Am. Chem. Soc. 1992, 114, 317). Our study is the first of its kind that shows the feasibility of a thorough kinetic analysis of a template-directed oligomerization and provides a detailed mechanistic model of these reactions.

Non-NASA Center↗

Search for projectile fragments with fractional charge in relativistic heavy ion collisions

We measured the charge of about 35000 projectile fragments with Z > or = 5e produced by 14.5 GeV/nucleon and 200 GeV/nucleon 16O beams in a Pb target using CR39 plastic nuclear track detectors. A minimum track length of 3 mm in the detector without nuclear interaction was required. No evidence for fragments carrying a fractional charge was found.

NASA Discipline Number 22-70↗

Viewpoint 9--molecular structure of aqueous interfaces

In this review we summarize recent progress in our understanding of the structure of aqueous interfaces emerging from molecular level computer simulations. It is emphasized that the presence of the interface induces specific structural effects which, in turn, influence a wide variety of phenomena occurring near the phase boundaries. At the liquid-vapor interface, the most probable orientations of a water molecule is such that its dipole moment lies parallel to the interface, one O-H bond points toward the vapor and the other O-H bond is directed toward the liquid. The orientational distributions are broad and slightly asymmetric, resulting in an excess dipole moment pointing toward the liquid. These structural preferences persist at interfaces between water and nonpolar liquids, indicating that the interactions between the two liquids in contact are weak. It was found that liquid-liquid interfaces are locally sharp but broadened by capillary waves. One consequence of anisotropic orientations of interfacial water molecules is asymmetric interactions, with respect to the sign of the charge, of ions with the water surface. It was found that even very close to the surface ions retain their hydration shells. New features of aqueous interfaces have been revealed in studies of water-membrane and water-monolayer systems. In particular, water molecules are strongly oriented by the polar head groups of the amphiphilic phase, and they penetrate the hydrophilic head-group region, but not the hydrophobic core. At infinite dilution near interfaces, amphiphilic molecules exhibit behavior different from that in the gas phase or in bulk water. This result sheds new light on the nature of hydrophobic effect in the interfacial regions. The presence of interfaces was also shown to affect both equilibrium and dynamic components of rates of chemical reactions. Applications of continuum models to interfacial problems have been, so far, unsuccessful. This, again, underscores the importance of molecular-level information about interfaces.

Review, Tutorial↗

Exclusive measurements of mean pion multiplicities in 4He-nucleus reactions from 200 to 800 MeV/nucleon

Mean multiplicities of pi+ and pi- in 4He collisions with C, Cu, and Pb at 200, 600, and 800 MeV/u, and with C and Pb at 400 MeV/u have been measured using the large solid angle detector Diogene. The independence of pion multiplicity on projectile incident energy, target mass and proton multiplicity is studied in comparison with intra-nuclear cascade predictions. The discrepancy between experimental results and theory is pointed out and discussed.

NASA Discipline Radiation Health↗

Sensitive detection of narrowband pulses

Highly monochromatic signals, such as TV carriers, can be detected sensitively by using a narrow filter (b < or = 1 Hz) and performing power accumulation on its output. If instead a low-duty-cycle pulsed signal of the same total energy is present, the sensitivity of a square law device, followed by a thresholding operation (to eliminate most samples containing only noise), followed by the algorithm to be described, is greater by about 7 dB in typical cases. This is particularly interesting to SETI because such a pulsed signal is exactly what is sent by a rotating beacon with a directional antenna. Such a pulsed signal is, therefore, a good candidate for an extraterrestrial beacon. Software for detecting this signal type is now ready for field testing with the NASA Multichannel Spectrum Analyzer (MCSA).

NASA Discipline Number 52-60↗

A kill curve for Phanerozoic marine species

A kill curve for Phanerozoic species is developed from an analysis of the stratigraphic ranges of 17,621 genera, as compiled by Sepkoski. The kill curve shows that a typical species' risk of extinction varies greatly, with most time intervals being characterized by very low risk. The mean extinction rate of 0.25/m.y. is thus a mixture of long periods of negligible extinction and occasional pulses of much higher rate. Because the kill curve is merely a description of the fossil record, it does not speak directly to the causes of extinction. The kill curve may be useful, however, to li inverted question markmit choices of extinction mechanisms.

NASA Discipline Number 52-40↗