SEARCH · Engineering Papers

Results for “Input Sensitivity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Expendable vs reusable propulsion systems cost sensitivity

One of the key trade studies that must be considered when studying any new space transportation hardware is whether to go reusable or expendable. An analysis is presented here for such a trade relative to a proposed Liquid Rocket Booster which is being studied at MSFC. The assumptions or inputs to the trade were developed and integrated into a model that compares the Life-Cycle Costs of both a reusable LRB and an expendable LRB. Sensitivities were run by varying the input variables to see their effect on total cost. In addition a Monte-Carlo simulation was run to determine the amount of cost risk that may be involved in a decision to reuse or expend.

Hamaker, Joseph W.↗

NASA NTRS2015

Simulating the X-Ray Image Contrast to Set-Up Techniques with Desired Flaw Detectability

The paper provides simulation data of previous work by the author in developing a model for estimating detectability of crack-like flaws in radiography. The methodology is being developed to help in implementation of NASA Special x-ray radiography qualification, but is generically applicable to radiography. The paper describes a method for characterizing X-ray detector resolution for crack detection. Applicability of ASTM E 2737 resolution requirements to the model are also discussed. The paper describes a model for simulating the detector resolution. A computer calculator application, discussed here, also performs predicted contrast and signal-to-noise ratio calculations. Results of various simulation runs in calculating x-ray flaw size parameter and image contrast for varying input parameters such as crack depth, crack width, part thickness, x-ray angle, part-to-detector distance, part-to-source distance, source sizes, and detector sensitivity and resolution are given as 3D surfaces. These results demonstrate effect of the input parameters on the flaw size parameter and the simulated image contrast of the crack. These simulations demonstrate utility of the flaw size parameter model in setting up x-ray techniques that provide desired flaw detectability in radiography. The method is applicable to film radiography, computed radiography, and digital radiography.

Koshti, Ajay M.↗

NASA NTRS2012

Prevention of Over-Pressurization During Combustion in a Sealed Chamber

The combustion of flammable material in a sealed chamber invariably leads to an initial pressure rise in the volume. The pressure rise is due to the increase in the total number of gaseous moles (condensed fuel plus chamber oxygen combining to form gaseous carbon dioxide and water vapor) and, most importantly, the temperature rise of the gas in the chamber. Though the rise in temperature and pressure would reduce with time after flame extinguishment due to the absorption of heat by the walls and contents of the sealed spacecraft, the initial pressure rise from a fire, if large enough, could lead to a vehicle over-pressure and the release of gas through the pressure relief valve. This paper presents a simple lumped-parameter model of the pressure rise in a sealed chamber resulting from the heat release during combustion. The transient model considers the increase in gaseous moles due to combustion, and heat transfer to the chamber walls by convection and radiation and to the fuel-sample holder by conduction, as a function of the burning rate of the material. The results of the model are compared to the pressure rise in an experimental chamber during flame spread tests as well as to the pressure falloff after flame extinguishment. The experiments involve flame spread over thin solid fuel samples. Estimates of the heat release rate profiles for input to the model come from the assumed stoichiometric burning of the fuel along with the observed flame spread behavior. The sensitivity of the model to predict maximum chamber pressure is determined with respect to the uncertainties in input parameters. Model predictions are also presented for the pressure profile anticipated in the Fire Safety-1 experiment, a material flammability and fire safety experiment proposed for the European Space Agency (ESA) Automated Transfer Vehicle (ATV). Computations are done for a range of scenarios including various initial pressures and sample sizes. Based on these results, various mitigation approaches are suggested to prevent vehicle over-pressurization and help guide the definition of the space experiment.

Gokoglu, Suleyman A.↗

DOE OSTI2017

Annual climate variation modifies nitrogen induced carbon accumulation of Pinus sylvestris forests

