Search NASA⌕ Search

SEARCH · Search NASA

Results for “Deterministic particle 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.

80 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Sensitivity Analysis of Particle-In-Cell Modeling Parameters in Settling Bed, Bubbling Fluidized Bed and Circulating Fluidized Bed

The objective of the work presented is to perform a preliminary sensitivity analysis of particle-in-cell (PIC) model parameters when applied to settling bed, bubbling fluidized bed, and circulating fluidized bed simulations. These examples correspond to widely different flow conditions commonly seen in chemical engineering applications. Simulations were performed using the PIC method in the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). As part of the non-intrusive uncertainty quantification (UQ) analysis, simulation campaigns were generated using Nodeworks. Sampling locations or settings for PIC model parameters were determined using the Latin Hypercube method. Response surfaces were created using radial basis functions (RBF), and Sobol’ indices were estimated to quantify the influence of model parameters on the quantities of interest (QoI). This study marks a first step towards systematically determining optimal ranges for model parameters used in MFiX-PIC. Based on limited experience, it is expected that these values would depend strongly on flow conditions. Given the complexity of the multiphase flow systems under analysis, a non-intrusive UQ based approach is used to identify the most influential parameters in each case. This prior knowledge will help in proposing an effective design of experiments (DoE) and determine optimal parameters through techniques such as deterministic or Bayesian calibration, which will be pursued in the future.

42 ENGINEERING↗
DOE OSTI2020

Methods and apparatus for making magnetic skyrmions

A system and a method for the deterministic generation of magnetic skyrmions includes a magnetic strip configured to store and transport skyrmions. The magnetic strip includes one or more spatial inhomogeneities configured to generate a skyrmion at known locations when excited by a current pulse. A current pulse generator is used to inject current pulses into the magnetic strip via contact pads electrically coupled to both the current pulse generator and the magnetic strip. The system also includes a magnetic field source to apply an out-of-plane magnetic field across the magnetic strip to facilitate generation of skyrmions. Skyrmions can be generated by applying an out-of-plane magnetic field to the magnetic strip and injecting a current pulse with sufficient current density towards the spatial inhomogeneities. Once a skyrmion is generated, another current pulse with sufficient current density can be injected to move the skyrmion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Some Notes on Neutron Up-Scattering and the Doppler-Broadening of High-Z Scattering Resonances

When neutrons are scattered by target nuclei at elevated temperatures, it is entirely possible that the neutron will actually gain energy (i.e., up-scatter) from the interaction. This phenomenon is in addition to the more usual case of the neutron losing energy (i.e., down-scatter). Furthermore, the motion of the target nuclei can also cause extended neutron down-scattering, i.e., the neutrons can and do scatter to energies lower than predicted by the simple asymptotic models. In recent years, more attention has been given to temperature-dependent scattering cross sections for materials in neutron multiplying systems. This has led to the inclusion of neutron up-scatter in deterministic codes like Partisn and to free gas scattering models for material temperature effects in Monte Carlo codes like MCNP and cross section processing codes like NJOY. The free gas scattering models have the effect of Doppler Broadening the scattering cross section output spectra in energy and angle. The current state of Doppler-Broadening numerical techniques used at Los Alamos for scattering resonances will be reviewed, and suggestions will be made for further developments. The focus will be on the free gas scattering models currently in use and the development of new models to include high-Z resonance scattering effects. These models change the neutron up-scattering behavior.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1994

Advantages of estimating rate corrections during dynamic propagation of spacecraft rates: Applications to real-time attitude determination of SAMPEX

This paper describes real-time attitude determination results for the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX), a gyroless spacecraft, using a Kalman filter/Euler equation approach denoted the real-time sequential filter (RTSF). The RTSF is an extended Kalman filter whose state vector includes the attitude quaternion and corrections to the rates, which are modeled as Markov processes with small time constants. The rate corrections impart a significant robustness to the RTSF against errors in modeling the environmental and control torques, as well as errors in the initial attitude and rates, while maintaining a small state vector. SAMPLEX flight data from various mission phases are used to demonstrate the robustness of the RTSF against a priori attitude and rate errors of up to 90 deg and 0.5 deg/sec, respectively, as well as a sensitivity of 0.0003 deg/sec in estimating rate corrections in torque computations. In contrast, it is shown that the RTSF attitude estimates without the rate corrections can degrade rapidly. RTSF advantages over single-frame attitude determination algorithms are also demonstrated through (1) substantial improvements in attitude solutions during sun-magnetic field coalignment and (2) magnetic-field-only attitude and rate estimation during the spacecraft's sun-acquisition mode. A robust magnetometer-only attitude-and-rate determination method is also developed to provide for the contingency when both sun data as well as a priori knowledge of the spacecraft state are unavailable. This method includes a deterministic algorithm used to initialize the RTSF with coarse estimates of the spacecraft attitude and rates. The combined algorithm has been found effective, yielding accuracies of 1.5 deg in attitude and 0.01 deg/sec in the rates and convergence times as little as 400 sec.

