Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Geometry”

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 955 records · Page 53

Angular-resolution and material-characterization measurements for a dual-particle imaging system with mixed-oxide fuel

A dual-particle imaging (DPI) system, capable of simultaneously imaging fast neutrons and gamma rays, has been operated in the presence of mixed-oxide (MOX) fuel to assess the system's angular resolution and material-characterization capabilities. The detection principle is based on the scattering physics of neutrons (elastic scattering) and gamma rays (Compton scattering) in organic and inorganic scintillators. The detection system is designed as a combination of a two-plane Compton camera and a neutron-scatter camera. The front plane consists of EJ-309 liquid scintillators and the back plane consists of interleaved EJ-309 and NaI(Tl) scintillators. MCNPX-PoliMi was used to optimize the geometry of the system and the resulting prototype was built and tested using a Cf-252 source as an SNM surrogate. A software package was developed to acquire and process data in real time. The software was used for a measurement campaign to assess the angular resolution of the imaging system with MOX samples. Measurements of two MOX canisters of similar isotopics and intensity were performed for 6 different canister separations (from 5° to 30°, corresponding to distances of 21 cm and 131 cm, respectively). The measurements yielded a minimum separation of 20° at 2.5 m (86-cm separation) required to see 2 separate hot spots. Additionally, the results displayed good agreement with MCNPX-PoliMi simulations. These results indicate an angular resolution between 15° and 20°, given the 5° step. Coupled with its large field of view, and its capability to differentiate between spontaneous fission and (α,n) sources, the DPI system shows its potential for nuclear-nonproliferation applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A 3D-printed microbial cell culture platform with in situ PEGDA hydrogel barriers for differential substrate delivery

Additive manufacturing, or 3D-printing techniques have recently begun to enable simpler, faster, and cheaper production of millifluidic devices at resolutions approaching 100–200 μ m. At this resolution, cell culture devices can be constructed that more accurately replicate natural environments compared with conventional culturing techniques. A number of microfluidics researchers have begun incorporating additive manufacturing into their work, using 3D-printed devices in a wide array of chemical, fluidic, and even some biological applications. In this work, we describe a 3D-printed cell culture platform and demonstrate its use in culturing Pseudomonas putida KT2440 bacteria for 44 h under a differential substrate gradient. Polyethylene glycol diacrylate (PEGDA) hydrogel barriers are patterned in situ within a 3D-printed channel. Transport of the toluidine blue tracer dye through the hydrogel barriers is characterized. Nutrients and oxygen were delivered to cells in the culture region by diffusion through the PEGDA hydrogel barriers from adjacent media or saline perfusion channels. Expression of green fluorescent protein by P. putida KT2440 enabled real time visualization of cell density within the 3D-printed channel, and demonstrated cells were actively expressing protein over the course of the experiment. Cells were observed clustering near hydrogel barrier boundaries where fresh substrate and oxygen were being delivered via diffusive transport, but cells were unable to penetrate the barrier. The device described here provides a versatile and easy to implement platform for cell culture in readily controlled gradient microenvironments. By adjusting device geometry and hydrogel properties, this platform could be further customized for a wide variety of biological applications.

42 ENGINEERING↗

FLUX EMERGENCE IN A MAGNETIZED CONVECTION ZONE

We study the influence of a dynamo magnetic field on the buoyant rise and emergence of twisted magnetic flux ropes and their influence on the global external magnetic field. We ran three-dimensional MHD numerical simulations using the ASH code (anelastic spherical harmonics) and analyzed the dynamical evolution of such buoyant flux ropes from the bottom of the convection zone until the post-emergence phases. The global nature of this model can only very crudely and inaccurately represent the local dynamics of the buoyant rise of the implanted magnetic structure, but nonetheless allows us to study the influence of global effects, such as self-consistently generated differential rotation and meridional circulation, and of Coriolis forces. Although motivated by the solar context, this model cannot be thought of as a realistic model of the rise of magnetic structures and their emergence in the Sun, where the local dynamics are completely different. The properties of initial phases of the buoyant rise are determined essentially by the flux-rope's properties and the convective flows and consequently are in good agreement with previous studies. However, the effects of the interaction of the background dynamo field become increasingly strong as the flux ropes evolve. During the buoyant rise across the convection zone, the flux-rope's magnetic field strength scales as B{proportional_to}{rho}{sup {alpha}}, with {alpha} {approx}< 1. An increase of radial velocity, density, and current density is observed to precede flux emergence at all longitudes. The geometry, latitude, and relative orientation of the flux ropes with respect to the background magnetic field influences the resulting rise speeds, zonal flow amplitudes (which develop within the flux ropes), and the corresponding surface signatures. This influences the morphology, duration and amplitude of the surface shearing, and the Poynting flux associated with magnetic flux-rope emergence. The emerged magnetic flux influences the system's global polarity, leading in some cases to a polarity reversal while inhibiting the background dynamo from doing so in others. The emerged magnetic flux is slowly advected poleward while being diffused and assimilated by the background dynamo field.