We report results from long-term simulated external nitrogen (N) input experiments in three northern Pinus sylvestris forests, two of moderately high and one of moderately low productivity, assessing effects on annual net primary production (NPP) of woody mass and its interannual variation in response to variability in weather conditions. A sigmoidal response of wood NPP to external N inputs was observed in the both higher and lower productivity stands, reaching a maximum of ~65% enhancement regardless of the native site productivity, saturating at an external N input of 4–5 g N·m -2 ·yr -1 . The rate of increase in wood NPP and the N response efficiency (REN, increase in wood NPP per external N input) were maximized at an external N input of ~3 g N·m -2 ·yr -1 , regardless of site productivity. The maximum REN was greater in the higher productivity than the lower productivity stand (~20 vs. ~14 g C/g N). The N-induced enhancement of wood NPP and its REN were, however, markedly contingent on climatic variables. In both of the higher and lower productivity stands, wood NPP increased with growing season precipitation (P), but only up to ~400 mm. The sensitivity of the response to P increased with increasing external N inputs. Increasing growing season temperature (T) somewhat increased the N-induced drought effect, whereas decreasing T reduced the drought effect. These responses of wood NPP infused a large temporal variation to REN, making the use of a fixed value unadvisable. Based on these results, we suggest that regional climate conditions and future climate scenarios should be considered when modeling carbon sequestration in response to N deposition in boreal P. sylvestris, and possibly other forests.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Graph-based Reversible Evaluation and Tangents Library

GRETL is a C++ library for evaluation, re-evaluation and algorithmic differentiation of functional operations on an arbitrary computational graph with limited memory usage. Similar to popular machine learning frameworks in Python, like PyTorch and JAX, it tracks and stores both operations and output data as functions are evaluated. Once this composition of functions is built up, the entire chain of operations can be back propagated to compute sensitivities of the final result with respect to any number of inputs. In contrast to most machine learning applications, memory usage becomes the bottleneck for back propagation in many physics applications, especially for time-dependent PDEs. Dynamic check pointing becomes essential. An important distinguishing feature of GRETL is its ability to limit the maximum memory usage by automatically dynamic checkpointing the data output for each graph operation (see Wang, Moin, Iaccarino, 2009). During backpropagation, parts of the graph that are no longer in memory are automatically re-evaluated from upstream checkpointed states as needed for derivative sensitivity calculations (or more precisely, for vector-Jacobian products). GRETL is particularly beneficial for applications, such as coupled multi-physics, where deriving adjoint-based sensitivities and managing checkpoint memory across modules becomes onerous. Cases which can be readily handled by the GRETL library include: different time-integration algorithms per physics (e.g., coupled predictor-corrector algorithms, IMEX, etc.), sub-cycling, asynchronous integrators, state dependent timestep sizes, iterative solvers and coupling algorithms, controller algorithms, and more.

