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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

Electron Anisotropies in Magnetotail Dipolarization Events

Anisotropic electron distributions can be a major source of free energy for the generation of microinstabilities and waves. Here we investigate specifically anisotropies of the suprathermal electrons associated with dipolarization events in the geomagnetic tail. The investigation is based on an MHD simulation of magnetotail reconnection, flow bursts and dipolarization. This simulation is used to trace test particles in the electromagnetic fields of the MHD simulation and investigate their acceleration and contributions to flux increases. The simulated velocity distributions yield anisotropies which can be dominantly parallel (“cigar”-shaped), perpendicular (“pancake”-shaped) to the magnetic field, or a combination thereof (e.g., “rolling-pin”-shaped), depending on location, relative timing, and energy.

electron anisotropies↗
DOE OSTI2011

Effect of Particle Size Distribution on Slurry Rheology: Nuclear Waste Simulant Slurries

Controlling the rheological properties of slurries has been of great interest in various industries such as cosmetics, ceramic processing, and nuclear waste treatment. Many physicochemical parameters, such as particle size, pH, ionic strength, and mass/volume fraction of particles, can influence the rheological properties of slurry. Among such parameters, the particle size distribution of slurry would be especially important for nuclear waste treatment because most nuclear waste slurries show a broad particle size distribution. We studied the rheological properties of several different low activity waste nuclear simulant slurries having different particle size distributions under high salt and high pH conditions. Using rheological and particle size analysis, it was found that the percentage of colloid-sized particles in slurry appears to be a key factor for rheological characteristics and the efficiency of rheological modifiers. This behavior was shown to be coupled with an existing electrostatic interaction between particles under a low salt concentration. Our study suggests that one may need to implement the particle size distribution as a critical factor to understand and control rheological properties in nuclear waste treatment plants, such as the U.S. Department of Energy’s Hanford and Savannah River sites, because the particle size distributions significantly vary over different types of nuclear waste slurries.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

ANALYZING METEOROID FLIGHTS USING PARTICLE FILTERS

Fireball observations from camera networks provide position and time information along the trajectory of a meteoroid that is transiting our atmosphere. The complete dynamical state of the meteoroid at each measured time can be estimated using Bayesian filtering techniques. A particle filter is a novel approach to modeling the uncertainty in meteoroid trajectories and incorporates errors in initial parameters, the dynamical model used, and observed position measurements. Unlike other stochastic approaches, a particle filter does not require predefined values for initial conditions or unobservable trajectory parameters. The Bunburra Rockhole fireball, observed by the Australian Desert Fireball Network (DFN) in 2007, is used to determine the effectiveness of a particle filter for use in fireball trajectory modeling. The final mass is determined to be 2.16±1.33 kg with a final velocity of 6030±216 m s{sup −1}, similar to previously calculated values. The full automatability of this approach will allow an unbiased evaluation of all events observed by the DFN and lead to a better understanding of the dynamical state and size frequency distribution of asteroid and cometary debris in the inner solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

SU-E-I-28: Evaluating the Organ Dose From Computed Tomography Using Monte Carlo Calculations

Purpose: To evaluate organ doses from computed tomography (CT) using Monte Carlo (MC) calculations. Methods: A Philips Brilliance CT scanner (64 slice) was simulated using the GMctdospp (IMPS, Germany) based on the EGSnrc user code. The X-ray spectra and a bowtie filter for MC simulations were determined to coincide with measurements of half-value layer (HVL) and off-center ratio (OCR) profile in air. The MC dose was calibrated from absorbed dose measurements using a Farmer chamber and a cylindrical water phantom. The dose distribution from CT was calculated using patient CT images and organ doses were evaluated from dose volume histograms. Results: The HVLs of Al at 80, 100, and 120 kV were 6.3, 7.7, and 8.7 mm, respectively. The calculated HVLs agreed with measurements within 0.3%. The calculated and measured OCR profiles agreed within 3%. For adult head scans (CTDIvol) =51.4 mGy), mean doses for brain stem, eye, and eye lens were 23.2, 34.2, and 37.6 mGy, respectively. For pediatric head scans (CTDIvol =35.6 mGy), mean doses for brain stem, eye, and eye lens were 19.3, 24.5, and 26.8 mGy, respectively. For adult chest scans (CTDIvol=19.0 mGy), mean doses for lung, heart, and spinal cord were 21.1, 22.0, and 15.5 mGy, respectively. For adult abdominal scans (CTDIvol=14.4 mGy), the mean doses for kidney, liver, pancreas, spleen, and spinal cord were 17.4, 16.5, 16.8, 16.8, and 13.1 mGy, respectively. For pediatric abdominal scans (CTDIvol=6.76 mGy), mean doses for kidney, liver, pancreas, spleen, and spinal cord were 8.24, 8.90, 8.17, 8.31, and 6.73 mGy, respectively. In head scan, organ doses were considerably different from CTDIvol values. Conclusion: MC dose distributions calculated by using patient CT images are useful to evaluate organ doses absorbed to individual patients.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2003

