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 1,477 records · Page 82

Which Combinations of Environmental Conditions and Microphysical Parameter Values Produce a Given Orographic Precipitation Distribution?

This study applies an idealized modeling framework, alongside a Bayesian Markov chain Monte Carlo (MCMC) algorithm, to explore which combinations of upstream environmental conditions and cloud microphysical parameter values can produce a particular precipitation distribution over an idealized two-dimensional, bell-shaped mountain. Simulations focus on orographic precipitation produced when an atmospheric river interacts with topography. MCMC-based analysis reveals that different combinations of parameter values produce a similar precipitation distribution, with the most influential parameters being relative humidity (RH), horizontal wind speed (U), surface potential temperature (θ sfc ), and the snow fall speed coefficient (A s ). RH, U, and A s exhibit interdependence: changes in one or more of these factors can be mitigated by compensating changes in the other(s) to produce similar orographic precipitation rates. The results also indicate that the parameter sensitivities and relationships can vary for spatial subregions and given different environmental conditions. In particular, high θ sfc values are more likely to produce the target precipitation rate and spatial distribution, and thus the ensemble of simulations shows a preference for liquid precipitation at the surface. The results presented here highlight the complexity of orographic precipitation controls, and have implications for flood and water management, observational efforts, and climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

All-optical nanoscopic spatial control of molecular reaction yields on nanoparticles

Molecular adsorbate reactions on nanoparticles play a fundamental role in areas such as nano-photocatalysis, atmospheric, and astrochemistry. They can be induced, enhanced, and controlled by field localization and enhancement on the nanoparticle surface. In particular, the ability to perform highly controlled near-field-mediated reactions is key to deepening our understanding of surface photoactivity on nanosystems. Here, using reaction nanoscopy, we experimentally demonstrate all-optical nanoscopic control of surface reaction yields by tailoring the near fields on nanoparticles with waveform-controlled linear and bicircular two-color laser pulses, respectively. We observe site-selective proton emission from the dissociative ionization of adsorbate molecules on Si O 2 nanoparticles as a function of the polarization and relative phase of the two-color pulses. The angularly resolved close-to-uniform mapping between the surface reaction yields and the measured ion momentum enables the observation and spatial control of molecular reactions on the nanoparticle surface with nanoscopic resolution. The experimental results are modeled and reproduced qualitatively by classical trajectory Monte Carlo simulations. Our work paves the way toward reliable all-optical control of photocatalytic chemical reactions on nanoscale surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulations of short-term exposure to NO 2 and PM 2.5 to inform capture efficiency standards

The California building code requires all new or renovated residential dwelling units to have kitchen exhaust ventilation to manage air pollutants and moisture generated during cooking. Current performance requirements are specified as a minimum airflow and maximum sound level. This report presents an analysis to support consideration of adding a capture efficiency requirement to the code. The analysis uses a physics-based simulation model to calculate air pollutant concentrations in homes that result from cooking and the inflow of outdoor air, mediated by loss and removal process including deposition, dwelling unit ventilation, and the use of a range hood with varying levels of capture efficiency. Calculated pollutant concentrations are compared to relevant, health-based guidelines. The analysis considers the highest 1-hour concentration of nitrogen dioxide (NO 2 ), which is a product of natural gas combustion and the highest 24-hour concentrations of fine particulate matter (PM 2.5 ) which is emitted in substantial quantities from frying, broiling and grilling among other cooking activities. For NO 2 , the analysis considers cooking of a dinner for 3-4 persons consisting of pasta, meat sauce, a par-boiled vegetable, and baked garlic bread. For PM 2.5 , the analysis considers a day in which breakfast emitted particles at the 80th percentile and lunch and dinner entailed particle emissions at the 50th percentile of cooking emission events reported in the literature. Model simulations were performed to determine the level of range hood capture efficiency that will allow these cooking scenarios to occur in the vast majority (>99%) of new homes being built in California while maintaining pollutant concentrations below the health-based guidelines, if the range hood is used throughout cooking. All homes were assumed to have dwelling unit ventilation at the rate required in the building code. Simulation model input parameters were specified using a Monte Carlo approach to represent a range of housing characteristics, outdoor conditions, and indoor pollutant dynamics. Simulation results suggest that requiring a minimum capture efficiency of at least 70% is needed to avoid unacceptably high NO 2 (1-h average concentration of 100 ppb or higher) and at least 60% to avoid unacceptably high PM 2.5 (24-h average of 25 μg/m³ or higher). These results were driven by multi-family homes, which have a smaller volume of air to dilute any pollutants not captured and removed at the cooking area, resulting in higher concentrations.