Tupek, MichaelR [Lawrence Livermore National Labor↗

NASA NTRS2017

Evaluation of HZETRN on the Martian Surface: Sensitivity Tests and Model Results

The Mars Science Laboratory Radiation Assessment Detector (MSLRAD) is providing continuous measurements of dose, dose equivalent, and particle flux on the surface of Mars. These measurements have been highly useful in validating environmental and radiation transport models that will be heavily relied upon for future deep space missions. In this work, the HZETRN code is utilized to estimate radiation quantities of interest on the Martian surface. A description of the modeling approach used with HZETRN is given along with the various input models and parameters used to define the galactic cosmic ray (GCR) environment and Martian geometry. Sensitivity tests are performed to gauge the impact of varying several input factors on quantities being compared to MSLRAD data. Results from these tests provide context for inter-code comparisons presented in a companion paper within this issue. It is found that details of the regolith and atmospheric composition have a minimal impact on surface flux, dose, and dose equivalent. Details of the density variation within the atmosphere and uncertainties associated with specifying the vertical atmospheric thickness are also found to have minimal impact. Two widely used GCR models are used as input into HZETRN and it is found that the associated surface quantities are within several percent of each other.

Tony C Slaba↗

DOE OSTI2021

Input-response space-filling designs

Traditional space-filling designs are a convenient way to explore throughout an input space of flexible dimension and have design points close to any region where future predictions might be of interest. In some applications, there may be a model connecting the input factors to the response(s), which provides an opportunity to consider the spacing not only in the input space but also in the response space. In this paper, we present an approach for leveraging current understanding of the relationship between inputs and responses to generate designs that allow the experimenter to flexibly balance the spacing in these two regions to find an appropriate design for the experimental goals. Applications where good spacing of the observed response values include calibration problems where the goal is to demonstrate the adequacy of the model across the range of the responses, sensitivity studies where the outputs from a submodel may be used as inputs for subsequent models, and inverse problems where the outputs of a process will be used in the inverse prediction for the unknown inputs. Here, we use the multi-objective optimization method of Pareto fronts to generate multiple non-dominated designs with different emphases on the input and response space-filling criteria from which the experimenter can choose. The methods are illustrated through several examples and a chemical engineering case study.

42 ENGINEERING↗

DOE OSTI2016

Adaptive surrogate modeling by ANOVA and sparse polynomial dimensional decomposition for global sensitivity analysis in fluid simulation

The Polynomial Dimensional Decomposition (PDD) is employed in this work for the global sensitivity analysis and uncertainty quantification (UQ) of stochastic systems subject to a moderate to large number of input random variables. Due to the intimate connection between the PDD and the Analysis of Variance (ANOVA) approaches, PDD is able to provide a simpler and more direct evaluation of the Sobol' sensitivity indices, when compared to the Polynomial Chaos expansion (PC). Unfortunately, the number of PDD terms grows exponentially with respect to the size of the input random vector, which makes the computational cost of standard methods unaffordable for real engineering applications. In order to address the problem of the curse of dimensionality, this work proposes essentially variance-based adaptive strategies aiming to build a cheap meta-model (i.e. surrogate model) by employing the sparse PDD approach with its coefficients computed by regression. Three levels of adaptivity are carried out in this paper: 1) the truncated dimensionality for ANOVA component functions, 2) the active dimension technique especially for second- and higher-order parameter interactions, and 3) the stepwise regression approach designed to retain only the most influential polynomials in the PDD expansion. During this adaptive procedure featuring stepwise regressions, the surrogate model representation keeps containing few terms, so that the cost to resolve repeatedly the linear systems of the least-squares regression problem is negligible. The size of the finally obtained sparse PDD representation is much smaller than the one of the full expansion, since only significant terms are eventually retained. Consequently, a much smaller number of calls to the deterministic model is required to compute the final PDD coefficients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

HED Hydrodynamics in the Common Modeling Framework

Within the Common Modeling Framework (CMF), the HED Hydrodynamics Project has implemented a common methodology, or "Authority", for modeling high energy density (HED) experiments by extending development of the Eulerian Applications Project code Authority, "EAPA," to include relevant physics models, customization tools and templates, and in situ and post-processing capabilities. Additionally, several models for small-scale HED experiments were successfully added to the CMF and simulations reproduce the results of preexisting, benchmarked input decks. We expect use of the CMF to improve our effectiveness in developing predictive capabilities for experiments ranging from small-scale planar single-interface, single-shock to multi-interface, multi-shock configurations and up through multi-shell inertial confinement fusion (ICF) implosions. Additionally, this will strengthen our ability to design future targeted experiments. Using the CMF infrastructure and Authorities for HED modeling has three key impacts. First, it enables a version-controlled, pedigreed and archived base model for each experiment with tools that allow easy setup, execution and analysis of simulations to assess and/or expose sensitivities to various parameters. Second, having a common base set of inputs and physics definitions helps identify systematic differences between models in the codes and experimental data to support code validation, and build confidence in our predictive capabilities. Finally, designed to work in concert with multiple programs and projects including PEM, IC, OES, V&V and DSW, the CMF supports improved collaboration and integration through shared model definitions, data, simulations and results, enhancing the workflow required to underwrite the knowledge base and capabilities ultimately required for stockpile stewardship

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Scan strategies in EBM-printed IN718 and the physics of bulk 3D microstructure development

Three-dimensional (3D) characterization provides opportunities for understanding processing-structure relationships in additively manufactured (AM) materials. Bulk samples of Inconel 718 were fabricated via electron beam melting (EBM) in order to study microstructural development as a function of energy input and beam scan strategy. TriBeam tomography of bulk Inconel 718 microstructures built under steady-state growth conditions reveals the sensitivity of microstructure formation and evolution to machine process parameters. Here, samples manufactured using a narrow range of energy input per unit build area result in varied grain morphologies and crystallographic textures. Using TRUCHAS, a thermal simulation software, the thermal history of bulk scan strategies was predicted, and combined with a calibrated microstructure-processing map to accurately predict bulk grain morphologies. The solidification parameters and the 3D measured nucleation density are used to predict the transition between columnar and equiaxed grain morphologies, providing a process map to guide AM parameter choices to locally control as-printed microstructure. A two-dimensional metric for characterizing bulk grain morphology was also found to agree well with predictions from the process map calibrated by 3D data. Combined with 3D tomography and thermal modelling, the physics of structure development were understood at a new level of detail with respect to the competing processes of grain nucleation and epitaxial growth.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Investigating the age distribution of fracture discharge using multiple environmental tracers, Bedrichov Tunnel, Czech Republic

We investigated the transit time distribution (TTD) of discharge collected from fractures in the Bedrichov Tunnel, Czech Republic, using lumped parameter models and multiple environmental tracers. We then utilize time series of δ 18 O, δ 2 H and 3 H along with CFC measurements from individual fractures in the Bedrichov Tunnel of the Czech Republic to investigate the TTD, and the uncertainty in estimated mean travel time in several fracture networks of varying length and discharge. We also compare several TTDs, including the dispersion distribution, the exponential distribution, and a developed TTD which includes the effects of matrix diffusion. The effect of seasonal recharge is explored by comparing several seasonal weighting functions to derive the historical recharge concentration. We identify best fit mean ages for each TTD by minimizing the error-weighted, multi-tracer χ2 residual for each seasonal weighting function. We use this methodology to test the ability of each TTD and seasonal input function to fit the observed tracer concentrations, and the effect of choosing different TTD and seasonal recharge functions on the mean age estimation. We find that the estimated mean transit time is a function of both the assumed TTD and seasonal weighting function. Best fits as measured by the χ2 value were achieved for the dispersion model using the seasonal input function developed here for two of the three modeled sites, while at the third site, equally good fits were achieved with the exponential model and the dispersion model and our seasonal input function. The average mean transit time for all TTDs and seasonal input functions converged to similar values at each location. The sensitivity of the estimated mean transit time to the seasonal weighting function was equal to that of the TTD. These results indicated that understanding seasonality of recharge is at least as important as the uncertainty in the flow path distribution in fracture networks and that unique identification of the TTD and mean transit time is difficult given the uncertainty in the recharge function. But, the mean transit time appears to be relatively robust to the structural model uncertainty. The results presented here should be applicable to other studies using environmental tracers to constrain flow and transport properties in fractured rock systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Life cycle assessment of novel heat exchanger for dry cooling of power plants based on encapsulated phase change materials

Cooling systems in power plants account for approximately 40% of total freshwater withdrawals in the U.S. Due to dwindling access to freshwater resources worldwide, continued operation of wet cooling systems poses a significant engineering challenge. To reduce water consumption, a novel air-cooled heat exchanger has been developed using encapsulated phase change material (EPCM) for dry cooling of power plants. Compared to traditional finned-tube air-cooled condensers, this novel EPCM heat exchanger improves the heat transfer coefficient and power plant efficiency while reducing the pressure drop and cooling system cost. Life cycle assessment (LCA) and techno-economic analysis (TEA) are used to evaluate the environmental and economic performance of EPCM heat exchangers from cradle-to-grave and to compare them to wet cooling and traditional air-cooled condensers. A thermodynamic model is developed to predict the EPCM heat exchanger performance for plant-scale operations. Equipment and construction costs for heat exchangers are estimated based on design parameters obtained from the thermodynamic model. Both process-LCA and economic-input–output LCA are used to simulate and test the sensitivity of EPCM alternatives with commercial wet and dry cooling technologies. We investigate options for EPCM end-of-life management upon retiring the heat exchanger and construct a process-based LCA model to estimate a greenhouse gas (GHG) emissions credit for recycling the EPCM. The life cycle GHG emission of the novel dry cooling technology is 1.16 kg CO 2 eq. /MWh compared with the 1.1–4.3 kg CO 2 eq. /MWh reported for commercial dry cooling technologies and consumes 9.5 L/MWh e of water for cradle-to-gate life cycle, which is significantly lower than that of wet cooling systems. The TEA shows many advantages of EPCM cooling technology over the state-of-art dry cooling solutions. Overall, the EPCM heat exchanger provides a better alternative compared to existing dry cooling and wet cooling technologies.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2016

iMatTOUGH: An open-source Matlab-based graphical user interface for pre- and post-processing of TOUGH2 and iTOUGH2 models

TOUGH2 and iTOUGH2 are powerful models that simulate the heat and fluid flows in porous and fracture media, and perform parameter estimation, sensitivity analysis and uncertainty propagation analysis. However, setting up the input files is not only tedious, but error prone, and processing output files is time consuming. Here, we present an open source Matlab-based tool (iMatTOUGH) that supports the generation of all necessary inputs for both TOUGH2 and iTOUGH2 and visualize their outputs. The tool links the inputs of TOUGH2 and iTOUGH2, making sure the two input files are consistent. It supports the generation of rectangular computational mesh, i.e., it automatically generates the elements and connections as well as their properties as required by TOUGH2. The tool also allows the specification of initial and time-dependent boundary conditions for better subsurface heat and water flow simulations. The effectiveness of the tool is illustrated by an example that uses TOUGH2 and iTOUGH2 to estimate soil hydrological and thermal properties from soil temperature data and simulate the heat and water flows at the Rifle site in Colorado.

58 GEOSCIENCES↗

DOE OSTI2021

Effects of different vegetation drag parameterizations on the tidal propagation in coastal marshlands

Vegetation drag is a fundamental quantity directly affecting results for both long- and short-term coastal marsh and geomorphological studies. The vegetation drag in coastal marshland has been modeled by various two-dimensional (2D) and threedimensional (3D) numerical parameterizations. 2D parameterizations treat coastal marshes as bottom roughness elements, while 3D parameterizations resolve the vertically-variable vegetation drag through the water column. However, differences in tidal propagation arising from different drag parameterizations within a single model are largely unknown, and clear guidance on parameterization selection is still missing. In this study, we implemented four vegetation drag parameterizations into the Model for Prediction Across Scales-Ocean (MPAS-O), which include 1) a 2D parameterization using land-cover type-determined Manning’s n (2DLM); 2) a 2D parameterization using vegetation-determined Manning’s n (2DVM); 3) a 3D parameterization for stiff vegetation (3DSV); and 4) a 3D parameterization for flexible vegetation (3DFV). Estimates of the flow resistance effects from these parameterizations were compared using a series of idealized tidal propagation simulations. Given the same tidal condition, flooding depth and flooding distance are the largest in the 2DLM simulations and the smallest in the 3DSV simulations. 2DVM results are the closest to the 2DLM results. 3DFV results are the closest to the average of 2DVM, 3DSV, and 3DFV results. 2DVM and 3DSV results are the least and most sensitive to the vegetation aboveground biomass, respectively. Based on the input data requirement and computational efficiency of each parameterization, a comparison summary is provided to help inform parameterization selection for specific applications. Here, the effects of these parameterizations on coastal geomorphology are further discussed, and the results demonstrate that estimates of the long-term evolution of coastal marshes and coastal morphology depend upon the selection of the vegetation drag parameterization

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Transient uncertainty quantification and Global Sensitivity Analysis of the open-source Molten Chloride Reactor Experiment (MCRE) using GP-PCA surrogate models

