Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear forces”

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 505 records · Page 28

Robotic Planetary Science Missions Enabled with Small NTR Engine/Stage Technologies

The high specific impulse (Isp) and engine thrust-to-weight ratio of liquid hydrogen (LH2)-cooled nuclear thermal rocket (NTR) engines makes them ideal for upper stage applications to difficult robotic planetary science missions. A small 15 thousand pound force (klbf) NTR engine using a uranium-zirconium-niobium 'ternary carbide' fuel (Isp approximately 960 seconds at approximately 3025K) developed in the Commonwealth of Independent States (CIS) is examined and its use on an expendable injection stage is shown to provide major increases in payload delivered to the outer planets (Saturn, Uranus, Neptune and Pluto). Using a single 'Titan IV-class' launch vehicle, with a lift capability to low Earth orbit (LEO) of approximately 20 metric tons (t), an expendable NTR upper stage can inject two Pluto 'Fast Flyby' spacecraft (PFF/SC) plus support equipment-combined mass of approximately 508 kg--on high energy, '6.5-9.2 year' direct trajectory missions to Pluto. A conventional chemical propulsion mission would use a liquid oxygen (LOX)/LH2 'Centaur' upper stage and two solid rocket 'kick motors' to inject a single PFF/SC on the same Titan IV launch vehicle. For follow on Pluto missions, the NTR injection stage would utilize a Jupiter 'gravity assist' (JGA) maneuver to launch a LOX/liquid methane (CH4) capture stage (Isp approximately 375 seconds) and a Pluto 'orbiter' spacecraft weighing between approximately 167-312 kg. With chemical propulsion, a Pluto orbiter mission is not a viable option because c inadequate delivered mass. Using a 'standardized' NTR injection stage and the same single Titan IV launch scenario, 'direct flight' (no gravity assist) orbiter missions to Saturn, Uranus and Neptune are also enabled with transit times of 2.3, 6.6, and 12.6 years, respectively. Injected mass includes a storable, nitrogen tetroxide/monomethyl hydrazine (N2O4/MMH) capture stage (Isp approximately 330 seconds) and orbiter payloads 340 to 820% larger than that achievable using a LOX/LH2-fueled injection stage. The paper discusses NTR technology and mission characteristics, shows NTR stage and payload accommodations within the 26.2 m long Titan IV payload fairing, and discusses NTR stage performance as a function of assumed cryogenic tank technology.

Borowski, Stanley K.↗

Measurement of the Shear Lift Force on a Bubble in a Channel Flow

Two-phase flow systems play vital roles in the design of some current and anticipated space applications of two-phase systems which include: thermal management systems, transfer line flow in cryogenic storage, space nuclear power facilities, design and operation of thermal bus, life support systems, propulsion systems, In Situ Resource Utilization (ISRU), and space processes for pharmaceutical applications. The design of two-phase flow systems for space applications requires a clear knowledge of the behaviors of the dispersed phase (bubble), its interaction with the continuous phase (liquid) and its effect on heat and mass transfer processes, The need to understand the bubble generation process arises from the fact that for all space applications, the size and distribution of bubbles are extremely crucial for heat and mass transfer control. One important force in two-phase flow systems is the lift force on a bubble or particle in a liquid shear flow. The shear lift is usually overwhelmed by buoyancy in normal gravity, but it becomes an important force in reduced gravity. Since the liquid flow is usually sheared because of the confining wall, the trajectories of bubbles and particles injected into the liquid flow are affected by the shear lift in reduced gravity. A series of experiments are performed to investigate the lift force on a bubble in a liquid shear flow and its effect on the detachment of a bubble from a wall under low gravity conditions. Experiments are executed in a Poiseuille flow in a channel. An air-water system is used in these experiments that are performed in the 2.2 second drop tower. A bubble is injected into the shear flow from a small injector and the shear lift is measured while the bubble is held stationary relative to the fluid. The trajectory of the bubble prior, during and after its detachment from the injector is investigated. The measured shear lift force is calculated from the trajectory of the bubble at the detachment point. These values for the shear lift are then compared with the theoretical predictions from various published works on shear lift in the open literature, which include asymptotic solutions at low bubble Reynolds number, potential flow predictions and numerical studies that deal with intermediate bubble Reynolds numbers.