54 ENVIRONMENTAL SCIENCES↗

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Trivializing Flow in 2D-O(3) model

The two dimensional O(3) sigma model, just as quantum chromodynamics, is an asymptotically free theory with a mass gap. Therefore, it is an interesting and simple toy model to investigate algorithms for Markov Chain Monte Carlo simulations of quantum chromodynamics. In this talk, we discuss the construction of a trivializing map, a field transformation from a given theory to a trivial one, through a suitably chosen gradient flow. An analytic solution for the generating functional of this trivializing flow has been obtained by a perturbative expansion in the flow time. Utilizing this solution allows for new approaches to be considered when proposing updates for a Markov Chain Monte Carlo algorithm.

Chamness, Christopher↗

A Review of the Multiple-Readout Concept and Its Application in an Integrally Active Calorimeter

A comprehensive multi-jet physics program is anticipated for experiments at future colliders. Key physics processes necessitate detectors that can distinguish signals from W and Z bosons and the Higgs boson. Typical examples include channels with or pairs and processes involving new physics in those cases where neutral particles must be disentangled from charged ones due to the presence of W or Z bosons in their final states. Such a physics program demands calorimetric energy resolution at or beyond the limits of traditional calorimetric techniques. Multiple-readout calorimetry, which aims to reduce fluctuations in energy measurements of hadronic showers, is a promising approach. The first part of this article reviews dual- and triple-readout calorimetry within a mathematical framework describing the underlying compensating mechanism. The second part proposes a potential implementation using an integrally active and total absorption detector. This model serves as the basis for several Monte Carlo studies, illustrating how the response of a multiple-readout calorimeter depends on construction parameters. Among the layouts considered, one configuration operating in triple-readout mode shows the potential to achieve an energy resolution approaching .

47 OTHER INSTRUMENTATION↗

Dark Energy Survey Year 3 Results: Measuring the Survey Transfer Function with Balrog

Abstract We describe an updated calibration and diagnostic framework, Balrog , used to directly sample the selection and photometric biases of the Dark Energy Survey (DES) Year 3 (Y3) data set. We systematically inject onto the single-epoch images of a random 20% subset of the DES footprint an ensemble of nearly 30 million realistic galaxy models derived from DES Deep Field observations. These augmented images are analyzed in parallel with the original data to automatically inherit measurement systematics that are often too difficult to capture with generative models. The resulting object catalog is a Monte Carlo sampling of the DES transfer function and is used as a powerful diagnostic and calibration tool for a variety of DES Y3 science, particularly for the calibration of the photometric redshifts of distant “source” galaxies and magnification biases of nearer “lens” galaxies. The recovered Balrog injections are shown to closely match the photometric property distributions of the Y3 GOLD catalog, particularly in color, and capture the number density fluctuations from observing conditions of the real data within 1% for a typical galaxy sample. We find that Y3 colors are extremely well calibrated, typically within ∼1–8 mmag, but for a small subset of objects, we detect significant magnitude biases correlated with large overestimates of the injected object size due to proximity effects and blending. We discuss approaches to extend the current methodology to capture more aspects of the transfer function and reach full coverage of the survey footprint for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Nuclear electromagnetic cascades from alpha particles incident on an iron absorber

Mean values and fluctuations of the nuclear-electromagnetic cascade development resulting from alpha particles incident on an iron absorber was determined for the primary energy range 10 to 300 GeV/nucleon. This has been accomplished by using a three-dimensional Monte Carlo simulation of the cascade process. The model used was first adjusted until it gave predictions for a small ionization spectrometer which agreed with measurements obtained by exposing the same spectrometer to 10, 20.5, and 28 GeV/c protons at an accelerator. The calcuations were then modified to apply to high energy alpha particles. The results show that measurements of the energies of alpha particles incident on an iron spectrometer with six interaction lengths total depth can be made with accuracies ranging from approximately 5% at 300 GeV/nucleon to approximately 16% at 10 GeV/nucleon.

