Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bundle method”

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 199 records · Page 11

Distributed Interplanetary Delay/Disruption Tolerant Network (DTN) Monitor and Control System

The main purpose of Distributed interplanetary Delay Tolerant Network Monitor and Control System as a DTN system network management implementation in JPL is defined to provide methods and tools that can monitor the DTN operation status, detect and resolve DTN operation failures in some automated style while either space network or some heterogeneous network is infused with DTN capability. In this paper, "DTN Monitor and Control system in Deep Space Network (DSN)" exemplifies a case how DTN Monitor and Control system can be adapted into a space network as it is DTN enabled.

Bundle Status Report↗

Analysis of the progressive failure of brittle matrix composites

This report investigates two of the most common modes of localized failures, namely, periodic fiber-bridged matrix cracks and transverse matrix cracks. A modification of Daniels' bundle theory is combined with Weibull's weakest link theory to model the statistical distribution of the periodic matrix cracking strength for an individual layer. Results of the model predictions are compared with experimental data from the open literature. Extensions to the model are made to account for possible imperfections within the layer (i.e., nonuniform fiber lengths, irregular crack spacing, and degraded in-situ fiber properties), and the results of these studies are presented. A generalized shear-lag analysis is derived which is capable of modeling the development of transverse matrix cracks in material systems having a general multilayer configuration and under states of full in-plane load. A method for computing the effective elastic properties for the damaged layer at the global level is detailed based upon the solution for the effects of the damage at the local level. This methodology is general in nature and is therefore also applicable to (0(sub m)/90(sub n))(sub s) systems. The characteristic stress-strain response for more general cases is shown to be qualitatively correct (experimental data is not available for a quantitative evaluation), and the damage evolution is recorded in terms of the matrix crack density as a function of the applied strain. Probabilistic effects are introduced to account for the statistical nature of the material strengths, thus allowing cumulative distribution curves for the probability of failure to be generated for each of the example laminates. Additionally, Oh and Finney's classic work on fracture location in brittle materials is extended and combined with the shear-lag analysis. The result is an analytical form for predicting the probability density function for the location of the next transverse crack occurrence within a crack bounded region. The results of this study verified qualitatively the validity of assuming a uniform crack spacing (as was done in the shear-lag model).

Thomas, David J.↗

Integral Field Probe for Mapping of Curved Magnets

The Single Stretched Wire (SSW) method allows highly precise integral field measurements by recording voltage across a tensioned wire mounted to 2-axis linear stages at either end of the magnet aperture. However, traditional SSW probes are not well suited for curved accelerator magnets, which are essential for steering charged particles along arced trajectories in storage rings or beamlines. The tension required to eliminate sag demands a purely straight path, making them incompatible with non-linear magnet geometries. To address this limitation for curved magnets, a modified approach was developed using a segmented, 3D-printed support structure that incorporates a pre-shaped anti-sag curve. Under its own weight and that of the wire bundle, the structure deforms to lie flat while conforming to the curvature of the magnet in the horizontal plane. The optimal geometry of the probe was derived using an iterative process combining FEA simulations in Ansys Mechanical with testing of various carbon fiber-reinforced filaments. The printed and assembled probe was successfully used to measure the SDD-059 magnet at Fermilab, yielding promising results.

Jakopin, Alexander [Northern Illinois U.]↗

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗

Validation of Ocean Color Satellite Data Products in Under Sampled Marine Areas

The planktonic marine cyanobacterium, Trichodesmium sp., is broadly distributed throughout the oligotrophic marine tropical and sub-tropical oceans. Trichodesmium, which typically occurs in macroscopic bundles or colonies, is noteworthy for its ability to form large surface aggregations and to fix dinitrogen gas. The latter is important because primary production supported by N2 fixation can result in a net export of carbon from the surface waters to deep ocean and may therefore play a significant role in the global carbon cycle. However, information on the distribution and density of Trichodesmium from shipboard measurements through the oligotrophic oceans is very sparse. Such estimates are required to quantitatively estimate total global rates of N2 fixation. As a result current global rate estimates are highly uncertain. Thus in order to understand the broader biogeochemical importance of Trichodesmium and N2 fixation in the oceans, we need better methods to estimate the global temporal and spatial variability of this organism. One approach that holds great promise is satellite remote sensing. Satellite ocean color sensors are ideal instruments for estimating global phytoplankton biomass, especially that due to episodic blooms, because they provide relatively high frequency synoptic information over large areas. Trichodesmium has a combination of specific ultrastructural and biochemical features that lend themselves to identification of this organism by remote sensing. Specifically, these features are high backscatter due to the presence of gas vesicles, and absorption and fluorescence of phycoerythrin. The resulting optical signature is relatively unique and should be detectable with satellite ocean color sensors such as the Sea-Viewing Wide Field-of-view Sensor (SeaWiFS).