Challa, M. S.↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

Spin–orbit torques in structures with asymmetric dusting layers

Current-induced spin–orbit torques (SOTs) in heavy metal/ferromagnet heterostructures have emerged as an efficient method for magnetization switching with applications in nonvolatile magnetic memory and logic devices. However, experimental realization of SOT switching of perpendicular magnetization requires an additional inversion symmetry breaking, calling for modifications of the conventional SOT heterostructures. In this work, we study SOTs and deterministic switching of perpendicular magnetization by inserting different asymmetric dusting layers at the heavy metal/ferromagnet interface. Similar to the previous works with lateral structural asymmetry, we study the emergence of current-induced perpendicular effective magnetic fields (H z eff ). By examining three different material combinations of heavy metal/dusting layers (W/IrMn, Pt/IrMn, and W/Ta), we shed light on the origins of H z eff ; we show that H z eff is generically created in all the studied asymmetric structures, has a close correlation with the interfacial magnetic anisotropy, and is independent of the signs of spin Hall angles of the materials. Furthermore, we show that the induction of H z eff enables field-free deterministic SOT switching of perpendicular magnetization. Our results can be used in designing SOT heterostructures for practical applications in nonvolatile technologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

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

Testing Nesting DXTRAN Spheres for Neutron Transport in Air

The deterministic transport sphere (DXTRAN sphere) is a technology that can be utilized in MCNP for improving sampling in the vicinity of a particular geometry, when the solid angle, between the source and particular geometry, is small. At a collision event, a DXTRAN particle is created and deterministically transported to the DXTRAN sphere, while the colliding particle (NONDXTRAN particle) continues transporting as if no DXTRAN event had occurred. Weight is balanced by killing the NONDXTRAN particle when traversing the DXTRAN sphere. If set up correctly, the DXTRAN particle placed on the DXTRAN sphere has a high probability of interacting with the geometry within the DXTRAN sphere. The Central Limit Theorem (CLT) states that the estimated mean approaches a normal distribution with known variance as the number of random samples grows sufficiently large. In reality the variance is seldom known and therefore must be estimated. In a Monte Carlo calculation, one must determine when the number of random samples is sufficiently large for the CLT to be satisfied. For the CLT to be satisfied, the first two moments (mean and variance) must exist. The second moment, the variance of a normal distribution is estimated. The slope test, which uses the largest 201 history scores of each tally fluctuation chart (TFC) bin fit to a Pareto function, is used to determine if the largest history scores decrease faster than 1/x{sup 3}. If the slope is less than 3, then enough histories have not yet been sampled such to meet the second moment existence requirement of the CLT. Usually, a slope of less than 3 indicates that the effects of high weight scores have not been satisfactorily sampled; such a consequence can also be manifested in examining trends in the 4. moment, the variance of the variance (VOV). For example, when examining the TFC, if the VOV is monotonically decreasing for a few sets of histories and then drastically increases for the next set of histories, chances are a rare high weight score, significantly greater than the average score, was encountered. As a result of the 1/R{sup 2} term, a point detector can achieve an extremely high weight, much greater than the average, as the distance from collision to point detector approaches zero. To remedy this difficulty, a 'radius of exclusion' is applied, where within the radius of exclusion 'the point detector estimate is assumed to be the average flux uniformly distributed'. A DXTRAN sphere can suffer from the same weight fluctuation as collisions in the vicinity of the DXTRAN sphere produce much larger weights than collisions happening far from the DXTRAN sphere. If the rare large weight put on the DXTRAN sphere then scores in a tally, there will be a significant increase in the VOV and the slope may drop below 3. To compliment that analysis, here we utilize an air composition and density that varies as a function of altitude, detailing the trials and tribulations of attempting to pass the 10 statistical checks for a variety of single and nested DXTRAN sphere approaches. The objective here is to illuminate important considerations for radiation transport problems with small solid angles in a scatter medium of few mean free paths (MFPs)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