Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coding”

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

Physical implementation of a Majorana fermion surface code for fault-tolerant quantum computation

We propose a physical realization of a commuting Hamiltonian of interacting Majorana fermions realizing Z 2 topological order, using an array of Josephson-coupled topological superconductor islands. The required multi-body interaction Hamiltonian is naturally generated by a combination of charging energy induced quantum phase-slips on the superconducting islands and electron tunneling between islands. Our setup improves on a recent proposal for implementing a Majorana fermion surface code (Vijay et al 2015 Phys. Rev. X 5 041038), a 'hybrid' approach to fault-tolerant quantum computation that combines (1) the engineering of a stabilizer Hamiltonian with a topologically ordered ground state with (2) projective stabilizer measurements to implement error correction and a universal set of logical gates. Our hybrid strategy has advantages over the traditional surface code architecture in error suppression and single-step stabilizer measurements, and is widely applicable to implementing stabilizer codes for quantum computation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

HPC-driven computational reproducibility in numerical relativity codes: a use case study with IllinoisGRMHD

Abstract Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produceconsistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using theIllinoisGRMHDcode that is part of an open-source community software for simulation in relativistic astrophysics, theEinstein Toolkit. We verify that numerical results of simulating a single isolated neutron star withIllinoisGRMHDcan be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find thatIllinoisGRMHDreproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version ofIllinoisGRMHDreproduces these results once we account for bug fixes which have occurred since the original publication.

Astronomy & Astrophysics↗
DOE OSTI2024

The fastVFP code for solution of the Vlasov–Fokker–Planck equation

Abstract We describe the fast VFP code for solution of the Vlasov–Fokker–Planck equation for non-local electron transport and the generation of magnetic field, especially for application to laser-produced plasmas. We describe the essential features of the code that make it fast and robust and suitable for inclusion as a transport package in a fluid simulation. We present a few sample results that demonstrate the abilities of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The CLASSgal code for relativistic cosmological large scale structure

We present accurate and efficient computations of large scale structure observables, obtained with a modified version of the CLASS code which is made publicly available. This code includes all relativistic corrections and computes both the power spectrum C{sub ℓ}(z{sub 1},z{sub 2}) and the corresponding correlation function ξ(θ,z{sub 1},z{sub 2}) of the matter density and the galaxy number fluctuations in linear perturbation theory. For Gaussian initial perturbations, these quantities contain the full information encoded in the large scale matter distribution at the level of linear perturbation theory. We illustrate the usefulness of our code for cosmological parameter estimation through a few simple examples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Investigation of neutral particle dynamics in Aditya tokamak plasma with DEGAS2 code

Neutral particle behavior in Aditya tokamak, which has a circular poloidal ring limiter at one particular toroidal location, has been investigated using DEGAS2 code. The code is based on the calculation using Monte Carlo algorithms and is mainly used in tokamaks with divertor configuration. This code has been successfully implemented in Aditya tokamak with limiter configuration. The penetration of neutral hydrogen atom is studied with various atomic and molecular contributions and it is found that the maximum contribution comes from the dissociation processes. For the same, H α spectrum is also simulated which was matched with the experimental one. The dominant contribution around 64% comes from molecular dissociation processes and neutral particle is generated by those processes have energy of ~ 2.0 eV. Furthermore, the variation of neutral hydrogen density and H α emissivity profile are analysed for various edge temperature profiles and found that there is not much changes in H α emission at the plasma edge with the variation of edge temperature (7 to 40 eV).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

AR reactor neutronics multi-code validation*