Jones, W. V.↗

Production of high-energy gamma rays by cosmic ray interactions in the atmosphere and lunar surface

The production of gamma rays above 10 MeV by the interaction of cosmic rays in the atmosphere and the lunar surface is simulated using Monte Carlo methods. The calculation incorporates a new model of high-energy nucleon-nucleus interactions based on empirical fits to inclusive cross sections for the production of pions and nucleons. The atmospheric gamma ray flux is calculated as a function of direction, energy, and atmospheric depth. These calculations are compared with observations from balloons and from the SAS 2 satellite. Estimates of the flux of earth albedo electrons produced by cosmic ray interactions are presented. The lunar gamma ray albedo is calculated and compared with an upper limit based on SAS 2 measurements.

Morris, D. J.↗

Impact of flying qualities on mission effectiveness for helicopter air combat, volume 1

A computer simulation to investigate the impact of flying qualities on mission effectiveness is described. The objective of the study was to relate the effects of flying qualities, such as precision of flight path control and pilot workload, to the ability of a single Scout helicopter, or helicopter team, to accomplish a specified anti-armor mission successfully. The model of the actual engagement is a Monte Carlo simulation that has the capability to assess the effects of helicopter characteristics, numbers, tactics and weaponization on the force's ability to accomplish a specific mission against a specified threat as a function of realistic tactical factors. A key feature of this program is a simulation of micro-terrain features and their effects on detection, exposure, and masking for nap-of-the-earth (NOE) flight.

Harris, T. M.↗

Direct simulation of rotational and vibrational nonequilibrium

The ways in which energy transfer is calculated in the Direct Simulation Monte Carlo method are presented. An energy exchange model that deals with translational and rotational modes is described. A model for simulating the transfer of energy between the translational and vibrational modes is presented as well.

Boyd, Iain D.↗

A correction algorithm for particle size distribution measurements made with the forward-scattering spectrometer probe

A correction algorithm for evaluating the particle size distribution measurements of atmospheric aerosols obtained with a forward-scattering spectrometer probe (FSSP) is examined. A model based on Poisson statistics is employed to calculate the average diameter and rms width of the particle size distribution. The dead time and coincidence errors in the measured number density are estimated. The model generated data are compared with a Monte Carlo simulation of the FSSP operation. It is observed that the correlation between the actual and measured size distribution is nonlinear. It is noted that the algorithm permits more accurate calculation of the average diameter and rms width of the distribution compared to uncorrected measured quantities.

Lock, James A.↗

The solar ionisation rate deduced from Ulysses measurements and its implications to interplanetary Lyman alpha-intensity

We have developed a computer code which can be used to study 3-dimensional and time-dependent effects of the solar cycle on the interplanetary (IP) hydrogen distribution. The code is based on the inverted Monte Carlo simulation. In this work we have modelled the temporal behaviour of the solar ionisation rate. We have assumed that during the most of the time of the solar cycle there is an anisotopic latitudinal structure but right at the solar maximum the anisotropy disappears. The effects of this behaviour will be discussed both in regard to the IP hydrogen distribution and IP Lyman a a-intensity.

Summanen, T.↗

Prototype Conflict Alerting Logic for Free Flight

This paper discusses the development of a prototype alerting system for a conceptual Free Flight environment. The concept assumes that datalink between aircraft is available and that conflicts are primarily resolved on the flight deck. Four alert stages are generated depending on the likelihood of a conflict. If the conflict is not resolved by the flight crews, Air Traffic Control is notified to take over separation authority. The alerting logic is based on probabilistic analysis through modeling of aircraft sensor and trajectory uncertainties. Monte Carlo simulations were used over a range of encounter situations to determine conflict probability. The four alert stages were then defined based on probability of conflict and on the number of avoidance maneuvers available to the flight crew. Preliminary results from numerical evaluations and from a piloted simulator study at NASA Ames Research Center are summarized.

Yang, Lee C.↗

Evaluation of the Effect of Water Vapor On the Performance of NASA's NMRO Catalysts for Carbon Monoxide Oxidation

