Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation tools”

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 1,333 records · Page 74

Reinstated JET ICRF ILA: Overview and Results

The works undertaken to reinstate the JET ICRF ILA are reviewed. The vacuum matching capacitors were replaced, an extensive calibration of all the measurements in the RF circuit was carried out, new simulation tools were created and new control algorithms were implemented for the – toroidal and poloidal – phase control of the array as well as for the matching of the second stage. A review of the contribution of the reinstated ILA to the JET programme during the last campaigns is given showing namely that the new controls allowed extending the range of the operation to lower (29MHz) and higher (51MHz) frequencies than previously achieved and allowed more flexible and reliable operation. Operation with coupled power levels up to 2.8MW and voltages up to 40kV was achieved. ILA results on plasma are discussed and emphasis is given to the features of interest for ITER.

Dumortier, Pierre↗

Direct Measurement of ICRF-Enhanced Plasma Potentials on WEST Using Reciprocating Emissive Probes

An extensive documentation of ICRF-enhanced plasma potentials has been conducted over two experimental campaigns on the WEST tokamak using reciprocating emissive probes magnetically connected to two ICRF antennas. The collected data spans a wide range of antenna electrical settings (coupled power, toroidal phasing, left–right power balance) and plasma parameters (density at the antenna limiter above and below the lower hybrid resonance, plasma current, minority fraction). By scanning the edge safety factor across multiple probe plunges, the magnetic connection between the probe and the antenna varied, enabling the construction of a 2D map of the plasma and floating potentials around an active ICRF antenna. This dataset will be used to validate RF simulation tools equipped with the sheath boundary condition and used to predict RF rectified potentials and ICRF-induced impurity sputtering in future machines. This paper presents the diagnostic and some initial measurements, while the rest will be reported elsewhere.