The affordable, robust, compact (ARC) reactor is a tokamak fusion reactor concept currently under development by Commonwealth Fusion Systems and Massachusetts Institute of Technology. There are three important neutronics considerations for the operation of the ARC reactor: (1) breeding of enough tritium in the blanket to sustain the D–T reaction in the plasma; (2) ensuring low fluence on the superconducting toroidal field coils; and (3) assessing neutron volumetric heating in structural components. Here, this work aims to perform a validation of the neutronics analysis approach by code-to-code comparison. State-of-the-art software stacks are employed for the neutronics analysis of the ARC reactor, and a computer-aided design (CAD) model is used directly for Monte Carlo (MC) neutron transport calculations. Three software stacks, Attila-MCNP, OpenMC-DAGMC, and Shift-DAGMC, are used to perform neutronic analyses of a 90° sector CAD model of the ARC reactor. In conclusion, results show that the flux tallies calculated by the three software stacks are very close. Also, the volumetric heating and tritium breeding values have less than 0.6% relative difference between codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MAS: A versatile Landau-fluid eigenvalue code for plasma stability analysis in general geometry

We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number (n) and frequency (ω) in a broad range of interest in general tokamak geometry, based on a five-field Landau-fluid description of thermal plasmas. Beyond keeping the necessary plasma fluid response, we further retain the important kinetic effects including diamagnetic drift, ion finite Larmor radius, finite parallel electric field (E||), ion and electron Landau resonances in a self-consistent and non-perturbative manner without sacrificing the attractive efficiency in computation. The physical capabilities of the code are evaluated and examined in the aspects of both theory and simulation. In theory, the comprehensive Landau-fluid model implemented in MAS can be reduced to the well-known ideal MHD model, electrostatic ion-fluid model, and drift-kinetic model in various limits, which clearly delineates the physics validity regime. In simulation, MAS has been well benchmarked with theory and other gyrokinetic and kinetic-MHD hybrid codes in a manner of adopting the unified physical and numerical framework, which covers the kinetic Alfv\'en wave (KAW), ion sound wave (ISW), low-n kink, high-n ion temperature gradient mode (ITG) and kinetic ballooning mode (KBM). Moreover, MAS is successfully applied to model the Alfv\'en eigenmode (AE) activities in DIII-D discharge #159243, which faithfully captures the frequency sweeping of reversed shear Alfv\'en eigenmode (RSAE), the tunneling damping of toroidal Alfv\'en eigenmode (TAE), as well as the polarization characteristics of kinetic beta-induced Alfv\'en eigenmode (KBAE) and beta-induced Alfv\'en-acoustic eigenmode (BAAE) being consistent with former gyrokinetic theory and simulation. With respect to the key progress contributed to the community, MAS has the advantage of combining rich physics ingredients, realistic global geometry and high computation efficiency together for plasma stability analysis in linear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Analysis and design of fast flow liquid Li divertor for fusion nuclear science facility (FNSF) using coupled plasma boundary and LM MHD/heat transfer codes *

The SOLPS-ITER code is utilized to analyze the boundary plasma associated with a fast-flow lithium (Li) divertor configuration in the fusion nuclear science facility (FNSF) tokamak and identify operational regimes with acceptable divertor and core conditions. Plasma transport from the SOLPS-ITER code has been coupled with a liquid metal (LM) MHD/heat transfer code to model a Li open-surface divertor design and assess its impact on the scrape-off-layer (SOL) and core plasma performance. Simulations with only Neon (Ne) impurity seeding have been conducted to evaluate its impact on meeting FNSF design demands for the divertor and upstream plasma parameters. Simulation results indicate that Ne seeding significantly mitigates divertor heat flux but potentially reduces both upstream electron and main ion density due to fuel dilution. The combined application of Ne seeding and deuterium (D 2 ) puffing is required to satisfy the FNSF design requirements on upstream density ($n^{OMP}_{e,sep}$~1× 10 20 m -3 ) and divertor energy flux ($q^{Odiv}_{\bot,max}$<10 MW m -2 ). D 2 puffing plays a role in counteracting upstream density drops and augmenting energy and momentum losses through atomic and molecular processes. The inlet Li flow velocity is systematically varied across a wide range to identify acceptable flows and corresponding LM surface temperatures. This comprehensive analysis identifies the acceptable Li flow parameters, LM surface temperature, and emitted Li fluxes necessary to meet the major design constraints. The emitted Li fluxes exhibit minimal impact on the main plasma at surface temperatures up to approximately ~525 °C, corresponding emitted Li fluxes of up to φ Li ~2X10 23 atoms s -1 . Uncertainties in the Li emission processes from the surface are also investigated, primarily influencing Li loss in the lower surface temperature range (<525 °C), with simulation results indicating a minor impact on the divertor and upstream plasma. Conversely, evaporation predominantly drives the Li loss processes at higher surface temperature ranges (>525 °C), contaminating both the divertor and upstream plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Divertor heat load estimates on NSTX and DIII-D using new and open-source 2D inversion analysis code

