Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed simulation”

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 379 records · Page 21

Comparison of film measurements and Monte Carlo simulations of dose delivered with very high-energy electron beams in a polystyrene phantom

Purpose: To measure radiation dose in a water-equivalent medium from very high-energy electron (VHEE) beams and make comparisons to Monte Carlo (MC) simulation results. Methods: Dose in a polystyrene phantom delivered by an experimental VHEE beam line was measured with Gafchromic films for three 50 MeV and two 70 MeV Gaussian beams of 4.0–6.9 mm FWHM and compared to corresponding MC-simulated dose distributions. MC dose in the polystyrene phantom was calculated with the EGSnrc/BEAMnrc and DOSXYZnrc codes based on the experimental setup. Additionally, the effect of 2% beam energy measurement uncertainty and possible non-zero beam angular spread on MC dose distributions was evaluated. Results: MC simulated percentage depth dose (PDD) curves agreed with measurements within 4% for all beam sizes at both 50 and 70 MeV VHEE beams. Central axis PDD at 8 cm depth ranged from 14% to 19% for the 5.4–6.9 mm 50 MeV beams and it ranged from 14% to 18% for the 4.0–4.5 mm 70 MeV beams. MC simulated relative beam profiles of regularly shaped Gaussian beams evaluated at depths of 0.64 to 7.46 cm agreed with measurements to within 5%. A 2% beam energy uncertainty and 0.286° beam angular spread corresponded to a maximum 3.0% and 3.8% difference in depth dose curves of the 50 and 70 MeV electron beams, respectively. Absolute dose differences between MC simulations and film measurements of regularly shaped Gaussian beams were between 10% and 42%. Conclusions: The authors demonstrate that relative dose distributions for VHEE beams of 50–70 MeV can be measured with Gafchromic films and modeled with Monte Carlo simulations to an accuracy of 5%. The reported absolute dose differences likely caused by imperfect beam steering and subsequent charge loss revealed the importance of accurate VHEE beam control and diagnostics.

07 ISOTOPE AND RADIATION SOURCES↗

System Engineering Strategy for Distributed Multi-Purpose Simulation Architectures

This paper describes the system engineering approach used to develop distributed multi-purpose simulations. The multi-purpose simulation architecture focuses on user needs, operations, flexibility, cost and maintenance. This approach was used to develop an International Space Station (ISS) simulator, which is called the International Space Station Integrated Simulation (ISIS)1. The ISIS runs unmodified ISS flight software, system models, and the astronaut command and control interface in an open system design that allows for rapid integration of multiple ISS models. The initial intent of ISIS was to provide a distributed system that allows access to ISS flight software and models for the creation, test, and validation of crew and ground controller procedures. This capability reduces the cost and scheduling issues associated with utilizing standalone simulators in fixed locations, and facilitates discovering unknowns and errors earlier in the development lifecycle. Since its inception, the flexible architecture of the ISIS has allowed its purpose to evolve to include ground operator system and display training, flight software modification testing, and as a realistic test bed for Exploration automation technology research and development.

Bhula, Dlilpkumar↗

Benchmarking Model Variants in Development of a Hardware-in-the-Loop Simulation System

Distributed engine control architecture presents a significant increase in complexity over traditional implementations when viewed from the perspective of system simulation and hardware design and test. Even if the overall function of the control scheme remains the same, the hardware implementation can have a significant effect on the overall system performance due to differences in the creation and flow of data between control elements. A Hardware-in-the-Loop (HIL) simulation system is under development at NASA Glenn Research Center that enables the exploration of these hardware dependent issues. The system is based on, but not limited to, the Commercial Modular Aero-Propulsion System Simulation 40k (C-MAPSS40k). This paper describes the step-by-step conversion from the self-contained baseline model to the hardware in the loop model, and the validation of each step. As the control model hardware fidelity was improved during HIL system development, benchmarking simulations were performed to verify that engine system performance characteristics remained the same. The results demonstrate the goal of the effort; the new HIL configurations have similar functionality and performance compared to the baseline C-MAPSS40k system.

Networked Systems↗

