Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “V 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 757 records · Page 42

EMPIRE-PIC: A Performance Portable Unstructured Particle-in-Cell Code

In this study we introduce EMPIRE-PIC, a finite element method particle-in-cell (FEM-PIC) application developed at Sandia National Laboratories. The code has been developed in C++ using the Trilinos library and the Kokkos Performance Portability Framework to enable running on multiple modern compute architectures while only requiring maintenance of a single codebase. EMPIRE-PIC is capable of solving both electrostatic and electromagnetic problems in two- and three-dimensions to second-order accuracy in space and time. In this paper we validate the code against three benchmark problems — a simple electron orbit, an electrostatic Langmuir wave, and a transverse electromagnetic wave propagating through a plasma. We demonstrate the performance of EMPIRE-PIC on four different architectures: Intel Haswell CPUs, Intel's Xeon Phi Knights Landing, ARM Thunder-X2 CPUs, and NVIDIA Tesla V100 GPUs attached to IBM POWER9 processors. This analysis demonstrates scalability of the code up to more than two thousand GPUs, and greater than one hundred thousand CPUs.

97 MATHEMATICS AND COMPUTING↗

Stability of the Tilt Modes of an Actively Controlled Flywheel Analyzed

Applications of strongly gyroscopic rotors are becoming important, including flywheels for terrestrial and space energy storage and various attitude control devices for spacecraft. Some of these applications, especially the higher speed ones for energy storage, will have actively controlled magnetic bearings. These bearings will be required where speeds are too high for conventional bearings, where adequate lubrication is undesirable or impossible, or where bearing losses must be minimized for efficient energy storage. Flywheel rotors are highly gyroscopic, and above some speed that depends on the bandwidth of the feedback system, they always become unstable in an actively controlled magnetic bearing system. To assess ways to prevent instability until speeds well above the desired operating range, researchers at the NASA Lewis Research Center used a commercial controls code to calculate the eigenvalues of the tilt modes of a rigid gyroscopic rotor supported by active magnetic bearings. The real part of the eigenvalue is the negative of the damping of the mode, and the imaginary part is approximately equal to the mode s frequency.

Brown, Gerald V.↗

Large Scale Finite Element Modeling Using Scalable Parallel Processing

An iterative solver for use with finite element codes was developed for the Cray T3D massively parallel processor at the Jet Propulsion Laboratory. Finite element modeling is useful for simulating scattered or radiated electromagnetic fields from complex three-dimensional objects with geometry variations smaller than an electrical wavelength.

finite element modeling parallel processing iterat↗

Need for advanced research reactors for the next-generation reactor physics, analysis tools, and technology

Full text of publication follows. There is an urgent need for design and deployment of advanced research and test reactors in support of design, licensing and operation of advanced power reactors and education of next generation nuclear workforce. Existing research reactors mostly were designed and constructed decades ago with the main objectives of training operators, performing reactor physics experiments, and educating nuclear engineers and scientists. There are already gaps and significant concern about future capabilities for the existing research reactor facilities to address modern instrumentation and/or flexible environments for: performing reactor physics studies for advanced designs which have significantly different core materials forms and compositions, reactor shapes and size; validation of advanced high-fidelity software; development of machine learning algorithms for enhancement of human-machine collaboration in support of reactor monitoring, operation and safeguards; and, effective education of the next-generation workforce. The authors will focus on the need for advanced research reactors to improve and validate fast and accurate simulation tools for high-fidelity modeling and analysis of nuclear reactors in support of their design, optimization, licensing, operation, and monitoring. In the past, simulation tools were limited to relatively coarse models using approximate methodologies that benefited from two main factors: i) allowance for large margins and tolerances; ii) ability to construct prototype (e.g., zero power) reactors for adjustment of approximate methodologies. The next generation reactors have to be designed mainly by using novel high-fidelity computational tools that are accurate and fast, and therefore can be used for parametric studies and uncertainty quantification. To sufficiently demonstrate the accuracy of these tools, advanced research reactors are needed. The authors argue the need for new computational paradigms such as the MRT (Multistage, Response- function Transport) methodology which has resulted in the development of the novel high-fidelity RAPID (Real-time Analysis for Particle-transport and In-situ Detection) code system. Such code systems have to be robust in modeling any complex system, and should be fast and accurate, henceforth their uncertainties can be quantified at reasonable costs. Again, advanced research reactors are needed for the validation of the fidelity and accuracy of new computational tools. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Metallic fuel transient fuel-cladding interface liquefaction model assessment platform enabled by integrating BISON with databases