Radar Reflectivity Simulated by a 2-D Spectra Bin Model: Sensitivity of Cloud-aerosol Interaction

The Goddard Cumulus Ensemble (GCE) model with bin spectra microphysics is used to simulate mesoscale convective systems.The model uses explicit bins to represent size spectra of cloud nuclei, water drops, ice crystals, snow and graupel. Each hydrometeorite category is described by 33 mass bins. The simulations provide a unique data set of simulated raindrop size distribution in a realistic dynamic frame. Calculations of radar parameters using simulated drop size distribution serve as an evaluation of numerical model performance. In addition, the GCE bin spectra modes is a very useful tool to study uncertainties related to radar observations; all the environmental parameters are precisely known. In this presentation, we concentrate on the discussion of Z-R (ZDR-R) relation in the simulated systems. Due to computational limitations, the spectra bin model has been run in two dimensions with 31 stretched vertical layers and 1026 horizontal grid points (1 km resolution). Two different cases, one in midlatitude continent, the other in tropical ocean, have been simulated. The continental case is a strong convection which lasted for two hours. The oceanic case is a persistent system with more than 10 hours' life span. It is shown that the simulated Z-R (ZDR-R) relations generally agree with observations using radar and rain gauge data. The spatial and temporal variations of Z-R relation in different locations are also analyzed. Impact of aerosols on cloud formation and raindrop size distribution was studied. Both clean (low CCN) and dirty (high CCN) cases are simulated. The Z-R relation is shown to vary considerable in the initial CCN concentrations.

Li, Kiaowen↗
DOE OSTI2019

X-ray Scattering and Coarse-Grained Simulations for Clustering and Interactions of Monoclonal Antibodies at High Concentrations

Attractive protein-protein interactions (PPI) in concentrated monoclonal antibody (mAb) solutions may lead to reversible oligomers (clusters) that impact colloidal stability and viscosity. Herein, the PPI are tuned for two mAbs via the addition of arginine (Arg), NaCl, or ZnSO 4 as characterized by the structure factor ( S eff (q)) with small angle x-ray scattering (SAXS). The SAXS data are fit with molecular dynamics simulations by placing a physically relevant short-ranged attractive interaction on selected beads in coarse-grained 12 bead models of the mAb shape. The optimized 12-bead models are then used to differentiate key microstructural properties, including center of mass radial distribution functions ( g COM (r)), coordination numbers and cluster size distributions (CSD). The addition of co-solutes results in more attractive Seff(q) relative to the no co-solute control for all systems tested, with the most attractive systems showing an upturn at low q . Only the All1 model with an attractive site in each Fab and Fc region (possessing Fab-Fab, Fab-Fc and Fc-Fc interactions) can reproduce this upturn, and the corresponding CSDs show the presence of larger clusters compared to the control. In general, for models with similar net attractions, i.e., second osmotic virial coefficients, the size of the clusters increases as the attraction is concentrated on a smaller number of evenly distributed beads. Here, the cluster size distributions from simulations are used to improve the understanding and prediction of experimental viscosities. The ability to discriminate between models with bead interactions at particular Fab and Fc bead sites from SAXS simulations, and to provide real space properties (CSD and g COM (r)), will be of interest in engineering protein sequence and formulating protein solutions for weak PPI to minimize aggregation and viscosities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Decomposition of multi-particle azimuthal correlations in Q-cumulant analysis*

Abstract The method of Q -cumulants is a powerful tool for studying the fine details of azimuthal anisotropies in high energy nuclear collisions. This paper presents a new method, based on mathematical induction, to evaluate the analytical form of high-order Q -cumulants. The capability of this method is demonstrated via a toy model that uses the elliptic power distribution to simulate the anisotropic emission of particles, quantified in terms of Fourier flow harmonics . The method can help in studying the large amount of event statistics that can be collected in the future and allow measurements of the very high central moments of the distribution. This can, in turn, facilitate progress in understanding the initial geometry, the input to the hydrodynamic calculations of medium expansion in high energy nuclear collisions, and the constraints on it.

