Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code”

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,117 records · Page 62

A swapped genetic code prevents viral infections and gene transfer

Engineering the genetic code of an organism has been proposed to provide a firewall from natural ecosystems by preventing viral infections and gene transfer. However, numerous viruses and mobile genetic elements encode parts of the translational apparatus, potentially rendering a genetic-code-based firewall ineffective. Here we show that such mobile transfer RNAs (tRNAs) enable gene transfer and allow viral replication in Escherichia coli despite the genome-wide removal of 3 of the 64 codons and the previously essential cognate tRNA and release factor genes. We then establish a genetic firewall by discovering viral tRNAs that provide exceptionally efficient codon reassignment allowing us to develop cells bearing an amino acid-swapped genetic code that reassigns two of the six serine codons to leucine during translation. This amino acid-swapped genetic code renders cells resistant to viral infections by mistranslating viral proteomes and prevents the escape of synthetic genetic information by engineered reliance on serine codons to produce leucine-requiring proteins. As these cells may have a selective advantage over wild organisms due to virus resistance, we also repurpose a third codon to biocontain this virus-resistant host through dependence on an amino acid not found in nature. Furthermore, our results may provide the basis for a general strategy to make any organism safely resistant to all natural viruses and prevent genetic information flow into and out of genetically modified organisms.

59 BASIC BIOLOGICAL SCIENCES↗

Status of Mercury and Imp: Two Monte Carlo Transport Codes Developed Using Shared Infrastructure at Lawrence Livermore National Laboratory

The Monte Carlo Transport Project at Lawrence Livermore National Laboratory develops two Monte Carlo transport codes used in production by a sizable internal user community. Mercury is a Monte Carlo particle transport code used to model the interaction of neutrons, gammas, and light ions with a material. Imp is an implicit Monte Carlo thermal x-ray photon transport code used to model the interaction of x-ray photons with a material. This paper describes the two codes and highlights recent developments.

Pozulp, Michael↗

Experimental realization of the analogy of quantum dense coding in classical optics

We report on the experimental realization of the analogy of quantum dense coding in classical optical communication using classical optical correlations. Compared to quantum dense coding that uses pairs of photons entangled in polarization, we find that the proposed design exhibits many advantages. Considering that it is convenient to realize in optical communication, the attainable channel capacity in the experiment for dense coding can reach 2 bits, which is higher than that of the usual quantum coding capacity (1.585 bits). This increased channel capacity has been proven experimentally by transmitting ASCII characters in 12 quaternary digitals instead of the usual 24 bits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparison of four large-eddy simulation research codes and effects of model coefficient and inflow turbulence in actuator-line-based wind turbine modeling

Here, large-eddy simulation (LES) of a wind turbine under uniform inflow is performed using an actuator line model (ALM). Predictions from four LES research codes from the wind energy community are compared. The implementation of the ALM in all codes is similar and quantities along the blades are shown to match closely for all codes. The value of the Smagorinsky coefficient in the subgrid-scale turbulence model is shown to have a negligible effect on the time-averaged loads along the blades. Conversely, the breakdown location of the wake is strongly dependent on the Smagorinsky coefficient in uniform laminar inflow. Simulations are also performed using uniform mean velocity inflow with added homogeneous isotropic turbulence from a public database. The time-averaged loads along the blade do not depend on the inflow turbulence. Moreover, and in contrast to the uniform inflow cases, the Smagorinsky coefficient has a negligible effect on the wake profiles. It is concluded that for LES of wind turbines and wind farms using ALM, careful implementation and extensive cross-verification among codes can result in highly reproducible predictions. Moreover, the characteristics of the inflow turbulence appear to be more important than the details of the subgrid-scale modeling employed in the wake, at least for LES of wind energy applications at the resolutions tested in this work.

17 WIND ENERGY↗

Comparative collisionless alpha particle confinement in stellarator reactors with the XGC gyrokinetic code

Whole volume simulations of stellarators are necessary to address a number of important physics and engineering issues, including turbulent transport prediction and optimising divertor operation. In this work, advances in development of a whole volume stellarator capability for the global gyrokinetic particle-in-cell code XGC are detailed. A 3D interpolation of the equilibrium magnetic eld to the last close flux surface, calculated using the VMEC ideal magnetohydrodynamic equilibrum code, has been implemented. Orbit tracing with the XGC code is verified for the Wendelstein 7-X stellarator by comparison with the BEAMS3D and EUTERPE codes. The XGC stellarator version has then been used to consider collisionless alpha particle con nement in representative quasi-isodynamic and quasi-axisymmetric designs. It has been shown that, at least without further optimisation for energetic particle confinement, the Wendelstein 7-X-like quasi-isodynamic reactor design exhibits better alpha particle confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes

XGC (X-point Gyrokinetic Code) is a whole-volume, total-f gyrokinetic particle-in-cell code developed for modeling tokamaks. Recently, XGC has been extended to model more general 3D toroidal magnetic configurations, such as stellarators. These improvements have concluded in the XGC-S version. In this paper, XGC-S is benchmarked in the reduced delta-f limit for linear electrostatic ion temperature gradient-driven microinstabilities, which can underlie turbulent transport in stellarators. An initial benchmark of XGC-S in tokamak geometry shows good agreement with the XGC1, ORB5, and global GENE codes. A benchmark between XGC-S and the EUTERPE global gyrokinetic code for stellarators has also been performed, here in the geometry of the optimized stellarator Wendelstein 7-X. Good agreement has been discovered for the mode number spectrum, mode structure, and growth rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cross-code comparison of the impact of the fill tube on high yield implosions on the National Ignition Facility

Fill tubes are used to inject deuterium and tritium fuel into inertial confinement fusion capsules fielded on the National Ignition Facility. These fill tubes have been shown to have a detrimental effect on capsule performance, primarily by introducing a low-density pathway into the central fuel region that enables the jetting of ablation material into the hot spot. Due to the complexity of the highly nonlinear flow associated with the fill tube and the challenge of diagnosing the evolution of the fill tube jet late in the implosion experiments, the uncertainty in how this perturbation source evolves is great. In this work, we report on the results of a detailed code comparison performed to understand uncertainties in computational modeling of the impact of fill tubes on implosion performance. The study employed two radiation-hydrodynamics codes, HYDRA and xRAGE, which employ very different meshing strategies and hydrodynamics solvers, as well as two radiation transport methodologies, discrete ordinates and multi-group diffusion. Our results demonstrate generally good agreement between codes through most of the implosion although they indicate sensitivity to opacity averaging methods. Late in the implosion, differences arise in the distribution and amount of contaminant although these differences have a remarkably small impact on the amount of yield reduction due to the fill tube. While these results demonstrate sensitivity in fill tube modeling to algorithmic choices, the observed differences between codes are small relative to known sensitivities due to expected variations in the fill tube geometry. Finally, we have developed a methodology for performing multi-group diffusion simulations that show good agreement with the more accurate discrete ordinates method

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Simulation of an inductively coupled plasma with a two-dimensional Darwin particle-in-cell code

A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of the Subchannel Code CTF for Single- and Two-Phase Flows

As part of the Consortium for Advanced Simulation of Light Water Reactors, the subchannel code CTF is being used for single-phase and two-phase flow analysis under light water reactor operating conditions. Accurate determination of flow distribution, pressure drop, and void content is crucial for predicting margins to thermal crisis and ensuring more efficient plant performance. In preparation for the intended applications, CTF has been validated against data from experimental facilities comprising the General Electric (GE) 3 × 3 bundle, the boiling water reactor full-size fine-mesh bundle tests (BFBTs), the Risø tube, and the pressurized water reactor subchannel and bundle tests (PSBTs). Meanwhile, the licensed, well-recognized subchannel code VIPRE-01 was used to generate a baseline set of simulations for the targeted tests and solution parameters were compared to the CTF results.The flow split verification problem and single-phase GE 3 × 3 results are essentially in perfect agreement between the two codes. For the two-phase GE 3 × 3 cases, flow and quality discrepancies arise in the annular-mist flow regime, yet significant improvement is observed in CTF when void drift and two-phase turbulent mixing enhancement are considered. The BFBT pressure drop benchmark shows close agreement between predicted and measured results in general, although considerable overprediction by CTF is observed at relatively high void locations of the facility. This overestimation tendency is confirmed by the Risø cases. While overall statistics are satisfactory, both BFBT and PSBT bubbly-to-churn flow void contents are markedly overpredicted by CTF.The issues with two-phase closures such as turbulent mixing, interfacial and wall friction, and subcooled boiling heat transfer need to be addressed. As a result, preliminary sensitivity studies are presented herein, but more advanced models and code stability analysis require further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multigroup Cross-Section Generation with the OpenMC Monte Carlo Particle Transport Code

High-fidelity deterministic transport codes require highly accurate multigroup cross sections (MGXS). Monte Carlo is increasingly cited as a reactor-agnostic approach to MGXS generation since it is unconstrained by the engineering-based approximations that limit the applicability of deterministic MGXS generation tools. This paper introduces a new framework that uses the OpenMC Monte Carlo code to generate MGXS for use in multigroup transport codes. The openmc.mgxs module is built atop OpenMC's Python application programming interface to process tally data output by the OpenMC executable. This paper validates the module to generate MGXS that enable the multigroup OpenMOC transport code to compute eigenvalues to within 50 pcm and fission rates to within 1% of reference solutions for two heterogeneous pressurized water reactor benchmarks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structural Design and Modeling of MARVEL Primary Coolant System Using the ASME Section III, Division 5, Code

This paper presents the structural design and supporting analysis for the Microreactor Applications Research Validation and Evaluation (MARVEL) primary coolant system (PCS) using the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code Section III, Division 5, rules. MARVEL is a liquid metal–cooled microreactor intended to provide experimental capabilities for the rapid testing and development of microreactor technologies. The PCS utilizes high-temperature sodium-potassium liquid metal as the primary coolant and operates at a design temperature of 570°C, necessitating the consideration of creep-related failure mechanisms. The base metal for the PCS is 316H stainless steel, and the weldments are made with a 16-8-2 filler. The design approach incorporates the current base code rules along with ASME code cases N-924, N-861, and N-862 to address primary load, ratcheting, and creep-fatigue evaluations, respectively. The reactor’s operation involves complex thermal and mechanical interactions due to natural convective flow and differential thermal expansion between components. In conclusion, this paper discusses the structural engineering challenges encountered, such as managing thermal stresses in the distribution plenum and guard vessel, and outlines the strategies implemented to meet the code requirements, including design modifications and operational constraints.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

