Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “M CODES”

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,171 records · Page 65

Harmonic Generation and Inverse Cascade in the z-Pinch Driven, Preseeded Multimode, Magneto-Rayleigh-Taylor Instability

We report the magneto-Rayleigh-Taylor instability (MRTI) plays an essential role in astrophysical systems and in magneto-inertial fusion, where it is known to be an important degradation mechanism of confinement and target performance. In this Letter, we show for the first time experimental evidence of mode mixing and the onset of an inverse-cascade process resulting from the nonlinear coupling of two discrete preseeded axial modes (400- and 550 - μm wavelengths) on an Al liner that is magnetically imploded using the 20-MA, 100-ns rise-time Z Machine at Sandia National Laboratories. Four radiographs captured the temporal evolution of the MRTI. We introduce a novel unfold technique to analyze the experimental radiographs and compare the results to simulations and to a weakly nonlinear model. We find good quantitative agreement with simulations using the radiation magnetohydrodynamics code hydra. Spectral analysis of the MRTI time evolution obtained from the simulations shows evidence of harmonic generation, mode coupling, and the onset of an inverse-cascade process. The experiments provide a benchmark for future work on the MRTI and motivate the development of new analytical theories to better understand this instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-Energy Uncertainty Relation for Noisy Quantum Metrology

Detection of very weak forces and precise measurement of time are two of the many applications of quantum metrology to science and technology. To sense an unknown physical parameter, one prepares an initial state of a probe system, allows the probe to evolve as governed by a Hamiltonian 𝐻 for some time 𝑡, and then measures the probe. If 𝐻 is known, we can estimate 𝑡 by this method; if 𝑡 is known, we can estimate classical parameters on which 𝐻 depends. The accuracy of a quantum sensor can be limited by either intrinsic quantum noise or by noise arising from the interactions of the probe with its environment. In this work, we introduce and study a fundamental trade-off, which relates the amount by which noise reduces the accuracy of a quantum clock to the amount of information about the energy of the clock that leaks to the environment. Specifically, we consider an idealized scenario in which a party Alice prepares an initial pure state of the clock, allows the clock to evolve for a time that is not precisely known, and then transmits the clock through a noisy channel to a party Bob. Meanwhile, the environment (Eve) receives any information about the clock that is lost during transmission. We prove that Bob’s loss of quantum Fisher information about the elapsed time is equal to Eve’s gain of quantum Fisher information about a complementary energy parameter. We also prove a similar, but more general, trade-off that applies when Bob and Eve wish to estimate the values of parameters associated with two noncommuting observables. We derive the necessary and sufficient conditions for the accuracy of the clock to be unaffected by the noise, which form a subset of the Knill-Laflamme error-correction conditions. A state and its local time-evolution direction, if they satisfy these conditions, are said to form a metrological code. We provide a scheme to construct metrological codes in the stabilizer formalism. We show that there are metrological codes that cannot be written as a quantum error-correcting code with similar distance in which the Hamiltonian acts as a logical operator, potentially offering new schemes for constructing states that do not lose any sensitivity upon application of a noisy channel. We discuss applications of the trade-off relation to sensing using a quantum many-body probe subject to erasure or amplitude-damping noise.

metrology↗

First Experiments and Radiographs on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) Dense Plasma Focus

We report a dense plasma focus (DPF) is a relatively compact coaxial plasma gun, which completes its discharge as a Z-pinch. These devices are designed to operate at a variety of scales to produce short (<100 ns) pulses of ions, X-rays, and/or neutrons. LLNL recently constructed and brought into operation a new device, the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF, which is designed for radiography and high-yield operations. This device has been commissioned and has achieved neutron yields of up to 3.8E11 neutrons/pulse at 2.5-MA peak current while operating at up to 1 MJ of stored energy in its original pulsed power configuration. MJOLNIR is equipped with a wide range of diagnostics, including nuclear activation detectors, neutron time-of-flight (nToF) detectors, a fast-framing camera, optical light gates, and a time-gated neutron and X-ray imager. LLNL also runs unique particle-in-cell (PIC) simulations of DPF discharges in the Chicago code and has gained significant insight into the various physical factors that influence neutron yield. MJOLNIR is one of the first DPFs whose design and continual upgrades are heavily influenced by model predictions. In this article, we describe insights from modeling, device operation, and recent results. Comparisons between modeling predictions and measurements, as well as X-ray and neutron images are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Verification and Validation Activities for the Multi-Fidelity Toolkit

