Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled particle transport”

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 127 records · Page 7

Comparison of Organ Dose and Dose Equivalent Using Ray Tracing of Male and Female Voxel Phantoms to Space Flight Phantom Torso Data

Phantom torso experiments have been flown on the space shuttle and International Space Station (ISS) providing validation data for radiation transport models of organ dose and dose equivalents. We describe results for space radiation organ doses using a new human geometry model based on detailed Voxel phantoms models denoted for males and females as MAX (Male Adult voXel) and Fax (Female Adult voXel), respectively. These models represent the human body with much higher fidelity than the CAMERA model currently used at NASA. The MAX and FAX models were implemented for the evaluation of directional body shielding mass for over 1500 target points of major organs. Radiation exposure to solar particle events (SPE), trapped protons, and galactic cosmic rays (GCR) were assessed at each specific site in the human body by coupling space radiation transport models with the detailed body shielding mass of MAX/FAX phantom. The development of multiple-point body-shielding distributions at each organ site made it possible to estimate the mean and variance of space dose equivalents at the specific organ. For the estimate of doses to the blood forming organs (BFOs), active marrow distributions in adult were accounted at bone marrow sites over the human body. We compared the current model results to space shuttle and ISS phantom torso experiments and to calculations using the CAMERA model.

Kim, Myung-Hee Y.↗

Pair plasma instability in homogeneous magnetic guide fields

Pair plasmas, collections of both matter and antimatter particles of equal mass, represent a paradigm for the study of basic plasma science, and many open questions exist regarding these unique systems. They are found in many astrophysical settings, such as gamma-ray bursts, and have recently also been produced in carefully designed laboratory experiments. A central research topic in plasma physics is instability; however, unlike their more common ion–electron siblings, pair plasmas are generally thought to be stable to cross field pressure gradients in homogeneous magnetic fields. It is shown here by means of kinetic full-f simulations that, when a pressure gradient is first established, the Gradient-driven Drift Coupling mode is destabilized and becomes turbulent. Force balance is eventually achieved by a combination of flattened pressure profiles due to turbulent transport and establishment of a magnetic field gradient, saturating the growth. During the unstable phase, key physics can be captured by a δf gyrokinetic description, where it is shown analytically and numerically that parallel particle motion results in a coupling of all electromagnetic field components. A fluid model derived therefrom accurately predicts linear eigenmodes and is used to resolve global profile effects. For laser-based electron–positron plasma experiments, prompt instability is predicted with growth times much shorter than plasma lifetimes. Similarly, growth rates are calculated for the planned APEX experiment as well as gamma-ray burst scenarios, suggesting that the instability may contribute to the early evolution of these systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗

Interplanetary crew exposure estimates for the August 1972 and October 1989 solar particle events

Detailed exposure estimates for the bone marrow, ocular lens, and skin of astronauts on manned missions beyond the earth's magnetosphere have been made for the large solar particle events (SPE) of August 1972 and October 1989. The estimates were made using the coupled neutron-proton space radiation transport computer code BRYNTRN and the CAM model for the human body. It is found that at least 10 g/sq cm aluminum shielding is needed to keep estimated dose equivalent values below current space radiation exposure guidelines recommended for LEO missions. Space suits currently being designed will not provide adequate crew protection and may not even ensure crew survival in the event of exposure to a large SPE during extravehicular activity.

Townsend, Lawrence W.↗

Development and Validation of a Multidisciplinary Tool for Accurate and Efficient Rotorcraft Noise Prediction (MUTE)

A physics-based, systematically coupled, multidisciplinary prediction tool (MUTE) for rotorcraft noise was developed and validated with a wide range of flight configurations and conditions. MUTE is an aggregation of multidisciplinary computational tools that accurately and efficiently model the physics of the source of rotorcraft noise, and predict the noise at far-field observer locations. It uses systematic coupling approaches among multiple disciplines including Computational Fluid Dynamics (CFD), Computational Structural Dynamics (CSD), and high fidelity acoustics. Within MUTE, advanced high-order CFD tools are used around the rotor blade to predict the transonic flow (shock wave) effects, which generate the high-speed impulsive noise. Predictions of the blade-vortex interaction noise in low speed flight are also improved by using the Particle Vortex Transport Method (PVTM), which preserves the wake flow details required for blade/wake and fuselage/wake interactions. The accuracy of the source noise prediction is further improved by utilizing a coupling approach between CFD and CSD, so that the effects of key structural dynamics, elastic blade deformations, and trim solutions are correctly represented in the analysis. The blade loading information and/or the flow field parameters around the rotor blade predicted by the CFD/CSD coupling approach are used to predict the acoustic signatures at far-field observer locations with a high-fidelity noise propagation code (WOPWOP3). The predicted results from the MUTE tool for rotor blade aerodynamic loading and far-field acoustic signatures are compared and validated with a variation of experimental data sets, such as UH60-A data, DNW test data and HART II test data.