Nahra, Henry K.↗

Scaling Law for Cloud-Rise Velocity vs. Scaled Height of Burst (SHOB)

The shockwave produced by a nuclear detonation leaves a low-density region in the vicinity of air zero (i.e., the point of detonation). As the shockwave degenerates into an acoustic wave, the buoyant force takes over and causes the low-density region to begin to rise at speeds that can reach as high as ~200-300 mph. The objective of this project was to find out if this initial cloud rise followed any scaling law, by stimulating a cloud rise for a 1 kiloton (kt) detonation and determining if a correlation existed between the scaled-height-of-burst (SHOB) and the scaled-cloud-rise velocity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Security-by-Design: Light Water Small Modular Reactor

The growing demand for nuclear power is increasing pressure to find solutions to cost prohibitive requirements of both construction and security. Offsite response has been proposed as an option to reduce costs associated with training and maintaining an onsite response force. A previous report explored this option and revealed that security could be provided at the required level, but cost savings was not a result of this methodology. An offsite response strategy required costly active and passive delay barriers to provide sufficient time for responders to muster and deploy to a site in time to interrupt a determined and well-equipped adversary. Also, contrary to the hypothesis, the number of responders required for this strategy exceeded that needed for an onsite response force, as the adversaries could avail themselves of advantageous positions within the facility to repel arriving responders. This report builds upon the previous evaluation by using the same hypothetical light water small modular reactor (LWSMR) facility model, but this time an onsite response strategy was assessed. The goal of this analysis was to show that an onsite response strategy could be implemented effectively at a cost point that removes barriers within the industry at this critical time of growth and development. The assessment of the facility design and response strategy was completed through modeling using Scribe3D© and subsequent scenario analysis over the course of a two-day tabletop exercise. Subject matter experts in nuclear security, nuclear facility design, and response strategy and tactics contributed to the effort to ensure accurate representation of hypothetical scenarios. Several adaptations were made to the layout of the LWSMR based on lessons learned during the first day of scenario analysis. The subsequent design evaluated on the second day proved to provide a robust response posture against a large and well-trained adversary force. This report details the process of the analysis and compares the cost of the final facility design with that of the LWSMR model used for evaluation of offsite response. Ultimately, the results of this effort indicate that, when implemented correctly, an onsite response strategy is the best option from a security and cost perspective.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Earthquake Testing

During NASA's Apollo program, it was necessary to subject the mammoth Saturn V launch vehicle to extremely forceful vibrations to assure the moonbooster's structural integrity in flight. Marshall Space Flight Center assigned vibration testing to a contractor, the Scientific Services and Systems Group of Wyle Laboratories, Norco, California. Wyle-3S, as the group is known, built a large facility at Huntsville, Alabama, and equipped it with an enormously forceful shock and vibration system to simulate the liftoff stresses the Saturn V would encounter. Saturn V is no longer in service, but Wyle-3S has found spinoff utility for its vibration facility. It is now being used to simulate earthquake effects on various kinds of equipment, principally equipment intended for use in nuclear power generation. Government regulations require that such equipment demonstrate its ability to survive earthquake conditions. In upper left photo, Wyle3S is preparing to conduct an earthquake test on a 25ton diesel generator built by Atlas Polar Company, Ltd., Toronto, Canada, for emergency use in a Canadian nuclear power plant. Being readied for test in the lower left photo is a large circuit breaker to be used by Duke Power Company, Charlotte, North Carolina. Electro-hydraulic and electro-dynamic shakers in and around the pit simulate earthquake forces.

Source record↗

Scattering matrix pole expansions for complex wave numbers in R -matrix theory