The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.

42 ENGINEERING↗

Technical Bases for Consequence Analyses Using MACCS (MELCOR Accident Consequence Code System)

The use of MELCOR Accident Consequence Code System (MACCS) in regulatory applications requires that input parameters used in MACCS models have a defensible and transparent technical basis. The existing technical basis for MACCS parameters was identified and updates were developed based on current data and information. The transparent way these updates are documented provides a technical basis for updating MACCS parameters with consistent application among the MACCS user community.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Verification and Validation of the PLTEMP/ANL Code for Thermal-Hydraulic Analysis of Experimental and Test Reactors, Volume 1

This is Volume 1 of a two-volume document that collects the verification and validation (V&V) works done for the PLTEMP/ANL code during the years of its development and improvement. Volume 1 describes the V&V of sixteen capabilities of the PLTEMP/ANL code that were identified by research reactor analysts as frequently used in their thermal-hydraulic analysis. Volume 2 describes the V&V of developments and improvements since the release of PLTEMP/ANL Version 4.3. Each chapter of the document focuses on verifying or validating a specific part of the software that calculates a specific phenomenon, e.g., channel flow calculation, coolant property calculation, heat transfer calculation, and flow instability calculation. Software verification is performed by comparing the code with hand calculation, Microsoft spreadsheet calculation, Mathematica calculation, or MATLAB calculation. The software validation is done by comparing the code with experimental data or a widely tested code like the RELAP5 code. In addition, some PLTEMP/ANL V&V works that are available in the open literature are simply cited in Volume 1 of the document. PLTEMP/ANL has been used in the safety analysis reports of several US and foreign research reactors licensed and converted from highly enriched uranium fuel to low-enriched uranium fuel. A list of such reactors is given in Volume 1 of the document.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development and Demonstration of a Risk Assessment Approach for Approval of a Transportation Package of a Transportable Nuclear Power Plant for Domestic Highway Shipment

The U.S. Department of Defense (DoD) Strategic Capabilities Office (SCO) has tasked PNNL to address the regulatory challenges associated with confirming the safe transport of Transportable Nuclear Power Plants (TNPPs) containing irradiated nuclear fuel. A previous report—Proposed Risk-Informed Regulatory Framework for Approval of Microreactor Transportation Packages (PNNL-31867)—determined that the expected radioactive inventory in the irradiated fuel of a TNPP would likely require shipment in an NRC-approved Type B package (or spent nuclear fuel cask) but that a TNPP “package” is unlikely to meet the entire suite of NRC requirements set forth in Part 71 of Title 10 of the Code of Federal Regulations (CFR) for a Type B package. It was therefore concluded that shipment of this initial TNPP transportation package, as well as possibly others, under existing regulations would likely require NRC approval using the 10 CFR 71.12 (“Specific exemptions”) process that relies on risk-informed decision making supported by quantitative risk assessment (i.e., Probabilistic Risk Assessment).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Electric Field Data from the North Slope of Alaska

File contains raw as well as calibrated electric field values (V/m) from two CS110 electric field meters. The file also contains raw 2-m U and V winds at the same site. The sampling rate is 1Hz. All files contain raw electric field data, as well as calibrated values to a one week ground-based calibration setup. During this week, simultaneous upward-facing ground based electric field measurements were taken alongside the two downward facing instruments on the pole (CS110-3 at 2m and CS110-2 at 5m). Slopes and intercepts were calculated and applied to the values influenced by the pole to convert the raw data to calibrated values. These calibrated values are absolute electric field measurements. There are no missing data codes. The missong data is not included in the files, and thus any files with missing data will be shorter than the normal 86,400 values. In July 2023, the CS110-3 was replaced with an identical instrument named CS110-8, and the CS110-2 was replaced with an identical instrument called CS110-9. All the variable names remained consistent and the same.

54 ENVIRONMENTAL SCIENCES↗

