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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FEgen (v.1): Field Emission Distribution Generator Freeware Based on Fowler–Nordheim Equation

This work presents an open-source code, FEgen, that generates initial field emission distributions based on the time-dependent Fowler–Nordheim equation. Current in-demand and open-source beam tracking software—such as ASTRA, IMPACT-T, and GPT—do not contain native distribution generator functions suitable for field emitted electron tracking and phase space analysis. Furthermore, the intuitive graphical user interface of FEgen generates input files shown to be fully compatible with both ASTRA and GPT. Here, detailed application examples demonstrating a subset of FEgen’s capabilities are given by creating a customized emitter pattern with application to transverse beam shaping and longitudinal beam dynamic experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Approximate heating analysis for the windward-symmetry plane of Shuttle-like bodies at large angle of attack

An engineering method has been developed for computing the windward-symmetry plane convective heat-transfer rates on Shuttle-like vehicles at large angles of attack. The engineering code includes an approximate inviscid flowfield technique, laminar and turbulent heating-rate expressions, an approximation to account for the variable-entropy effects on the surface heating and the concept of an equivalent axisymmetric body to model the windward-ray flowfields of Shuttle-like vehicles at angles of attack from 25 to 45 degrees. The engineering method is validated by comparing computed heating results with corresponding experimental data measured on Shuttle and advanced transportation models over a wide range of flow conditions and angles of attack from 25 to 40 degrees and also with results of existing prediction techniques. The comparisons are in good agreement.

Zoby, E. V.↗
NASA NTRS1991

Numerical prediction of turbulent oscillating flow and associated heat transfer

A crucial point for further development of engines is the optimization of its heat exchangers which operate under oscillatory flow conditions. It has been found that the most important thermodynamic uncertainties in the Stirling engine designs for space power are in the heat transfer between gas and metal in all engine components and in the pressure drop across the heat exchanger components. So far, performance codes cannot predict the power output of a Stirling engine reasonably enough if used for a wide variety of engines. Thus, there is a strong need for better performance codes. However, a performance code is not concerned with the details of the flow. This information must be provided externally. While analytical relationships exist for laminar oscillating flow, there has been hardly any information about transitional and turbulent oscillating flow, which could be introduced into the performance codes. In 1986, a survey by Seume and Simon revealed that most Stirling engine heat exchangers operate in the transitional and turbulent regime. Consequently, research has since focused on the unresolved issue of transitional and turbulent oscillating flow and heat transfer. Since 1988, the University of Minnesota oscillating flow facility has obtained experimental data about transitional and turbulent oscillating flow. However, since the experiments in this field are extremely difficult, lengthy, and expensive, it is advantageous to numerically simulate the flow and heat transfer accurately from first principles. Work done at the University of Minnesota on the development of such a numerical simulation is summarized.

Koehler, W. J.↗
DOE OSTI2021

NSLS-II longitudinal impedance budget

From the beginning of the NSLS-II project, the vacuum chamber of the NSLS-II storage ring has been optimized to reduce its impedance. The optimization, performed in close collaboration with the Diagnostic, Vacuum, RF, and Mechanical Engineering groups, has been achieved via calculating the impedance of the vacuum chamber components with the use of both analytical and numerical tools, with the numerical computations carried out with the GdfidL code. Fore reliable NSLS-II operation with the total beam current of 500 mA, each component of the vacuum chamber has been individually optimized to minimize the beam-induced heating and to prevent single-bunch and coupled-bunch instabilities. In this study, the impedance calculation and optimization know-how are summarized for all the chamber components, and the total longitudinal impedance budget is analyzed. As a result of this extensive research and engineering effort, the design beam current of 500 mA has been recently accumulated in NSLS-II without any major problem caused by collective effects.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Investigation of compression of indirect-drive targets under conditions of the NIF facility using one-dimensional modelling

The indirect compression dynamics of targets containing capsules with ablators of a plastic, high-density carbon and beryllium is simulated in the framework of a one-dimensional model based on the 1D RADIAN code. Experiments with such targets are performed on the NIF facility in the Livermore Laboratory (USA) in 2014 – 2018. The 1D simulation data are consistent with the results of experiments and calculations made at this Laboratory. The effect of the hard part of hohlraum radiation on capsule compression parameters is confirmed. We demonstrate the feasibility of eliminating this influence not only by selecting the hohlraum material but also by introducing admixtures into the capsule that absorb this radiation. It is shown how varying the amount of admixture in the capsule ablator varies the spectrum of the radiation that heats the DT fuel. (interaction of laser radiation with matter. laser plasma)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigating the impact of intermediate-mode perturbations on diagnosing plasma conditions in DT cryogenic implosions via synthetic x-ray Thomson scattering