Physics↗
DOE OSTI2016

Finite-element analysis of transverse compressive and thermal loads on Nb 3 Sn wires with voids

High-field superconducting magnets play a very important role in many large-scale physics experiments, particularly particle colliders and fusion confinement devices such as Large Hadron Collider (LHC) and International Thermonuclear Experimental Reactor (ITER). The two most common superconductors used in these applications are NbTi and Nb 3 Sn. Nb 3 Sn wires are favored because of their significantly higher J c (critical current density) for higher field applications. The main disadvantage of Nb 3 Sn is that the superconducting performance of the wire is highly strain sensitive and it is very brittle. This strain sensitivity is strongly influenced by two factors: plasticity and cracked filaments. Cracks are induced by large stress concentrators that can be traced to the presence of voids in the wire. We develop detailed 2-D and 3-D finite-element models containing wire filaments and different possible distributions of voids in a bronze-route Nb 3 Sn wire. We apply compressive transverse loads for various cases of void distributions to simulate the stress and strain response of a Nb 3 Sn wire under the Lorentz force. Furthermore, this paper improves our understanding of the effect voids have on the Nb 3 Sn wire's mechanical properties, and in so, the connection between the distribution of voids and performance degradation such as the correlation between irreversible strain limit and the void-induced local stress concentrations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Design outline for a new multiman ATC simulation facility at NASA-Ames Research Center

A new and unique facility for studying human factors aspects in aeronautics is being planned for use in the Man-Vehicle Systems Research Division at the NASA-Ames Research Center. This facility will replace the existing three cockpit-single ground controller station and be expandable to include approximately seven cockpits and two ground controller stations. Unlike the previous system, each cockpit will be mini-computer centered and linked to a main CPU to effect a distributed computation facility. Each simulator will compute its own flight dynamic and flight path predictor. Mechanical flight instruments in each cockpit will be locally supported and CRT cockpit displays of (e.g.) traffic and or RNAV information will be centrally computed and distributed as a means of extending the existing computational and graphical resources. An outline of the total design is presented which addresses the technical design options and research possibilities of this unique man-machine facility and which may also serve as a model for other real time distributed simulation facilities.

Kreifeldt, J. G.↗
DOE OSTI2016

Simulation of local ion transport in lamellar block copolymer electrolytes based on electron micrographs

A method is presented in this paper to relate local morphology and ionic conductivity in a solid, lamellar block copolymer electrolyte for lithium batteries, by simulating conductivity through transmission electron micrographs. The electrolyte consists of polystyrene-block-poly(ethylene oxide) mixed with lithium bis(trifluoromethanesulfonyl) imide salt (SEO/LiTFSI), where the polystyrene phase is structural phase and the poly(ethylene oxide)/LiTFSI phase is ionically conductive. The electric potential distribution is simulated in binarized micrographs by solving the Laplace equation with constant potential boundary conditions. A morphology factor, f, is reported for each image by calculating the effective conductivity relative to a homogenous conductor. Images from two samples are examined, one annealed with large lamellar grains and one unannealed with small grains. The average value off is 0.45 ± 0.04 for the annealed sample, and 0.37 ± 0.03 for the unannealed sample, both close to the value predicted by effective medium theory, 1/2. Simulated conductivities are compared to published experimental conductivities. The value of f Unannealed /f Annealed is 0.82 for simulations and 6.2 for experiments. Simulation results correspond well to predictions by effective medium theory but do not explain the experimental measurements. Finally, observation of nanoscale morphology over length scales greater than the size of the micrographs (~1 μm) may be required to explain the experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2016

SU-F-R-36: Validating Quantitative Radiomic Texture Features for Oncologic PET: A Digital Phantom Study