A Test Problem for Codes Solving the Discretized Diffusion Equation in Cartesian Geometry Derived Via Discrete Green’s Functions

We obtain here a solution to a zone-centered discretization of the one dimensional time-dependent diffusion equation with arbitrary initial conditions and source, constant absorption and scattering opacity, and constant zone size and time step. The solution is obtained using the discrete Green’s functions of the discretized equation. The solution of the discretized equation is useful in the testing of computer codes, because the code can be expected to agree with the solution to the discrete equation, to within small errors caused by roundoff. This is in contrast to solutions of the differential equation, with which the code results only approximately agree. Finally, the usefulness of the solution for tests of an inertial confinement fusion code is demonstrated.

97 MATHEMATICS AND COMPUTING↗

AN OPEN-SOURCE NEUTRINO RADIATION HYDRODYNAMICS CODE FOR CORE-COLLAPSE SUPERNOVAE

We present an open-source update to the spherically symmetric, general-relativistic hydrodynamics, core-collapse supernova (CCSN) code GR1D. The source code is available at http://www.GR1Dcode.org. We extend its capabilities to include a general-relativistic treatment of neutrino transport based on the moment formalisms of Shibata et al. and Cardall et al. We pay special attention to implementing and testing numerical methods and approximations that lessen the computational demand of the transport scheme by removing the need to invert large matrices. This is especially important for the implementation and development of moment-like transport methods in two and three dimensions. A critical component of neutrino transport calculations is the neutrino–matter interaction coefficients that describe the production, absorption, scattering, and annihilation of neutrinos. In this article we also describe our open-source neutrino interaction library NuLib (available at http://www.nulib.org). We believe that an open-source approach to describing these interactions is one of the major steps needed to progress toward robust models of CCSNe and robust predictions of the neutrino signal. We show, via comparisons to full Boltzmann neutrino-transport simulations of CCSNe, that our neutrino transport code performs remarkably well. Furthermore, we show that the methods and approximations we employ to increase efficiency do not decrease the fidelity of our results. We also test the ability of our general-relativistic transport code to model failed CCSNe by evolving a 40-solar-mass progenitor to the onset of collapse to a black hole.

79 ASTRONOMY AND ASTROPHYSICS↗

Majorana subsystem qubit codes that also correct odd-weight errors

Abstract A potential platform for topological quantum computation is the Majorana-based tetron architecture. Its building blocks are superconducting islands called tetrons, which host four Majorana zero modes. Existing error correcting codes can correct even-weight errors on tetrons. In a previous proposal by us, we had shown that incorporating tetrons in the stabilizer group allows us to correct a combination of odd-weight errors and even-weight errors on tetrons. In this work, we show that inclusion of tetrons in the gauge group lets us create subsystem codes from conventional Pauli stabilizer codes, which can correct both kinds of errors. Compared to the previous approach, the current approach lets us construct codes with fewer stabilizer generators. This leads to shorter fault-tolerant sequence length, and improves the fault-tolerant pseudothreshold by as much as 84%.

Physics↗

MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas

Abstract MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MultiPhATE2: code for functional annotation and comparison of phage genomes

To address a need for improved tools for annotation and comparative genomics of bacteriophage genomes, we developed multiPhATE2. As an extension of multiPhATE, a functional annotation code released previously, multiPhATE2 performs gene finding using multiple algorithms, compares the results of the algorithms, performs functional annotation of coding sequences, and incorporates additional search algorithms and databases to extend the search space of the original code. MultiPhATE2 performs gene matching among sets of closely related bacteriophage genomes, and uses multiprocessing to speed computations. MultiPhATE2 can be re-started at multiple points within the workflow to allow the user to examine intermediate results and adjust the subsequent computations accordingly. In addition, multiPhATE2 accommodates custom gene calls and sequence databases, again adding flexibility. MultiPhATE2 was implemented in Python 3.7 and runs as a command-line code under Linux or MAC operating systems. Full documentation is provided as a README file and a Wiki website.

59 BASIC BIOLOGICAL SCIENCES↗

Covariant Quantum Error-Correcting Codes with Metrological Entanglement Advantage

Here, we show that a subset of the basis for the irreducible representations of a tensor-product SU(2) rotation forms a covariant approximate quantum error-correcting code with transversal U(1) logical gates. Generalizing previous work on “thermodynamic codes” to general local spin and different irreducible representations using only properties of the angular momentum algebra, we obtain bounds on the code inaccuracy under generic noise on any known 𝑑 sites, under independent and identically distributed noise, and under heralded 𝑑-local erasures. We demonstrate that this family of codes protects a probe state with quantum Fisher information surpassing the standard quantum limit when the sensing parameter couples to the generator of the U(1) logical gate.

quantum error correction↗