A thermography inversion algorithm has been developed in the open-source Python-based computer code, HYPERION, to calculate the heat flux incident on plasma-facing components (PFCs) in axisymmetric tokamaks. The chosen mesh size at the surface significantly affects the calculated transient heat flux results. The calculated transient heat flux will exceed the real value when the mesh size tends to zero but will underestimate the real value when the mesh size is large. A criterion for determining the appropriate mesh size for the transient heat flux calculation will be discussed. The numerical scheme for HYPERION uses a 2D fully implicit finite-difference approach, allowing temperature-dependent thermal properties of PFC materials. The inversion algorithm is benchmarked against established heat flux calculation codes, TACO and THEODOR, based on thermography data from NSTX and DIII-D respectively. The primary benefits of HYPERION compared to TACO and THEODOR are that it is open-source and it allows for the optimization of mesh thickness along the substrate. The algorithm also accounts for the thermal properties of thin surface layers that characteristically form on PFCs due to plasma-material interactions. The agreement between HYPERION and THEODOR is excellent, as the percent difference between the codes is ~5% on average in the case of the DIII-D data for moderate to high heat flux. Verification tests with TACO show slightly higher average percent differences of 8% and 12%. In using HYPERION to study filaments in heat flux, the initial results indicate that small ELMs filaments significantly broaden the divertor heat flux, and decrease divertor peak flux. Compared to the inter-ELM, the small ELM filaments decrease the divertor peak surface temperature. With intermittent divertor filaments, the divertor heat flux width is comparable with that found in L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Transport and losses of energetic particles in tokamaks in the presence of Alfvén activity using the new full orbit TAPaS code coupled to FAR3d

Recent developments and tools integrated into the TAPaS code are presented, enabling realistic scenario simulations of particle dynamics within experimental tokamak magnetic equilibria. In particular, the enhanced capabilities of TAPaS enable seamless coupling with external simulations, provided the metric and equilibrium magnetic field of the external code are known. Coupling TAPaS with the gyro-fluid code FAR3d, the transport and losses of energetic particles in the presence Alfvén eigenmodes (AEs) in DIII-D plasma discharge #159243 were investigated. Detailed analyses of prompt losses with and without collisions were performed. Then, further analysis was performed in the presence of electromagnetic perturbations resulting from AEs activity. The results indicate that, for the energies and the initial conditions considered here, the presence of AEs enhances the particle losses.

Alfven eigenmodes↗
DOE OSTI2024

Aeroelastic code comparison using the IEA 22MW reference turbine

Reference wind turbine designs and the associated aeroelastic models are widely used in both research and industry. Reference models representing future concepts are of particular interest. Current state of the art aeroelastic tools are relied upon to design the next generation of large wind turbines. However, modelling assumptions may be invalidated by upcoming very large turbines, and different aeroelastic tools may give inconsistent results. A 22MW turbine model has been defined as part of International Energy Agency (IEA) Wind Task 55 on Reference Wind Turbines and Farms to represent future turbines to be deployed in the 2030s. In this study, an aeroelastic model of this turbine has been created in four tools; Bladed, HAWC2, OpenFAST, and QBlade. Code comparisons are presented for steady state operation, linear stability analysis, and time domain power production simulations in steady and turbulent wind. Generally, the codes show a good agreement, but with some differences present in the linear stability analysis, periodic azimuthal variation, and time domain simulations. The models are a good basis for further study with the IEA 22MW turbine, and further code comparison exercises.

