Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transport models”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Enhanced production of strange baryons in high-energy nuclear collisions from a multiphase transport model

Here, we introduce additional coalescence factors for the production of strange baryons in a multiphase transport (AMPT) model in order to describe the enhanced production of multistrange hadrons observed in Pb + Pb collisions at $\sqrt{s_{NN}}$ = 2.76 TeV at the Large Hadron Collider (LHC) and Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV at Relativistic Heavy-Ion Collider (RHIC). This extended AMPT model is found to also give a reasonable description of the multiplicity dependence of the strangeness enhancement observed in high multiplicity events in $\textit{pp}$ collisions at $\sqrt{s}$ = 7 TeV and $\textit{p}$-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV. We find that the coalescence factors depend on the system size but not much on whether the system is produced from A + A or $\textit{p}$ + A collisions. Here, the extended AMPT model thus provides a convenient way to model the mechanism underlying the observed strangeness enhancement in collisions of both small and large systems at RHIC and LHC energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A review on charged-particle transport modeling for laser direct-drive fusion

Inertial confinement fusion (ICF) with the laser-indirect-drive scheme has recently made a tremendous breakthrough recently after decades of intensive research effort. Taking this success to the next step, the ICF community is coming to a general consensus that laser direct-drive (LDD) fusion might be the viable way for enabling inertial fusion energy (IFE) and high-gain targets for other applications. Designing and understanding LDD fusion targets heavily rely on radiation-hydrodynamic code simulations, in which charged-particle transport plays an essential role in modeling laser-target energy coupling and bootstrap heating of fusion-produced α-particles. To better simulate charged-particle transport in LDD targets, over the past four decades the plasma physics community has advanced transport calculations from simple plasma physics models to sophisticated computations based on first-principles methods. In this review, we give an overview of the current status of charged-particle transport modeling for LDD fusion, including what challenges we still face and the possible paths moving forward to advance transport modeling for ICF simulations. We hope this review will provide a summary of exciting challenges to stimulate young minds to enter the field, facilitate further progress in understanding warm-dense matter physics, and ultimately bridge toward the success of reliable LDD fusion designs for IFE and other high-gain ICF applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A priori analysis of a power-law sub-filter scale mixing model for transported PDF modeling of high Karlovitz turbulent premixed flames

Accurate modeling of mixing in large-eddy simulation (LES) / transported probability density function (PDF) modeling of turbulent combustion remains an outstanding issue. The issue is particularly salient in turbulent premixed combustion under extreme conditions such as high-Karlovitz number Ka. Here, the present study addresses this issue by conducting an a priori analysis of a power-law scaling based mixing timescale model for the transported PDF model. A recently produced DNS dataset of a high-Ka turbulent jet flame is used for the analysis. A power-law scaling is observed for a scaling factor used to model the sub-filter scale mixing timescale in this high-Ka turbulent premixed DNS flame when the LES filter size is much greater than the characteristic thermal thickness of a laminar premixed flame. The sensitivity of the observed power-law scaling to the different viewpoints (local or global) and to the different scalars for the data analysis is examined and the dependence of the model parameters on the dimensionless numbers Ka and Re (the Reynolds number) is investigated. Different model formulations for the mixing timescale are then constructed and assessed in the DNS flame. The proposed model is found to be able to reproduce the mixing timescale informed by the high-Ka DNS flame significantly better than a previous model.

42 ENGINEERING↗
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Utilizing Environmental Tracers to Reduce Groundwater Flow and Transport Model Parameter Uncertainties

Abstract Non‐uniqueness in groundwater model calibration is a primary source of uncertainty in groundwater flow and transport predictions. In this study, we investigate the ability of environmental tracer information to constrain groundwater model parameters. We utilize a pilot point calibration procedure conditioned to subsets of observed data including: liquid pressures, tritium ( 3 H), chlorofluorocarbon‐12 (CFC‐12), and sulfur hexafluoride (SF 6 ) concentrations; and groundwater apparent ages inferred from these environmental tracers, to quantify uncertainties in the heterogeneous permeability fields and infiltration rates of a steady‐state 2‐D synthetic aquifer and a transient 3‐D model of a field site located near Riverton, Wyoming (USA). To identify the relative data worth of each observation data type, the post‐calibration uncertainties of the optimal parameters for a given observation subset are compared to that from the full observation data set. Our results suggest that the calibration‐constrained permeability field uncertainties are largest when liquid pressures are used as the sole calibration data set. We find significant reduction in permeability uncertainty and increased predictive accuracy when the environmental tracer concentrations, rather than apparent groundwater ages, are used as calibration targets in the synthetic model. Calibration of the Riverton field site model using environmental tracer concentrations directly produces infiltration rate estimates with the lowest uncertainties, however; permeability field uncertainties remain similar between the environmental tracer concentration and apparent groundwater age calibration scenarios. This work provides insight on the data worth of environmental tracer information to calibrate groundwater models and highlights potential benefits of directly assimilating environmental tracer concentrations into model parameter estimation procedures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of electron transport models on capillary discharge plasmas

