Search NASA⌕ Search

SEARCH · Search NASA

Results for “Beam code development”

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 433 records · Page 24

Multi-Body Analysis of the 1/5 Scale Wind Tunnel Model of the V-22 Tiltrotor

The paper presents a multi-body analysis of the 1/5 scale wind tunnel model of the V-22 tiltrotor, the Wing and Rotor Aeroelastic Testing System (WRATS), currently tested at NASA Langley Research Center. An original multi-body formulation has been developed at the Dipartimento di Ingegneria Aerospaziale of the Politecnico di Milano, Italy. It is based on the direct writing of the equilibrium equations of independent rigid bodies, connected by kinematic constraints that result in the addition of algebraic constraint equations, and by dynamic constraints, that directly contribute to the equilibrium equations. The formulation has been extended to the simultaneous solution of interdisciplinary problems by modeling electric and hydraulic networks, for aeroservoelastic problems. The code has been tailored to the modeling of rotorcrafts while preserving a complete generality. A family of aerodynamic elements has been introduced to model high aspect aerodynamic surfaces, based on the strip theory, with quasi-steady aerodynamic coefficients, compressibility, post-stall interpolation of experimental data, dynamic stall modeling, and radial flow drag. Different models for the induced velocity of the rotor can be used, from uniform velocity to dynamic in flow. A complete dynamic and aeroelastic analysis of the model of the V-22 tiltrotor has been performed, to assess the validity of the formulation and to exploit the unique features of multi-body analysis with respect to conventional comprehensive rotorcraft codes; These are the ability to model the exact kinematics of mechanical systems, and the possibility to simulate unusual maneuvers and unusual flight conditions, that are particular to the tiltrotor, e.g. the conversion maneuver. A complete modal validation of the analytical model has been performed, to assess the ability to reproduce the correct dynamics of the system with a relatively coarse beam model of the semispan wing, pylon and rotor. Particular care has been used to model the kinematics of the gimbal joint, that characterizes the rotor hub, and of the control system, consisting in the entire swashplate mechanism. The kinematics of the fixed and the rotating plates have been modeled, with variable length control links used to input the controls, the rotating flexible links, the pitch horns and the pitch bearings. The investigations took advantage of concurring wind tunnel test runs, that were performed in August 1998, and allowed the acquisition of data specific to the multi-body analysis.

Ghiringhelli, G. L.↗

NASA Tech Briefs, August 2011

Topics covered include: Miniature, Variable-Speed Control Moment Gyroscope; NBL Pistol Grip Tool for Underwater Training of Astronauts; HEXPANDO Expanding Head for Fastener-Retention Hexagonal Wrench; Diagonal-Axes Stage for Pointing an Optical Communications Transceiver; Improvements in Speed and Functionality of a 670-GHz Imaging Radar; IONAC-Lite; Large Ka-Band Slot Array for Digital Beam-Forming Applications; Development of a 150-GHz MMIC Module Prototype for Large-Scale CMB Radiation; Coupling Between Waveguide-Fed Slot Arrays; PCB-Based Break-Out Box; Multiple-Beam Detection of Fast Transient Radio Sources; Router Agent Technology for Policy-Based Network Management; Remote Asynchronous Message Service Gateway; Automatic Tie Pointer for In-Situ Pointing Correction; Jitter Correction; MSLICE Sequencing; EOS MLS Level 2 Data Processing Software Version 3; DspaceOgre 3D Graphics Visualization Tool; Metallization for Yb14MnSb11-Based Thermoelectric Materials; Solvent/Non-Solvent Sintering To Make Microsphere Scaffolds; Enhanced Fuel-Optimal Trajectory-Generation Algorithm for Planetary Pinpoint Landing; Self-Cleaning Coatings and Materials for Decontaminating Field-Deployable Land and Water-Based Optical Systems; Separation of Single-Walled Carbon Nanotubes with DEP-FFF; Li Anode Technology for Improved Performance; Post-Fragmentation Whole Genome Amplification-Based Method; Microwave Tissue Soldering for Immediate Wound Closure; Principles, Techniques, and Applications of Tissue Microfluidics; Robotic Scaffolds for Tissue Engineering and Organ Growth; Stress-Driven Selection of Novel Phenotypes; Method for Accurately Calibrating a Spectrometer Using Broadband Light; Catalytic Microtube Rocket Igniter; Stage Cylindrical Immersive Display; Vacuum Camera Cooler; Atomic Oxygen Fluence Monitor; Thermal Management Tools for Propulsion System Trade Studies and Analysis; Introduction to Physical Intelligence; Technique for Solving Electrically Small to Large Structures for Broadband Applications; Accelerated Adaptive MGS Phase Retrieval; Large Eddy Simulation Study for Fluid Disintegration and Mixing; Tropospheric Correction for InSAR Using Interpolated ECMWF Data and GPS Zenith Total Delay; Technique for Calculating Solution Derivatives With Respect to Geometry Parameters in a CFD Code; Acute Radiation Risk and BRYNTRN Organ Dose Projection Graphical User Interface; Probabilistic Path Planning of Montgolfier Balloons in Strong, Uncertain Wind Fields; Flight Simulation of ARES in the Mars Environment; Low-Outgassing Photogrammetry Targets for Use in Outer Space; Planning the FUSE Mission Using the SOVA Algorithm; Monitoring Spacecraft Telemetry Via Optical or RF Link; and Robust Thermal Control of Propulsion Lines for Space Missions.