Subramaniam, Ajit↗

Miniaturized Hollow-Waveguide Gas Correlation Radiometer (GCR) for Trace Gas Detection in the Martian Atmosphere

Gas correlation radiometry (GCR) has been shown to be a sensitive and versatile method for detecting trace gases in Earth's atmosphere. Here, we present a miniaturized and simplified version of this instrument capable of mapping multiple trace gases and identifying active regions on the Mars surface. Reduction of the size and mass of the GCR instrument has been achieved by implementing a lightweight, 1 mm inner diameter hollow-core optical fiber (hollow waveguide) for the gas correlation cell. Based on a comparison with an Earth orbiting CO2 gas correlation instrument, replacement of the 10 meter mUltipass cell with hollow waveguide of equivalent pathlength reduces the cell mass from approx 150 kg to approx 0.5 kg, and reduces the volume from 1.9 m x 1.3 m x 0.86 m to a small bundle of fiber coils approximately I meter in diameter by 0.05 m in height (mass and volume reductions of >99%). This modular instrument technique can be expanded to include measurements of additional species of interest including nitrous oxide (N2O), hydrogen sulfide (H2S), methanol (CH3OH), and sulfur dioxide (SO2), as well as carbon dioxide (CO2) for a simultaneous measure of mass balance.

Wilson, Emily L.↗

A stochastic solver based on the residence time algorithm for crystal plasticity models

Abstract The deformation of crystalline materials by dislocation motion takes place in discrete amounts determined by the Burgers vector. Dislocations may move individually or in bundles, potentially giving rise to intermittent slip. This confers plastic deformation with a certain degree of variability that can be interpreted as being caused by stochastic fluctuations in dislocation behavior. However, crystal plasticity (CP) models are almost always formulated in a continuum sense, assuming that fluctuations average out over large material volumes and/or cancel out due to multi-slip contributions. Nevertheless, plastic fluctuations are known to be important in confined volumes at or below the micron scale, at high temperatures, and under low strain rate/stress deformation conditions. Here, we develop a stochastic solver for CP models based on the residence-time algorithm that naturally captures plastic fluctuations by sampling among the set of active slip systems in the crystal. The method solves the evolution equations of explicit CP formulations, which are recast as stochastic ordinary differential equations and integrated discretely in time. The stochastic CP model is numerically stable by design and naturally breaks the symmetry of plastic slip by sampling among the active plastic shear rates with the correct probability. This can lead to phenomena such as intermittent slip or plastic localization without adding external symmetry-breaking operations to the model. The method is applied to body-centered cubic tungsten single crystals under a variety of temperatures, loading orientations, and imposed strain rates.

42 ENGINEERING↗

Computational design of closely related proteins that adopt two well-defined but structurally divergent folds

The plasticity of naturally occurring protein structures, which can change shape considerably in response to changes in environmental conditions, is critical to biological function. While computational methods have been used for de novo design of proteins that fold to a single state with a deep free-energy minimum [P.-S. Huang, S. E. Boyken, D. Baker, Nature 537, 320–327 (2016)], and to reengineer natural proteins to alter their dynamics [J. A. Davey, A. M. Damry, N. K. Goto, R. A. Chica, Nat. Chem. Biol.13, 1280–1285 (2017)] or fold [P. A. Alexander, Y. He, Y. Chen, J. Orban, P. N. Bryan, Proc. Natl. Acad. Sci. U.S.A. 106, 21149–21154 (2009)], the de novo design of closely related sequences which adopt well-defined but structurally divergent structures remains an outstanding challenge. We designed closely related sequences (over 94% identity) that can adopt two very different homotrimeric helical bundle conformations—one short (~66 Å height) and the other long (~100 Å height)—reminiscent of the conformational transition of viral fusion proteins. Crystallographic and NMR spectroscopic characterization shows that both the short- and long-state sequences fold as designed. We sought to design bistable sequences for which both states are accessible, and obtained a single designed protein sequence that populates either the short state or the long state depending on the measurement conditions. The design of sequences which are poised to adopt two very different conformations sets the stage for creating large-scale conformational switches between structurally divergent forms.