Transport processes in the stratosphere: Model simulations and comparisons with satellite observations

A three dimensional atmospheric model was used to study transport processes and to simulate the distribution of chemically active species in the stratosphere. The present results are part of a long term simulation of the seasonally varying structure of stratospheric trace constituents. The occurrence of a midwinter stratospheric warming and the associated transport of O3 and HNO3 during the simulation are described. Comparison of the simulated distributions of O3 and HNO3 are made with data from the Limb Infrared Monitor of the Stratosphere (LIMS) experiment. In addition, distributions of Ertel's potential vorticity on isentropic surfaces (IPV) are evaluated as a diagnostic for interpreting transport processes. Comparisons are made with IPV distributions inferred from LIMS temperature data.

Grose, W. L.↗

Benchmarking Model Variants in Development of a Hardware-in-the-Loop Simulation System

Distributed engine control architecture presents a significant increase in complexity over traditional implementations when viewed from the perspective of system simulation and hardware design and test. Even if the overall function of the control scheme remains the same, the hardware implementation can have a significant effect on the overall system performance due to differences in the creation and flow of data between control elements. A Hardware-in-the-Loop (HIL) simulation system is under development at NASA Glenn Research Center that enables the exploration of these hardware dependent issues. The system is based on, but not limited to, the Commercial Modular Aero-Propulsion System Simulation 40k (C-MAPSS40k). This paper describes the step-by-step conversion from the self-contained baseline model to the hardware in the loop model, and the validation of each step. As the control model hardware fidelity was improved during HIL system development, benchmarking simulations were performed to verify that engine system performance characteristics remained the same. The results demonstrate the goal of the effort; the new HIL configurations have similar functionality and performance compared to the baseline C-MAPSS40k system.

Propulsion Simulation↗

Investigating Characteristic Droplet Size Distributions in Large Eddy Simulations of Stratocumulus Clouds

Cloud processes relevant to radiative and precipitation properties depend on the shape of the cloud droplet size distribution. Recent holographic observations revealed that cloud droplet populations do not have the same size distribution shapes throughout but form regions of characteristic distributions with similar microphysical properties. We investigate the existence and properties of these characteristic distributions within Large‐Eddy Simulations of stratocumulus clouds using Lagrangian and bin microphysics schemes. Distribution types are identified, revealing localized characteristic distributions that vary on the scale of the largest convective cell for simulations with bin microphysics. The results from the Lagrangian microphysics scheme hint at similar behavior. Compared to observations, the simulated clouds are much more uniform. Analysis of the LES results suggests a connection to the local entrainment rate, so the poorly resolved entrainment interface in LES may be a cause of the uniformity. The uniformity of the large‐scale forcing could also be a factor.

cloud droplet size distributions↗

ASSESSING ASTROPHYSICAL UNCERTAINTIES IN DIRECT DETECTION WITH GALAXY SIMULATIONS

We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions that are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Furthermore, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulations. We show that the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way’s estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models that are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulated Craters on Venus

The thick atmosphere of Venus prevents all but the largest impactors from cratering the surface. The number of small craters on Venus provides an interesting, and statistically significant test of models for the disruption and deceleration of impacting bodies. Here we compare Monte Carlo simulated crater distributions to the observed crater distribution on Venus. The simulation assumes: (1) a power law mass distribution for impactors of the form N(sub cum) alpha m (exp-b) where b=0.8; (2) isotropic incidence angles; (3) velocity at the top of the atmosphere of 20 kilometers per second (more realistic velocity distributions are also considered); (4) Schmidt-Housen crater scaling, modified such that only the normal component of the impact velocity contributes to cratering, and using crater slumping as parameterized (5) and modern populations (60% carbonaceous, 40% stone, 3% iron) and fluxes of asteroids. We use our previously developed model for the disruption and deceleration of large bodies striking thick planetary atmospheres to calculate the impact velocity at the surface as a function of impactor mass, incident velocity, and incident angle. We use a drag coefficient c(sub d) =1; other parameters are as described in Chyba et al. We set a low velocity cutoff of 500 meters per second on crater-forming impacts. Venus's craters are nicely matched by the simulated craters produced by 700 million years of striking asteroids. Shown for comparison are the simulated craters produced by incident comets over the same period, where for comets we have assumed b=0.7 and a flux at 10(exp 14) g 30% that of asteroids. Systematic uncertainties in crater scaling and crater slumping may make the surface age uncertain by a factor of two.