Source record↗

Synchrotron micro-computed tomography analysis of neutron-irradiated U-Mo fuel

The three-dimensional (3D) microstructure of neutron-irradiated uranium-10 wt.% molybdenum (U-10Mo) fuel with a burn-up of 9.8 × 10 21 fissions/cm 3 was characterized using a novel, multi-modal synchrotron micro-computed tomography approach combining propagation-based phase-contrast enhanced and absorption contrast techniques. The porosity development, porosity interconnectedness, swelling, composition, local thickness of the zirconium (Zr) diffusion barrier, and the influence of the fuel–cladding interaction on the local composition and pore morphology, were uniquely determined in 3D. Two cuboids were produced using a focused ion beam-scanning electron microscope at the Zr diffusion barrier–fuel interface and in the bulk fuel. The bulk fuel sample swelled by 53.3 [+9.7/−3.1]%, while the fuel near the Zr–fuel interface swelled by 63.3 [+14.7/−7.2]%. The average local thickness of the Zr diffusion barrier decreased by 53 %, compared to the expected pre-irradiated thickness. Four pore morphology regions were identified initiating parallel to the fuel–Zr interaction region: (1) an interaction layer of suppressed porosity, (2) a layer of elongated and interconnected porosity, (3) a transition zone of low porosity, and (4) a layer of unoriented porosity representative of the bulk fuel behavior. The increase in porosity near the diffusion barrier corresponded to a higher U concentration compared to that in the bulk fuel. The interconnected porosity in the fuel near the diffusion barrier was extensive and oriented parallel to the diffusion barrier, while the bulk fuel had more compact and isolated pore networks. The interaction layer, despite having suppressed porosity, was nearly 100 wt.% U. Porosity suppression at the diffusion barrier corresponds to the expected reduction in radiation-driven diffusion of Xe at the interface despite the anticipated increase in fission product nucleation originating from a higher U concentration. In conclusion, the novel 3D insights of the porosity, swelling, and compositional variations characterized herein can improve the fidelity of fuel performance codes for proliferation-resistant fuels for research and test reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Update On the Status of the FLUKA Monte Carlo Transport Code*

The FLUKA Monte Carlo transport code is a well-known simulation tool in High Energy Physics. FLUKA is a dynamic tool in the sense that it is being continually updated and improved by the authors. We review the progress achieved since the last CHEP Conference on the physics models, some technical improvements to the code and some recent applications. From the point of view of the physics, improvements have been made with the extension of PEANUT to higher energies for p, n, pi, pbar/nbar and for nbars down to the lowest energies, the addition of the online capability to evolve radioactive products and get subsequent dose rates, upgrading of the treatment of EM interactions with the elimination of the need to separately prepare preprocessed files. A new coherent photon scattering model, an updated treatment of the photo-electric effect, an improved pair production model, new photon cross sections from the LLNL Cullen database have been implemented. In the field of nucleus-- nucleus interactions the electromagnetic dissociation of heavy ions has been added along with the extension of the interaction models for some nuclide pairs to energies below 100 MeV/A using the BME approach, as well as the development of an improved QMD model for intermediate energies. Both DPMJET 2.53 and 3 remain available along with rQMD 2.4 for heavy ion interactions above 100 MeV/A. Technical improvements include the ability to use parentheses in setting up the combinatorial geometry, the introduction of pre-processor directives in the input stream. a new random number generator with full 64 bit randomness, new routines for mathematical special functions (adapted from SLATEC). Finally, work is progressing on the deployment of a user-friendly GUI input interface as well as a CAD-like geometry creation and visualization tool. On the application front, FLUKA has been used to extensively evaluate the potential space radiation effects on astronauts for future deep space missions, the activation dose for beam target areas, dose calculations for radiation therapy as well as being adapted for use in the simulation of events in the ALICE detector at the LHC.

Ferrari, A.↗

Theory and Simulation of Ultrafast Multidimensional Nonlinear X-ray Spectroscopy of Molecules (Final Report)