A novel platform has been developed within the BISON fuel performance code to assess models of fuel-cladding interface liquefaction for sodium-cooled fast reactor (SFR) metallic fuels. Here, this platform is crucial because liquefaction at the fuel-cladding interface significantly impacts fuel performance and may compromise fuel pin integrity during transient events. To ensure accurate predictions, the platform integrates data collected during the Integral Fast Reactor (IFR) program, now archived in metallic fuel databases. This integration supports verification and validation (V&V) of the models in BISON. Leveraging the extensive US experience with metallic fuel liquefaction and the collections of preserved legacy data, the platform serves as a powerful tool for evaluating existing models and advancing the development of new ones.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing topological spin liquids on a programmable quantum simulator

Synthesizing topological order Topologically ordered matter exhibits long-range quantum entanglement. However, measuring this entanglement in real materials is extremely tricky. Now, two groups take a different approach and turn to synthetic systems to engineer the topological order of the so-called toric code type (see the Perspective by Bartlett). Satzinger et al . used a quantum processor to study the ground state and excitations of the toric code. Semeghini et al . detected signatures of a toric code–type quantum spin liquid in a two-dimensional array of Rydberg atoms held in optical tweezers. —JS

Science & Technology - Other Topics↗

Corrosion of SS 316, Inconel 617, Inconel 625 in molten chloride salt with emphasis on the speciation of Cr, Ni, and Co in molten salt matrix by XAFS

Dr. Hartmann’ report summarizes the corrosion behavior of materials used in the construction and operation of molten salt reactors which is is of great interest to the molten salt reactor (MSR) community. In the present study we intent to decipher the very elemental phenomena on the corrosion of austenitic chromium-containing stainless steel SS316, as well as Inconel 625 and ASME-coded high temperature Alloy (Inconel) 617. We initiated and are in continuation of conducting ampoule type corrosion testing at 550 °C and 650 °C to determine chromium diffusivity and the related activation energy as well as the speciation of the corrosion products.

36 MATERIALS SCIENCE↗

Negative streamer development in FEP teflon

A computational model is developed which describes the evolution and propagation of an ionizing front (negative streamer) in solid materials. The ionization front consists of drifting avalanching electrons moving self-consistently under the influence of their own space-charge field together with an applied external field. The required input information for the model consists of the functional dependence of the macroscopic transport coefficients on the local electric field, the initial conditions for beginning the calculation, and the strength of the applied field. A computational approach for specifying the transport coefficients and initional conditions is also described. The approach has been implemented by constructing three computer codes which sequentially interface, beginning with single electron scattering, and ending with streamer development. Computational results are presented for model calculations in Teflon. The overall model is perceived to provide a picture of the initiation phase of a propagating discharge in electron-irradiated dielectrics.

Beers, B. L.↗

Turbofan forced mixer-nozzle internal flowfield. Volume 2: Computational fluid dynamic predictions

A general program was conducted to develop and assess a computational method for predicting the flow properties in a turbofan forced mixed duct. The detail assessment of the resulting computer code is presented. It was found that the code provided excellent predictions of the kinematics of the mixing process throughout the entire length of the mixer nozzle. The thermal mixing process between the hot core and cold fan flows was found to be well represented in the low speed portion of the flowfield.

Werle, M. J.↗