Microreactor Optimization Using Simulation And Economics (mouse)

Microreactor Optimization Using Simulation and Economics (MOUSE) is a tool that integrates both nuclear microreactor design and reactor economics to provide comprehensive evaluations and optimizations. This tool enables stakeholders to explore the interplay between technical and economic variables, guiding them towards effective and competitive microreactor solutions. For the reactor core simulations, MOUSE leverages the OpenMC Monte Carlo Particle Transport Code to perform detailed core simulations for various microreactor designs. The included OpenMC models are 2D core designs of a Liquid Metal Thermal Microreactor (LMTR), a Gas-Cooled TRISO-Fueled Microreactor (GCMR), and a Heat Pipe Microreactor. Beyond core design, MOUSE includes simplified calculations for: - Calculating the masses of heat exchangers within the system. - Mechanical power of pumps. - Estimating the area occupied by various buildings within the nuclear plant. For the economic analysis, MOUSE provides detailed bottom-up cost estimates, encompassing a wide range of costs including preconstruction costs, direct costs, indirect costs, training costs, financial costs, operation & maintenance (O&M) costs, and fuel costs. These cost estimations are developed using data from the MARVEL project and additional literature sources, enabling the calculation of total capital costs and levelized cost of energy for both first-of-a-kind and nth-of-a-kind microreactors. MOUSE also enables analysis of the cost drivers and competitiveness in the electricity market. MOUSE allows users to modify a wide array of technical and economic parameters to evaluate different scenarios and their impacts. Examples of these parameters include: Fuels, coolants, or reflector materials Enrichment levels Control drum materials and geometry Fuel pin geometry and materials Moderator pin geometry and materials Reactor core and reflector dimensions Packing factor for the TRISO particles Nuclear reactor power and reactor burnup Number of sensors Shielding thickness Reactor vessel and guard vessel dimensions Operational staff requirements Number of emergency shutdowns Levelization period Interest rate Construction duration Since MOUSE is powered by the WATTS toolkit, it supports optimization studies, parametric analyses, and uncertainty calculations/propagation. The optimization techniques enable users to identify optimal design and economic configurations. The parametric analysis tools allow users to explore the sensitivity of various parameters, while uncertainty propagation helps quantify the impact of uncertainties on overall performance and cost. User Interface and Workflow: Currently, MOUSE is a command-line-based tool. Users can input various reactor design or economic parameters, modify the designs, run simulations, and visualize results through comprehensive data visualization and reporting capabilities. The typical workflow involves setting up the reactor model, defining economic parameters, running simulations, and analyzing the results to make informed decisions. By combining advanced design calculations with detailed economic modeling, MOUSE provides a robust framework for optimizing nuclear microreactor technologies, enhancing their competitiveness, and guiding stakeholders towards innovative and cost-effective solutions.