The pursuit of inertial confinement fusion ignition target designs requires precise experimental validation of the conditions within imploding capsules, in particular the density and temperature of the compressed shell. Previous work has identified x-ray Thomson scattering (XRTS) as a viable diagnostic tool for inferring the in-flight compressed deuterium-tritium shell conditions during capsule implosions (Poole et al 2022 Phys. Plasmas 29 072703). However, this study focused on one-dimensional simulations, which do not account for the growth of hydrodynamic instabilities. In this work, two-dimensional DRACO simulations incorporating intermediate-mode perturbations up to Legendre mode l = 50 were used to generate synthetic XRTS spectra with the SPECT3D code. The analysis employed Markov-Chain Monte Carlo techniques to infer plasma conditions from these spectra. The results demonstrate that the XRTS diagnostic platform can effectively discern the in-flight compressed shell conditions for targets with varying adiabats, even in the presence of intermediate-mode perturbations. This work underscores the potential of XRTS for realistic inertial confinement fusion experiments, providing a robust method for probing the complex dynamics of fusion implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Development of a three-dimensional Navier-Stokes code on CDC star-100 computer

A three-dimensional code in body-fitted coordinates was developed using MacCormack's algorithm. The code is structured to be compatible with any general configuration, provided that the metric coefficients for the transformation are available. The governing equations are developed in primitive variables in order to facilitate the incorporation of physical boundary conditions and turbulence-closure models. MacCormack's two-step, unsplit, time-marching algorithm is used to solve the unsteady Navier-Stokes equations until steady-state solution is achieved. Cases discussed include (1) flat plate in supersonic free stream; (2) supersonic flow along an axial corner; (3) subsonic flow in an axial corner at M infinity = 0.95; and (4) supersonic flow in an axial corner at M infinity 1.5.

Vatsa, V. N.↗
NASA NTRS1984

Scalar and momentum turbulent transport experiments with swirling and nonswirling flows

Combustor models for the aircraft gas turbine industry have been obtained because of the need to reduce the costs of developing improved performance and more durable engines. A few years ago, it became apparent that the mass concentration and velocity predictions provided by the computer codes were not representing the data measured in some confined recirculating flows. It is pointed out that errors in the mass concentration distribution are an especially serious problem because of their influence on the heat release, temperature, and reactant distributions. Combined mass and momentum turbulent transport experiments with swirling and nonswirling flow have been conducted with the objective to obtain an experimental data base which can be used to evaluate and improve the turbulent transport submodes employed in the aerothermal models. The present paper is mainly concerned with the overall characteristics of the mass turbulent transport processes in complex flows with recirculation and the deficiencies of the conventional models.

Johnson, B. V.↗
NASA NTRS2003

Stormtime Particle Energization with AMIE Potentials

Simulations were conducted to investigate the influence of rapid electric field fluctuations on electron energization in the inner magnetosphere based on the assimilative mapping of ionospheric electrodynamics (AMIE) technique. Simulations for four different magnetic storms were run, namely those that occurred on May 15,1997, May 4, 1998, September 25, 1998, and October 19, 1998. Here, we have examined the formation of high energy electrons in the inner magnetosphere during these storm events with our recently-developed relativistic radiation belt transport code. The point of this numerical experiment is to show that a simulation of a real event must have the high time resolution electric field input files in order to produce the seed population for the radiation belts, which are often observed to increase in the days following a magnetic storm. Specifically, a cadence of the global electric field pattern of 5 minutes or less produces inner magnetospheric fluxes that are larger (by up to 5 orders of magnitude) than fluxes produced with a longer cadence. Differences were particularly large relative to simulation results with a 3-hour time cadence, analogous to a Kp-driven electric field model.

Khazanov, George V.↗
NASA NTRS2004

Stormtime Particle Energization with High Temporal Resolution AMIE Potentials

Simulations were conducted to investigate the influence of rapid electric field fluctuations on electron energization in the inner magnetosphere based on the assimilative mapping of ionospheric electrodynamics (AMIE) technique. Simulations for four different magnetic storms were run, namely those that occurred on May 15, 1997, May 4, 1998, September 25, 1998, and October 19, 1998. Here, we have examined the formation of high energy (10-1000 keV) electrons in the inner magnetosphere during these storm events with our recently-developed relativistic radiation belt transport code. The point of this numerical experiment is to show that a simulation of a real event must have the high time resolution electric field input files in order to produce the seed population for the radiation belts, which are often observed to increase in the days following a magnetic storm. Specifically, a cadence of the global electric field pattern of 5 minutes or less produces inner magnetospheric fluxes that are larger (by up to several orders of magnitude) than fluxes produced with a longer cadence. Differences were particularly large relative to simulation results with a 3-hour time cadence, analogous to a Kp-driven electric field model.

Khazanov, George V.↗
NASA NTRS2005

