Search NASA⌕ Search

SEARCH · Search NASA

Results for “process simulation”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Visualization and Analysis Requirements for In Situ Processing for a Large-Scale Fusion Simulation Code

In situ techniques have become a very active research area since they have been shown to be an effective way to combat the issues associated with the ever growing gap between computation and I/O bandwidth. In order to take full advantage of in situ techniques with a large-scale simulation code, it is critical to understand the breadth and depth of its analysis requirements. In this paper, we present the results of a survey done with members of the XGC1 fusion simulation code team in order to gather their requirements for analysis and visualization. We look at these requirements from the perspective of in situ processing and present a list of XGC1 analysis tasks performed by its physicists, engineers, and visualization specialists. This analysis of the specific needs and use cases of a single code is important in understanding the nature of the needs that simulations have in terms of data movement and usage for visualization and analysis, now and in the future.

Kress, James M.↗
DOE OSTI2017

AGN Heating in Simulated Cool-core Clusters

We analyze heating and cooling processes in an idealized simulation of a cool-core cluster, where momentum-driven AGN feedback balances radiative cooling in a time-averaged sense. We find that, on average, energy dissipation via shock waves is almost an order of magnitude higher than via turbulence. Most of the shock waves in the simulation are very weak shocks with Mach numbers smaller than 1.5, but the stronger shocks, although rare, dissipate energy more effectively. We find that shock dissipation is a steep function of radius, with most of the energy dissipated within 30 kpc, more spatially concentrated than radiative cooling loss. However, adiabatic processes and mixing (of post-shock materials and the surrounding gas) are able to redistribute the heat throughout the core. A considerable fraction of the AGN energy also escapes the core region. The cluster goes through cycles of AGN outbursts accompanied by periods of enhanced precipitation and star formation, over gigayear timescales. The cluster core is under-heated at the end of each cycle, but over-heated at the peak of the AGN outburst. During the heating-dominant phase, turbulent dissipation alone is often able to balance radiative cooling at every radius but, when this is occurs, shock waves inevitably dissipate even more energy. Our simulation explains why some clusters, such as Abell 2029, are cooling dominated, while in some other clusters, such as Perseus, various heating mechanisms including shock heating, turbulent dissipation and bubble mixing can all individually balance cooling, and together, over-heat the core.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Summertime low clouds mediate the impact of the large-scale circulation on Arctic sea ice

The rapid Arctic sea ice retreat in the early 21 st century is believed to be driven by several dynamic and thermodynamic feedbacks, such as ice-albedo feedback and water vapor feedback. However, the role of clouds in these feedbacks remains unclear since the causality between clouds and these processes is complex. Here, we use NASA CERES satellite products and NCAR CESM model simulations to suggest that summertime low clouds have played an important role in driving sea ice melt by amplifying the adiabatic warming induced by a stronger anticyclonic circulation aloft. The upper-level high pressure regulates low clouds through stronger downward motion and increasing lower troposphere relative humidity. The increased low clouds favor more sea ice melt via emitting stronger longwave radiation. Then decreased surface albedo triggers a positive ice-albedo feedback, which further enhances sea ice melt. Considering the importance of summertime low clouds, accurate simulation of this process is a prerequisite for climate models to produce reliable future projections of Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Exploring a Multi-resolution Approach Using AMIP Simulations

This study presents a diagnosis of a multi-resolution approach using the Model for Prediction Across Scales - Atmosphere (MPAS-A) for simulating regional climate. Four AMIP experiments are conducted for 1999-2009. In the first two experiments, MPAS-A is configured using global quasi-uniform grids at 120 km and 30 km grid spacing. In the other two experiments, MPAS-A is configured using variable-resolution (VR) mesh with local refinement at 30 km over North America and South America embedded inside a quasi-uniform domain at 120 km elsewhere. Precipitation and related fields in the four simulations are examined to determine how well the VR simulations reproduce the features simulated by the globally high-resolution model in the refined domain. In previous analyses of idealized aqua-planet simulations, the characteristics of the global high-resolution simulation in moist processes only developed near the boundary of the refined region. In contrast, the AMIP simulations with VR grids are able to reproduce the high-resolution characteristics across the refined domain, particularly in South America. This indicates the importance of finely resolved lower-boundary forcing such as topography and surface heterogeneity for the regional climate, and demonstrates the ability of the MPAS-A VR to replicate the large-scale moisture transport as simulated in the quasi-uniform high-resolution model. Outside of the refined domain, some upscale effects are detected through large-scale circulation but the overall climatic signals are not significant at regional scales. Our results provide support for the multi-resolution approach as a computationally efficient and physically consistent method for modeling regional climate.