The Noble Metal Reducible Oxide (NMRO) catalysts for the low temperature oxidation of carbon monoxide were developed by NASA for the reoxidation of carbon monoxide which forms by the dissociation of carbon dioxide during the operation of sealed carbon dioxide lasers. The NMRO catalyst, which consists of a noble metal in conjunction with a reducible metal oxide, was evaluated under conditions that will be encountered in a CO2 laser operation, namely temperatures in the range 298 to 373 K and no significant reaction gas components other than CO, CO2 and O2. The NMRO catalysts may have significant potential for spin-off applications such as the prevention of carbon monoxide build-up in closed spaces as in space vehicle cabins or submarines, and the elimination of the cold start-up problem of automobile exhaust catalysts. The most significant difference in the conditions of these possible future applications is the high moisture content of the gases to be processed. Lack of understanding of the effects of water vapor and high temperature on catalyst activity and operation for extended periods are currently the main stumbling blocks for the transfer of this NASA technology to be used for commercial purposes. In the original proposal the following objectives were stated: To obtain experimental data on the adsorption, desorption and reaction characteristics of CO and CO2 the catalysts under high moisture conditions; to collect evidence on the presence of carbonate and hydroxyl surface species and their involvement in the CO oxidation mechanism; and to model the reaction system using a Monte-Carlo simulation to gain insight on the various steps involved. After the work has commenced the NASA technical monitor Mr. David Scheyer informed us that there was increased interest in the possible use of the NMRO catalysts as automobile exhaust catalysts and therefore NASA wanted to know whether the catalysts can operate at high temperatures as well as with high moisture gases. At that meeting it was decided that investigation of the high temperature performance of the NMRO catalysts should be given priority and replace the Monte-Carlo simulation objective. As a result, the modified objectives of the investigation were taken as the investigation of the high-temperature activity of the NMRO catalysts, and the effect of water vapor on the performance of these catalysts.

Akyurtlu, Ates↗

Estimates of M2 Tidal Energy Dissipation from TOPEX/Poseidon Altimeter Data

Most of the tidal energy dissipation in the ocean occurs in shallow seas, as has long been recognized. However, recent work has suggested that a significant fraction of the dissipation, perhaps 1 TW or more, occurs in the deep ocean. This paper builds further evidence for that conclusion. More than 6 years of data from the TOPEX/Poseidon satellite altimeter are used to map the tidal dissipation rate throughout the world ocean. The dissipation rate is estimated as a balance between the rate of working by tidal forces and the energy flux divergence, computed using currents derived by least squares fitting of the altimeter data and the shallow water equations. Such calculations require dynamical assumptions, in particular about the nature of dissipation. To assess sensitivity of dissipation estimates to input assumptions, a large suite of tidal inversions based on a wide range of drag parameterizations and employing both real and synthetic altimeter data are compared. These experiments and Monte Carlo error fields from a generalized inverse model are used to establish error uncertainties for the dissipation estimates. Owing to the tight constraints on tidal elevation fields provided by the altimeter, area integrals of the energy balance are remarkably insensitive to required dynamical assumptions. Tidal energy dissipation is estimated for all major shallow seas (excluding individual polar seas) and compared with previous model and data-based estimates. Dissipation in the open ocean is significantly tnhanced around major bathymetric features, in a manner consistent with simple theories the generation of baroclinic tides.

Egbert, Gary D.↗

State Estimation of International Space Station Centrifuge Rotor With Incomplete Knowledge of Disturbance Inputs

This thesis develops a state estimation algorithm for the Centrifuge Rotor (CR) system where only relative measurements are available with limited knowledge of both rotor imbalance disturbances and International Space Station (ISS) thruster disturbances. A Kalman filter is applied to a plant model augmented with sinusoidal disturbance states used to model both the effect of the rotor imbalance and the 155 thrusters on the CR relative motion measurement. The sinusoidal disturbance states compensate for the lack of the availability of plant inputs for use in the Kalman filter. Testing confirms that complete disturbance modeling is necessary to ensure reliable estimation. Further testing goes on to show that increased estimator operational bandwidth can be achieved through the expansion of the disturbance model within the filter dynamics. In addition, Monte Carlo analysis shows the varying levels of robustness against defined plant/filter uncertainty variations.

SINUSOIDAL DISTURBANCE STATES↗