Liu, Yi↗

MESSENGER: Exploring Mercury's Magnetosphere

The MESSENGER mission to Mercury offers our first opportunity to explore this planet s miniature magnetosphere since the brief flybys of Mariner 10. Mercury s magnetosphere is unique in many respects. The magnetosphere of Mercury is among the smallest in the solar system; its magnetic field typically stands off the solar wind only - 1000 to 2000 km above the surface. For this reason there are no closed drift paths for energetic particles and, hence, no radiation belts. The characteristic time scales for wave propagation and convective transport are short and kinetic and fluid modes may be coupled. Magnetic reconnection at the dayside magnetopause may erode the subsolar magnetosphere allowing solar wind ions to impact directly the regolith. Inductive currents in Mercury s interior may act to modify the solar wind interaction by resisting changes due to solar wind pressure variations. Indeed, observations of these induction effects may be an important source of information on the state of Mercury s interior. In addition, Mercury s magnetosphere is the only one with its defining magnetic flux tubes rooted in a planetary regolith as opposed to an atmosphere with a conductive ionospheric layer. This lack of an ionosphere is probably the underlying reason for the brevity of the very intense, but short-lived, - 1-2 min, substorm-like energetic particle events observed by Mariner 10 during its first traversal of Mercury s magnetic tail. Because of Mercury s proximity to the sun, 0.3 - 0.5 AU, this magnetosphere experiences the most extreme driving forces in the solar system. All of these factors are expected to produce complicated interactions involving the exchange and re-cycling of neutrals and ions between the solar wind, magnetosphere, and regolith. The electrodynamics of Mercury s magnetosphere are expected to be equally complex, with strong forcing by the solar wind, magnetic reconnection at the magnetopause and in the tail, and the pick-up of planetary ions all driving field-aligned electric currents. However, these field-aligned currents do not close in an ionosphere, but in some other manner. In addition to the insights- into magnetospheric physics offered by study of the solar wind - Mercury system, quantitative specification of the "external" magnetic field generated by magnetospheric currents is necessary for accurate determination of the strength and multi-polar decomposition of Mercury s intrinsic magnetic field. MESSENGER S highly capable instrumentation and broad orbital coverage will greatly advance our understanding of both the origin of Mercury s magnetic field and the acceleration of charged particles in small magnetospheres. In. this article, we review what is known about Mercury s magnetosphere and describe the MESSENGER science team s strategy for obtaining answers to the outstanding science questions surrounding the interaction of the solar wind with Mercury and its small, but dynamic, magnetosphere.

Slavin, James A.↗

Apparatus, Method and Program Storage Device for Determining High-Energy Neutron/Ion Transport to a Target of Interest

An apparatus, method and program storage device for determining high-energy neutron/ion transport to a target of interest. Boundaries are defined for calculation of a high-energy neutron/ion transport to a target of interest; the high-energy neutron/ion transport to the target of interest is calculated using numerical procedures selected to reduce local truncation error by including higher order terms and to allow absolute control of propagated error by ensuring truncation error is third order in step size, and using scaling procedures for flux coupling terms modified to improve computed results by adding a scaling factor to terms describing production of j-particles from collisions of k-particles; and the calculated high-energy neutron/ion transport is provided to modeling modules to control an effective radiation dose at the target of interest.

Wilson, John W.↗

An open-source hybrid unstructured mesh - CAD fusion multiphysics analysis workflow in SALAMANDER

Plasma facing components in fusion devices will endure extreme neutron and heat fluxes. To facilitate their design using simulation tools, the open-source Fusion Module, Fusion ENergy Integrated multiphys-X (FENIX) framework is being developed to model these components with a high-fidelity multi-physics multi-dimensional approach. It can iteratively resolve couplings between all the physics at play, from neutron radiation, to thermomechanics, to near-wall plasma dynamics. This framework is based on the Multiphysics Object Oriented Simulation Environment (MOOSE), which is developed by a collaboration of US National Laboratories since 2008, for advanced nuclear, geomechanics simulations and other applications. FENIX couples numerous simulation tools, including OpenMC, the Tritium Migration Analysis Program v8, the NekRS CFD software, and most MOOSE modules. For the coupling of radiation transport and other physics, FENIX supports a hybrid workflow between Computer Assisted Design (CAD) and unstructured mesh geometries. The CAD can be generated from skinning the unstructured mesh, to enable a coarse geometry for efficient particle transport, but still resolving the local material compositions and temperature gradients. Neutron transport is performed using DAGMC on the CAD, and Cardinal, integrated in FENIX, maps tallied quantities, such as the heat deposition or the tritium generation rates, from a tally volumetric mesh to the other physics’ unstructured mesh. This coupling was exercised on a simplified tokamak geometry, coupling neutron transport with the heat conduction equation, and on a monoblock divertor problem, coupling additionally with tritium migration. Mesh convergence studies highlight the importance of the mapping conservativeness. Coupling with thermo-mechanics is further enabled by the generalization of the approach to moving meshes. The presentation will include these coupled analysis as well as an update on status of the FENIX framework.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Correlated particle transport enables biological free energy transduction