Sakaguchi, Koichi↗
DOE OSTI2026

Bayesian Exploration and Surrogate Emulation of Nonlinear Beam-Response Geometry in the LBNF Beamline

Next-generation long-baseline neutrino experiments aim to achieve multi-MW proton beam power while reducing accelerator-induced systematic uncertainties. At Fermilab, the LBNF beamline is designed for 1.2 MW operation with PIP-II and is upgradeable to 2.4 MW. DUNE will probe the three-flavor neutrino paradigm and search for CP violation, requiring precise neutrino-flux normalization and improved control of accelerator-related uncertainties. Within the LBNF beamline, the System for On-Axis Neutrino Detection (SAND) will constrain flux uncertainties using precision near-detector measurements, while the Muon Monitor System (MuMS) will provide beamline diagnostics sensitive to the proton beam, target, and horn configuration. However, the pion phase space relevant for DUNE depends simultaneously on many correlated parameters, including beam centroid, beam width, horn current and alignment, target position, optics shifts, and radiation-induced changes. Consequently, MuMS observables exhibit nonlinear and coupled responses that are difficult to characterize using traditional one-parameter scans. To address this challenge, we are developing a Bayesian Exploration framework coupled to physics-informed surrogate emulators trained on Geant4 beamline simulations. Gaussian-process emulators provide both fast predictions and uncertainty estimates, enabling adaptive selection of new simulation points in beam-parameter space. As an initial demonstration, we construct surrogate emulators for MuMS response observables using a verified simulation campaign spanning proton-beam steering conditions. The emulators reproduce the simulated dependence of MuMS centroid and gradient observables while providing predictive uncertainties, and serve as the foundation for future multidimensional exploration including beam width, horn current, and additional beamline parameters. This work establishes a framework for uncertainty-aware beam monitoring, adaptive simulation campaigns, and rapid beam-response inference for future DUNE operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482↗
DOE OSTI2021

High-level fuel fabrication facility designs from discrete-event simulation

Like other industrial processes, the production of metallic nuclear fuels (MNF) requires that fabrication facilities be able to reliably meet production demands, operate efficiently, and adhere to federal and local safety regulations. In turn, the set of design variables employed by such a facility, such as operations policies, infrastructure and machinery purchased, and the type and number of staff hired, directly impact a facility’s ability to satisfy these goals. Therefore, facility designers must carefully determine which set of design variable values optimally satisfies these constraints. In this paper, we explore how values for these high-level design variables, namely hiring requirements, can be determined in the context of nuclear fuel manufacturing through the coupling of physics-based and discrete-event simulation technologies. Using the Versatile Test Reactor (VTR) program as a case study, we demonstrate how SCALE and MCNP nuclear physics model outputs can be integrated into ExtendSim discrete-event simulation (DES) models of the fuel fabrication process to determine the optimal number of staff hired to ensure fuel production goals are met, operations comply with effective dose limit regulations, and overall project costs are reduced.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

WE-DE-202-01: Connecting Nanoscale Physics to Initial DNA Damage Through Track Structure Simulations