An algebraic turbulence model for three-dimensional viscous flows

An algebraic turbulence model is proposed for use with three-dimensional Navier-Stokes analyses. It incorporates features of both the Baldwin-Lomax and Cebeci-Smith models. The Baldwin-Lomax model uses the maximum of a function f(y) to determine length and velocity scales. An analysis of the Baldwin-Lomax model shows that f(y) can have a spurious maximum close to the wall, causing numerical problems and non-physical results. The proposed model uses integral relations to determine delta(*) u(sub e) and delta used in the Cebeci-Smith mode. It eliminates a constant in the Baldwin-Lomax model and determines the two remaining constants by comparison to the Cebeci-Smith formulation. Pressure gradient effects, a new wake model, and the implementation of these features in a three-dimensional Navier-Stokes code are also described. Results are shown for a flat plate boundary layer, an annular turbine cascade, and endwall heat transfer in a linear turbine cascade. The heat transfer results agree well with experimental data which shows large variations in endwall Stanton number contours with Reynolds number.

Chima, R. V.↗

Effect of convective transport in edge and scrape-off-layer plasmas of ADITYA-U tokamak

The 2-D edge plasma fluid transport code, UEDGE has been used to simulate the edge and scrape-off layer regions of circular limiter plasmas of ADITYA-U. The computational mesh defining the limiter geometry of ADITYA-U has been introduced in the UEDGE code through an in-house developed routine. It has successfully been integrated with the UEDGE code and the measured radial profile of electron density (n e ) in the edge region is modelled using the simulations. The plasma discharge for the present study has plasma current I P ~ 130 kA, duration ~ 150 ms, chord-averaged electron density ~ 1.3 x 10 19 m -3 , and central electron temperature ~ 270 eV has been observed. The maximum toroidal magnetic field is ~ 1.2 T. It has been found that an inward convective velocity, v conv ~ 1.5 m/s is required in addition to the constant perpendicular diffusion coefficient, D ⟂ ~ 0.2 m 2 /s to adequately match the measured n e profile in typical discharges of ADITYA-U. The value of D ⟂ ~ 0.2 m 2 /s is found to lie between the estimated neoclassical diffusivity and Bohm diffusivity.

Physics - Plasma physics↗

Ballistic limit regression analysis for Space Station Freedom meteoroid and space debris protection system