Emerging X-ray free electron laser (XFEL) beam sources offer new types of probes of matter with unprecedented spatial and temporal resolutions. These experimental advances must be met by robust theoretical and computational tools that provide predictive modeling capacity of the underlining electronic and structural dynamics. The latter will be essential for the design of sophisticated multi-pulse experiments and for their interpretation. The proposed research effort will focus on developing cutting- edge simulation tools for nonlinear multidimensional X-ray/optical spectroscopies and aims to address key questions in Priority Research Opportunities 1 (Probing and controlling electron motion within a molecule) and 3 (Capturing rare events and intermediate states in the transformation of matter) as mentioned in the BES Roundtable Report “Opportunities for Basic Research at the Frontiers of XFEL Ultrafast Science”. XFEL multidimensional nonlinear techniques, which combine sequences of X-ray and possibly optical pulses, provide a unique experimental toolbox for probing the dynamics of core and valence electronic excitations, as well as material structure. Predictive modeling of these dynamical processes requires the combination of analytical theory for nonlinear interactions of light and matter, robust quantum-chemical methodologies for the accurate description of electronic structure of various materials, and multiscale ab initio electron and nuclear dynamics techniques operating beyond Born-Oppenheimer approximation. These challenges will be addressed with three research thrusts (i) Develop and implement theoretical apparatus for modeling a broad range of multidimensional spectroscopic techniques enabled by present and upcoming XFEL facilities. This thrust also includes the incorporation of a computational module in the DOE supported open-source NWChem computational chemistry package as well as the development of other open-source codes ready for dissemination across a broad user base; (ii) Propose and design new multi-pulse experiments that make use of the capabilities of the incoming LCLS-II facility; (iii) Perform selected applications to specific molecular systems that can be carried out at LCLS-II and demonstrate how these X-ray sources may be used to study nonadiabatic dynamics through conical interactions, electronic correlations in multi-core excitons, and charge transfer/energy transfer processes. The proposed research will be carried out by a multi-disciplinary four-institution research team which combines academia and national laboratories and spans the broad and necessary expertise in theoretical spectroscopy, nonlinear optics, quantum chemistry, molecular non-adiabatic dynamics and code development. The work will be performed in a highly interactive team environment with junior researchers shared between institutions thus cementing cross- disciplinary interactions. The developed simulation tools will be immediately deployed for XFEL facility users, both experimentalists and theorists, via freely distributed codes and databases. Altogether, this project will facilitate establishing XFEL-based multidimensional spectroscopies as a novel diagnostic tool for monitoring electronic and structural dynamics in molecular materials.

74 ATOMIC AND MOLECULAR PHYSICS↗

Instruments for Reading Direct-Marked Data-Matrix Symbols

Improved optoelectronic instruments (specially configured digital cameras) for reading direct-marked data-matrix symbols on the surfaces of optically reflective objects (including specularly reflective ones) are undergoing development. Data-matrix symbols are two-dimensional binary patterns that are used, like common bar codes, for automated identification of objects. The first data-matrix symbols were checkerboard-like patterns of black-and-white rectangles, typically existing in the forms of paint, ink, or detachable labels. The major advantage of direct marking (the marks are more durable than are painted or printed symbols or detachable labels) is offset by a major disadvantage (the marks generated by some marking methods do not provide sufficient contrast to be readable by optoelectronic instruments designed to read black-and-white data-matrix symbols). Heretofore, elaborate lighting, lensing, and software schemes have been tried in efforts to solve the contrast problem in direct-mark matrix- symbol readers. In comparison with prior readers based on those schemes, the readers now undergoing development are expected to be more effective while costing less. All of the prior direct-mark matrix-symbol readers are designed to be aimed perpendicularly to marked target surfaces, and they tolerate very little angular offset. However, the reader now undergoing development not only tolerates angular offset but depends on angular offset as a means of obtaining the needed contrast, as described below. The prototype reader (see Figure 1) includes an electronic camera in the form of a charge-coupled-device (CCD) image detector equipped with a telecentric lens. It also includes a source of collimated visible light and a source of collimated infrared light for illuminating a target. The visible and infrared illumination complement each other: the visible illumination is more useful for aiming the reader toward a target, while the infrared illumination is more useful for reading symbols on highly reflective surfaces. By use of beam splitters, the visible and infrared collimated lights are introduced along the optical path of the telecentric lens, so that the target is illuminated and viewed from the same direction.

Schramm, Harry F.↗

Hyper-Velocity Nanoparticle Plasma Jet as Fast Probe for Runaway Electrons in Tokamak Disruptions