Radiation therapy for the treatment of cancer has been established as a highly precise and effective way to eradicate a localized region of diseased tissue. To achieve further significant gains in the therapeutic ratio, we need to move towards biologically optimized treatment planning. To achieve this goal, we need to understand how the radiation-type dependent patterns of induced energy depositions within the cell (physics) connect via molecular, cellular and tissue reactions to treatment outcome such as tumor control and undesirable effects on normal tissue. Several computational biology approaches have been developed connecting physics to biology. Monte Carlo simulations are the most accurate method to calculate physical dose distributions at the nanometer scale, however simulations at the DNA scale are slow and repair processes are generally not simulated. Alternative models that rely on the random formation of individual DNA lesions within one or two turns of the DNA have been shown to reproduce the clusters of DNA lesions, including single strand breaks (SSBs), double strand breaks (DSBs) without the need for detailed track structure simulations. Efficient computational simulations of initial DNA damage induction facilitate computational modeling of DNA repair and other molecular and cellular processes. Mechanistic, multiscale models provide a useful conceptual framework to test biological hypotheses and help connect fundamental information about track structure and dosimetry at the sub-cellular level to dose-response effects on larger scales. In this symposium we will learn about the current state of the art of computational approaches estimating radiation damage at the cellular and sub-cellular scale. How can understanding the physics interactions at the DNA level be used to predict biological outcome? We will discuss if and how such calculations are relevant to advance our understanding of radiation damage and its repair, or, if the underlying biological processes are too complex for a mechanistic approach. Can computer simulations be used to guide future biological research? We will debate the feasibility of explaining biology from a physicists’ perspective. Learning Objectives: Understand the potential applications and limitations of computational methods for dose-response modeling at the molecular, cellular and tissue levels Learn about mechanism of action underlying the induction, repair and biological processing of damage to DNA and other constituents Understand how effects and processes at one biological scale impact on biological processes and outcomes on other scales J. Schuemann, NCI/NIH grantsS. McMahon, Funding: European Commission FP7 (grant EC FP7 MC-IOF-623630)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Investigating the effect of oxy-fuel combustion and light coal volatiles interaction: A mass spectrometric study

Given the multi-physical nature of coal combustion, the development and validation of detailed chemical models reproducing coal volatiles combustion under oxy-fuel conditions is a crucial step towards the advancement of predictive full-scale simulations. During the devolatilization process, a large variety of gases is released and undergoes secondary pyrolysis and oxidation reactions. Therefore, the ability to capture their interactions is a prerequisite for each chemical model used in its detailed or reduced form to simulate these processes. In this work, a high-resolution time-of-flight molecular-beam mass spectrometer was employed to enable fast and simultaneous detection of stable and unstable species in counterflow flames of typical light volatiles. Following an approach of increasing complexity, carbon dioxide and methane were progressively added to an argon diluted acetylene base flame. For the three flames investigated here, results showed a significant increase in the concentration of C 2 and C 3 hydrocarbons and oxygenated compounds caused by methane addition to the acetylene flame. By hindering the production of the butadienyl radical, the addition of methane induces the reduction of benzene which triggers the decrease of aromatic species. Conversely, CO 2 addition did not have significant effects on intermediates. To guide and interpret the measurements, numerical simulations with two existing chemical models were performed and the results were found to be consistent with the experimental data for small hydrocarbons. Some discrepancies were found between the two model predictions and between simulations and experiments for C 4 and C 5 species. In addition, numerical simulations were found to overestimate the role of the methyl radical in aromatics formation.

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

Active Interrogation using Photofission Technique for Nuclear Materials Control and Accountability

Innovative systems with increased sensitivity and resolution are in great demand to detect diversion and to prevent misuse in support of nuclear materials management for the U.S. fuel cycle. Nuclear fission is the most important multiplicative process involved in non-destructive active interrogation. This process produces the most easily recognizable signature for nuclear materials. In addition to thermal or high-energy neutrons, high-energy gamma rays can also excite a nucleus and cause fission through a process known as photofission. Electron linear accelerators (linacs) are widely used as the interrogating photon sources for inspection methods involving photofission technique. After photofission reactions, prompt signals are much stronger than the delayed signals, but it is difficult to quantify them in practical measurements. Delayed signals are easily distinguishable from the interrogating radiation. Linac-based, advanced inspection techniques utilizing the delayed signals after photofission have been extensively studied for homeland security applications. Previous research also showed that a unique delayed gamma ray energy spectrum exists for each fissionable isotope. In this work, high-energy delayed γ-rays were demonstrated to be signatures for detection, identification, and quantification of special nuclear materials. Such γ-rays were measured in between linac pulses using independent data acquisition systems. A list-mode system was developed to measure low-energy delayed γ-rays after irradiation. Photofission product yields of 238 U and 239 Pu were determined based on the measured delayed γ-ray spectra. The differential yields of delayed γ-rays were also proven to be able to discriminate nuclear from non-nuclear materials. The measurement outcomes were compared with Monte Carlo simulation results. It was demonstrated that the current available codes have capabilities and limitations in the simulation of photofission process. A two-fold approach was used to address the high-rate challenge in used nuclear fuel assay based on photofission technique. First, a standard HPGe preamplifier was modified to improve its capabilities in high-rate pulsed photofission environment. Second, advanced pulse processing algorithms were shown to greatly improve throughput rate without large sacrifice in energy resolution at ultra-high input count rate. Two customized gamma spectroscopy systems were also developed in real-time on FPGAs. They were shown to have promising performance matching available commercial units.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Success of digital adiabatic simulation with large Trotter step