59 BASIC BIOLOGICAL SCIENCES↗

Design, Development, and Testing of a Water Vapor Exchanger for Spacecraft Life Support Systems

Thermal and environmental control systems for future exploration spacecraft must meet challenging requirements for efficient operation and conservation of resources. Maximizing the use of regenerative systems and conserving water are critical considerations. This paper describes the design, development, and testing of an innovative water vapor exchanger (WVX) that can minimize the amount of water absorbed in, and vented from, regenerative CO2 removal systems. Key design requirements for the WVX are high air flow capacity (suitable for a crew of six), very high water recovery, and very low pressure losses. We developed fabrication and assembly methods that enable high-efficiency mass transfer in a uniform and stable array of Nafion tubes. We also developed analysis and design methods to compute mass transfer and pressure losses. We built and tested subscale units sized for flow rates of 2 and 5 cu ft/min (3.4–8.5 cu m/hr). Durability testing demonstrated that a stable core geometry was sustained over many humid/dry cycles. Pressure losses were very low (less than 0.5 in. H2O (125 Pa) total) and met requirements at prototypical flow rates. We measured water recovery efficiency across a range of flow rates and humidity levels that simulate the range of possible cabin conditions. We measured water recovery efficiencies in the range of 80 to 90%, with the best efficiency at lower flow rates and higher cabin humidity levels. We compared performance of the WVX with similar units built using an unstructured Nafion tube bundle. The WVX achieves higher water recovery efficiency with nearly an order of magnitude lower pressure drop than unstructured tube bundles. These results show that the WVX provides uniform flow through flow channels for both the humid and dry streams and can meet requirements for service on future exploration spacecraft. The WVX technology will be best suited for long-duration exploration vehicles that require regenerative CO2 removal systems while needing to conserve water.

Izenson, Michael G.↗

Classifying photonic topology using the spectral localizer and numerical K -theory

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

97 MATHEMATICS AND COMPUTING↗

Noninterference Systems Developed for Measuring and Monitoring Rotor Blade Vibrations

In the noninterference measurement of blade vibrations, a laser light beam is transmitted to the rotor blade tips through a single optical fiber, and the reflected light from the blade tips is collected by a receiving fiber-optic bundle and conducted to a photodetector. Transmitting and receiving fibers are integrated in an optical probe that is enclosed in a metal tube which also houses a miniature lens that focuses light on the blade tips. Vibratory blade amplitudes can be deduced from the measurement of the instantaneous time of arrival of the blades and the knowledge of the rotor speed. The in-house noninterference blade-vibration measurement system was developed in response to requirements to monitor blade vibrations in several tests where conventional strain gauges could not be installed or where there was a need to back up strain gauges should critical gauges fail during the test. These types of measurements are also performed in the aircraft engine industry using proprietary in-house technology. Two methods of measurement were developed for vibrations that are synchronous with a rotor shaft. One method requires only one sensor; however, it is necessary to continuously record the data while the rotor is being swept through the resonance. In the other method, typically four sensors are employed and the vibratory amplitude is deduced from the data by performing a least square fit to a harmonic function. This method does not require continuous recording of data through the resonance and, therefore, is better suited for monitoring. The single-probe method was tested in the Carl facility at the Wright- Patterson Air Force Base, and the multiple-probe method was tested in NASA Glenn Research Center's Spin Rig facility, which uses permanent magnets to excite synchronous vibrations. Representative results from this test are illustrated in the bar chart. Nonsynchronous vibrations were measured online during testing of the Quiet High Speed Fan in Glenn s 9- by 15-Foot Low-Speed Wind Tunnel. Three sensors were employed, enabling a reconstruction of the vibratory patterns at the leading and trailing edges at the tip span, as well as a determination of vibratory amplitudes for every blade.

Kurkov, Anatole P.↗

Development of Core Cladding Fabrication Techniques for Phase I Fission Propulsion Systems