The large energy focalized deposition of the runaway electron (RE) beam on the first wall can potentially yield to destructive effects on the reactor vessel during the explosive instability of a major disruption in a large tokamak (e.g. RE beam estimation is multi-MA of ~10 MeV energy in ITER.) Impurity injection is seen as an effective way to suppress REs after thermal quench (TQ). Injected material in solid state is advantageous, but it must be in small size grain form, with large surface-to-mass ratio to facilitate its rapid ablation within plasma. However, the results using shattered pellet injection (SPI), a technique aiming to satisfy these condition, indicated that the relatively unexplored RE evolution phases and RE beam-plasma interaction are far from being properly diagnosed and fully understood. We proposed the idea to probe REs with a pulsed power hyper-velocity large-mass nanoparticle plasma jet (NPPJ). NPPJ provides a high specific surface area and swift injection to core, hence fast increase of three basic parameters: 1) the electron density ne (bound and free through ionization, if the cold post-TQ plasma has a still high enough T e ), 2) the effective charge (Z eff ) which rises both the dynamical friction force and plasma resistivity, and, finally, 3) the value of the (Dreicer and critical) threshold E-field above which electrons are accelerated and become REs. But the fast transient phenomena of RE ”seed” creation, RE ”avalanche” multiplication, and RE beam formation and interaction with the background plasma all take place dominantly in the core plasma, where the induction E-field has its maximum. Thus, the RE probe success requires a rapid (1 - 2 ms) injection of sufficient mass able to penetrate through the strong toroidal magnetic field (2 - 5 T) over a long distance (~1 - 2 m) with a large assimilation fraction in the core plasma. As a RE probe, NPPJ is a unique combination of fast response with trigger-to-delivery time of few milliseconds, large mass-velocity (>75 mg at several km/s) and deep injection into post-TQ central plasma. The Objectives of this Category #1 project (concept development) were as follows: 1: We investigated the essential basic physics aspects of creating a plasma with ~10 15 to 10 16 cm -3 boron nitride (BN) NP density suited for a pulsed power plasma gun acceleration to several km/s. We concluded that, for RE NPPJ probe, BN (diameter Ø10 to 100 nm NPs) are in all respects usable. 2: We advanced solutions to get BN NPPJ by adapting FAR-TECH’s technology for C 60 NPPJ. 3: We developed self-consistent models for RE current density j RE (R; t) evolution in the TQ aftermath, with fast and simultaneous increase of n e (R; t) and Z eff (R; t) > 1 upon the specific injection of C 60 (i.e. C ions) and BN NP (i.e. B and N ions), respectively. We showed that prompt post-TQ injection and high velocity are of essence in RE probe diagnostics and suppression. 4: We used the in-house HEM-2D code to simulate the C$_{60}^{q+}$ NPPJ transverse penetration through tokamak-level B-field (BT ~2 to 5 T) and distance (~1 to 2m), with fragmentation of C$_{60}^{q+}$ and developed the algorithm for BN Q+ NPPJ ablative sublimation of BN Q+ NP by the collisions with post-TQ plasma electrons. HEM-2D results indicated that fragmenting C$_{60}^{1+}$ NPPJ penetrates up to ~60 cm in 100 us when injecting at 1.6 ms after TQ and gradually released C 2 fragments and C ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Exploring Deuterium-Tritium Gas-Filled Hohlraums as a High-Yield Neutron Source

Inverted-corona targets as neutron sources for high-energy-density applications were tested experimentally using the OMEGA laser. In a corona fusion target, laser beams directly ablate a layer of fusion fuel on the inner surface of a capsule. As the ablative flows converge, the plasma particles transition to collisional stagnation, heating the ions and generating fusion reactions. Based on a recent publication these targets show promise as a neutron source, with extrapolated yields at NIF scale of up to 1e17. Additionally, these targets are robust to low-mode drive asymmetries, enabling the possibility of single-sided laser drive, making them an ideal backlighter source for neutron-radiography applications. A number of different targets were fielded in experiments on the OMEGA laser. The 1.8-mm diameter plastic (CH) capsules featured either one or two laser-entrance holes (LEH’s), and were either lined with deuterated plastic (CD), or filled with 1.5-atm deuterium gas. Contrary to predictions by the hydrodynamics code HYDRA, the experiments exhibit a strong correlation between total neutron yield and CD liner thickness, indicating substantial mix of CH wall material into the ablated CD coronal plasma. The thickest liner targets with approximately 10 kJ of incident laser energy performed consistent with previously published data, while the data at 18 kJ laser energy may also be subject to mix from wall material. The yield from the gas-filled capsules was less than from 10-μm CD lined targets at the same laser energy, which may als be explained by wall material mixing into the fusion fuel. Finally, the experiments demonstrated similar performance between a 1-LEH and a 2-LEH target irradiated with comparable laser-energies, a key finding for demonstrating the utility of the inverted-corona neutron source for radiography applications. Future work will focus on developing the platform for NIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗

Second Report of the Nuclear Data Subcommittee of the Nuclear Science Advisory Committee