Hanna, Botros [Idaho National Laboratory (INL), Id↗

Differential Gene Expression in Host Ubiquitination Processes in Childhood Malarial Anemia

Background: Malaria remains one of the leading global causes of childhood morbidity and mortality. In holoendemic Plasmodium falciparum transmission regions, such as western Kenya, severe malarial anemia [SMA, hemoglobin (Hb) < 6.0 g/dl] is the primary form of severe disease. Ubiquitination is essential for regulating intracellular processes involved in innate and adaptive immunity. Although dysregulation in ubiquitin molecular processes is central to the pathogenesis of multiple human diseases, the expression patterns of ubiquitination genes in SMA remain unexplored. Methods: To examine the role of the ubiquitination processes in pathogenesis of SMA, differential gene expression profiles were determined in Kenyan children ( n = 44, aged <48 mos) with either mild malarial anemia (M l MA; Hb ≥9.0 g/dl; n = 23) or SMA (Hb <6.0 g/dl; n = 21) using the Qiagen Human Ubiquitination Pathway RT 2 Profiler PCR Array containing a set of 84 human ubiquitination genes. Results: In children with SMA, 10 genes were down-regulated ( BRCC3 , FBXO3 , MARCH5 , RFWD2 , SMURF2 , UBA6 , UBE2A , UBE2D1 , UBE2L3 , UBR1 ), and five genes were up-regulated ( MDM2 , PARK2 , STUB1 , UBE2E3 , UBE2M ). Enrichment analyses revealed Ubiquitin-Proteasomal Proteolysis as the top disrupted process, along with altered sub-networks involved in proteasomal, protein, and ubiquitin-dependent catabolic processes. Conclusion: Collectively, these novel results show that protein coding genes of the ubiquitination processes are involved in the pathogenesis of SMA.

60 APPLIED LIFE SCIENCES↗

Non-local slicing approaches for NNLO QCD in MCFM

We present the implementation of several processes at Next-to-Next-to Leading Order (NNLO) accuracy in QCD in the parton-level Monte Carlo program MCFM. The processes treated are pp → H, W ± , Z, W ± H, ZH, W ± γ, Zγ and γγ and, for the first time in the code, W + W - , W ± Z and ZZ. Decays of the unstable bosons are fully included, resulting in a flexible fully differential Monte Carlo code. The NNLO corrections have been calculated using two non-local slicing approaches, isolating the doubly unresolved region by cutting on the zero-jettiness, $\mathscr{T}_0$, or on q T , the transverse momentum of the colour singlet final-state particles. We find that for most, but not all processes the q T slicing method leads to smaller power corrections for equal computational burden.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

hPIC2: A hardware-accelerated, hybrid particle-in-cell code for dynamic plasma-material interactions

The exascale era of high performance computing promises to bring the field of computational plasma physics ever closer to the goal of accurate multiscale modeling. Such computers will rely on hardware acceleration to offload work to dedicated components, notably general-purpose graphics processing units (GPUs). However, devices from different manufacturers require software to be written with different parallel programming models, greatly increasing the code maintenance burden of applications designed to perform on more than one such device. hPIC2 is a hybrid plasma simulation code developed with the Kokkos performance portability framework to target the architectures that will drive exascale computing for the foreseeable future. As a hybrid simulation code, hPIC2 investigates the simultaneous use of various plasma models on the same domain, at the same time. hPIC2 also optionally couples to RustBCA, which accurately models ion-material interactions using the binary collision approximation (BCA) method. In conclusion, hPIC2 therefore achieves scalable performance on a variety of computing architectures when simulating complex and diverse plasmas, particularly near plasma-material interfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature

Abstract Mountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.

54 ENVIRONMENTAL SCIENCES↗

Isotope impact on Alfvén eigenmodes and fast ion transport in DIII-D

Abstract Measurements of beam driven Alfvén Eigenmode (AE) activity in matched deuterium (D) and hydrogen (H) DIII-D plasmas show a dramatic difference in unstable mode activity and fast ion transport for a given injected beam power. The dependence of the unstable AE spectrum in reversed magnetic shear plasmas on beam and thermal species is investigated in the current ramp by varying beam power in a sequence of discharges for fixed thermal and beam species at fixed density. In general, a spectrum of Reversed Shear Alfvén Eigenmodes (RSAEs) and Toroidal Alfvén Eigenmodes (TAEs) are driven unstable with sub-Alfvénic D beam injection while primarily only RSAEs are driven unstable for the H beam cases investigated. Further, for a given beam power, the driven AE amplitude is always reduced with H beams relative to D and for H thermal plasma relative to pure D or mixed D/H plasmas. Estimates of the fast ion stored energy combined with modeling using the hybrid kinetic-MHD code MEGA indicate that the dominant mechanism contributing to the difference between H and D beam drive is the faster classical slowing down of H beam ions relative to D and the resultant lower beam ion pressure. Calculations of the AE induced stored energy deficits using the reduced critical gradient model TGLFEP show quantitative agreement with the observed dependencies on injected power, isotope and minimum safety factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A machine learning approach to galaxy properties: joint redshift–stellar mass probability distributions with Random Forest

We demonstrate that highly accurate joint redshift–stellar mass probability distribution functions (PDFs) can be obtained using the Random Forest (RF) machine learning (ML) algorithm, even with few photometric bands available. As an example, we use the Dark Energy Survey (DES), combined with the COSMOS2015 catalogue for redshifts and stellar masses. We build two ML models: one containing deep photometry in the griz bands, and the second reflecting the photometric scatter present in the main DES survey, with carefully constructed representative training data in each case. We validate our joint PDFs for 10 699 test galaxies by utilizing the copula probability integral transform and the Kendall distribution function, and their univariate counterparts to validate the marginals. Benchmarked against a basic set-up of the template-fitting code bagpipes, our ML-based method outperforms template fitting on all of our predefined performance metrics. In addition to accuracy, the RF is extremely fast, able to compute joint PDFs for a million galaxies in just under 6 min with consumer computer hardware. Such speed enables PDFs to be derived in real time within analysis codes, solving potential storage issues. As part of this work we have developed galpro 1, a highly intuitive and efficient python package to rapidly generate multivariate PDFs on-the-fly. galpro is documented and available for researchers to use in their cosmology and galaxy evolution studies.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensor-Assisted Fault Mitigation in Quantum Computation

We propose a method for assisting fault mitigation in quantum computation through the use of sensors, co-located near physical qubits. Specifically, we consider using transition edge sensors co-located on silicon substrates hosting superconducting qubits to monitor for energy injection from ionizing radiation that has been demonstrated to increase decoherence in transmon qubits. We generalize from these two physical device concepts and explore the potential advantages of co-located sensors to assist fault mitigation in quantum computation. In the simplest scheme, co-located sensors beneficially assist rejection of calculations potentially affected by environmental disturbances. Investigating the potential computational advantage further required development of an extension to the standard formulation of quantum error correction. In a specific case the standard three-qubit, bit flip quantum correction code, we show that given a 20% overall error probability per qubit, approximately 90% of repeated calculation attempts are correctable. However, when sensor-detectable errors account for 45% of overall error probability, the use of colocated sensors uniquely associated with independent qubits, boosts the fraction of correct final-state calculations to 96% at the cost of rejecting 7% of repeated calculation attempts.

Orrell, John L.↗

Encoding an Oscillator into Many Oscillators

An outstanding challenge for quantum information processing using bosonic systems is Gaussian errors such as excitation loss and added thermal noise errors. Thus, bosonic quantum error correction is essential. Most bosonic quantum error correction schemes encode a finite-dimensional logical qubit or qudit into noisy bosonic oscillator modes. In this case, however, the infinite-dimensional bosonic nature of the physical system is lost at the error-corrected logical level. On the other hand, there are several proposals for encoding an oscillator mode into many noisy oscillator modes. However, these oscillator-into-oscillators encoding schemes are in the class of Gaussian quantum error correction. Therefore, these codes cannot correct practically relevant Gaussian errors due to the established no-go theorems that state that Gaussian errors cannot be corrected by using only Gaussian resources. Here, we circumvent these no-go results and show that it is possible to correct Gaussian errors by using Gottesman-Kitaev-Preskill (GKP) states as non-Gaussian resources. In particular, we propose a non-Gaussian oscillator-into-oscillators code, namely the GKP two-mode squeezing code, and demonstrate that it can quadratically suppress additive Gaussian noise errors in both the position and momentum quadratures except for a small sublogarithmic correction. Furthermore, we demonstrate that our GKP two-mode squeezing code is near optimal in the weak noise limit by proving via quantum information theoretic tools that quadratic noise suppression is optimal when we use two physical oscillator modes. Lastly, in this work we show that our non-Gaussian oscillator encoding scheme can also be used to correct excitation loss and thermal noise errors, which are dominant error sources in many realistic bosonic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutronics Enhanced Meshing Operations (NEMO) Users Manual V.1.0

The is a Python framework developed at Idaho National Laboratory that uses the Cubit Python to build Cubit models and export them in EXODUS II format. These EXODUS mesh files can be used by any -based application, but NEMO is specialized for coupled neutronics calculations. Therefore, the primary purpose of NEMO is to prepare mesh files for reactor geometries and generate a homogenization overlay file for the Serpent Monte Carlo code. This homogenization overlay file ensures the consistency of the mesh with the generation of homogenized cross sections and reference fluxes for the Superhomogenization method. NEMO also supports the preparation of mesh files for thermal and thermo-mechanical models, particularly the discontinuous finite element gap heat transfer model in Griffin. This manual introduces the user to the basic capabilities in NEMO for prismatic geometries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