In this followup article to Ducru et al., we establish new results on scattering matrix pole expansions for complex wave numbers in R-matrix theory. In the past, two branches of theoretical formalisms emerged to describe the scattering matrix in nuclear physics: R-matrix theory and pole expansions. The two have been quite isolated from one another. Recently, our study of Brune's alternative parametrization of R-matrix theory has shown the need to extend the scattering matrix (and the underlying R-matrix operators) to complex wave numbers. Two competing ways of doing so have emerged from a historical ambiguity in the definitions of the shift S and penetration P functions: the legacy Lane and Thomas's “force closure” approach versus analytic continuation (which is the standard in mathematical physics). The R-matrix community has not yet come to a consensus as to which to adopt for evaluations in standard nuclear data libraries, such as ENDF. Here, in this article, we argue in favor of analytic continuation of R-matrix operators. We bridge R-matrix theory with the Humblet-Rosenfeld pole expansions, and discover new properties of the Siegert-Humblet radioactive poles and widths, including their invariance properties to changes in channel radii a c . We then show that analytic continuation of R-matrix operators preserves important physical and mathematical properties of the scattering matrix—canceling spurious poles and guaranteeing generalized unitarity—while still being able to close channels below thresholds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Toward Exascale: Overview of Large Eddy Simulations and Direct Numerical Simulations of Nuclear Reactor Flows with the Spectral Element Method in Nek5000

At the beginning of the last decade, Petascale supercomputers (i.e., computers capable of more than 1 petaFLOP) emerged. Now, at the dawn of exascale supercomputing, we provide a review of recent landmark simulations of portions of reactor components with turbulence-resolving techniques that this computational power has made possible. In fact, these simulations have provided invaluable insight into flow dynamics, which is difficult or often impossible to obtain with experiments alone. We focus on simulations performed with the spectral element method, as this method has emerged as a powerful tool to deliver massively parallel calculations at high fidelity by using large eddy simulation or direct numerical simulation. We also limit this paper to constant-property incompressible flow of a Newtonian fluid in the absence of other body or external forces, although the method is by no means limited to this class of flows. We briefly review the fundamentals of the method and the reasons it is compelling for the simulation of nuclear engineering flows. We review in detail a series of Petascale simulations, including the simulations of helical coil steam generators, fuel assemblies, and pebble beds. Even with Petascale computing, however, limitations for nuclear modeling and simulation tools remain. In particular, the size and scope of turbulence-resolving simulations are still limited by computing power and resolution requirements, which scale with the Reynolds number. In the final part of this paper, we discuss the future of the field, including recent advancements in emerging architectures such as GPUbased supercomputers, which are expected to power the next generation of high-performance computers.

computational fluid dynamics↗

Mesoscale Modeling of the Effects of Accelerated Burnup on UO2 Microstructural Evolution

Accelerating the nuclear fuel qualification process will rely on some combination of advanced modeling and simulation techniques with accelerated irradiation testing and separate effects experiments to enable the development of new fuel concepts in a shorter time frame. One of the key challenges to successfully leveraging accelerated irradiation tests will be understanding the artifacts that may be introduced with accelerated accumulation of dose and/or burnup. This work presents phase field (MARMOT) simulations of the evolution of representative 2D UO2 microstructures up to 40 MWd/kgU. Simulations were performed under both commercial light water reactor fuel conditions as well as those that would be expected for highly accelerated (~10x) burnup conditions similar to those used in the MiniFuel irradiations in Oak Ridge National Laboratory’s High Flux Isotope Reactor. The phase field model was coupled with a discrete nucleation algorithm to model re- structuring at high burnup. The effect of the different fission rates in both microstructures was investigated at two temperatures: 650?C and 800?C. The lower temperature simulations both showed an onset of restructuring at nearly 60 MWd/kgU. More extensive restructuring was obtained in the MiniFuel microstructure compared with that of the PWR fuel. At 800?C, no restructuring was obtained as a result of the thermally activated diffusion of Xe atoms and U vacancies to fission gas bubbles, which reduces the nucleation driving force. These results highlight the importance of using modeling and simulation tools to inform the environmental conditions during targeted accelerated irradiation tests to extract the most useful fuel performance data.

accelerated fuel qualification, Phase Field, Restr↗

Fundamental Symmetries, Neutrons, and Neutrinos (FSNN): Whitepaper for the 2023 NSAC Long Range Plan

Through the exploration of fundamental symmetries, and by using nuclei, neutrons, and neutrinos, nuclear physics addresses some of the most profound questions in science. Why does the universe contain so much more matter than antimatter? Are neutrinos their own antiparticles and where do their masses come from? What objects make up the dark matter that is responsible for most of the universe’s mass? Does nature contain more forces than the four we know about? Our Standard Model of nature’s particles and forces is incomplete because it does not answer these questions; new physics, from beyond the Standard Model (BSM) is needed. With that physics not appearing at the high energy frontier, it has become imperative to realize the potential of the burgeoning program of precision nuclear-physics measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Model