Zahnle, Kevin↗

Penumbra: A spatially distributed, mechanistic model for simulating ground-level incident solar energy across heterogeneous landscapes

Landscape solar energy is a significant environmental driver, yet it remains complicated to model well. Several solar radiation models simplify the complexity of light by estimating it at discrete point locations or by averaging values over larger areas. These modeling approaches may be useful in certain cases, but they are unable to provide spatially distributed and temporally dynamic representations of solar energy across entire landscapes. We created a landscape-scale ground-level shade and solar energy model called Penumbra to address this deficiency. Penumbra simulates spatially distributed ground-level shade and incident solar energy at user-defined timescales by modeling local and distant topographic shading and vegetative shading. Spatially resolved inputs of a digital elevation model, a normalized digital surface model, and landscape object transmittance are used to estimate spatial variations in solar energy at user-defined temporal timesteps. The research goals for Penumbra included: 1) simulations of spatiotemporal variations of shade and solar energy caused by both objects and topographic features, 2) minimal user burden and parameterization, 3) flexible user defined temporal parameters, and 4) flexible external model coupling. We test Penumbra’s predictive skill by comparing the model’s predictions with monitored open and forested sites, and achieve calibrated mean errors ranging from -17.3 to 148.1 μmoles/m 2 /s. Penumbra is a dynamic model that can produce spatial and temporal representations of shade percentage and ground-level solar energy. Outputs from Penumbra can be used with other ecological models to better understand the health and resilience of aquatic, near stream terrestrial, and upland ecosystems.

54 ENVIRONMENTAL SCIENCES↗

Transonic pressure and load distributions for a group of simulated launch vehicles

Pressure and load distributions for a related group of simulated launch vehicle configurations are presented. The configurations were selected so that the nose cone and interstage transition flare components were relatively close to one another and subject to mutual interference effects. Tests extended over a Mach number range from 0.40 to 1.20 at angles of attack from 0 deg to about 10 deg. The test Reynolds numbers, based on main stage diameter, were of the order of 0.00000098.

Kelly, T. C.↗

Caustic-induced features in microlensing magnification probability distributions

Numerical simulations have uncovered a previously unrecognized 'bump' in the macroimage magnification probabilities produced by a planar distribution of point masses. The result could be relevant to cases of microlensing by star fields in single galaxies, for which this lensing geometry is an excellent approximation. The bump is produced by bright pairs of microimages formed by sources lying near the caustics of the lens. The numerically calculated probabilities for the magnifications in the range between 3 and 30 are significantly higher than those given by the asymptotic relation derived by Schneider. The bump present in the two-dimensional lenses appears not to exist in the magnification probability distribution produced by a fully three-dimensional lens.

Rauch, Kevin P.↗

Margin Evaluation in the Presence of Deformation, Rotation, and Translation in Prostate and Entire Seminal Vesicle Irradiation With Daily Marker-Based Setup Corrections