The central importance of the nuclear data curated by the US Nuclear Data Program (USNDP) for clean energy generation, national security, nonproliferation, medical applications, and space exploration as well as basic science was described in a prior report issued by the DOE/NSF Nuclear Science Advisory Committee subcommittee on Nuclear Data (NSAC-ND) in September 2022. In this report, we present a set of fourteen (14) recommendations that will enhance and advance DOE-NP's stewardship of nuclear data. The first three recommendations focus on the existing core USNDP capabilities, namely: 1) Support the nuclear structure evaluation workforce to improve the currency, consistency, and accessibility of the Evaluated Nuclear Structure Data File (ENSDF); 2) Enhance nuclear reaction evaluation within the USNDP in support of the Evaluated Nuclear Data File (ENDF) through expansion of the workforce and integration of high-performance computing, automation, and machine learning and; 3) Continue atomic mass evaluation in support AME and NUBASE databases. This is followed by eight (8) recommendations representing new cross-cutting initiatives involving both measurement and evaluation to address outstanding nuclear data needs. These new initiatives require a highly trained, diverse workforce that includes personnel with expertise from both inside and outside the nuclear physics community from which evaluators have traditionally been recruited. As such, many of these initiatives are accomplished via a Topical Nuclear Data Collaborations (TNDC). A TNDC is made up of domestic and international stakeholders, subject matter and nuclear data experts, and nuclear data evaluators and features a workforce development plan to ensure that nuclear data evaluators maintain currency in the relevant applications and are seen as equity partners in the endeavor. These include: 1) Establish a coordinated effort to improve evaluation and modeling in nuclear astrophysics for stellar dynamics, multi-messenger astronomy and nucleosynthesis; 2) Initiate a TNDC to develop and maintain nuclear structure evaluation beyond discrete states, including nuclear level densities, photon strength functions and photonuclear data for improved reaction modeling, and exploring nuclear structure at finite temperature; 3) Create a TNDC to perform correlated fission data evaluation, including cross sections, fragment yields, v(A), v(E n ) for nuclear energy, national security, nonproliferation and basic science; 4) From a panel of subject matter experts to establish and annually update a roster of key decay data to nurture its accelerated dissemination including both measurement and evaluation for targeted high-value nuclides for national security, nonproliferation and medical applications; 5) Comprehensive, consistent neutron-induced structure and reaction data for nuclear energy, national security, nonproliferation and planetary nuclear spectroscopy; 6) Charged-particle stopping powers for detector design, space effects and ion beam therapy; 7) High-energy reactions for space exploration and medical nuclide production, and; 8) The creation of an infrastructure for open data and data preservation for use by the entire nuclear physics community. All told, these initiatives require approximately $6.5M increase in NP support of the USNDP in fiscal year 2023 dollars and would require at least 3-5 years to carry out due to the length of time needed to recruit and train new nuclear data researchers. This relatively modest investment would help ensure that the fruits of the nuclear data research carried out by DOE-NP and its collaborators would be brought to bear to address some of the most important needs of our nation and the world. To ensure effective execution of this plan, we present an overview of recruitment, training, and retention goals for the USNDP, the centerpiece of which is a mutually agreed upon code of conduct. Finally, we identify the facility and instrumentation needed to perform the recommended experimental activities. This includes a short review of target fabrication capabilities, reactors, neutron beam, light- and heavy-stable ion, gamma-ray, high-energy and radioactive ion beam facilities. Lastly, a more complete appendix of experimental facilities previously compiled is included with new input provided for 6 facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Porosity, swelling, and composition evolution in high-burnup monolithic U-Mo fuel

The microstructural progression of very high-burnup (>8 × 10 21 fissions/cm 3 ) monolithic uranium-molybdenum (U-10wt.%Mo) was analyzed, providing crucial insights into the behavior of post-recrystallized nuclear fuel, where scant data exists. Three focused ion beam cuboids sourced from a fuel plate with varying local burnups of 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 were characterized. The porosity and composition of the samples were evaluated to characterize the evolution of the microstructure as a function of fission density and different locations on the fuel plate, while simultaneously isolating plate-specific parameters such as Zr diffusion barrier thickness, hot-isostatic-pressing conditions, enrichment, and reactor conditions. The porosity was segmented, and the three-dimensional distribution of the porosity was extrapolated from the two-dimensional segmentation. The composition was assessed and quantified using energy-dispersive X-ray spectroscopy areal mapping. The porosity fraction increased as a function of the burnup from 27.77±0.51, 35.12±1.54, and 37.71±0.44 % for 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 , respectively. When compared to literature, the porosity volume fraction plateaus at burnups greater than 6 × 10 21 fissions/cm 3 , while the pore size grows linearly as a function of fission density. The number of large pores increased in number density as a function of burnup, while the smallest pores (<0.3 µm) increased up to 9.05 × 10 21 fissions/cm 3 , followed by a decrease at 9.36 × 10 21 fissions/cm 3 . The delamination and cracking in the fuel plate propagated through an interconnected porosity sublayer identified ∼5 µm from the diffusion barrier. The local swelling of the specimens was within or near the prediction bounds of the Robinson-Williams model for local swelling. The fission products, strontium, barium, cerium, and cesium, precipitated into the pores, while neodymium accumulated adjacent to the pores. Furthermore, these findings have direct implications for the development of fuel performance codes and the accurate documentation of the microstructure evolution in high burnup U-Mo, thus enhancing the safety and efficiency of nuclear fuel usage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics and Multiscale Modeling of Coupled Transport of Chloride Ions in Concrete