The simulation of adiabatic evolution has deep connections with adiabatic quantum computation, the quantum approximate optimization algorithm, and adiabatic state preparation. Here we address the error analysis problem in quantum simulation of adiabatic process using Trotter formulas. Here we show that with additional conditions, the circuit depth can be linear in simulation time T. The improvement comes from the observation that the fidelity error here can't be estimated by the norm distance between evolution operators. This phenomenon is termed the robustness of discretization in digital adiabatic simulation. It can be explained in three steps, from analytical and numerical evidence: (1) The fidelity error should be estimated by applying adiabatic theorem on the effective Hamiltonian instead. (2) Because of the specialty of Riemann-Lebesgue lemma, most adiabatic process is naturally robust against discretization. (3) As the Trotter step gets larger, the spectral gap of effective Hamiltonian tends to close, which results in the failure of digital adiabatic simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Life-cycle Analysis of Bioproducts and Their Conventional Counterparts in GREET

To further expand upon the literature in this field and to develop a platform for bioproduct LCA, we developed LCA results for ten bioproducts produced either from algal glycerol or from corn stover-derived sugars. We used Argonne National Laboratory’s Greenhouse gases, Regulated Emissions, and Energy use in Transportation (GREETTM) model as the platform for this study. The data and calculations reported herein are available to GREET users in a bioproducts module included in the fall 2015 GREET release. This report documents our approach to this analysis and the results. In Chapter 2, we review the process we underwent to select the bioproducts for analysis based on market and technology readiness criteria. In Chapter 3, we review key parameters for production of the two feedstocks we considered: corn stover and algae. Given the lack of publicly available information about the production of bioproducts, which is caused in large part by the emerging nature of the industry, we developed Aspen Plus® simulations of the processes that could be used to produce each bioproduct. From these simulations, we extracted the energy and material flows of these processes, which were important inputs to the GREET bioproducts module. Chapter 4 provides the details of these Aspen Plus simulations. It is important to compare the LCA results for bioproducts to those for their petroleum counterparts. We therefore also developed material and energy flow data for conventional products based mostly on the literature. These data are described in Chapter 5 and are also included in the GREET bioproducts module. In Chapter 6, we present results from this analysis and examine areas for refinement and future research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

SU-F-J-187: The Statistical NTCP and TCP Models in the Proton Therapy

Purpose: The statistical models (SM) are typically used as a subjective description of a population for which there is only limited sample data, and especially in cases where the relationship between variables is known. The normal tissue complications and tumor control are frequently stochastic effects in the Radiotherapy (RT). Based on probabilistic treatments, it recently has been formulated new NTCP and TCP models for the RT. Investigating the particular requirements for their clinical use in the proton therapy (PT) is the goal of this work. Methods: The SM can be used as phenomenological or mechanistic models. The former way allows fitting real data and getting theirparameters. In the latter one, we should do efforts for determining the parameters through the acceptable estimations, measurements, and/or simulation experiments. Experimental methodologies for determination of the parameters have been developed from the fraction cells surviving the proton irradiation curves in tumor and OAR, and precise RBE models are used for calculating the variable of effective dose. As the executions of these methodologies have a high costs, so we have developed computer tools enable to perform simulation experiments as complement to limitations of the real ones. Results: The requirements for the use of the SM in the PT, such as validation and improvement of the elaborated and existent methodologies for determining the SM parameters and effective dose respectively, were determined. Conclusion: The SM realistically simulates the main processes in the PT, and for this reason these can be implemented in this therapy, which are simples, computable and they have other advantages over some current models. It has been determined some negative aspects for some currently used probabilistic models in the RT, like the LKB NTCP and others derived from logistic functions; which can be improved with the proposed methods in this study.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Force fields for monovalent and divalent metal cations in TIP3P water based on thermodynamic and kinetic properties