79 ASTRONOMY AND ASTROPHYSICS↗

Design and Construction of a High-Resolution Hodoscope for the GlueX Experiment with High-Statistics Analysis of the p0, ¿, and ¿1 Photoproduction Cross Sections from the RadPhi Experiment

Differential cross sections for forward-angle photoproduction of p0, ¿, and ¿ 1 pseudoscalar mesons were measured using data from the RadPhi experiment conducted in Hall B at Jef ferson Lab. RadPhi utilized a tagged bremsstrahlung photon beam incident on a stationary 9Be target, with a detector system configured to trigger on a recoil proton in coincidence with multiple neutral showers in the calorimeter. Events were reconstructed and subjected to kinematic constraints, with background suppressed via sideband subtraction guided by Monte Carlo modeling of background contributions. Cross sections were extracted over the photon energy range 4.4– 5.4 GeV and binned in invariant momentum transfer t, providing measurements from one of the first high-statistics experiments of forward ¿ and ¿1 pro duction from a nuclear target at these energies. Acceptance corrections were applied using a detailed GEANT-based simulation of the detector geometry and response. The resulting cross sections are consistent with 2020 CLAS results, when scaled by the number of protons in beryllium, and show broad agreement with other data and theoretical models. In parallel, a high-resolution photon tagger detector, the Tagger Microscope (TAGM), was designed, constructed, and commissioned for the GlueX experiment in Hall D at Jefferson Lab. The TAGM was developed to provide high-rate tagging capability in the coherent bremsstrahlung peak by detecting post-bremsstrahlung electrons across a one GeV range along the focal plane of the tagging spectrometer. The detector consists of a 5ˆ102 array of 2ˆ2 mm2 square BCF-20 plastic scintillating fibers thermally fused to BCF-98 light guide fibers optically coupled to silicon photomultipliers. These fibers are mounted in a precision machined framework enabling fine positional adjustments to maintain precise alignment with post-bremsstrahlung electron trajectories, while ensuring mechanical rigidity, thermal stability, optical isolation, minimal inactive area, and radiation shielding for electronics. The construction effort involved extensive testing of fiber quality, light transmission, thermal fusing, radiation hardness, and defect analysis using SEM and EDX techniques. Following its installation and commissioning, the TAGM became a critical component of the GlueX beamline, enabling high-rate tagging essential for studies of hybrid mesons and gluonic ex citations.