Chloride ions (Cl−)-induced corrosion is one of the main degradation mechanisms in reinforced concrete (RC) structures. In most situations, the degradation initiates with the transport of Cl− from the surface of the concrete towards the reinforcing steel. The accumulation of Cl− at the steel-concrete interface could initiate reinforcement corrosion once a threshold Cl− concentration is achieved. An accurate numerical model of the Cl− transport in concrete is required to predict the corrosion initiation in RC structures. However, existing numerical models lack a representation of the heterogenous concrete microstructure resulting from the varying environmental conditions and the indirect effect of time dependent temperature and relative humidity (RH) on the water adsorption and Cl− binding isotherms. In this study, a numerical model is developed to study the coupled transport of Cl− with heat, RH and oxygen (O2) into the concrete. The modeling of the concrete microstructure is performed using the Virtual Cement and Concrete Testing Laboratory (VCCTL) code developed by the U.S. National Institute of Standards and Technology (NIST). The concept of equivalent maturation time is utilized to eliminate the limitation of simulating concrete microstructure using VCCTL in specific environmental conditions such as adiabatic. Thus, a time-dependent concrete microstructure, which depends on the hydration reactions coupled with the temperature and RH of the environment, is achieved to study the Cl− transport. Additionally, Cl− binding isotherms, which are a function of the pH of the concrete pore solution, Cl− concentration, and weight fraction of mono-sulfate aluminate (AFm) and calcium-silicate-hydrate (C-S-H), obtained from an experimental study by the same authors are utilized to account for the Cl− binding of cement hydration products. The temperature dependent RH diffusion was considered to account for the transport of Cl− with moisture transport. The temperature and RH diffusion in the concrete domain, composite theory, and Cl− binding and water adsorption isotherms are used in combination, to estimate the ensuing Cl− diffusion field within the concrete. The coupled transport process of heat, RH, Cl−, and O2 is implemented in the Multiphysics Object-Oriented Simulation Environment (MOOSE) developed by the U.S. Idaho National Laboratory (INL). The model was verified and validated using data from multiple experimental studies with different concrete mixture proportions, curing durations, and environmental conditions. Additionally, a sensitivity analysis was performed to identify that the water-to-cement (w/c) ratio, the exposure duration, the boundary conditions: temperature, RH, surface Cl− concentration, Cl− diffusion coefficient in the capillary water, and the critical RH are the important parameters that govern the Cl− transport in RC structures. In a case study, the capabilities of the developed numerical model are demonstrated by studying the complex 2D diffusion of Cl− in a RC beam located in two different climatic regions: warm and humid weather in Galveston, Texas, and cold and dry weather in North Minnesota, Minnesota, subjected to time varying temperature, RH, and surface Cl− concentrations.

composite theory↗

Free-Flying Magnetometer Data System