Metal cations are essential in many vital processes. In order to capture the role of different cations in all-atom molecular dynamics simulations of biological processes, an accurate parametrization is crucial. Here, we develop force field parameters for the metal cations Li + , Na + , K + , Cs + , Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ in combination with the TIP3P water model that is frequently used in biomolecular simulations. In progressing toward improved force fields, the approach presented here is an extension of previous efforts and allows us to simultaneously reproduce thermodynamic and kinetic properties of aqueous solutions. We systematically derive the parameters of the 12-6 Lennard-Jones potential which accurately reproduces the experimental solvation free energy, the activity derivative, and the characteristics of water exchange from the first hydration shell of the metal cations. In order to reproduce all experimental properties, a modification of the Lorentz-Berthelot combination rule is required for Mg 2+ . As a result, using a balanced set of solution properties, the optimized force field parameters aim to capture the fine differences between distinct metal cations including specific ion binding affinities and the kinetics of cation binding to biologically important anionic groups.

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

Application of Generic Disposal System Models

This report describes specific GDSA activities in fiscal year 2015 (FY2015) toward the development of the enhanced disposal system modeling and analysis capability for geologic disposal of nuclear waste. The GDSA framework employs the PFLOTRAN thermal-hydrologic-chemical multi-physics code (Hammond et al., 2011) and the Dakota uncertainty sampling and propagation code (Adams et al., 2013). Each code is designed for massively-parallel processing in a high-performance computing (HPC) environment. Multi-physics representations in PFLOTRAN are used to simulate various coupled processes including heat flow, fluid flow, waste dissolution, radionuclide release, radionuclide decay and ingrowth, precipitation and dissolution of secondary phases, and radionuclide transport through the engineered barriers and natural geologic barriers to a well location in an overlying or underlying aquifer. Dakota is used to generate sets of representative realizations and to analyze parameter sensitivity.

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

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗
DOE OSTI2022

Prediction of grain structure after thermomechanical processing of U-10Mo alloy using sensitivity analysis and machine learning surrogate model

Abstract Hot rolling and annealing are critical intermediate steps for controlling microstructures and thickness variations when fabricating uranium alloyed with 10% molybdenum (U-10Mo), which is highly relevant to worldwide nuclear non-proliferation efforts. This work proposes a machine-learning surrogate model combined with sensitivity analysis to identify and predict U-10Mo microstructure development during thermomechanical processing. Over 200 simulations were collected using physics-based microstructure models covering a wide range of thermomechanical processing routes and initial alloy grain features. Based on the sensitivity analysis, we determined that an increase in rolling reduction percentage at each processing pass has the strongest effect in reducing the grain size. Multi-pass rolling and annealing can significantly improve recrystallization regardless of the reduction percentage. With a volume fraction below 2%, uranium carbide particles were found to have marginal effects on the average grain size and distribution. The proposed stratified stacking ensemble surrogate predicts the U-10Mo grain size with a mean square error four times smaller than a standard single deep neural network. At the same time, with a significant speedup (1000×) compared to the physics-based model, the machine learning surrogate shows good potential for U-10Mo fabrication process optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Simulator for SUPO, a Benchmark Aqueous Homogeneous Reactor (AHR)

A simulator has been developed for SUPO (Super Power) an aqueous homogeneous reactor (AHR) that operated at Los Alamos National Laboratory (LANL) from 1951 to 1974. During that period SUPO accumulated approximately 600,000 kWh of operation. It is considered the benchmark for steady-state operation of an AHR. The SUPO simulator was developed using the process that resulted in a simulator for an accelerator-driven subcritical system, which has been previously reported.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

1D EDIPIC

EDIPIC codes were developed for simulation of kinetic effects in low-temperature plasmas. They can be used for simulation of plasmas in electrostatic approximation and include many atomistic and plasma-surface interaction processes needed for simulations of low temperature plasmas. These codes include a detailed model of secondary electron emission induced by electrons and ions. The Monte-Carlo model of electron-neutral collisions and the Langevin model of Coulomb collisions for electrons are also implemented.

Sydorenko, Dmytro↗