Diab, Raymond [Massachusetts Inst. of Technology (↗

Visibility of two-dimensional layered materials on various substrates

For the investigation of 2D layered materials such as graphene, transition-metal dichalcogenides, boron nitride, and their heterostructures, dedicated substrates are required to enable unambiguous identification through optical microscopy. A systematic study is conducted, focusing on various 2D layered materials and substrates. The simulated colors are displayed and compared with microscopy images. Additionally, the issue of defining an appropriate index for measuring the degree of visibility is discussed. For a wide range of substrate stacks, layer thicknesses for optimum visibility are given along with the resulting sRGB colors. Further simulations of customized stacks can be conducted using our simulation tool, which is available for download and contains a database featuring a wide range of materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Precise parameterization of the recombination velocity at passivated phosphorus doped surfaces

We investigate the surface recombination velocity S{sub p} at the silicon-dielectric interface of phosphorus-doped surfaces for two industrially relevant passivation schemes for crystalline silicon solar cells. A broad range of surface dopant concentrations together with a high accuracy of evaluating the latter is achieved by incremental back-etching of the surface. The analysis of lifetime measurements and the simulation of the surface recombination consistently apply a set of well accepted models, namely, the Auger recombination by Richter et al. [Phys. Rev. B 86, 1–14 (2012)], the carrier mobility by Klaassen [Solid-State Electron. 35, 953–959 (1992); 35, 961–967 (1992)], the intrinsic carrier concentration for undoped silicon by Altermatt et al. [J. Appl. Phys. 93, 1598–1604 (2003)], and the band-gap narrowing by Schenk [J. Appl. Phys. 84, 3684–3695 (1998)]. The results show an increased S{sub p} at textured in respect to planar surfaces. The obtained parameterizations are applicable in modern simulation tools such as EDNA [K. R. McIntosh and P. P. Altermatt, in Proceedings of the 35th IEEE Photovoltaic Specialists Conference, Honolulu, Hawaii, USA (2010), pp. 1–6], PC1Dmod [Haug et al., Sol. Energy Mater. Sol. Cells 131, 30–36 (2014)], and Sentaurus Device [Synopsys, Sentaurus TCAD, Zürich, Switzerland] as well as in the analytical solution under the assumption of local charge neutrality by Cuevas et al. [IEEE Trans. Electron Devices 40, 1181–1183 (1993)].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

MPRAD: A Monte Carlo and ray-tracing code for the proton radiography in high-energy-density plasma experiments

Proton radiography is used in various high-energy-density (HED) plasma experiments. In this paper, we describe a Monte Carlo and raytracing simulation tool called multimegaelectronvolt proton radiography (MPRAD) that can be used for modeling the deflection of proton beams in arbitrary three dimensional electromagnetic fields as well as the diffusion of the proton beams by Coulomb scattering and stopping power. Specifically, the Coulomb scattering and stopping power models in cold matter and fully ionized plasma are combined using interpolation. We discuss the application of MPRAD in a few setups relevant to HED plasma experiments where the plasma density can play a role in diffusing the proton beams and affecting the prediction and interpretation of the proton images. It is shown how the diffusion due to plasma density can affect the resolution and dynamical range of the proton radiography.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effects of initial conditions and transport on ram pressure, Mach number, and uniformity for plasma liner formation and implosion

Plasma jet driven magneto-inertial fusion involves the merging of high Mach number plasma jets in a spherically converging geometry to form an imploding plasma liner. A three-dimensional smooth particle hydrodynamic simulation tool was used to study the number of plasma jets, jet density, and implosion velocity on the effects of formation and subsequent implosion of a plasma liner. The effects of different transport physics are presented. Thermal radiation and two-temperature modeling were found to have the most significant impact on the ram pressure. Solid angle averaged ram pressure and Mach number degradation were of particular interest in this study. Higher jet numbers were found to result in higher peak ram pressure. A spherical harmonic analysis was conducted for each case, providing a quantitative measure of the non-uniformity of the liner during the implosion process and how non-uniformities are affected by the different parameters explored in this study.

Schillo, Kevin↗

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmark between antenna code TOPICA, RAPLICASOL and Petra-M for the ICRH ITER antenna

ITER will be equipped with three plasma heating systems: neutral beam (NB), electron cyclotron (EC), and ion cy-clotron resonance heating (ICRH). The latter consists of two identical ICRH antennas to deliver 20 MW to the plasma (baseline, upgradable to 40 MW). ICRH will play a crucial role in the ignition and sustainment of burning plasmas in ITER. A high fidelity and robust modeling effort to understand the interaction of the IC waves with the scrape-off-layer (SOL) plasma is a very important aspect. Among the main important research topics, we have the assessment of the antenna loading for different plasma scenarios, the role of the lower hybrid resonance in front of the antenna and how to include it in our models, and the RF sheath boundary conditions to evaluate the antenna impurity generation. In this work, we tackle the first of these by reporting on ICRF simulations employing the Petra-M code, which is an electromagnetic simulation tool for modeling RF wave propagation based on MFEM [http://mfem.org] for the ITER ICRH antenna. Moreover, a benchmark between the well tested antenna codes TOPICA, RAPLI-CASOL, which is based on COMSOL [www.comsol.com], and the Petra-M code is also presented. In conclusion, S- and Z-matrices and wave electric field are compared showing an excellent agreement among these codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mixed-flow design for microfluidic printing of two-component polymer semiconductor systems

The rational creation of two-component conjugated polymer systems with high levels of phase purity in each component is challenging but crucial for realizing printed soft-matter electronics. Here, we report a mixed-flow microfluidic printing (MFMP) approach for two-component π -polymer systems that significantly elevates phase purity in bulk-heterojunction solar cells and thin-film transistors. MFMP integrates laminar and extensional flows using a specially microstructured shear blade, designed with fluid flow simulation tools to tune the flow patterns and induce shear, stretch, and pushout effects. This optimizes polymer conformation and semiconducting blend order as assessed by atomic force microscopy (AFM), transmission electron microscopy (TEM), grazing incidence wide-angle X-ray scattering (GIWAXS), resonant soft X-ray scattering (R-SoXS), photovoltaic response, and field effect mobility. For printed all-polymer (poly[(5,6-difluoro-2-octyl-2H-benzotriazole-4,7-diyl)-2,5-thiophenediyl[4,8-bis[5-(2-hexyldecyl)-2-thienyl]benzo[1,2-b:4,5-b′]dithiophene-2,6-diyl]-2,5-thiophenediyl]) [J51]:(poly{[N,N′-bis(2-octyldodecyl)naphthalene-1,4,5,8-bis(dicarboximide)-2,6-diyl]-alt-5,5′-(2,2′-bithiophene)}) [N2200]) solar cells, this approach enhances short-circuit currents and fill factors, with power conversion efficiency increasing from 5.20% for conventional blade coating to 7.80% for MFMP. Moreover, the performance of mixed polymer ambipolar [poly(3-hexylthiophene-2,5-diyl) (P3HT):N2200] and semiconducting:insulating polymer unipolar (N2200:polystyrene) transistors is similarly enhanced, underscoring versatility for two-component π -polymer systems. Mixed-flow designs offer modalities for achieving high-performance organic optoelectronics via innovative printing methodologies.

42 ENGINEERING↗

Multiphase Species Transport Modeling for Molten Salt Reactors in the System Analysis Module: Generation, Decay, Deposition, and Extraction of Insoluble Fission Products

With the increase of interests in the design and deployment of advanced reactor systems, a desire for simulation tools supporting system analysis of reactor operation and safety is rising. Molten salt reactors (MSRs), one of the advanced reactor systems, utilize liquid fused salt fuel as both coolant and fuel. During operation, MSR generates insoluble fission products, including noble metals and gases. The buildup of these species in fuel salt presents safety concerns as they may deposit on surfaces of critical components and produce excessive decay heat, causing the failure of system components. Timely removal of these noble metals and gases would ensure the safe operation of the reactor system. The dynamic nature of salt fuel system, involving the generation, decay, deposition, and extraction of noble metals and gases, calls for robust species transport models to facilitate system analysis and monitoring, and design of efficient species removal components. This paper concentrates on the development of a computational framework for species transport, consisting of multiphase transport model formulation, mass transfer between phases, numerical implementation in MOOSE environment, verification through Method of Manufacture Solutions (MMS) and validation against experimental data from the Molten Salt Reactor Experiment (MSRE). Integrating this framework into the System Analysis Module (SAM) code further enhances SAM’s capabilities for advanced reactor analysis in the future.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Enhancing Monte Carlo Workflows for Nuclear Reactor Analysis with Metamodel-Driven Modeling

Monte Carlo codes are essential components of many reactor physics simulation workflows as high-fidelity continuous-energy neutron transport solvers. Among Monte Carlo radiation transport codes, MCNP is particularly notable due to its diverse simulation capabilities, large user base, and long validation history. Despite being a powerful simulation tool, MCNP provides limited capabilities to allow automated execution, model transformation, or support for user-defined logic and abstractions that limit its compatibility with modern workflows. Here, to better integrate MCNP into a modern scientific workflow, we have developed an intuitive yet full-featured MCNP Application Program Interface (API) in Python, named MCNPy, which provides a specialized set of classes for MCNP input development. Moreover, to guarantee that our reading, writing, and modeling capabilities remain self-consistent (and to render the huge scope of the MCNP API manageable), we have adopted a strategy of model-driven software development in which a generalized model of the MCNP input format has been created. From this generalized model, or “metamodel,” problem-specific implementations such as an engine for input validation or a codebase for programmatic operations may be automatically generated. Since MCNPy primarily acts as a Python front-end to the underlying Java API that directly interfaces with the metamodel, it is intrinsically linked to the metamodel and thus remains maintainable. With MCNPy, users can programmatically read, write, and modify any syntactically valid MCNP input file regardless of its origin. These capabilities allow users to automate complicated tasks like design optimization and model translation for nuclear systems. As examples, this work demonstrates the use of MCNPy to find the critical radius of a plutonium sphere and to translate a 9000+ line MCNP input file into a corresponding OpenMC model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Bayesian Measurment Error Model for Misaligned Radiographic Data

An understanding of the inherent variability in micro-computed tomography (micro-CT) data is essential to tasks such as statistical process control and the validation of radiographic simulation tools. The data present unique challenges to variability analysis due to the relatively low resolution of radiographs, and also due to minor variations from run to run which can result in misalignment or magnification changes between repeated measurements of a sample. Positioning changes artificially inflate the variability of the data in ways that mask true physical phenomena. We present a novel Bayesian nonparametric regression model that incorporates both additive and multiplicative measurement error in addition to heteroscedasticity to address this problem. We also use this model to assess the effects of sample thickness and sample position on measurement variability for an aluminum specimen. Supplementary materials for this article are available online.

97 MATHEMATICS AND COMPUTING↗

Guest Editorial

The design of X-ray optical systems is very different from that of conventional optical systems, whether it be for the peculiarities of synchrotron radiation or the grazing incidence operation of most X-ray optical components. Further, specific simulation tools are thus required to ensure that designs work as expected, and there is now a vibrant landscape of open-source, freely available software that can be used to that end, and explore new design concepts.

47 OTHER INSTRUMENTATION↗

Sparse-Data Deep Learning Strategies for Radiographic Non-Destructive Testing

Radiography is an imaging technique used in a variety of applications, such as medical diagnosis, airport security, and nondestructive testing. We present a deep learning system for extracting information from radiographic images. We perform various prediction tasks using our system, including material classification and regression on the dimensions of a given object that is being radiographed. Our system is designed to address the sparse-data issue for radiographic nondestructive testing applications. It uses a radiographic simulation tool for synthetic data augmentation, and it uses transfer learning with a pre-trained convolutional neural network model. Using this system, our preliminary results indicate that the object geometry regression task saw an improvement of 70% in the R-squared value when using a multi-regime model. In addition, we increase the performance of the object material classification tasks by utilizing data from different imaging systems. In particular, using neutron imaging improved the material classification accuracy by 20% when compared to x-ray imaging.

convolutional neural networks↗

Calculating Measurement Uncertainty of the “Conventional Value of the Result of Weighing in Air”

The conventional value of the result of weighing in air is frequently used in commercial calibrations of balances. The guidance in OIML D-028 for reporting uncertainty of the conventional value is too terse. When calibrating mass standards at low measurement uncertainties, it is necessary to perform a buoyancy correction before reporting the result. When calculating the conventional result after calibrating true mass, the uncertainty due to calculating the conventional result is correlated with the buoyancy correction. We show through Monte Carlo simulations that the measurement uncertainty of the conventional result is less than the measurement uncertainty when reporting true mass. The Monte Carlo simulation tool is available in the online version of this article.

42 ENGINEERING↗

THREE-DIMENSIONAL RADIO AND X-RAY MODELING AND DATA ANALYSIS SOFTWARE: REVEALING FLARE COMPLEXITY

Many problems in solar physics require analysis of imaging data obtained in multiple wavelength domains with differing spatial resolution in a framework supplied by advanced three-dimensional (3D) physical models. To facilitate this goal, we have undertaken a major enhancement of our IDL-based simulation tools developed earlier for modeling microwave and X-ray emission. The enhanced software architecture allows the user to (1) import photospheric magnetic field maps and perform magnetic field extrapolations to generate 3D magnetic field models; (2) investigate the magnetic topology by interactively creating field lines and associated flux tubes; (3) populate the flux tubes with user-defined nonuniform thermal plasma and anisotropic, nonuniform, nonthermal electron distributions; (4) investigate the spatial and spectral properties of radio and X-ray emission calculated from the model; and (5) compare the model-derived images and spectra with observational data. The package integrates shared-object libraries containing fast gyrosynchrotron emission codes, IDL-based soft and hard X-ray codes, and potential and linear force-free field extrapolation routines. The package accepts user-defined radiation and magnetic field extrapolation plug-ins. We use this tool to analyze a relatively simple single-loop flare and use the model to constrain the magnetic 3D structure and spatial distribution of the fast electrons inside this loop. We iteratively compute multi-frequency microwave and multi-energy X-ray images from realistic magnetic flux tubes obtained from pre-flare extrapolations, and compare them with imaging data obtained by SDO, NoRH, and RHESSI. We use this event to illustrate the tool's use for the general interpretation of solar flares to address disparate problems in solar physics.

79 ASTRONOMY AND ASTROPHYSICS↗

LHC Forward Physics

The goal of this report is to give a comprehensive overview of the rich field of forward physics, with a special attention to the topics that can be studied at the LHC. The report starts presenting a selection of the Monte Carlo simulation tools currently available, chapter 2, then enters the rich phenomenology of QCD at low, chapter 3, and high, chapter 4, momentum transfer, while the unique scattering conditions of central exclusive production are analyzed in chapter 5. The last two experimental topics, Cosmic Ray and Heavy Ion physics are presented in the chapter 6 and 7 respectively. Chapter 8 is dedicated to the BFKL dynamics, multiparton interactions, and saturation. The report ends with an overview of the forward detectors at LHC. Each chapter is correlated with a comprehensive bibliography, attempting to provide to the interested reader with a wide opportunity for further studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

PIV/BOS synthetic image generation in variable density environments for error analysis and experiment design

In this work, we present an image generation methodology based on ray tracing that can be used to render realistic images of Particle Image Velocimetry (PIV) and Background Oriented Schlieren (BOS) experiments in the presence of density/refractive index gradients. This methodology enables the simulation of aero-thermodynamics experiments for experiment design, error, and uncertainty analysis. Images are generated by emanating light rays from the particles or dot pattern, and propagating them through the density gradient field and the optical elements, up to the camera sensor. The rendered images are realistic, and can replicate the features of a given experimental setup, like optical aberrations and perspective effects, which can be deliberately introduced for error analysis. We demonstrate this methodology by simulating a BOS experiment with a known density field obtained from direct numerical simulations (DNS) of homogeneous buoyancy driven turbulence, and comparing the light ray displacements from ray tracing to results from BOS theory. The light ray displacements show good agreement with the reference data. This methodology provides a framework for further development of simulation tools for use in experiment design and development of image analysis tools for PIV and BOS applications. Finally, an implementation of the proposed methodology in a Python-CUDA program is made available as an open source software for researchers.

47 OTHER INSTRUMENTATION↗