The Free-Flying Magnetometer (FFM) is an autonomous "sensorcraft" developed at the Jet Propulsion Laboratory (JPL) for the Enstrophy sounding rocket mission. This mission was a collaborative project between the University of New Hampshire, Cornell University and JPL. The science goal of the mission was the study of current filamentation phenomena in the northern auroral region through multipoint measurements of magnetic field. The technical objective of the mission was the proof of concept of the JPL FFM design and the demonstration of an in-situ multipoint measurement technique employing many free-flying spacecraft. Four FFMs were successfully deployed from a sounding rocket launched from Poker Flats, Alaska on February 11, 1999. These hockey-puck-sized (80 mm diameter, 38 mm. height, 250 gram mass) free flyers each carry a miniature 3-axis flux-gate magnetometer that output +/- 2 V signals corresponding to a +/- 60,000 nT measurement range for each axis. The FFM uses a synchronized four-channel Sigma(Delta) Analog-to-Digital Converter (ADC) having a dynamic range of +/- 2.5V and converting at a rate of 279 samples/second/channel. Three channels are used to digitize the magnetometer signals to 17-bit (1.144 nT/bit) resolution. The fourth ADC channel is multiplexed for system monitoring of four temperature sensors and two battery voltages. The FFM also contains two sun sensors, a laser diode which emits a fan-shaped beam, a miniature S-band transmitter for direct communication to the ground station antennas, an ultra-stable Temperature Compensated Crystal Oscillator (TCXO) clock, an integrated data subsystem implemented in a Field-Programmable Gate Array (FPGA), a 4 Mbit Static Random Access Memory (SRAM) for data storage and Lithium Thionyl Chloride batteries for power. Communicating commands to the FFM prior to deployment is achieved with an infrared (IR) link. The FFM IR receiver responds to 9-bit pulse coded signals that are generated by an IR Light Emitting Diode (LED) in the payload for turning FFM power on or off and placing the FFM in a test mode or flight mode. The IR links are also used to synchronize (zero) the clocks onboard all the FFMs through a reset pulse originating from the payload GPS receiver that is issued when the FFMs are in flight mode. The FPGA based data subsystem manages continuous data collection from the four ADC channels and sun sensors, formatting and storing the data to SRAM, and controlling downlink transmission. The transmitter is powered only after a 2547 frame SRAM buffer has been filled (approx. 5 minutes of data). The data is Viterbi encoded and sent to the S-band transmitter via a First-In-First-Out (FIFO) buffer who's output is clocked at 100 bits/second. After the 26-second transmission, the transmitter is turned off to reduce noise coupling to the sensitive magnetometer. The data subsystem control consists of a master state machine that performs data flow management and is interfaced through a prioritized interrupt scheme to state machines that service the ADC, sun sensors and transmitter FIFO. Continuous data collection prevents the missing of data during transmission and provides implicit time tagging of the data acquired by the ADC because of synchronization with the TCXO clock.

Blaes, B.↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

Computational Model Prediction and Biological Validation Using Simplified Mixed Field Exposures for the Development of a GCR Reference Field

The yield of chromosomal aberrations has been shown to increase in the lymphocytes of astronauts after long-duration missions of several months in space. Chromosome exchanges, especially translocations, are positively correlated with many cancers and are therefore a potential biomarker of cancer risk associated with radiation exposure. Although extensive studies have been carried out on the induction of chromosomal aberrations by low- and high-LET radiation in human lymphocytes, fibroblasts, and epithelial cells exposed in vitro, there is a lack of data on chromosome aberrations induced by low dose-rate chronic exposure and mixed field beams such as those expected in space. Chromosome aberration studies at NSRL will provide the biological validation needed to extend the computational models over a broader range of experimental conditions (more complicated mixed fields leading up to the galactic cosmic rays (GCR) simulator), helping to reduce uncertainties in radiation quality effects and dose-rate dependence in cancer risk models. These models can then be used to answer some of the open questions regarding requirements for a full GCR reference field, including particle type and number, energy, dose rate, and delivery order. In this study, we designed a simplified mixed field beam with a combination of proton, helium, oxygen, and iron ions with shielding or proton, helium, oxygen, and titanium without shielding. Human fibroblasts cells were irradiated with these mixed field beam as well as each single beam with acute and chronic dose rate, and chromosome aberrations (CA) were measured with 3-color fluorescent in situ hybridization (FISH) chromosome painting methods. Frequency and type of CA induced with acute dose rate and chronic dose rates with single and mixed field beam will be discussed. A computational chromosome and radiation-induced DNA damage model, BDSTRACKS (Biological Damage by Stochastic Tracks), was updated to simulate various types of CA induced by acute exposures of the mixed field beams used for the experiments. The chromosomes were simulated by a polymer random walk algorithm with restrictions to their respective domains in the nucleus [1]. The stochastic dose to the nucleus was calculated with the code RITRACKS [2]. Irradiation of a target volume by a mixed field of ions was implemented within RITRACKs, and the fields of ions can be delivered over specific periods of time, allowing the simulation of dose-rate effects. Similarly, particles of various types and energies extracted from a pre-calculated spectra of galactic cosmic rays (GCR) can be used in RITRACKS. The number and spatial location of DSBs (DNA double-strand breaks) were calculated in BDSTRACKS using the simulated chromosomes and local (voxel) dose. Assuming that DSBs led to chromosome breaks, and simulating the rejoining of damaged chromosomes occurring during repair, BDSTRACKS produces the yield of various types of chromosome aberrations as a function of time (only final yields are presented). A comparison between experimental and simulation results will be shown.

Hada, M.↗

A New Statistical Estimate of the Radar Coverage of the Low Earth Orbit Debris Environment