Low Noise Exhaust Nozzle Technology Development

NASA and the U.S. aerospace industry have been assessing the economic viability and environmental acceptability of a second-generation supersonic civil transport, or High Speed Civil Transport (HSCT). Development of a propulsion system that satisfies strict airport noise regulations and provides high levels of cruise and transonic performance with adequate takeoff performance, at an acceptable weight, is critical to the success of any HSCT program. The principal objectives were to: 1. Develop a preliminary design of an innovative 2-D exhaust nozzle with the goal of meeting FAR36 Stage III noise levels and providing high levels of cruise performance with a high specific thrust for Mach 2.4 HSCT with a range of 5000 nmi and a payload of 51,900 lbm, 2. Employ advanced acoustic and aerodynamic codes during preliminary design, 3. Develop a comprehensive acoustic and aerodynamic database through scale-model testing of low-noise, high-performance, 2-D nozzle configurations, based on the preliminary design, and 4. Verify acoustic and aerodynamic predictions by means of scale-model testing. The results were: 1. The preliminary design of a 2-D, convergent/divergent suppressor ejector nozzle for a variable-cycle engine powered, Mach 2.4 HSCT was evolved, 2. Noise goals were predicted to be achievable for three takeoff scenarios, and 3. Impact of noise suppression, nozzle aerodynamic performance, and nozzle weight on HSCT takeoff gross weight were assessed.

Majjigi, R. K.↗
NASA NTRS2010

A Model for the Sources of the Slow Solar Wind

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed-field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind has large angular width, up to approximately 60 degrees, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far front the heliospheric current sheet. We then use an MHD code and MIDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

Antiochos, Spiro K.↗
NASA NTRS1990

Computation of turbine flowfields with a Navier-Stokes code

A new technique has been developed for the solution of the incompressible Navier-Stokes equations. The numerical technique, derived from a pressure substitution method (PSM), overcomes many of the deficiencies of the pressure crrection method. This technique allows for the direct solution of the actual pressure in the form of a Poisson equation which is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. In two-dimensions a turbine flowfield, including heat transfer, has been computed with this method and the prediction of the cascade performance is presented. The extension of the pressure correction method for the solution of three-dimensional flows is also presented for laminar flow in an S-shaped duct and turbulent flow in the end-wall region of a turbine cascade.

Hobson, G. V.↗
DOE OSTI2020

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Results from a synthetic model of the ITER XRCS-Core diagnostic based on high-fidelity x-ray ray tracing

A high-fidelity synthetic diagnostic has been developed for the ITER core x-ray crystal spectrometer diagnostic based on x-ray ray tracing. This synthetic diagnostic has been used to model expected performance of the diagnostic, to aid in diagnostic design, and to develop engineering tolerances. The synthetic model is based on x-ray ray tracing using the recently developed xicsrt ray tracing code and includes a fully three-dimensional representation of the diagnostic based on the computer aided design. The modeled components are: plasma geometry and emission profiles, highly oriented pyrolytic graphite pre-reflectors, spherically bent crystals, and pixelated x-ray detectors. Plasma emission profiles have been calculated for Xe 44+ , Xe 47+ , and Xe 51+ , based on an ITER operational scenario available through the Integrated Modelling & Analysis Suite database, and modeled within the ray tracing code as a volumetric x-ray source; the shape of the plasma source is determined by equilibrium geometry and an appropriate wavelength distribution to match the expected ion temperature profile. All individual components of the x-ray optical system have been modeled with high-fidelity producing a synthetic detector image that is expected to closely match what will be seen in the final as-built system. Particular care is taken to maintain preservation of photon statistics throughout the ray tracing allowing for quantitative estimates of diagnostic performance.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Polarization–frequency hyperentangled photons: generation, characterization, and manipulation

Frequency-bin encoding is massively parallelizable and robust for optical fiber transmission. When coupled with an additional degree of freedom (DoF), the expansion of the Hilbert space allows for deterministic controlled operations between two DoFs within a single photon. Such capabilities, when combined with photonic hyperentanglement, are of great value for quantum communication protocols, including dense coding and single-copy entanglement distillation. In this talk, we present an all-fiber-coupled, ultrabroadband polarization–frequency hyperentangled source and conduct comprehensive quantum state tomography across multiple dense wavelength division multiplexing channels spanning the optical C+L-band (1530–1625 nm). In addition, we design and implement a high-fidelity controlled-NOT (cnot) operation between polarization and frequency DoFs by exploiting electro-optic phase modulation within a fiber Sagnac loop. Collectively, our hyperentangled source and two-qubit gate should unlock new opportunities for harnessing polarization–frequency resources in established telecommunication fiber networks for future quantum applications.

Lu, Hsuan-Hao↗
DOE OSTI2023

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