Studies of biological transport frequently neglect the explicit statistical correlations among particle site occupancies (i.e., they use a mean-field approximation). Neglecting correlations sometimes captures biological function, even for out-of-equilibrium and interacting systems. Here, we show that neglecting correlations fails to describe free energy transduction, mistakenly predicting an abundance of slippage and energy dissipation, even for networks that are near reversible and lack interactions among particle sites. Interestingly, linear charge transport chains are well described without including correlations, even for networks that are driven and include site-site interactions typical of biological electron transfer chains. We examine three specific bioenergetic networks: a linear electron transfer chain (as found in bacterial nanowires), a near-reversible electron bifurcation network (as in complex III of respiration and other recently discovered structures), and a redox-coupled proton pump (as in complex IV of respiration).

bacterial nanowires↗

NAIRAS Model Nowcasting and Forecasting of the Aviation Radiation Environment

The Nowcast of Aerospace Ionizing RAdiation System (NAIRAS) predicts dosimetric quantities for quantifying human radiation exposure and differential/integral flux/fluence quantities for assessing single event effects (SEE) in avionic systems from galactic cosmic rays (GCR), trapped inner belt protons, and solar energetic particle (SEP) events from the Earth’s surface to the space environment. Real-time predictions of the aviation radiation environment are available at NASA Goddard Space Flight Center’s Community Coordinated Modeling Center (CCMC) integrated Space Weather Analysis (iSWA) data feeds and cygnets, and a run-on-request (RoR) capability has also been deployed at CCMC. Recent model improvements include a more accurate atmospheric ionizing radiation transport methodology and more robust and reliable SEP nowcast dose predictions. In addition, preliminary results of SEP dose forecasts are shown by coupling the University of Malaga Solar Energetic Particle (UMASEP) model of integral proton flux forecasts with the NAIRAS model. These model updates and improvements are presented, and results are shown for aircraft, high-latitude balloon, and low-Earth orbit flights during quiescent and solar-geomagnetic disturbed conditions. Model comparisons with flight measurements are also shown.

NAIRAS↗

Chemical vapor infiltration of additively manufactured preforms: Pore-resolved simulations and experimental validation

The densification of additively manufactured porous preforms by chemical vapor infiltration (CVI) is studied using pore-resolved simulations and experiments. Experimentally, 3D printed silicon carbide (SiC) preforms are subject to CVI synthesis using methyltrichlorosilane (MTS) precursor to obtain high purity SiC/SiC composites. Optical images of the cross sections of the processed preforms are analyzed to obtain the spatial porosity distribution. The numerical method is based on a level set formulation to capture the spatial distribution and time evolution of the pore scale microstructural characteristics. The coupled transport and kinetic effects are represented using a dimensionless Thiele modulus. Simulations are initialized using representative synthetic preform geometries comprising of packed particles based on the size distribution of the powder used for 3D printing. The simulation results are validated against the experimental observations in terms of total density and the distribution of residual porosity. The densification characteristics, porosity classification, concentration profiles, and structure functions are analyzed as functions of processing temperature and Thiele modulus.

36 MATERIALS SCIENCE↗

Implications of reduced-complexity aerosol thermodynamics on organic aerosol mass concentration and composition over North America

Atmospheric organic aerosol (OA) mass concentrations can be affected by water uptake through its impact on the gas–particle partitioning of semi-volatile compounds. Current chemical transport models (CTMs) neglect this process. We have implemented the Binary Activity Thermodynamics model coupled to a volatility basis set partitioning scheme in the GEOS-Chem CTM, providing an efficient reduced-complexity OA model that predicts relative-humidity-dependent mixing and partitioning thermodynamics, while limiting the impact on computational efficiency. We provide a quantitative assessment of this water-sensitive OA treatment, focusing on a subdomain over North America. The updated OA scheme predicts a spatiotemporal mean enhancement in surface-level OA mass concentration of 145 % for January 2019 and 76 % for July 2019 compared to GEOS-Chem's most advanced OA scheme. The temporal mean surface-level OA organic mass concentration can increase by up to ∼590 % for January 2019 and ∼280 % for July 2019, with the greatest enhancements occurring over the ocean. The updated OA scheme also quantifies the OA-associated water content. The simulations show how different OA precursors and related OA surrogates contribute and respond to water uptake, including that due to changes in temperature and relative humidity over the diurnal cycle in selected winter and summer months. These results are independent of future CTM improvements involving updates to chemical reaction schemes and emission inventories. Our water-sensitive OA scheme allows for a better representation of the seasonal and regional variations in OA mass concentration in CTMs.