Purpose: To develop a method for margin evaluation accounting for all measured displacements during treatment of prostate cancer. Methods and Materials: For 21 patients treated with stereographic targeting marker-based online translation corrections, dose distributions with varying margins and gradients were created. Sets of possible cumulative delivered dose distributions were simulated by moving voxels and accumulating dose per voxel. Voxel motion was simulated consistent with measured distributions of systematic and random displacements due to stereographic targeting inaccuracies, deformation, rotation, and intrafraction motion. The method of simulation maintained measured correlation of voxel motions due to organ deformation. Results: For the clinical target volume including prostate and seminal vesicles (SV), the probability that some part receives <95% of the prescribed dose, the changes in minimum dose, and volume receiving 95% of prescription dose compared with planning were 80.5% {+-} 19.2%, 9.0 {+-} 6.8 Gy, and 3.0% {+-} 3.7%, respectively, for the smallest studied margins (3 mm prostate, 5 mm SV) and steepest dose gradients. Corresponding values for largest margins (5 mm prostate, 8 mm SV) with a clinical intensity-modulated radiotherapy dose distribution were 46.5% {+-} 34.7%, 6.7 {+-} 5.8 Gy, and 1.6% {+-} 2.3%. For prostate-only clinical target volume, the values were 51.8% {+-} 17.7%, 3.3 {+-} 1.6 Gy, and 0.6% {+-} 0.5% with the smallest margins and 5.2% {+-} 7.4%, 1.8 {+-} 0.9 Gy, and 0.1% {+-} 0.1% for the largest margins. Addition of three-dimensional rotation corrections only improved these values slightly. All rectal planning constraints were met in the actual reconstructed doses for all studied margins. Conclusion: We developed a system for margin validation in the presence of deformations. In our population, a 5-mm margin provided sufficient dosimetric coverage for the prostate. In contrast, an 8-mm SV margin was still insufficient owing to deformations. Addition of three-dimensional rotation corrections was of minor influence.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Study of In-Trap Ion Clouds by Ion Trajectory Simulations

Gaussian distribution has been utilized to describe the global number density distribution of ion cloud in the Paul trap, which is known as the thermal equilibrium theory and widely used in theoretical modeling of ion clouds in the ion traps. Using ion trajectory simulations, however, the ion clouds can now also be treated as a dynamic ion flow field and the location-dependent features could now be characterized. This study was carried out to better understand the in-trap ion cloud properties, such as the local particle velocity and temperature. The local ion number densities were found to be heterogeneously distributed in terms of mean and distribution width; the velocity and temperature of the ion flow varied with pressure depending on the flow type of the neutral molecules; and the “quasi-static” equilibrium status can only be achieved after a certain number of collisions, for which the time period is pressure-dependent. This work provides new insights of the ion clouds that are globally stable but subjected to local rf heating and collisional cooling. .

74 ATOMIC AND MOLECULAR PHYSICS↗

Impact of DER Communication Delay in AGC: Cyber-Physical Dynamic Co-Simulation: Preprint

Distributed energy resource (DER) frequency regulations are promising technologies for future grid operation. Unlike conventional generators, DERs might require open communication networks to exchange signals with control centers, possibly through DER aggregators; therefore, the impacts of the communication variations including the latency on the system frequency stability need to be investigated. This paper develops a cyber-physical dynamic simulation (CPDS) model based on the Hierarchical Engine for Large-Scale Co-Simulation (HELICS) to evaluate the impact of the communication variations, such as delays in DER frequency regulations. The feasible delay range can be obtained under different AGC parameter settings. The results show that the risk of instability generally increases with the communication delay and that the difference of AGC provided from DER and conventional generation units should be recognized when designing frequency regulation controllers.

61 RADIATION PROTECTION AND DOSIMETRY↗

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 2: Stratospheric and free-tropospheric response

The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Simulation to Study Speed Distributions in a Solar Plasma

We have carried out a numerical simulation of a plasma with characteristics similar to those found in the core of the Sun. Particular emphasis is placed on the Coulomb interaction between the ions and electrons, which could result in a relative velocity distribution different from the Maxwell-Boltzmann (MB) distribution generally assumed for a plasma. The fact that the distribution may not exactly follow the MB distribution could have very important consequences for a variety of problems in solar physics, especially the neutrino problem. Very briefly, the neutrino problem is that the observed neutrino detections from the Sun are smaller than what the standard solar theory predicts. In Section I we introduce the problem and in section II we discuss the approach to try to solve the problem: i.e., a molecular dynamics approach. In section III we provide details about the integration method, and any simplifications that can be applied to the problem. In section IV (the core of this report) we state our results. First for the specific case of 1000 particles and then for other cases with different number of particles. In section V we summarize our findings and state our conclusions. Sections VI VII and VIII provide the list of figures, reference material and acknowledgements respectively.

Cheeseman, Peter↗