Purpose: Radiomic texture features derived from the oncologic PET have recently been brought under intense investigation within the context of patient stratification and treatment outcome prediction in a variety of cancer types; however, their validity has not yet been examined. This work is aimed to validate radiomic PET texture metrics through the use of realistic simulations in the ground truth setting. Methods: Simulation of FDG-PET was conducted by applying the Zubal phantom as an attenuation map to the SimSET software package that employs Monte Carlo techniques to model the physical process of emission imaging. A total of 15 irregularly-shaped lesions featuring heterogeneous activity distribution were simulated. For each simulated lesion, 28 texture features in relation to the intensity histograms (GLIH), grey-level co-occurrence matrices (GLCOM), neighborhood difference matrices (GLNDM), and zone size matrices (GLZSM) were evaluated and compared with their respective values extracted from the ground truth activity map. Results: In reference to the values from the ground truth images, texture parameters appearing on the simulated data varied with a range of 0.73–3026.2% for GLIH-based, 0.02–100.1% for GLCOM-based, 1.11–173.8% for GLNDM-based, and 0.35–66.3% for GLZSM-based. For majority of the examined texture metrics (16/28), their values on the simulated data differed significantly from those from the ground truth images (P-value ranges from <0.0001 to 0.04). Features not exhibiting significant difference comprised of GLIH-based standard deviation, GLCO-based energy and entropy, GLND-based coarseness and contrast, and GLZS-based low gray-level zone emphasis, high gray-level zone emphasis, short zone low gray-level emphasis, long zone low gray-level emphasis, long zone high gray-level emphasis, and zone size nonuniformity. Conclusion: The extent to which PET imaging disturbs texture appearance is feature-dependent and could be substantial. It is thus advised that use of PET texture parameters for predictive and prognostic measurements in oncologic setting awaits further systematic and critical evaluation.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2000

A Simulation of Bromoform's Contribution to Stratospheric Bromine

Many chlorinated and brominated compounds that are inert in the troposphere are destroyed in the stratosphere and act as an in-situ source of stratospheric reactive chlorine and bromine. Other halogenated compounds that are reactive in the troposphere might contribute to the stratosphere's halogen budget in two ways. First, like their unreactive companions, rapid convective transport might carry them to the upper troposphere and make them available for subsequent advection by the mean circulation into the stratosphere before they are oxidized or photolyzed. Second, it is more likely that they are destroyed in the troposphere, and the chlorine and bromine that is released might then be transported to the stratosphere. We evaluate the relative influence of these processes on stratospheric bromine in a three-dimensional chemistry and transport model which simulates the distribution of bromoform (CHBr3). CHBr3 is parameterized as a short-lived, ocean-surface source gas whose destruction by photolysis and reaction with hydroxyl (OH) in the troposphere and stratosphere yields inorganic bromine (Br(sub y)). Many of the observed features of CHBr3 are simulated well, and comparisons with observations are used to show that the model represents aspects of transport in the upper troposphere and lower stratosphere that are critical to the evaluation. In particular, the model maintains the observed troposphere-stratosphere distinctness in transport pathways and reproduces the observed seasonal dependence of the mixture of air in the middle- and high-latitude lowermost stratosphere. We estimate that adding CHBr3 to models which already include the long-lived organic brominated compounds (halons and methyl bromide) will increase the simulated stratospheric mass of Br(sub y) by about 15 percent. In-situ stratospheric destruction of CHBr3 produces Br(sub y) in amounts which are comparable to that transported into the stratosphere after photolysis and oxidation of CHBr3 in the troposphere. In our simulations the mass of Br(sub y) produced from the destruction of CHBr3 does not exceed the mass of Br(sub y) produced from the destruction of the long-lived compounds at any level in the stratosphere. However, Br(sub y) from the loss of CHBr3 accounts for approximately one-third of the total Br(sub y) in the lowest kilometer of the stratosphere.

Nielsen, J. Eric↗
DOE OSTI2013

Evaluation of the use of color-set geometry during lattice physics constants generation for boiling water reactor simulation

Current methods for BWR nuclear design and analysis consist of using lattice physics neutron transport methods to generate the two-group homogenized cross-sections that are then used in a nodal diffusion theory code. The lattice transport solutions are performed for a single assembly with reflective boundary conditions, which is a practical approximation. A method is developed to account for assembly exposure distributions (environment) in the core within the lattice transport calculations with the use of color-sets (2x2) geometry. The loading pattern is examined and an appropriate number of characteristic color-set cells are selected for analysis. Treatment of the co-resident exposed fuel within this method is also presented. The calculation process was followed for a recent BWR cycle design with comparisons being performed on both a lattice and core-wide basis to evaluate the proposed method. The lattice based comparisons show noticeable differences in the pin power distribution predictions, which require further investigation to see how this translates into core performance calculations. The core-wide comparisons show minor differences and are generally in a good agreement, which is expected with this small perturbation. A slight improvement was noticed in the reduction of the power distribution uncertainty. However, given the additional amount of work and computer run time increase, further evaluation, especially of core pin power predictions, is needed to consider this method for production level design and safety analysis calculations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2011