Magnetohydrodynamics (MHD) can be used to model capillary discharge waveguides in laser-wakefield accelerators. However, the predictive capability of MHD can suffer due to poor microscopic closure models. Here, we study the impact of electron heating and thermal conduction on the capillary waveguide performance as part of an effort to understand and quantify uncertainties in modeling and designing next-generation plasma accelerators. To do so, we perform two-dimensional high-resolution MHD simulations using an argon-filled capillary discharge waveguide with three different electron transport coefficients models. The models tested include (i) Davies et al., (ii) Spitzer, and (iii) Epperlein–Haines (EH). We found that the EH model overestimates the electron temperature inside the channel by over 20% while predicting a lower azimuthal magnetic field. Furthermore, the Spitzer model, often used in MHD simulations for plasma-based accelerators, predicts a significantly higher electron temperature than the other models suggest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Permeate fluxes from desalination of brines and produced waters: A reactive transport modeling study

The increasing interest in the use of membrane systems to desalinate inland brackish water, agricultural drainage, and industrially produced wastewater demands improved means of predicting desalination system performance under variable feedwater compositions. The interaction among water flow, solute transport, and chemical composition in these systems impacts permeate flux evolution. Here, an established multicomponent reactive transport simulator that accounts for these coupled processes is applied to compute osmotic pressure and permeate fluxes in reverse osmosis (RO) systems. The model is first validated by predicting permeate fluxes for a set of benchtop crossflow experiments subject to a range of feed flow rates and compositions, under fouling and non-fouling conditions. Results compare favorably with measured data that show that solutions with similar total dissolved solids concentrations but different compositions result in different permeate fluxes. The model is then applied to predict permeate fluxes from the desalination of produced waters using a commercial spiral wound RO module. For NaCl-dominant brines, at total dissolved salt concentrations (TDS) below about 70 g/L, permeate fluxes are inversely proportional to water mole fraction as the latter is a reasonable approximation of water activity (i.e. ideal mixing). In the case of Ca–Cl-, Na–CO3- and Na–SO4-dominant brines below about 70 g/L TDS, this relationship does not hold as well and tends to overpredict osmotic pressure and thus underpredict permeate fluxes. However, the opposite becomes true at higher TDS values for typical produced waters. The scaling potential of these waters is also computed by allowing the precipitation of minerals above their saturation limit on the RO membrane. This work demonstrates how reactive transport models developed for the analysis of waters from geological systems can be extended to improve process design, optimization, and control in desalination systems from produced waters and beyond.

Molins, Sergi↗
DOE OSTI2022

Comparison of AGR-3/4 Fission Product Transport Model to Measurements and Extraction of Diffusivities via Analytical Fits

A one-dimensional finite-element model based on the MOOSE framework has been constructed, and overall transport estimates are predicted based on available diffusivities and sorption isotherms in the extant literature. The predicted transport profiles are compared to experimental results, where such results are available, including total inventory analysis of gamma emitters, liquid scintillation results for Sr-90 concentrations, radial deconsolidation leach-burn-leach within compacts, tomographic reconstructions based on detailed gamma scanning and destructive physical sampling of several rings. The magnitude of expected transport to the sink ring is generally predicted by the model for silver and cesium, however, modeled detailed concentration profiles are not consistent with experimental results. Strontium and europium transport is not well characterized by model results.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Availability of Critical Benchmark Experiments for the Pebble Tanker Transportation Model for Nuclear Criticality Safety Validation of TRISO Pebbles

This study addresses the need for comprehensive investigations into TRi-structural ISOtropic (TRISO) fuel pebble transportation validation. In this work, an exploratory model, the pebble tanker(PT), was developed with the aim of facilitating the validation of nuclear criticality safety calculations in the context of industrial-scale transportation of TRISO fuel. The PT model was designed to investigate the availability and applicability of critical benchmark experiments crucial for assessing the transportation of these pebbles. This work incorporated sensitivity/uncertainty (S/U) similarity studies to quantify the applicability of critical benchmark experiments and to address nuclear data uncertainties in the context of TRISO transportation. Two container models were investigated: one for the Hermes-type pebble and one for the Pebble Bed Modular Reactor (PBMR)–type pebble. The models were simplified, considering fuel, containment, and either water or air, to enable a focus on the underlying physics of applications involving TRISO fuel pebbles using the PT model. A crucial aspect under consideration was the capacity of the transport package to hold pebbles while ensuring subcriticality in the flooded state. An approach in the criticality validation process involves assessing the similarity between systems through an integral index parameter evaluation. This involves calculating a correlation coefficient (referred to as c k ) based on shared nuclear data–induced uncertainty between a benchmark experiment and the application of the PT model. To facilitate this analysis, the SCALE tools, particularly the CSAS6-Shift, TSUNAMI-3D-Shift, and TSUNAMI-IP sequences, were employed for comprehensive studies in neutronics and S/U analysis. Our findings showed that there are sufficient critical experimental benchmarks to perform this validation of the PT model in the most reactive state, i.e. when the tanker is flooded. This paper provides valuable insights into validating a transport package for Generation IV TRISO fuel pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Validation of heat transport modeling using directly driven beryllium spheres