Uncertainties in the thermophysical properties of molten salts impact both the steady-state and transient behavior of Molten Salt Reactors (MSRs). In this work, we aim to quantify the influence of such uncertainties on the transient operation of the Molten Chloride Reactor Experiment (MCRE), utilizing the open-source specifications provided for this reactor. Seven representative transient scenarios are considered. For each scenario, we evaluate the impact of thermophysical property uncertainties on four key multiphysics model output variables of interest (VoIs): maximum power density, maximum fuel temperature, maximum reflector temperature, and average fuel velocity magnitude. In addition, we perform a Global Sensitivity Analysis (GSA) by computing Sobol’ indices for the uncertain input parameters to determine their contribution to the variability of each VoI. Conducting GSA is computationally intensive due to the large number of required evaluations of the high-fidelity multiphysics model. To mitigate this cost, we develop a surrogate modeling framework that combines Gaussian Process (GP) regression with Principal Component Analysis (PCA), enabling efficient sample generation for the GSA. Our results show that for energy-related VoIs, thermal conductivity is the dominant contributor to uncertainty. In contrast, for flow-related VoIs, density and dynamic viscosity are the primary sources of uncertainty. The specific heat of the fuel salt was found to play a secondary role in the transient analyses.

42 - ENGINEERING↗

DOE OSTI2020

