Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Numerical studies on saturated kink and sawtooth induced fast ion transport in JET ITER-like plasmas

Here, this presentation examines the energetic particle transport induced by saturated kink modes and sawtooth crashes in JET deuterium plasmas. It is known that kink mode-resonant transport and phase-space redistribution from sawtooth crashes can drive strong fast ion transport with dependencies on particle pitch and energy. Measurements with JET's Faraday cup fast ion loss detector array have shown that the internal kink growth phase preceding sawtooth crashes produces substantial fast ion losses. This report will numerically investigate the dominant energetic particle transport mechanism with a detailed examination of the fast ion phase-space dependencies, resonances, orbit topology changes, induced losses, and redistribution associated with the long-lived, resonant, kink mode and non-resonant sawtooth crash. The ORBIT-kick model forms the basis of the transport studies with realistic fast ion distributions produced from TRANSP. A recently created reduced model for sawtooth induced transport is used while the saturated kink modes are modeled with ideal MHD codes. The simulations were further validated against experiment with a newly developed synthetic Faraday cup fast ion loss detector in addition to scintillator probe and neutron measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Modeling heat transport processes in enhanced geothermal systems: Validation study from EGS Collab Experiment 1

Heat recovery from enhanced geothermal systems (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Wu, Hui↗
NASA NTRS1986

Evaluation of 3 numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finitie volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, S. F.↗
NASA NTRS1986

Evaluation of three numerical methods for propulsion integration studies on transonic transport configurations

An effort has been undertaken at the NASA Langley Research Center to assess the capabilities of available computational methods for use in propulsion integration design studies of transonic transport aircraft, particularly of pylon/nacelle combinations which exhibit essentially no interference drag. The three computer codes selected represent state-of-the-art computational methods for analyzing complex configurations at subsonic and transonic flight conditions. These are: EULER, a finite volume solution of the Euler equation; VSAERO, a panel solution of the Laplace equation; and PPW, a finite difference solution of the small disturbance transonic equations. In general, all three codes have certain capabilities that allow them to be of some value in predicting the flows about transport configurations, but all have limitations. Until more accurate methods are available, careful application and interpretation of the results of these codes are needed.

Yaros, Steven F.↗
NASA NTRS2018

Numerical Studies of Hall Thruster Acceleration Region Electron Transport

Azimuthal direction 1-D and azimuthal and axial direction 2-D Particle in Cell (PIC) simulations of the Hall thruster acceleration and near-plume regions are presented. The axial extent of the 2-D computational domain is based on the ExB drift parameter, 𝑫𝑬𝒙𝑩, the number of times around the channel an electron ExB drifts before being scattered. With only high energy main beam ions, the simulations produced a broad potential profile in the acceleration region. The addition of a small amount of ionization in the acceleration domain at a rate consistent with Hall thruster parameters causes the potential profile to develop two distinct regions. In the upstream portion of the computational domain, the potential profile is very steep, typical of Hall thruster acceleration regions. In the downstream portion, the potential profile is almost flat, like potentials in Hall thruster near plume regions. This bifurcation happens naturally; there is no other change in code parameters between the two results. Without the addition of a small amount of ionization in the acceleration region, the axial electron transport is very low. With ionization, the calculated ratio of electron current to ion current is the same order of those measured in Hall thrusters. This suggests that in regions with lower energy electrons and ions, plasma waves are effective in allowing crossfield transport, and axial electric fields are small. In the acceleration region, the plasma waves are less effective at scattering the energetic electrons, and high axial electric fields are needed to produce the observed cross-field electron transport.

Ortega, Alejandro Lopez↗
DOE OSTI2021

A Poisson–Nernst–Planck Model of Ion Transport and Interface Segregation in Metal–Insulator–Semiconductor Structures and Solar Cells

A numerical model that describes the transport of mobile ionic species in metal–insulator–semiconductor (MIS) and photovoltaic (PV) devices subject to temperature and voltage stress is presented. The finite element method (FEM) is used to solve the Nernst–Planck equation while imposing Poisson's equation self-consistently as a restriction for the electrostatic potential. This allows the contribution of the ionic species to the potential to be taken into account. Using a variational formulation eases the implementation of diverse boundary conditions, including the incorporation of segregation kinetics at the device interfaces. Segregation across the dielectric–semiconductor interface is relevant to modeling the electronic device degradation in systems where contamination reaches the semiconductor. The model in closed systems with no-flux boundary conditions is validated first. In the limiting case of low contamination levels with respect to the gate bias, the FEM solution matches analytically derived approximations. Then, the implementation is broadened to include an open boundary at the dielectric–semiconductor interface to account for leakage of ions. The predicted time dependence of the flatband voltage in Na-contaminated MIS test structures agrees well with measurements. This model successfully captures the role of long-range ion transport at concentrations of relevance to electronic and PV device instability and neuromorphic computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Conditional Point Sampling: A Monte Carlo Method for Radiation Transport in Stochastic Media.

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2002

Characterizing Thermal Properties of Melting Te Semiconductor: Thermal Diffusivity Measurements and Simulation

Tellurium is an element for many II-VI and I-III-VI(sub 2) compounds that are useful materials for fabricating many devices. In the melt growth techniques, the thermal properties of the molten phase are important parameter for controlling growth process to improve semiconducting crystal quality. In this study, thermal diffusivity of molten tellurium has been measured by a laser flash method in the temperature range from 500 C to 900 C. A pulsed laser with 1064 nm wavelength is focused on one side of the measured sample. The thermal diffusivity can be estimated from the temperature transient at the other side of the sample. A numerical simulation based on the thermal transport process has been also performed. By numerically fitting the experimental results, both the thermal conductivity and heat capacity can be derived. A relaxation phenomenon, which shows a slow drift of the measured thermal conductivity toward the equilibrium value after cooling of the sample, was observed for the first time. The error analysis and the comparison of the results to published data measured by other techniques will be discussed.

Zhu, Shen↗
NASA NTRS2002

Characterizing Thermal Properties of Melting Te Semiconductor: Thermal Diffusivity Measurements and Simulation

Tellurium is an element for many II-VI and I-III-VI(sub 2) compounds that are useful materials for fabricating many devises. In the melt growth techniques, the thermal properties of the molten phase are important parameter for controlling growth process to improve semiconducting crystal quality. In this study, thermal diffusivity of molten tellurium has been measured by a laser flash method in the temperature range from 500 C to 900 C. A pulsed laser with 1064 nm wavelength is focused on one side of the measured sample. The thermal diffusivity can be estimated from the temperature transient at the other side of the sample. A numerical simulation based on the thermal transport process has been also performed. By numerically fitting the experimental results, both the thermal conductivity and heat capacity can be derived. A relaxation phenomenon, which shows a slow drift of the measured thermal conductivity toward the equilibrium value after cooling of the sample, was observed for the first time. The error analysis and the comparison of the results to published data measured by other techniques will be discussed in the presentation.

Zhu, Shen↗
DOE OSTI2021

Building Blocks of Electric Vehicle Deployment: A Guide for Developing Countries

Countries can use electric transportation to help fulfill numerous goals, including greenhouse gas (GHG) emissions targets, local air quality goals, mobility objectives, energy security, and transportation resiliency. Vehicle electrification is a promising pathway to achieving clean energy transitions in the transport sector at scale. As vehicles electrify, the traditionally siloed electricity and transport sectors increasingly converge to create technical, institutional, and economic opportunities and challenges. To navigate this transition effectively, we propose the following foundational pillars or “building blocks” that undergird effective electric vehicle (EV) deployment. Jurisdictions with experience and mature EV markets offer useful lessons learned that may enable developing countries to leapfrog over common roadblocks. Nonetheless, developing countries face distinct challenges collectively, and individually, from developed countries that require careful consideration. From our work around the world, the U.S. Agency for International Development (USAID) and National Renewable Energy Laboratory (NREL) have witnessed interest growing in EVs for a variety of reasons. Lao People’s Democratic Republic (PDR), for example, looks toward EVs as a way to use their surplus of hydropower to displace expensive oil. Thailand and Pakistan seek economic development opportunities in EV manufacturing. Cities including Mexico City, Surat, India, and Kingston, Jamaica see the potential of bus electrification to improve local air quality and reduce traffic congestion. Many USAID partner countries in Southeast Asia are trying to understand how to reach ambitious EV deployment targets and implement EV and transportation plans. This pursuit is prompting questions about EV supply equipment (EVSE) standards, tariff design, and business models that affirm the importance of the building blocks outlined in this report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1981

Modeling of the float zone process

The interaction of heat, mass, and momentum transport in the floating zone method for growing single crystals from the melt is examined. Methods for detailed numerical simulation of the transport phenomena in a floating zone are developed. Results of the calculations are combined with experiments to determine the effects of solidification induced, surface tension driven, and buoyancy driven convection in establishing dopant redistribution in the melt and the roles of heat transfer in crystal and melt and melt/solid interface shape in determining crystal quality. State of the art finite element techniques were developed for calculating the influence of natural convection in the melt on the shape of a melt/crystal interface and dopant segregation in the crystal. These techniques are demonstrated for solidification by the Bridgman technique. Numerical techniques are developed that calculate the shapes of both the melt/solid and melt/gas interfaces simultaneously with the thermal fields in melt and solid. Models for the fluid flows due to the rotation of the feed and crystal rods are completed and the effects of these flows on dopant segregation are studied, especially in the case of zones longer than can be achieved on Earth.

Brown, R. A.↗
NASA NTRS1999

Numerical Modeling of Crystal of ZnSe by Physical Vapor Transport - Towards a more Comprehensive Formulations

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of PVT growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided. The effect of gravity on the process from numerical computations shows subtle effects although experimental evidence from vertically and horizontally grown samples show dramatic evidence of gravitational effects. The shortcomings of the problem formulation will be discussed and a framework will be provided leading up towards a more comprehensive model of PVT systems.

Ramachandran, N.↗
NASA NTRS1993

Growth of zinc selenide single crystals by physical vapor transport in microgravity

The goals of this research were the optimization of growth parameters for large (20 mm diameter and length) zinc selenide single crystals with low structural defect density, and the development of a 3-D numerical model for the transport rates to be expected in physical vapor transport under a given set of thermal and geometrical boundary conditions, in order to provide guidance for an advantageous conduct of the growth experiments. In the crystal growth studies, it was decided to exclusively apply the Effusive Ampoule PVT technique (EAPVT) to the growth of ZnSe. In this technique, the accumulation of transport-limiting gaseous components at the growing crystal is suppressed by continuous effusion to vacuum of part of the vapor contents. This is achieved through calibrated leaks in one of the ground joints of the ampoule. Regarding the PVT transport rates, a 3-D spectral code was modified. After introduction of the proper boundary conditions and subroutines for the composition-dependent transport properties, the code reproduced the experimentally determined transport rates for the two cases with strongest convective flux contributions to within the experimental and numerical error.

Rosenberger, Franz↗
NASA NTRS1996

Investigation of Thermal Creep and Thermal Stress Effects in Microgravity Physical Vapor Transport

Reported here are the results of our numerical investigation into the mechanisms which affect the transport and growth processes in physical vapor transport (PVT) crystal growth ampoules. The first part of the report consists of a brief summary of the major accomplishments and conclusions of our work. The second part consists of two manuscripts, submitted to the Journal of Crystal Growth, which provided a detailed description of the findings in our investigation.

Mackowski, D. W.↗
DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗
DOE OSTI2021

Exponential Time Differencing Schemes for Fuel Depletion and Transport in Molten Salt Reactors: Theory and Implementation

A numerical framework for modeling depletion and mass transport in liquid-fueled molten salt reactions is presented based on exponential time differencing. The solution method involves using the finite volume method to transform the system of partial differential equations (PDEs) into a much larger system of ordinary differential equations. The key part of this method involves solving for the exponential of a matrix. We explore six different algorithms to compute the exponential in a series of progression problems that explore physical transport phenomena in molten salt reactors. This framework shows good results for solving linear parabolic PDEs with each of the six matrix exponential algorithms. For large problems, the series solvers such as Padé and Taylor have large run times, which can be mitigated by using the Krylov subspace.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1995

Applications of direct numerical simulation of turbulence in second order closures

This paper discusses two methods of developing models for the rapid pressure-strain correlation term in the Reynolds stress transport equation using direct numerical simulation (DNS) data. One is a perturbation about isotropic turbulence, the other is a perturbation about two-component turbulence -- an extremely anisotropic turbulence. A model based on the latter method is proposed and is found to be very promising when compared with DNS data and other models.

Shih, Tsan-Hsing↗