Phase I fission propulsion systems focus on safety, timely development, and affordability. Prototype and flight units can be tested at full thrust, using resistance heaters to closely simulate heat from a fission reaction. In Phase I ground testing, one goal is to establish a reliable and affordable manufacturing technique for fabricating a flight-like core. A refractory metal (Mo) has been suggested for the core substrate, primarily due to the existence of a significant database for Mo/LJ02 fuel. The core can be fabricated by bundling Mo tubes with a bonding system that meets preliminary test goals. These criteria include materials compatibility, ability to maintain thermal and structural integrity during 10,000 hours of operation, and fabrication with existing facilities. This paper describes an effort to investigate several fabrication techniques in a cost-effective manner. First, inexpensive materials were tested at low temperatures to determine the relative effectiveness of such techniques as welding, brazing, plating, and vacuum plasma spraying (VPSing). Promising techniques were chosen for further evaluation, including thermal and structural studies, using ceramic tubing at intermediate temperatures. The most desirable technique will be tested on actual Mo tubing at anticipated operating temperatures. This work is being performed by the National Aeronautics & Space Administration (NASA) at George C. Marshall Space Flight Center (MSFC), Los Alamos National Laboratory (LANL), and Advanced Methods & Materials (AMM), Inc.

Salvail, Patrick G.↗

Analysis of Nd3+:glass, solar-pumped, high-powr laser systems

The operating characteristics of Nd(3+):glass lasers energized by a solar concentrator were analyzed for the hosts YAG, silicate glass, and phosphate glass. The modeling is based on the slab zigzag laser geometry and assumes that chemical hardening methods for glass are successful in increasing glass hardness by a factor of 4. On this basis, it was found that a realistic 1-MW solar-pumped laser might be constructed from phosphate glass 4 sq m in area and 2 mm thick. If YAG were the host medium, a 1-MW solar-pumped laser need only be 0.5 sq m in area and 0.5 cm thick, which is already possible. In addition, Nd(3+) doped glass fibers were found to be excellent solar-pumped laser candidates. The small diameter of fibers eliminates thermal stress problems, and if their diameter is kept small (10 microns), they propagate a Gaussian single mode which can be expanded and transmitted long distances in space. Fiber lasers could then be used for communications in space or could be bundled and the individual beams summed or phase-matched for high-power operation.

Zapata, L. E.↗

Modeling the Formation and Eruption of Coronal Structures by Linking Data-Driven Magnetofrictional and MHD Simulations for AR 12673

Context .The data-driven and time-dependent modeling of coronal magnetic fields is crucial for understanding solar eruptions. These efforts are complicated by the challenges of finding a balance between physical realism and computing efficiency. One possible technique is to couple two modeling approaches. Aims . Our aim here is to showcase our progress in using time-dependent magnetofrictional model (TMFM) results as input to dynamical magnetohydrodynamic (MHD) simulations. However, due to the different evolution processes in these two models, using TMFM snapshots in an MHD simulation is nontrivial. We address these issues, both physically and numerically, discuss the incompatibility of the TMFM output to serve as the initial condition in MHD simulations, and show our methods of mitigating this. The evolution of the flux systems and the cause of the eruption are investigated. Methods . TMFM is a prevalent approach that has proven to be a very useful tool in the study of the formation of unstable structures in the solar corona. In particular, it is capable of incorporating observational data as initial and boundary conditions and requires shorter computational time compared to MHD simulations. To leverage the efficiency of data-driven TMFM and also to simulate eruptive events in the MHD framework, one can apply TMFM up to a certain time before the expected eruption(s) and then proceed with the simulation in the full or ideal MHD regime in order to more accurately capture the eruption process. Results . We show the results of a benchmark test case with a linked TMFM and MHD simulation to study the evolution of NOAA active region 12673. A rise of a twisted flux bundle through the MHD simulation domain is observed, but we find that the rate of the rise and the altitude reached depends on the time of the TMFM snapshot that was used to initialize the MHD simulation and the helicity injected into the system. The analysis suggested that torus instability and slip-running reconnection could play an important role in the eruption. Conclusions . The results show that the linkage of TMFM and zero-β MHD models can be successfully used to model the eruptive coronal magnetic fields.

F. Daei↗

Methods and Usability Enhancements in Shift for Non-LWR Applications