Uncertainty evaluation for twist drilling stability model

This paper describes the first uncertainty analysis for drilling stability using a frequency-domain drilling stability model. The stability model inputs include: the modal parameters for the torsional-axial vibration mode from the twist drill-holder-spindle axial frequency response function; and the mechanistic coefficients that relate the torque and thrust force to chip area for the selected drill-workpiece material combination. Furthermore, Monte Carlo simulation is applied to propagate the input uncertainties to output uncertainty in the predicted stability map, which separates stable from unstable (chatter) zones in the spindle speed-chip width parameter space. Additionally, the mean stability boundary and its 95% confidence intervals are determined for five cases: varying all four inputs simultaneously and varying them individually. This enables the individual sensitivities to be compared. Experimental results from drilling tests are included for comparison to the prediction. Additionally, Matlab code is provided to implement the stability model and Monte Carlo uncertainty analysis.

42 ENGINEERING↗

DOE OSTI2025

Life-cycle greenhouse gas emissions analysis of battery-grade lithium production in Finland

Various countries are undertaking initiatives to domestically produce battery-related critical materials. Within Finland, Keliber Technology Oy is developing capabilities for battery-grade lithium hydroxide monohydrate (LHM) production from spodumene ores. A detailed life-cycle assessment (LCA) of this pathway is conducted to determine its life-cycle GHG impacts using Argonne's R&D GREET (Research and Development Greenhouse gases, Regulated Emissions, and Energy use in Technologies) model. The analysis shows life-cycle GHG emissions of similar to 9.2 kg CO 2-eq /kg LHM, dominated by contributions from three energy sources - diesel, natural gas, and electricity - and two material inputs - lime (CaO) and soda ash (Na 2 CO 3 ). Sensitivity analyses highlight the potential to reduce these impacts using low-carbon electricity, sequestration of process CO 2 emissions generated during CaO and Na 2 CO 3 production, and bio-based energy for LHM production (by similar to 15 % each). A comparative analysis shows lower impacts for Keliber's LHM than for existing LHM production from Australian spodumene ores processed in China (by similar to 40 %).