Modeling Impacts of Climate Change Mitigation Technologies on Power Grids

This paper describes a modeling approach that simulate the impacts of different climate change mitigation technologies on power grids for power system planning purposes. Because the historical data is less credible when new technologies are being deployed to the system, it is then critical to model them to address their impacts. This paper illustrated how to integrate modeling results obtained from different modeling tools to give a reasonable forecast of the future. Building simulation tools, distribution power grid modeling tools, and power system planning tools are used to model and aggregate the impacts from the end-use to the system level. Electricity generation, production cost, emission, and transmission congestions are used to quantify the influence of different mitigation technologies. Modeling results have shown that the cross-discipline modeling approach provided the modeler with the necessary time resolution and input details to address the variables that influence the modeling results. Different modeling issues are also addressed in the paper.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

Propulsion Electric Grid Simulator (PEGS) for Future Turboelectric Distributed Propulsion Aircraft

NASA Glenn Research Center, in collaboration with the aerospace industry and academia, has begun the development of technology for a future hybrid-wing body electric airplane with a turboelectric distributed propulsion (TeDP) system. It is essential to design a subscale system to emulate the TeDP power grid, which would enable rapid analysis and demonstration of the proof-of-concept of the TeDP electrical system. This paper describes how small electrical machines with their controllers can emulate all the components in a TeDP power train. The whole system model in Matlab/Simulink was first developed and tested in simulation, and the simulation results showed that system dynamic characteristics could be implemented by using the closed-loop control of the electric motor drive systems. Then we designed a subscale experimental system to emulate the entire power system from the turbine engine to the propulsive fans. Firstly, we built a system to emulate a gas turbine engine driving a generator, consisting of two permanent magnet (PM) motors with brushless motor drives, coupled by a shaft. We programmed the first motor and its drive to mimic the speed-torque characteristic of the gas turbine engine, while the second motor and drive act as a generator and produce a torque load on the first motor. Secondly, we built another system of two PM motors and drives to emulate a motor driving a propulsive fan. We programmed the first motor and drive to emulate a wound-rotor synchronous motor. The propulsive fan was emulated by implementing fan maps and flight conditions into the fourth motor and drive, which produce a torque load on the driving motor. The stator of each PM motor is designed to travel axially to change the coupling between rotor and stator. This feature allows the PM motor to more closely emulate a wound-rotor synchronous machine. These techniques can convert the plain motor system into a unique TeDP power grid emulator that enables real-time simulation performance using hardware-in-the-loop (HIL).

electric power grid emulator↗
DOE OSTI2014

New approach for absolute fluence distribution calculations in Monte Carlo simulations of light propagation in turbid media

A novel way to attain three dimensional fluence rate maps from Monte-Carlo simulations of photon propagation is presented in this work. The propagation of light in a turbid medium is described by the radiative transfer equation and formulated in terms of radiance. For many applications, particularly in biomedical optics, the fluence rate is a more useful quantity and directly derived from the radiance by integrating over all directions. Contrary to the usual way which calculates the fluence rate from absorbed photon power, the fluence rate in this work is directly calculated from the photon packet trajectory. The voxel based algorithm works in arbitrary geometries and material distributions. It is shown that the new algorithm is more efficient and also works in materials with a low or even zero absorption coefficient. The capabilities of the new algorithm are demonstrated on a curved layered structure, where a non-scattering, non-absorbing layer is sandwiched between two highly scattering layers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Online Robot Dead Reckoning Localization Using Maximum Relative Entropy Optimization With Model Constraints

The principle of Maximum relative Entropy optimization was analyzed for dead reckoning localization of a rigid body when observation data of two attached accelerometers was collected. Model constraints were derived from the relationships between the sensors. The experiment's results confirmed that accelerometers each axis' noise can be successfully filtered utilizing dependency between channels and the dependency between time series data. Dependency between channels was used for a priori calculation, and a posteriori distribution was derived utilizing dependency between time series data. There was revisited data of autocalibration experiment by removing the initial assumption that instantaneous rotation axis of a rigid body was known. Performance results confirmed that such an approach could be used for online dead reckoning localization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