The ITER electron cyclotron emission (ECE) diagnostic system is located at diagnostics shielding module 2 (DSM2), equatorial port 9 (EP9), to measure electron temperature profile and electron temperature fluctuations, and also assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two hot sources for calibration and two couples of mirrors for two different optical views. It shall receive that the large electromagnetic (EM) loads up to 100-MN/m3 force density due to the eddy currents on the copper mirrors generated by the short transient plasma disruption. The simplified EM analysis model using magnetic field ( B ) data and flux variations (dB/dt) method based on the worst case of plasma disruption, MD_DW_EXP16MS_CATIII, has been developed to calculate the force and moments for each ECE component and support structure for EP9 DSM2 Bay2 and Bay3. The B and dB/dt methods use the constant B and dB/dt assumed for each local bay space, but the real B and dB/dt are basically varied along the radial direction of each bay. The complex global model with the local ECE components and support structure using Maxwell transient has been developed to do EM analysis to compare the results of the B and dB/dt models. The volumetric EM force density of the whole ECE components and support structure can be used for the subsequent structural analysis. In conclusion, combined with the thermal and nuclear loads and seismic and initial loads, the ECE integration analysis can be finalized for the operation case.

Fang, J. [Princeton Plasma Physics Laboratory (PPP↗

Heat transfer from glass melt to cold cap: Computational fluid dynamics study of cavities beneath cold cap

Efficient glass production depends on the continuous supply of heat from the glass melt to the floating layer of batch, or cold cap. Computational fluid dynamics (CFD) are employed to investigate the formation and behavior of gas cavities that form beneath the batch by gases released from the collapsing primary foam bubbles, ascending secondary bubbles, and in the case of forced bubbling, from the rising bubbling gas. The gas phase fraction, temperature, and velocity distributions below the cold cap are used to calculate local and average heat transfer rates as a function of the bubbling rate. It is shown that the thickness of the cavities is nearly independent of the cold cap shape and the amount of foam evolved during batch conversion. Furthermore, it is ~7 mm and up to ~15 mm for the cases without and with forced bubbling used to promote circulation within the melt, respectively. Using computed velocity and temperature profiles, the melting rate of the simulated high-level nuclear waste glass batch was estimated to increase with the bubbling rate to the power of ~0.3 to 0.9, depending on the flow pattern. The simulation results are in good agreement with experimental data from laboratory- and pilot-scale melter tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Combined Radiation and Forced Convection for Flow of a Transparent Gas in a Tube with Sinusoidal Axial Wall Heat Flux Distribution

An analysis is made of the influence that radiation exchanges between elements on the inside surface of a tube have on the wall and gas temperature distributions for forced-convection flow. The wall heat generation has a chopped sine distribution with length along the tube, which is a distribution often encountered in nuclear-reactor channels. The flowing gas is assumed transparent to thermal radiation and hence does not participate directly in the radiative exchange process. Axial heat conduction is neglected in the gas and tube wall, and the convective-heat-transfer coefficient and fluid properties are assumed constant. Several numerical examples are given to illustrate the effects of eight independent parameters such as wall emissivity, Stanton number, and length-diameter ratio. In some instances the radiation exchanges reduced the peak wall temperature or caused a reduction in the exit gas temperature because of radiation losses from the tube end openings.

TUBE↗

Flight test techniques for validating simulated nuclear electromagnetic pulse aircraft responses

An attempt has been made to determine the effects of nuclear EM pulses (NEMPs) on aircraft systems, using a highly instrumented NASA F-106B to document the simulated NEMP environment at the Kirtland Air Force Base's Vertically Polarized Dipole test facility. Several test positions were selected so that aircraft orientation relative to the test facility would be the same in flight as when on the stationary dielectric stand, in order to validate the dielectric stand's use in flight configuration simulations. Attention is given to the flight test portions of the documentation program.

Winebarger, R. M.↗

Simplifying silicon burning: Application of quasi-equilibrium to (alpha) network nucleosynthesis

While the need for accurate calculation of nucleosynthesis and the resulting rate of thermonuclear energy release within hydrodynamic models of stars and supernovae is clear, the computational expense of these nucleosynthesis calculations often force a compromise in accuracy to reduce the computational cost. To redress this trade-off of accuracy for speed, the authors present an improved nuclear network which takes advantage of quasi- equilibrium in order to reduce the number of independent nuclei, and hence the computational cost of nucleosynthesis, without significant reduction in accuracy. In this paper they will discuss the first application of this method, the further reduction in size of the minimal alpha network. The resultant QSE- reduced alpha network is twice as fast as the conventional alpha network it replaces and requires the tracking of half as many abundance variables, while accurately estimating the rate of energy generation. Such reduction in cost is particularly necessary for future generation of multi-dimensional models for supernovae.

Hix, W. R.↗

Final Report - Hypergravity Effects on Chromatin Conformation and Nuclear Structure in Cultured Cells using High-throughput Sequencing

The conformation of the human genome is known to play an important role in transcriptional control of gene expression. Our aim is to assess whether exposure of cultured human cardiomyocytes to hypergravity can induce changes in chromatin organization as assayed by Hi-C proximity ligation technique. Hi-C is used to analyze chromatin interactions by using formaldehyde to crosslink regions of chromatin that are in close proximity. The DNA is then fragmented using a restriction enzyme and ligated under dilute conditions to favor intramolecular ligation of cross-linked fragments. Finally, the DNA is sequenced, allowing reconstruction of genomic structure. iCell cardiomyocyte cultures will be exposed to 40 g via an engineered device (incu-fuge) that allows for the chemical fixation of cells on a spinning centrifuge. Additionally, we plan on using fluorescence microscopy to analyze the nuclear and actin cytoskeletal conformation of cardiomyocytes post hypergravity exposure. We hypothesize that chromatin will rearrange in hypergravity conditions, both rapidly due to direct mechanical forces, as well as over longer time frames due to changes to structures necessary for transcriptional responses such as the de novo formation of promoter-enhancer loops. A negative result–no change between 1 g and hypergravity conditions–would nevertheless be an important data point in our understanding of the rheology of the nucleus, and how cells and cellular structures respond to different gravity fields. The cellular physiology of hypergravity has clinical relevance for spaceflight, and can further inform our understanding of microgravity physiology.

chromatin↗

Caustic Neutralization and Precipitation of Acidic Dissolved Simulated Stainless Steel–Clad Plutonium and Plutonium/Uranium Nuclear Fuel

Simulated dissolved stainless steel (SS) clad Pu and Pu/U nuclear fuel in HNO 3 was neutralized to a free hydroxide (OH - ) concentration of 0.6 M. A thermal neutron poison, Gd, was added to the simulants at concentrations of either ~3 - 6 g/L or ~37 - 38 g/L. The supernate Pu concentration the day of neutralization ranged from 0.48 to 8.75 mg/L. The supernate Pu concentration of a simplified simulant neutralized to 0.6 M OH - above precipitated solids containing Pu was demonstrated to decrease over 18 days. A significant portion of precipitated Pu was found to be insoluble in 8 M HNO 3 at ambient temperature, but essentially quantitative Pu dissolution was achieved in 11.5 M HNO 3 /0.1 M KF at 100 °C. The difficulty in dissolving the Pu precipitate is believed to be due to the formation of refractory PuO 2 •xH 2 O during the neutralization process. Initial Gd concentrations of ~37 – 38 g/L were found to result in a greater Al precipitation when neutralized to 0.6 M OH - than initial Gd concentrations of ~3 – 6 g/L. Physical properties of the resultant slurries were measured and used to calculate limiting flowrates and slurry velocities by gravity only in transfer piping between the Savannah River Site’s H-Canyon Facility and the Concentration, Storage, and Transfer Facility (CSTF). These results were compared to calculated deposition velocities to predict if solids would settle during the transfer. The Newtonian model was found to be reasonable for each diluted slurry evaluated. Deposition velocities of Pu containing slurries are lower than nuclear fuel slurries primarily composed of U due to the high density of Pu solids. In conclusion, dilution of slurries reduces the margin between the slurry and deposition velocities due to the reduction in viscosity because higher viscous forces on the particles promote maintained suspension.

Actinide Neutralization↗