CCS↗

DOE OSTI2016

The effects of electron thermal radiation on laser ablative shock waves from aluminum plasma into ambient air

The effect of electron thermal radiation on 7 ns laser ablative shock waves from aluminum (Al) plasma into an ambient atmospheric air has been numerically investigated using a one-dimensional, three-temperature (electron, ion, and radiation) radiation hydrodynamic code MULTI. The governing equations in Lagrangian form are solved using an implicit scheme for planar, cylindrical, and spherical geometries. The shockwave velocities (V{sub sw}) obtained numerically are compared with our experimental values obtained over the intensity range of 2.0 × 10{sup 10} to 1.4 × 10{sup 11 }W/cm{sup 2}. It is observed that the numerically obtained V{sub sw} is significantly influenced by the thermal radiation effects which are found to be dominant in the initial stage up to 2 μs depending on the input laser energy. Also, the results are found to be sensitive to the co-ordinate geometry used in the simulation (planar, cylindrical, and spherical). Moreover, it is revealed that shock wave undergoes geometrical transitions from planar to cylindrical nature and from cylindrical to spherical nature with time during its propagation into an ambient atmospheric air. It is also observed that the spatio-temporal evolution of plasma electron and ion parameters such as temperature, specific energy, pressure, electron number density, and mass density were found to be modified significantly due to the effects of electron thermal radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