54 ENVIRONMENTAL SCIENCES↗

Investigation of CAD-based Geometry Workflows for Multiphysics Fusion Problems Using OpenMC and MOOSE

Fusion system designs are complex and require intricate and accurate meshes to be properly modeled. In this study, we investigate the use of CAD-based geometry workflows in fusion systems multiphysics problems. A simplified tokamak was introduced and modeled in CAD using a multiphysics coupling of OpenMC Monte Carlo transport and MOOSE heat conduction. The meshed geometry was prepared using direct accelerated geometry Monte Carlo (DAGMC) for particle transport, and a volumetric mesh was also prepared to be used in MOOSE and to tally OpenMC results. Cardinal was used to run OpenMC Monte Carlo particle transport within MOOSE framework. The heat source distribution and tritium production were calculated in OpenMC. The data transfer system was used to transfer heat source and temperature distribution between OpenMC and MOOSE. Two computational studies related to mesh refinement were performed: (1) refining the DAGMC and volumetric meshes used for tallying results and solving heat conduction and (2) only refining the DAGMC particle transport mesh. The refinement of the tally mesh has a much larger effect on the runtime compared to the refinement of the DAGMC particle transport surface mesh.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the Ohmic-dominant heating mode of capacitively coupled plasma inverted by boundary electron emission

Electron emission from the boundary is ubiquitous in a capacitively coupled plasma (CCP) and precipitates nonnegligible influence on the discharge properties. Here, we present Particle-in-Cell/Monte Carlo Collision simulation of an Ohmic-dominant heating mode of the capacitively coupled plasma, where the stochastic heating vanishes and only Ohmic heating sustains the discharge due to sheath inversion by boundary electron emission. The inverted CCP features negative sheath potential without Bohm presheath, hence excluding plasma heating due to sheath edge oscillation. The particle and energy transport of the proposed heating mode is analyzed. The influence of boundary electron emission flux, source voltage, and neutral pressure on the transition between classic and Ohmic-dominant CCP heating modes is shown with designated simulation scans. Further, a modified inverse sheath–plasma coupling due to excessive ionization is discovered. In the end, key indicators of the proposed heating mode in plasma diagnostics are provided for future experimental verifications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗

High-order diamond differencing schemes for the Boltzmann Fokker-Planck equation in 3D Cartesian geometries

The Boltzmann Fokker-Planck, an approximate form of the linear Boltzmann equation is commonly used to treat efficiently the transport of charged particles in matter. This paper introduces the application of high-order diamond differencing schemes (HODD), specifically the DD1 and DD2 schemes which are 4- and 6-order accurate respectively, to handle the spatial discretization of that equation in 3D Cartesian geometries. The energy deposition solutions for the coupled transport of electrons and photons presented in this work shows that HODD, compared to classical DD scheme, provides correction to the oscillations and a reduced propensity to yield negative fluxes. They are useful tools to minimize local error, notably in regions with abrupt variations of the flux solution. They also can be used to reduced execution time by decreasing the needed number of voxels to obtain a fixed accuracy. On the tested benchmarks, the DD1 scheme is 87%- 92%-91% more accurate than the classical DD scheme for total, mean per-voxels and maximum deviation of energy deposition values respectively. For comparison, a calculation with 8 times more voxels, requiring roughly 2.5 times more time to execute, is 92%-90%-77% more accurate. (authors)

97 MATHEMATICS AND COMPUTING↗

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity of Simulated Warm Rain Formation to Collision and Coalescence Efficiencies, Breakup, and Turbulence: Comparison of Two Bin-Resolved Numerical Models

Numerical models that resolve cloud particles into discrete mass size distributions on an Eulerian grid provide a uniquely powerful means of studying the closely coupled interaction of aerosols, cloud microphysics, and transport that determine cloud properties and evolution. However, such models require many experimentally derived paramaterizations in order to properly represent the complex interactions of droplets within turbulent flow. Many of these parameterizations remain poorly quantified, and the numerical methods of solving the equations for temporal evolution of the mass size distribution can also vary considerably in terms of efficiency and accuracy. In this work, we compare results from two size-resolved microphysics models that employ various widely-used parameterizations and numerical solution methods for several aspects of stochastic collection.

Fridlind, Ann↗