Recent experiments involving directly driven beryllium spheres are reported. Plasma conditions are measured using Thomson scattering with the probe beam pointed 200, 300, and 400 μ m from the surface of the sphere. Laser coupling is assessed using calorimeters that collect scattered light placed at various locations within the target chamber. Laser intensities of 1 0 14 W / c m 2 and 2.5 × 1 0 14 W / c m 2 are chosen to minimize unmodeled laser-plasma interactions (LPIs) that lead to laser-target decoupling. Two-dimensional simulations are compared to the interpreted data using the radiation-hydrodynamics code Lasnex. Heat transport is simulated using flux-limited Spitzer–Harm with both high ( f = 0.15) and low ( f = 0.03) flux limiters and the Schurtz–Nicolai-Busquet (SNB) model. At 1 0 14 W / c m 2 , all three heat transport models agree well with the measurement, demonstrating that the heat flux is local at low intensities near the measurement locations. At 2.5 × 1 0 14 W / c m 2 , the SNB and high flux model roughly match the plasma conditions but predict 2% uncoupled light compared to 10% measured. Additionally, the use of drive multipliers to match the measured coupled light does not alter the agreement between measured and simulated plasma conditions, suggesting that decoupling due to LPI is unlikely to alter this agreement. The low flux model cannot match the plasma conditions and results in 19% scattered light. The use of a resonant absorption model can be used to bring the simulated scattered light into agreement, but the simulated plasma conditions are still in disagreement with the measurement. For this reason, the low flux model is rejected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX part 2: Impurity transport modeling and experimental comparison

The Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is a pulsed linear plasma device that uses a radio frequency (RF) helicon antenna and a helicon ceramic window as the plasma ionization source. Modeling performed within this work on the helicon operations predicts hot spots on the helicon ceramic window due to an RF-induced sheath. The RF sheath potential leads to physical sputtering of the ceramic window predominantly where the helicon antenna sits around the helicon window. The role of oxygen in sputtering within Proto-MPEX is also investigated and is needed to match the experimental data. Experiments on Proto-MPEX show reasonable agreement with the modeled erosion/redeposition pattern on the helicon window, and the modeled impurity flux to the target shows radial transport within the device as expected from classical collisional impurity transport. Temperature screening of impurities is shown to be needed for Proto-MPEX to reduce impurity flux to the target. The modeling approach used in this work can be applied to other plasma RF applications for the determination of impurity production and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Machine Learning Emulation of Spatial Deposition from a Multi-Physics Ensemble of Weather and Atmospheric Transport Models

In the event of an accidental or intentional hazardous material release in the atmosphere, researchers often run physics-based atmospheric transport and dispersion models to predict the extent and variation of the contaminant spread. These predictions are imperfect due to propagated uncertainty from atmospheric model physics (or parameterizations) and weather data initial conditions. Ensembles of simulations can be used to estimate uncertainty, but running large ensembles is often very time consuming and resource intensive, even using large supercomputers. In this paper, we present a machine-learning-based method which can be used to quickly emulate spatial deposition patterns from a multi-physics ensemble of dispersion simulations. We use a hybrid linear and logistic regression method that can predict deposition in more than 100,000 grid cells with as few as fifty training examples. Logistic regression provides probabilistic predictions of the presence or absence of hazardous materials, while linear regression predicts the quantity of hazardous materials. The coefficients of the linear regressions also open avenues of exploration regarding interpretability—the presented model can be used to find which physics schemes are most important over different spatial areas. A single regression prediction is on the order of 10,000 times faster than running a weather and dispersion simulation. However, considering the number of weather and dispersion simulations needed to train the regressions, the speed-up achieved when considering the whole ensemble is about 24 times. Ultimately, this work will allow atmospheric researchers to produce potential contamination scenarios with uncertainty estimates faster than previously possible, aiding public servants and first responders.

97 MATHEMATICS AND COMPUTING↗