17 WIND ENERGY↗
DOE OSTI2015

Large Eddy Simulation of Wind Turbine Wakes. Detailed Comparisons of Two Codes Focusing on Effects of Numerics and Subgrid Modeling

In this work we report on results from a detailed comparative numerical study from two Large Eddy Simulation (LES) codes using the Actuator Line Model (ALM). The study focuses on prediction of wind turbine wakes and their breakdown when subject to uniform inflow. Previous studies have shown relative insensitivity to subgrid modeling in the context of a finite-volume code. The present study uses the low dissipation pseudo-spectral LES code from Johns Hopkins University (LESGO) and the second-order, finite-volume OpenFOAMcode (SOWFA) from the National Renewable Energy Laboratory. When subject to uniform inflow, the loads on the blades are found to be unaffected by subgrid models or numerics, as expected. The turbulence in the wake and the location of transition to a turbulent state are affected by the subgrid-scale model and the numerics.

17 WIND ENERGY↗
DOE OSTI2024

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based translated sequence annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long-read sequencing data and in the context of pseudogenes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

A new moment-based general-relativistic neutrino-radiation transport code: Methods and first applications to neutron star mergers

ABSTRACT We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The cis-regulatory codes of response to combined heat and drought stress in Arabidopsis thaliana

Abstract Plants respond to their environment by dynamically modulating gene expression. A powerful approach for understanding how these responses are regulated is to integrate information about cis-regulatory elements (CREs) into models called cis-regulatory codes. Transcriptional response to combined stress is typically not the sum of the responses to the individual stresses. However, cis-regulatory codes underlying combined stress response have not been established. Here we modeled transcriptional response to single and combined heat and drought stress in Arabidopsis thaliana. We grouped genes by their pattern of response (independent, antagonistic and synergistic) and trained machine learning models to predict their response using putative CREs (pCREs) as features (median F-measure = 0.64). We then developed a deep learning approach to integrate additional omics information (sequence conservation, chromatin accessibility and histone modification) into our models, improving performance by 6.2%. While pCREs important for predicting independent and antagonistic responses tended to resemble binding motifs of transcription factors associated with heat and/or drought stress, important synergistic pCREs resembled binding motifs of transcription factors not known to be associated with stress. These findings demonstrate how in silico approaches can improve our understanding of the complex codes regulating response to combined stress and help us identify prime targets for future characterization.

Azodi, Christina B.↗
DOE OSTI2022

Why don't we share data and code? Perceived barriers and benefits to public archiving practices

The biological sciences community is increasingly recognizing the value of open, reproducible and transparent research practices for science and society at large. Despite this recognition, many researchers fail to share their data and code publicly. This pattern may arise from knowledge barriers about how to archive data and code, concerns about its reuse, and misaligned career incentives. Here, we define, categorize and discuss barriers to data and code sharing that are relevant to many research fields. We explore how real and perceived barriers might be overcome or reframed in the light of the benefits relative to costs. By elucidating these barriers and the contexts in which they arise, we can take steps to mitigate them and align our actions with the goals of open science, both as individual scientists and as a scientific community.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2023

Multispin Clifford codes for angular momentum errors in spin systems

The physical symmetries of a system play a central role in quantum error correction. Here, in this work we encode a qubit in a collection of systems with angular momentum symmetry (spins), extending the tools developed by Gross for single large spins. By considering large spins present in atomic systems and focusing on their collective symmetric subspace, we develop codes with octahedral symmetry capable of correcting errors up to second order in angular momentum operators. These errors include the most physically relevant noise sources such as microwave control errors and optical pumping. We additionally explore qubit codes that exhibit distance scaling commensurate with the surface code while permitting transversal single-qubit Clifford operations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