McIntyre, James [Univ. of Connecticut, Storrs, CT ↗

Geometric treatment of conduction electron scattering by crystal lattice strains and dislocations

The problem of conduction electron scattering by inhomogeneous crystal lattice strains is addressed using a tight-binding formalism and the differential geometric treatment of deformations in solids. In this approach, the relative positions of neighboring atoms in a strained lattice are naturally taken into account, even in the presence of crystal dislocations, resulting in a fully covariant Schrödinger equation in the continuum limit. Unlike previous work, the developed formalism is applicable to cases involving purely elastic strains as well as discrete and continuous distributions of dislocations—in the latter two cases, it clearly demarcates the effects of the dislocation strain field and core. It also differentiates between elastic and plastic strain contributions, respectively. The electrical resistivity due to the strain field of edge dislocations is then evaluated and the resulting numerical estimate for Cu shows good agreement with reported experimental values. This indicates that the electrical resistivity of edge dislocations in metals is not entirely due to the core, contrary to current models. Application to the study of strain effects in constrained quantum systems is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

ME8373 Spring 2015 ICME Proposal ICME Analysis of Fatigue Crack Growth Through a Weld in SA-516 Grade 70 Plate

This paper describes work accomplished to predict the service life of a flexure joint design which is a component of a diffuser duct in the A3 Test Stand, an altitude simulation rocket engine test facility at NASA's Stennis Space Center. The duct has two pressure shells separated by cooling water passages and connected by stiffening ribs and flexure joints. Rocket exhaust flows within the duct and heats the inner pressure shell while the outer pressure shell remains at ambient temperature. The flexure joints allow for differential thermal expansion of the inner and outer pressure shells and are subject to in-service loading by this thermal expansion along with water pressure in the cooling water passage, atmospheric pressure outside the duct, near vacuum conditions within the duct, and vibrational loads from operation of the facility and rocket engine. Figure 1 shows a schematic axisymmetric cross section of the diffuser pressure shells and flexure joints with a zoomed in view of the flexure joint. The flexure joints are expected to eventually fail by fatigue cracking leading to leaks from the cooling water passages to the outside. The zoomed in view in Figure 1 indicates where cracking is expected to occur, namely through a weld bead between two plates of SA-516 Grade 70 steel. This weld bead acts as the fulcrum of the flexure joint and it is clear from inspection of the geometry and loading represented in the zoomed in portion of Figure 1 that inherent in the design there is a severe notch formed between the flexure plate, weld bead, and stiffening ring that will be the site of crack initiation and location from which the crack grows to the outer surface of the weld bead.

Woods, Jody L.↗

Theoretical studies of the outer envelopes of young stellar objects

With the Monte Carlo code developed by Whitney and Hartmann, a series of models was computed of scattering in disks around young stellar objects. The code calculates scattering by dust, including polarization, in arbitrary geometries. By computing model images, it was found that disk, by themselves, around young stellar objects would be very difficult to detect with present day imaging techniques. In comparing these images to observations of young stellar objects which show diffuse structure, little resemblance was found. A flared disk system will only give high polarization when viewed edge-on, and the position angle is always oriented perpendicular to the disk plane. This suggests that an envelope, perhaps the remnant infalling envelope, must be present to scatter more stellar light than a disk can, and obscure the star at many inclinations. A grid was computed of models of scattering in a disk+envelope system. Evidence is presented that the wind of the pre-main sequence object FU Orionis arises from the surface of the luminous prostellar accretion disk. A disk wind model calculated assuming radiative equilibrium explains the differential behavior of the observed asymmetrical absorption line profiles. The model predicts that strong lines should be asymmetric and blueshifted, while weak lines should be symmetric and doubled peaked due to disk rotation, in agreement with observations.

Hartmann, Lee↗

High-resolution sub-Doppler infrared spectroscopy of atmospherically relevant Criegee precursor CH 2 I radicals: CH 2 stretch vibrations and “charge-sloshing” dynamics

The combination of a pulsed supersonic slit-discharge source and single-mode difference frequency direct absorption infrared spectroscopy permit first high resolution infrared study of the iodomethyl (CH 2 I) radical, with the CH 2 I radical species generated in a slit jet Ne/He discharge and cooled to 16 K in the supersonic expansion. Dual laser beam detection and collisional collimation in the slit expansion yield sub-Doppler linewidths (60 MHz), an absolute frequency calibration of 13 MHz, and absorbance sensitivities within a factor of two of the shot-noise limit. Fully rovibrationally resolved direct absorption spectra of the CH 2 symmetric stretch mode (ν 2 ) are obtained and fitted to a Watson asymmetric top Hamiltonian with electron spin-rotation coupling, providing precision rotational constants and spin-rotation tensor elements for the vibrationally excited state. Assessment of the asymmetric top rotational constants confirms a vibrationally averaged planar geometry in both the ground- and first-excited vibrational levels. Sub-Doppler resolution permits additional nuclear spin hyperfine structures to be observed, with splittings in excellent agreement with microwave measurements on the ground state. Spectroscopic data on CH 2 I facilitate systematic comparison with previous studies of halogen-substituted methyl radicals, with the periodic trends strongly correlated with the electronegativity of the halogen atom. Surprisingly, we do not observe any asymmetric CH 2 stretch transitions, despite S/N ≈ 25:1 on strongest lines in the corresponding symmetric CH 2 stretch manifold. This dramatic reversal of the more typical 3:1 antisymmetric/symmetric CH 2 stretch intensity ratio signals a vibrational transition moment poorly described by simple “bond-dipole” models. Instead, the data suggest that this anomalous intensity ratio arises from “charge sloshing” dynamics in the highly polar carbon-iodine bond, as supported by ab initio electron differential density plots and indeed consistent with observations in other halomethyl radicals and protonated cluster ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Algebraic Structures in the Coupling of Gravity to Gauge Theories

Highlights: • Geometric foundation for Quantum General Relativity coupled to Quantum Electrodynamics. • Rigorous introduction to Hopf algebraic renormalization. • Appropriate modifications in case renormalization Hopf algebra is ill-defined. • Generalization of Furry’s Theorem to include gravitons and ghosts. • Discussion of obstructions to multiplicative renormalization. This article is an extension of the author’s second master thesis, Prinz (2017). It aims to introduce to the theory of perturbatively quantized General Relativity coupled to Spinor Electrodynamics, provide the results thereof and set the notation to serve as a starting point for further research in this direction. It includes the differential geometric and Hopf algebraic background, as well as the corresponding Lagrange density and some renormalization theory. Then, a particular problem in the renormalization of Quantum General Relativity coupled to Quantum Electrodynamics is addressed and solved by a generalization of Furry’s Theorem. Next, the restricted combinatorial Green’s functions for all two-loop propagators and all one-loop divergent subgraphs thereof are presented. Finally, relations between these one-loop restricted combinatorial Green’s functions necessary for multiplicative renormalization are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Parametric studies of stitching effectiveness for preventing substructure disbond

A methodology is desired that will allow a designer to select appropriate amounts of through-thickness reinforcement needed to meet design requirements. The goal is to use a relatively simple analysis to minimize the amount of testing that needs to be performed, and to make test results from simple configurations applicable to more general structures. Using this methodology one should be able to optimize the selection of stitching materials, the weight of the yarn, and the stitching density. The analysis approach is to treat substructure disbond as a crack propagation problem. In this approach, the stitches have little influence until a delamination begins to grow. Once the delamination reaches, or extends beyond a stitch, the stitch serves to reduce the strain-energy-release-rate (G) at the crack tip for a given applied load. The reduced G can then be compared to the unstitched materials toughness to predict the load required to further extend the crack. The current model treats the stitch as a simple spring which responds to displacements in the vertical (through-thickness) direction. In concept, this approach is similar to that proposed by other authors. Test results indicate that the model should be refined to include the shearing stiffness of the stitch. The strain-energy-release-rate calculations are performed using a code which uses interconnected higher-order plates to model built-up composite cross-sections. When plates are stacked vertically, the interfacial tractions between the plates can be computed. The plate differential equations are solved in closed-form. The code, called SUBLAM, was developed as part of this section in one dimension. Because of this limitation, rows of stitches are treated as a two-dimensional sheet. The spring stiffness of a row of stitches can be estimated from the stitch material, weight, and density. As a practical and conservative approach, we can assume that the stitch is bonded until a crack passes the stitch location. After the crack passes, it is fully bonded. A series of tests were performed to exercise this methodology and incorporated an attached flange such that the sudden change in thickness initiated a delamination. The analysis was used to estimate the materials' critical G from that of the unstitched specimens. With this data, a prediction was made for the load required to delaminate the stitched specimens. Using the methodology, design charts have been created for simplified geometries. These charts give stitch force, along with G(sub 1) and G(sub 2) as as function of the stitch spring stiffness. Using the charts, it should be possible to determine the stitch spring stiffness and strength required to reduce the G to a desired level. From these parameters, the actual stitching material, weight, and density can be computed.

Flanagan, Gerry↗

An Investigation of Differential Deposition for Figure Corrections in Full-Shell Grazing-Incidents X-Ray Optics

We are investigating differential deposition as a way of correcting small figure errors inside full-shell grazing-incidence x-ray optics. The optics in our study are fabricated using the electroformed-nickel-replication technique, and the figure errors arise from fabrication errors in the mandrel, from which the shells are replicated, as well as errors induced during the electroforming process. Combined, these give sub-micron-scale figure deviations which limit the angular resolution of the optics to approx. 10 arcsec. Sub-micron figure errors can be corrected by selectively depositing (physical vapor deposition) material inside the shell. The requirements for this filler material are that it must not degrade the ultra-smooth surface finish necessary for efficient x-ray reflection (approx. 5 A rms), and must not be highly stressed. In addition, a technique must be found to produce well controlled and defined beams within highly constrained geometries, as some of our mirror shells are less than 3 cm in diameter.

Gubarev, Mikhail V.↗

SU-F-T-261: Reconstruction of Initial Photon Fluence Based On EPID Images

Purpose: Verifying an algorithm to reconstruct relative initial photon fluence for clinical use. Clinical EPID and CT images were acquired to reconstruct an external photon radiation treatment field. The reconstructed initial photon fluence could be used to verify the treatment or calculate the applied dose to the patient. Methods: The acquired EPID images were corrected for scatter caused by the patient and the EPID with an iterative reconstruction algorithm. The transmitted photon fluence behind the patient was calculated subsequently. Based on the transmitted fluence the initial photon fluence was calculated using a back-projection algorithm which takes the patient geometry and its energy dependent linear attenuation into account. This attenuation was gained from the acquired cone-beam CT or the planning CT by calculating a water-equivalent radiological thickness for each irradiation direction. To verify the algorithm an inhomogeneous phantom consisting of three inhomogeneities was irradiated by a static 6 MV photon field and compared to a reference flood field image. Results: The mean deviation between the reconstructed relative photon fluence for the inhomogeneous phantom and the flood field EPID image was 3% rising up to 7% for off-axis fluence. This was probably caused by the used clinical EPID calibration, which flattens the inhomogeneous fluence profile of the beam. Conclusion: In this clinical experiment the algorithm achieved good results in the center of the field while it showed high deviation of the lateral fluence. This could be reduced by optimizing the EPID calibration, considering the off-axis differential energy response. In further progress this and other aspects of the EPID, eg. field size dependency, CT and dose calibration have to be studied to realize a clinical acceptable accuracy of 2%.

60 APPLIED LIFE SCIENCES↗

GQLink: an implementation of Quantized State Systems (QSS) methods in Geant4

Simulations in high energy physics (HEP) often require the numerical solution of ordinary differential equations (ODE) to determine the trajectories of charged particles in a magnetic field when particles move throughout detector volumes. Each crossing of a volume interrupts the underlying numerical method that solves the equations of motion, triggering iterative algorithms to estimate the intersection point within a given accuracy. The computational cost of this procedure can grow significantly depending on the application at hand. Quantized State System (QSS) is a recent family of discrete-event driven numerical methods exhibiting attractive features for this type of problems, such as native dense output (sequences of polynomial segments updated only by accuracy-driven events) and lightweight detection and handling of volume crossings. In this work we present GQLink, a proof-of-concept integration of QSS with the Geant4 simulation toolkit which stands as an interface for co-simulation that orchestrates robustly and transparently the interaction between the QSS simulation engine and aspects such as geometry definition and physics processes controlled by Geant4. In conclusion, we validate the accuracy and study the performance of the method in simple geometries (subject to intense volume crossing activity) and then in a realistic HEP application using a full CMS detector configuration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Examining Graphite Degradation in Molten Salt Environments: A Chemical, Physical, and Material Analysis

Molten-salt reactors (MSRs) are Generation IV nuclear reactors that use liquid salt as a coolant and/or fuel. In several MSR designs, graphite serves as a moderator and/or reflector. However, due to limited experimental data and operational experience, our understanding of graphite behavior in molten salt environments remains incomplete. This report aims to identify the degradation mechanisms of nuclear graphite in MSRs, detail the mechanisms of each factor, and provide an initial assessment of their impact on the structural integrity of graphite components. This assessment is based on an extensive literature review and insights from subject matter experts. Furthermore, given the limited data, a modeling strategy using existing Grizzly software is proposed for a more thorough analysis where appropriate. Additionally, it presents mitigation strategies where applicable. The report covers physical degradation mechanisms such as infiltration, erosion, and abrasion, as well as chemical degradation mechanisms including fluorination, intercalation, corrosion, and oxidation. Molten salt can infiltrate the porous structure of graphite, leading to several detrimental effects. Entrapment of fissile products within the graphite pores can cause radiation damage and could pose challenges in the handling and disposal of contaminated components. The differential thermal expansion between the infiltrated salt and graphite, along with internal stress from pressurized molten salt and volumetric heating, can compromise the structural integrity of graphite. To mitigate these effects, employing ultra-fine graphite grades and applying sealants and coatings are effective strategies. A computational model based on coupled solid mechanics and heat transfer phenomena could be used to predict the internal stresses using Grizzly software. In pebble-bed MSRs, graphite fuel pebbles can cause abrasion against reactor components due to friction and wear. The severity of wear is influenced by various factors such as temperature, environment, and the presence of lubricants. Tribological studies reveal that higher temperatures and molten salt environments, such as FLiBe, significantly reduce wear rates compared to dry conditions. Additionally, the chemical composition of the salt can further optimize graphite's tribological performance. Long-term wear effects can be modeled by incorporating surface defects into the geometry and predict stresses under thermal and radiation effects using Grizzly software. Chemical degradation of graphite in a molten salt environment can occur through fluorination and intercalation. Fluorination can occur via replacement of hydrogen or oxygen atoms, or at the active sites, but does not cause structural degradation. Intercalation, on the other hand, can lead to exfoliation, where layers of graphite separate and peel away, damaging the graphite. Protective coatings can enhance graphite's resistance to intercalation. Graphite generally exhibits good chemical stability in molten salt environments, though it can corrode under specific conditions, particularly in the presence of impurities or oxidants. Studies have shown that protective coatings, such as plasma-sprayed partially stabilized zirconia (PSZ), can effectively prevent such degradation. Corrosion behavior varies significantly with different graphite grades and coating applications, underscoring the need for detailed studies on uncoated and coated graphite to understand and mitigate corrosion mechanisms in MSRs. Research indicates that the presence of oxidants and impurities can accelerate graphite degradation in molten salts, making it essential to explore acceptable impurity limits. Oxidation is another critical degradation mechanism, leading to weight loss and structural damage due to the formation of CO and CO 2 from the reaction of carbon atoms with oxygen. This process creates new porosity and compromises graphite's integrity. While extensive research on graphite oxidation has been conducted for gas-cooled reactors, studies specific to MSRs are limited. Findings from the coal industry suggest that molten alkali metal salts can significantly accelerate graphite oxidation, a hypothesis worth exploring for fluoride salts in MSRs. Understanding oxidation behavior in MSRs is vital for developing protective measures. The analysis of post-irradiated graphite from the MSRE experiment demonstrated exceptional chemical compatibility with molten fluoride salt, suggesting that the extent of chemical attack on graphite largely depends on the salt's infiltration capability. Therefore, the use of ultra-fine grade graphite could help mitigate chemical degradation effects. Existing oxidation modeling capabilities in Grizzly, which use reaction-diffusion equations to model graphite-air interactions, could be adapted to simulate the chemical degradation effects of graphite in molten salt environments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational modeling of graphite degradation in molten salt reactors: Role of infiltration

Molten salt reactors (MSRs) often employ graphite as a moderator and reflector. An important challenge for deploying graphite in these reactors is that, due to limited experimental data, our understanding of graphite’s structural integrity in molten salt environments remains incomplete. Here, this study addresses heat generation from fuel-bearing salt that has infiltrated open pores in the graphite, driven primarily by pressure differentials. This is one of multiple identified physical and chemical mechanisms through which molten salt could potentially degrade graphite. Thermally driven stresses are quantified using the Molten-Salt Reactor Experiment (MSRE) graphite moderator elements as a case study. Finite element simulations predict stress distributions at varying infiltration levels, indicating that thermal stresses increase with higher infiltration. Rare-event simulations using the parallel subset simulation framework identify the combinations and corresponding ranges of input parameters that lead to stresses above a specified threshold. In particular, combinations involving high infiltration amounts, high power density, and low thermal conductivity tend to induce the highest stresses. Under the inputs and assumptions considered in this work, the magnitudes of the thermally driven stresses are quite low, with a very low likelihood of causing failure due to exceeding the graphite’s tensile strength. Additionally, rare-event simulations were performed for two more scenarios: a scaled-up moderator geometry and a localized hotspot in the original geometry. Both cases resulted in increased susceptibility to failure, though not to a detrimental extent. Furthermore, the combined effects of irradiation and infiltration-induced thermal stresses were evaluated. The results showed that thermal stresses from infiltration were negligible compared to those caused by irradiation. The findings of such a study are inherently component-specific, but the methodology presented here could be used for similar assessments of salt-infiltration effects in other graphite components.

36 - MATERIALS SCIENCE↗

Dynamic System Simulation of Fissile Solution Systems

Dynamic System Simulation (DSS) utilizes state-variables described by differential or difference equations to model system evolution. DSS techniques have been applied to develop a family of models to examine the time-dependent operational behavior of fissile solution systems. The base, or generic model, predicts the behavior of an aqueous homogeneous reactor (AHR), configured with a single cooling loop, operating at atmospheric pressure. A model specific for SUPO (Super Power), a uranium solution fueled AHR that operated at Los Alamos National Laboratory from 1951 to 1974, tailored from the generic model demonstrated close correlation with experimental data for steady-state operation1. SUPO is considered the benchmark for steady-state AHR operation. Subsequently specific models were tailored from the generic model for a variety of historical reactors including KEWB (Kinetics Experiments Water Boiler) “A-2” and “B-5” cores and Silene, which is considered the benchmark for AHR pulse operations. These models demonstrated DSS techniques could reliably be extended to different core geometries (SUPO and KEWB “A-2” were spherical; KEWB “B-5” and Silene cylindrical)2. This family of models also shows close correlation with experimental data in all modes of fissile solution system operation including pulse, free evolution, and steady-state. System response due to rate and amplitude of reactivity insertion closely matches experimental data. This report describes additional extensions of this family of models to include a wider design space for fissile solution systems. Pressurized cores, cooling schemes involving multiple cooling loops of various geometries, and accelerator-driven sub-critical system concepts are modeled.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Application of the core shell model for strengthening polymer filament interfaces

While surface segregation has been intensively studied in thin films, relatively little has been investigated in filament geometries. In contrast to thin films, the filament surface forms a continuum, with interactions both at the free air and internal filament interface, which makes it difficult to achieve an equilibrium configuration. This asymmetry has become particularly apparent with the increasing use of additive manufacturing where entire structures are constructed via layered deposition of filamentous layers from a moving thermal nozzle. In addition, this process is highly non-equilibrated with poor thermal retention and is responsible for insufficient filament–filament interfacial fusion and internal pore formation along the printing direction. These features all contribute to poor mechanical and structural properties, which are exacerbated in semi-crystalline polymers, such as Polylactic Acid (PLA), where the melting point, T m , is much higher than the glass transition, T g . Even though system is liquid the molecules become ordered within the crystal structure and this ordering interferes with interdiffusion, which in turn prevents proper interfacial fusion. Capitalizing on the differential in surface tension, we show that addition of less than 1%, or trace amounts, of longer chain polymers, an oriented core shell structure, will be formed with a self-limiting thickness of three time the radius of gyration, R g , of the tracer chains. The lower T g of the tracers and their orientation makes them poised for interdiffusion when subsequent layers are deposited producing structures with significant improvement of the mechanical properties and internal structural integrity.

36 MATERIALS SCIENCE↗

SU-F-T-97: Outlier Identification in Radiation Therapy Knowledge Modeling

Purpose: To investigate the impact of outliers on knowledge modeling in radiation therapy, and develop a systematic workflow for identifying and analyzing geometric and dosimetric outliers using pelvic cases. Methods: Four groups (G1-G4) of pelvic plans were included: G1 (37 prostate cases), G2 (37 prostate plus lymph node cases), and G3 (37 prostate bed cases) are all clinical IMRT cases. G4 are 10 plans outside G1 re-planned with dynamic-arc to simulate dosimetric outliers. The workflow involves 2 steps: 1. identify geometric outliers, assess impact and clean up; 2. identify dosimetric outliers, assess impact and clean up.1. A baseline model was trained with all G1 cases. G2/G3 cases were then individually added to the baseline model as geometric outliers. The impact on the model was assessed by comparing leverage statistic of inliers (G1) and outliers (G2/G3). Receiver-operating-characteristics (ROC) analysis was performed to determine optimal threshold. 2. A separate baseline model was trained with 32 G1 cases. Each G4 case (dosimetric outliers) was then progressively added to perturb this model. DVH predictions were performed using these perturbed models for remaining 5 G1 cases. Normal tissue complication probability (NTCP) calculated from predicted DVH were used to evaluate dosimetric outliers’ impact. Results: The leverage of inliers and outliers was significantly different. The Area-Under-Curve (AUC) for differentiating G2 from G1 was 0.94 (threshold: 0.22) for bladder; and 0.80 (threshold: 0.10) for rectum. For differentiating G3 from G1, the AUC (threshold) was 0.68 (0.09) for bladder, 0.76 (0.08) for rectum. Significant increase in NTCP started from models with 4 dosimetric outliers for bladder (p<0.05), and with only 1 dosimetric outlier for rectum (p<0.05). Conclusion: We established a systematic workflow for identifying and analyzing geometric and dosimetric outliers, and investigated statistical metrics for detecting. Results validated the necessity for outlier detection and clean-up to enhance model quality in clinical practice. Research Grant: Varian master research grant.

60 APPLIED LIFE SCIENCES↗