Several development and analysis tasks were undertaken in FY21 under the Nuclear Energy Advanced Modeling and Simulation program to enhance modeling of non light–water reactors (LWRs) with Shift. Specifically, these efforts targeted enhancements for tristructural isotropic (TRISO) fuel modeling. A new Shift user interface was developed that allows for much better usability and ease of modeling for non-LWR problems and TRISO fuel. Performance studies were conducted using an HTR-10 fuel pebble model by comparing different geometry packages in Shift, KENO-VI, and Serpent. These studies showed that the new geometry package in Shift performs well compared to Serpent for TRISO fuel modeling with consistent tracking options between both packages. The studies also identified the most critical areas of improvement for more efficiently performing Monte Carlo transport on TRISO fuel models with Shift. Tally calculations in Shift were optimized for non-LWR cross section generation and depletion calculations, and areas for further optimization and accuracy improvements were identified. Finally, initial collaboration efforts were formed between Idaho National Laboratory, Argonne National Laboratory, and the Nuclear Regulatory Commission to use Shift for Comprehensive Reactor Analysis Bundle support.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Biosensors for EVA: Improved Instrumentation for Ground-based Studies

During lunar excursions in the EVA suit, real-time measurement of metabolic rate is required to manage consumables and guide activities to ensure safe return to the base. Metabolic rate, or oxygen consumption (VO2), is normally measured from pulmonary parameters but cannot be determined with standard techniques in the oxygen-rich environment of a spacesuit. Our group has developed novel near infrared spectroscopic (NIRS) methods to calculate muscle oxygen saturation (SmO 2), hematocrit, and pH, and we recently demonstrated that we can use our NIRS sensor to measure VO 2 on the leg during cycling. Our NSBRI project has 4 objectives: (1) increase the accuracy of the metabolic rate calculation through improved prediction of stroke volume; (2) investigate the relative contributions of calf and thigh oxygen consumption to metabolic rate calculation for walking and running; (3) demonstrate that the NIRS-based noninvasive metabolic rate methodology is sensitive enough to detect decrement in VO 2 in a space analog; and (4) improve instrumentation to allow testing within a spacesuit. Over the past year we have made progress on all four objectives, but the most significant progress was made in improving the instrumentation. The NIRS system currently in use at JSC is based on fiber optics technology. Optical fiber bundles are used to deliver light from a light source in the monitor to the patient, and light reflected back from the patient s muscle to the monitor for spectroscopic analysis. The fiber optic cables are large and fragile, and there is no way to get them in and out of the test spacesuit used for ground-based studies. With complimentary funding from the US Army, we undertook a complete redesign of the sensor and control electronics to build a novel system small enough to be used within the spacesuit and portable enough to be used by a combat medic. In the new system the filament lamp used in the fiber optic system was replaced with a novel broadband near infrared LED light source. The compact grating spectrometer was replaced with a chip-scale spectrometer. With this new design, the sensor is 4 in 2 in 0.5 in, weighs 60 g, and no fiber optic cables are needed. The sensor, which contains the light source and the spectrometer, is adhered directly to the skin with medical grade adhesive. The sensor can be powered via the USB port of the laptop computer that controls the sensor operation. Alternatively, for studies in the spacesuit, the sensor can be powered by a small battery pack and operated by an ultra-portable hand-held computer. Both the handheld computer and battery pack will easily fit within the PLSS of the test spacesuit. System automation was significantly improved, to add features suggested by our colleagues in the Cardiovascular Laboratory and the NASA JSC Exercise Physiology and Countermeasures Project. The functionality and portability of this system were demonstrated in our UMass laboratory.

Soller, B.↗

Multiscale thermal hydraulic coupling methods for boiling water reactor simulation

During the last two years, the Virtual Environment for Reactor Applications (VERA) has been extended to simulate boiling water reactors (BWRs). The thermal hydraulic effects present in BWRs are far more complex than those in pressurized water reactors (PWRs). Therefore, the runtime is significantly increase and convergence behavior worse compared to PWRs. Most of the additional runtime is spent during the first few coupled iterations, when the power shape is still rapidly evolving, dramatically affecting the thermal hydraulics (TH). To alleviate the increased computational expense, a multiscale TH coupling approach was developed in VERA in which a highly efficient, simplified TH model is solved for several coupled iterations until the power shape is partially converged. The simplified solution is used to inform the assembly-wise flow distribution in the high-fidelity TH solver CTF, reducing the amount of work required to properly balance the pressure drop in each channel and reducing runtime for the couple calculation. Because the final TH calculations are performed with CTF, there is ultimately no impact on the accuracy of the converged solution. This paper presents the details of this multiscale TH coupling approach, along with results for a 4x4 array of GE-14 fuel bundles and whole-core coupled simulations of the Hatch core. These two cases show that the multiscale approach dramatically reduces the runtime of coupled BWR simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