Relationships defining the ballistic limit of Space Station Freedom's (SSF) dual wall protection systems have been determined. These functions were regressed from empirical data found in Marshall Space Flight Center's (MSFC) Hypervelocity Impact Testing Summary (HITS) for the velocity range between three and seven kilometers per second. A stepwise linear least squares regression was used to determine the coefficients of several expressions that define a ballistic limit surface. Using statistical significance indicators and graphical comparisons to other limit curves, a final set of expressions is recommended for potential use in Probability of No Critical Flaw (PNCF) calculations for Space Station. The three equations listed below represent the mean curves for normal, 45 degree, and 65 degree obliquity ballistic limits, respectively, for a dual wall protection system consisting of a thin 6061-T6 aluminum bumper spaced 4.0 inches from a .125 inches thick 2219-T87 rear wall with multiple layer thermal insulation installed between the two walls. Normal obliquity is d(sub c) = 1.0514 v(exp 0.2983 t(sub 1)(exp 0.5228). Forty-five degree obliquity is d(sub c) = 0.8591 v(exp 0.0428) t(sub 1)(exp 0.2063). Sixty-five degree obliquity is d(sub c) = 0.2824 v(exp 0.1986) t(sub 1)(exp -0.3874). Plots of these curves are provided. A sensitivity study on the effects of using these new equations in the probability of no critical flaw analysis indicated a negligible increase in the performance of the dual wall protection system for SSF over the current baseline. The magnitude of the increase was 0.17 percent over 25 years on the MB-7 configuration run with the Bumper II program code.

Jolly, William H.↗

Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave

Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.

Cambier, Jean-Luc↗

NuGraph2 with context-aware inputs: physics-inspired improvements in semantic segmentation

Graph neural networks have recently shown strong promise for event reconstruction tasks in Liquid Argon Time Projection Chambers, yet their performance remains limited for underrepresented classes of particles, such as Michel electrons. In this work, we investigate physics-informed strategies to improve semantic segmentation within the NuGraph2 architecture. We explore three complementary approaches: (i) enriching the input representation with context-aware features derived from detector geometry and track continuity, (ii) introducing auxiliary decoders to capture class-level correlations, and (iii) incorporating energy-based regularization terms motivated by Michel electron energy distributions. Experiments on MicroBooNE public datasets show that physics-inspired feature augmentation yields the largest gains, particularly boosting Michel electron precision and recall by disentangling overlapping latent space regions. In contrast, auxiliary decoders and energy-regularization terms provided limited improvements, partly due to the hit-level nature of NuGraph2, which lacks explicit particle- or event-level representations. Our findings highlight that embedding physics context directly into node-level inputs is more effective than imposing task-specific auxiliary losses, and suggest that future hierarchical architectures such as NuGraph3, with explicit particle- and event-level reasoning, will provide a more natural setting for advanced decoders and physics-based regularization. The code for this work is publicly available on Github at https://github.com/vitorgrizzi/nugraph_phys/tree/main_phys.

Other Experiments↗

Metastable states from multinucleon excitations in 202 Tl and 203 Pb

The excited level structures of 202 Tl and 203 Pb, above the 7 + and 29/2 – isomers, respectively, have been studied. An isomer with I π = 20 + and T 1/2 = 215(10) μs has been established in 202 Tl, and the level scheme extended from I = 10 to 20ℏ with the placement of fifteen new transitions. In 203 Pb, the I π = 37/2 + state is established to be metastable, with T 1/2 = 2.5(3) ns. Levels in both nuclei arise from intrinsic excitations, with likely particle-hole character for the higher-lying states in 203 Pb. The 20 + isomer in 202 Tl is most likely associated with a $πh$ $^{–1}_{11/2}$ Ⓧ ν($i$ $^{–2}_{13/2}$, $f$ $^{–1}_{5/2}$) configuration, while the 37/2 + state in 203 Pb results from the excitation of five neutrons. Furthermore, calculations, using both an empirical approach and the oxbash code, have been performed to aid in the description of the excited level structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Effect of upstream sidewall boundary layer removal on an airfoil test

Sidewall boundary layer effects were investigated by applying partial upstream sidewall boundary layer removal in the Langley 0.3-m transonic cryogenic tunnel. Over the range of sidewall boundary layer displacement thickness of these tests the influence on pressure distribution was found to be small for subcritical conditions; however, for supercritical conditions the shock position was affected by the sidewall boundary layer. For these tests (with and without boundary layer remove) comparisons with predictions of the GRUMFOIL computer code indicated that Mach number corrections due to the sidewall boundary layer improve the agreement for both subcritical and supercritical conditions. The results also show that sidewall boundary layer removal reduces the magnitude of the sidewall correction; however, a suitable correction must still be made.

Johnson, C. B.↗

Application of a full-potential solver to bending-torsion flutter in cascades

A two-dimensional, unsteady, full-potential cascade aerodynamic code has been used for flutter calculations. The code, which uses a time-marching implicit scheme with internal Newton-iterations at each time step, allows oscillatory blade motions with phase lag between adjacent blades to be modeled using a single inter-blade passage. This aerodynamic code has been used with a two-degrees-of-freedom typical section structural model to determine the aeroelastic stability of a bladed disk. The unsteady pressure results calculated from the aerodynamic code have been verified by comparison with results from the linear and linearized potential theories. Flutter results have been presented for two cases. A comparison with linear theory shows identical trends; however, the present calculations show a lower flutter frequency.

Bakhle, Milind A.↗