For over three decades, the NASA Orbital Debris Program Office (ODPO) has used the Goldstone Orbital Debris Radar, Haystack Ultrawideband Satellite Imaging Radar (HUSIR), and Haystack Auxiliary (HAX) radar assets to collect data on the low Earth orbit (LEO) debris environment. Each radar, with its unique beamwidth, altitude and inclination coverage, and limiting size threshold, operates in a beam park mode to statistically sample the orbital debris population in LEO. Provided that these assets are shared with other users, the orbital debris data collection is not continuous; rather, intermittent data collects are acquired and sent to the NASA ODPO. To understand the sampling process conducted by each radar over time and any related observational biases, a Statistical Confirmation of Radar Uniformity or Bias (SCRUB) code has been implemented to model the coverage of these assets for informing future operations, as well as usage of the data collected from these ground-based sensors. For this analysis, Right Ascension of the Ascending Node (RAAN) is used as a metric to measure statistical coverage. A complete survey of the LEO environment is understood to be measurements that sufficiently sample all values of RAAN for each altitude-inclination pair visible from the radar asset. Regions of incompleteness or statistical bias can help inform future observation campaigns. The SCRUB tool evaluates the coverage of the LEO environment, not by examining individual objects that may pass through a sensor’s field of view (FOV), but by determining which orbit planes pass through the FOV. Once a pointing geometry for a radar site, observation time, and range extent are configured by the user, SCRUB computes the inclinations and altitudes that are visible by the sensor. For each inclination-altitude pair, there is a distinct pair of possible RAAN values, corresponding to the ascending and descending orbit passing through that point in the sensor cone. Repeating this process for all points in the beam, and for multiple time periods during an observation window, creates a matrix of all inclination-altitude-RAAN combinations that are visible during a sensor run. This process can then be repeated for all observations within a year (for an annual survey), propagating all the RAAN values to a common epoch, typically the start of a calendar year, and combining the observations to assemble a full estimate of the RAAN coverage of the LEO environment. This RAAN coverage can then be analyzed for uniformity of sampling, within a certain inclination-altitude pair, or between larger regions of the space environment. This paper provides a general overview of the radar assets utilized by ODPO, typical analysis data products assuming circular orbits, and a discussion of the algorithms that feed between modeling and measurement operations. Following the description of the algorithm, estimates of coverage using HUSIR radar data collected from multiple years are developed and compared.

Chris L Ostrom↗

Simple Laser Communications Terminal for Downlink from Earth Orbit at Rates Exceeding 10 Gb/s

A compact, low-cost laser communications transceiver was prototyped for downlinking data at 10 Gb/s from Earth-orbiting spacecraft. The design can be implemented using flight-grade parts. With emphasis on simplicity, compactness, and light weight of the flight transceiver, the reduced-complexity design and development approach involves: 1. A high-bandwidth coarse wavelength division multiplexed (CWDM) (4 2.5 or 10-Gb/s data-rate) downlink transmitter. To simplify the system, emphasis is on the downlink. Optical uplink data rate is modest (due to existing and adequate RF uplink capability). 2. Highly simplified and compact 5-cm diameter clear aperture optics assembly is configured to single transmit and receive aperture laser signals. About 2 W of 4-channel multiplexed (1,540 to 1,555 nm) optically amplified laser power is coupled to the optical assembly through a fiber optic cable. It contains a highly compact, precision-pointing capability two-axis gimbal assembly to coarse point the optics assembly. A fast steering mirror, built into the optical path of the optical assembly, is used to remove residual pointing disturbances from the gimbal. Acquisition, pointing, and tracking are assisted by a beacon laser transmitted from the ground and received by the optical assembly, which will allow transmission of a laser beam. 3. Shifting the link burden to the ground by relying on direct detection optical receivers retrofitted to 1-m-diameter ground telescopes. 4. Favored mass and volume reduction over power-consumption reduction. The two major variables that are available include laser transmit power at either end of the link, and telescope aperture diameter at each end of the link. Increased laser power is traded for smaller-aperture diameters. 5. Use of commercially available spacequalified or qualifiable components with traceability to flight qualification (i.e., a flight-qualified version is commercially available). An example is use of Telecordia-qualified fiber optic communication components including active components (lasers, amplifiers, photodetectors) that, except for vacuum and radiation, meet most of the qualifications required for space. 6. Use of CWDM technique at the flight transmitter for operation at four channels (each at 2.5 Gb/s or a total of 10 Gb/s data rate). Applying this technique allows utilization of larger active area photodetectors at the ground station. This minimizes atmospheric scintillation/turbulence induced losses on the received beam at the ground terminal. 7. Use of forward-error-correction and deep-interleaver codes to minimize atmospheric turbulence effects on the downlink beam. Target mass and power consumption for the flight data transmitter system is less than 10 kg and approximately 60 W for the 400-km orbit (900-km slant range), and 12 kg and 120 W for the 2,000-km orbit (6,000-km slant range). The higher mass and power for the latter are the result of employing a higher-power laser only.

Kovalik, Joseph M.↗