Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element problems”

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.

1,074 records · Page 60

Development of a Thin Film Primary Surface Heat Exchanger for Advanced Power Cycles

This project objective is to develop a high-temperature design upgrade for an existing primary surface heat exchanger so that the redesigned hardware is capable of operation in CO 2 at temperatures up to 1,510°F (821°C) and pressure differentials up to 130 psi (9 bar). The heat exchanger is proposed for use as a recuperator in an advanced low-pressure oxy-fuel Brayton cycle that is predicted to achieve over 50% thermodynamic efficiency, although the heat exchanger could also be used in other high-temperature, low-differential pressure cycles. This report describes the progress to date, which includes continuing work performed to select and test new candidate materials for the recuperator redesign, final mechanical and thermal performance analysis results of various redesign concepts, and the preliminary design of a test loop for the redesigned recuperator including a budgetary estimate for detailed test loop design, procurement, and test operation. A materials search was performed in order to investigate high-temperature properties of many candidate materials, including high-temperature strength and nickel content. These properties were used to rank the candidate materials, resulting in a reduced list of nine materials for corrosion testing. Multiple test rigs were considered and analyzed for short-term corrosion testing and Thermal Gravimetric Analysis (TGA) was selected as the most cost-effective option for evaluating corrosion resistance of the candidate materials. In addition, tantalum, niobium, and chromium coatings were identified as potential options for increased corrosion resistance. The test results show that many materials exhibit relatively low weight gain rates, and that niobium and tantalum coatings may improve corrosion resistance for many materials, while chromium coatings appear to oxidize and debond quickly. Metallurgical analysis of alloys was also performed, showing evidence of intergranular attack in 282 that may cause long-term reliability problems in CO 2 service at these temperatures. However, long-term testing in a flowing environment is recommended in order to understand accurately the severity of the attack. Detailed economic modeling of the existing air cycle recuperator and CO 2 cycle recuperator options was also completed, including costs for material, fabrication, fuel, maintenance, and operation. The analysis results show that the increased capital cost for high-temperature materials may be offset by higher cycle efficiencies, decreasing the overall lifetime cost of the system. The economic analysis also examines costs associated with increased pressure drop and material changes for two redesign options. These results show that, even with slightly reduced performance and/or higher material costs, the lifetime cost per energy production may still be reduced by over 12%. The existing recuperator design information was provided by Solar Turbines, Inc. via several models, drawings, and design handoff meetings. Multiple fluid/thermal and structural models were created in order to analyze critical recuperator performance and mechanical strength in critical areas throughout the redesign process. These models were analyzed for a baseline condition (consistent with current Mercury 50 operation) for validation purposes. Results are presented for heat transfer coefficients and pressure drops, matching well with the existing operational data. Simulation of higher-temperature CO 2 conditions was also performed, showing a slight expected increase in both heat transfer and pressure drop. Mechanical analysis results for critical areas on the cross-flow and counter-flow sheets have also been obtained for air and CO 2 cases. These results show similar stresses in both cases but significantly reduced safety factors for the CO 2 case due to reduced yield and creep rupture strengths of alloy 625 at the higher temperatures. A concept brainstorm session and initial down-selection were completed in order to identify promising redesign options for further analysis. Detailed analysis of all promising redesign options was performed via finite element and computational fluid dynamic simulations in order to characterize mechanical and thermal-fluid performance of each option. These options included material change, various sheet thickness configurations, pitch and phasing of cross-flow and counter-flow sheets, and separator sheets. The analysis results have identified two viable redesign options that maintain existing safety margins optimally through a material change to Haynes 282 and (A) sheet thickness increases of 40% on the counter-flow sheet and 75% on the hot side cross-flow corrugation sheet or (B) addition of a separator sheet in the counter-flow section while maintaining the original counter-flow sheet thickness and increasing the cross-flow corrugation sheet thickness by 90% to account for the increase in cell height. While both options satisfy mechanical stress constraints, the separator sheet design has a higher part count, slightly reduced heat transfer, and slightly higher pressure drop than the first option and is not preferred. Finally, several test loop concepts have been developed for different full-scale and reduced-scale recuperator testing options. For each option, various loop components, such as heat exchangers, valves, heaters, and compressors, were evaluated in an effort to maximize utilization of existing resources. All concepts utilize an existing 3-MW CO 2 compressor, heater, and loop coolers, but the concepts vary by incorporating different amounts of new equipment for achieving various flow rates (all concepts operate at design pressure and temperature). The third concept achieves a 1 kg/s test without purchasing any costly equipment (coolers, heaters, blowers, etc.). Since the stacked cell design of the recuperator results in the same flow conditions at each core cell (even for a reduced-scale test). Thus, test loop Concept #3 was selected for the preliminary design. This loop design is detailed within the report, culminating in a budgetary estimate of $1,013,000.00 for the detailed design, construction, commissioning, and operation of a high-temperature recuperator test loop.

03 NATURAL GAS↗

A non-uniform corrosion model and meso-scale fracture modelling of concrete

Corrosion-induced concrete cracking is a significant durability problem for reinforced concrete structures. Considerable research has been carried out in the last few decades to understand and model the expansion mechanism of the corrosion products around the reinforcing bar and simulate the cracking behaviour of the concrete cover. In this paper, a new corrosion model based on non-uniform corrosion expansion is formulated and validated against experimental data. A meso-scale fracture model, consisting of aggregates, cement paste/mortar and ITZ, is established for the cases of both middle and side reinforcing bars. Under the developed corrosion and concrete fracture model, the cracking phenomena of the concrete cover are accurately simulated. It has been found that the non-uniform corrosion model can be used to express the realistic corrosion rust progression around the reinforcing bar, with the best accuracy. It has also been found that some microcracks occur before they are connected to form the dominating discrete crack which usually appears on the concrete surface. Moreover, the effects of the corrosion variables, as well as other key material and geometric parameters, on surface cracking of concrete are investigated.

36 MATERIALS SCIENCE↗

MOOSE Framework Meshing Enhancements to Support Reactor Analysis

MOOSE-based physics codes require an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The use of difficult-touse, external licensed software is often required to generate high quality meshes for reactor geometries. High-fidelity geometry modeling also requires elaborate tracking of groups of elements for material property assignment and output reporting which can be considerably complex for the user to identify and maintain. Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, several meshingrelated enhancements have been developed for the MOOSE framework to address user challenges in creating finite element meshes for advanced reactor geometries. MOOSE mesh generators have been developed to mesh hexagonal geometries (pins, ducted assemblies, and cores) commonly found in liquid-metal cooled fast reactor concepts. The mesh generator used for hexagonal pin cells is generic for regular polygons and therefore may also be used for Cartesian pin cells. Hexagonal pin cells can be stitched into ducted assemblies, and assemblies can be stitched together into a core. The user may specify region ids, region names, and other preferences on the mesh. This control is useful for later material mapping in the MOOSE-based physics codes input. A capability was also developed for meshing rotating control drums including determination of material volume fractions in each mesh element as a function of time. Control drum meshes may be stitched to other hexagonal assemblies to create a core configuration. Additional mesh generators were developed that wrap around the hexagonal meshing capabilities and utilize “extra element integer” ID values on each element. In regular Cartesian or hexagonal assemblies or cores, the bookkeeping of element groups for both material assignment and output reporting can now be automated through assignment of pin, assembly, core, axial and depletion id values stored as extra element integers. The extra element tags on the mesh greatly speed the reactor analyst’s efforts to map materials to meshes, track depletion zones, and parse output such as axial pin power distributions. At the highest level, pin, assembly, and core mesh generators (with this reactor terminology) have also been developed to easily generate regular Cartesian and hexagonal cores, including axial extrusion. These reactor geometry builders call upon the previously mentioned capabilities to produce analysis-ready 3D meshes including material assignments. Open source mesh triangulation capabilities were also investigated for integration into the MOOSE framework to address the need for meshing the core periphery region which extends from the irregular outer assembly border to a cylindrical boundary. Options are limited due to licensing constraints, and the recommendation is pursue building a native MOOSE Delaunay triangulator routine with full functionality. Finally, a series of verification problems were performed with NEAMS physics tools. All developed capabilities will be available in the new open-source “Reactor” module of the MOOSE framework, which is accessible to any MOOSE-based NEAMS physics tool.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗

Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation

This Final Scientific and Technical Report summarizes work performed under the Phase IIA SBIR project “Enabling the Broader Use of MOOSE for Nuclear Energy and Other Simulation” (DE-SC0020906) from August 2023 through August 2025. The objective of the Phase IIA effort was to mature and harden capabilities developed during Phase II, with the goal of enabling practical interoperability between Coreform’s isogeometric analysis (IGA) technologies and the Multiphysics Object-Oriented Simulation Environment (MOOSE), while improving robustness, performance, and scalability for complex, nuclear-relevant geometries. Over the course of Phase IIA, the project established and validated an extraction-based interoperability pathway between Coreform tools and MOOSE. A combined mesh and matrix format was defined collaboratively with MOOSE developers and integrated into the solver, enabling standard MOOSE workflows to operate on data exported from Coreform’s IGA and Flex Representation Method (FRM) pipelines. Early demonstrations validated architectural compatibility using linear solid mechanics problems, while later efforts focused on benchmark testing and external use. By the end of the project period, engineers at BWXT were able to independently set up and execute a simulation using the Coreform–MOOSE workflow and provide direct feedback that informed further refinement. In parallel, substantial effort was devoted to improving the robustness of trimmed U-spline construction for complex CAD geometries. A growing test suite of nuclear-relevant models was compiled through collaboration with multiple stakeholders and used to drive extensive bug fixing and reliability improvements. These efforts resulted in improved robustness and performance, including the addition of fallback capabilities that enhance reliability when the underlying commercial CAD kernel fails. Performance-oriented work progressed later in the project, with the development and demonstration of methods to decompose complex geometries into structured subregions and updated data representations to support more efficient solver processing. Additionally, extensive enhancements to threadsafe parallel data structures and trimming operations established a foundation for scalable processing of large assemblies. Collaboration with Sandia National Laboratories on the SGM geometric modeling kernel advanced to a functioning interface test case, positioning the workflow for future kernel integration. Overall, the Phase IIA effort successfully transitioned the project from architectural proof-of-concept to externally exercised, solver-integrated capability, while clarifying remaining technical challenges related to standardization, performance optimization, and kernel integration.

42 ENGINEERING↗

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

A STUDY OF CHANGES IN DEFORMATION AND METABOLISM IN LEFT VENTRICLE AS A FUNCTION OF HYPERTROPHY IN SPONTANEOUS HYPERTENSIVE RATS USING MICROPET TECHNOLOGY

Problem: In the case of hypertrophy caused by pressure overload (hypertension) there is an increase in cardiac mass and modification cardiac metabolism. Aim: This study was designed to study the changes in glucose metabolism, ejection fraction, and deformation in the left ventricle with the progression of hypertrophy in spontaneous hypertensive rats (SHR). Methods: Dynamic PET data were acquired using the microPET II at UC Davis. Two rats were imaged at 10-week intervals for 18 months. Each time a dose of approximately 1- 1.5 mCi of F-18-FDG was injected into a normotensive Wistar Kyoto (WKY) rat and the same dose was injected into a SHR rat. Each rat was imaged using a gated dynamic acquisition for 80 minutes acquiring list mode data with cardiac gating of approximately 600-900 million total counts. For the analysis of glucose of metabolism, the list mode data were histogrammed into a dynamic sequence (42 frames over 80 mins). For each time frame, projection data of 1203 140x210 sinograms of 0.582 mm bins were formed by summing the last three gates before and one after the R-wave trigger to correspond to the diastolic phase of the cardiac cycle. Dynamic sequences of 128x128x83 matrices of 0.4x0.4x0.582 mm3 voxels in x, y, and z were reconstructed using an iterative MAP reconstruction which used a prior that penalized the high frequency components of the reconstruction using appropriate weighting between 26 nearest neighboring voxels. Time activity curves were generated from the dynamic reconstructed sequence for the blood and left ventricular tissue regions of interest which were fit to a 2-compartment model to obtain a least squares fit for the kinetic parameters. For the analysis of deformation, the list mode data were histogrammed into 8 gates of the cardiac cycle, each gate was the total sum of the later 60 mins of the 80 min acquisition. Images of 128x128x83 matrices for each gate were reconstructed using the same iterative MAP reconstruction used to reconstruct the dynamic sequence. The in-plane image dimensions were doubled to 256x256x83 in order to increase the resolution for the Warping analyses. These image data sets were then cropped to 128x128x83. The end-systolic image data sets were designated as the template images and the end-diastole image data sets were designated as the target images, thus providing an analysis of the diastolic relaxation and filling phases of the cardiac cycle. The template images were manually segmented to create surface definitions representing the epi- and endocardial surfaces. Finite element models of the left ventricles were created using the segmented surfaces and defining a transversely isotropic material with fiber angles varying from the epicardial surface to the endocardial surface. A Warping analyses was performed to obtain the LV strain tensor and fiber stretch distributions. Results: In one study, the average first principal Green-Lagrange strain, fiber stretch, ejection fraction, and metabolic rate of F-18-FDG was 0.22, 1.08, 80%, 0.1 for the WKY rat and 0.16, 1.06, 50%, 0.25 for the SHR rat, respectively. These same rats studied a year later presented with a metabolic rate of F-18-FDG of 0.11 and 0.25 for the WKY and SHR, respectively. A follow-up study the average strain (n=10) and ejection fraction (n=18) was 0.21, 72.7% for WKY and 0.17, 69.8% for the SHR, respectively. Conclusion: In the case of pressure overload there is an increased reliance on carbohydrate oxidation in an attempt to maintain contractile function.

Sitek, Arkadiusz↗

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiphysics modeling of two-phase film boiling within porous corrosion deposits

Porous corrosion deposits on nuclear fuel cladding, known as CRUD, can cause multiple operational problems in light water reactors (LWRs). CRUD can cause accelerated corrosion of the fuel cladding, increase radiation fields and hence greater exposure risk to plant workers once activated, and induce a downward axial power shift causing an imbalance in core power distribution. In order to facilitate a better understanding of CRUD's effects, such as localized high cladding surface temperatures related to accelerated corrosion rates, we describe an improved, fully-coupled, multiphysics model to simulate heat transfer, chemical reactions and transport, and two-phase fluid flow within these deposits. Our new model features a reformed assumption of 2D, two-phase film boiling within the CRUD, correcting earlier models' assumptions of single-phase coolant flow with wick boiling under high heat fluxes. This model helps to better explain observed experimental values of the effective CRUD thermal conductivity. Finally, we propose a more complete set of boiling regimes, or a more detailed mechanism, to explain recent CRUD deposition experiments by suggesting the new concept of double dryout specifically in thick porous media with boiling chimneys. - Highlights: • A two-phase model of CRUD's effects on fuel cladding is developed and improved. • This model eliminates the formerly erroneous assumption of wick boiling. • Higher fuel cladding temperatures are predicted when accounting for two-phase flow. • Double-peaks in thermal conductivity vs. heat flux in experiments are explained. • A “double dryout” mechanism in CRUD is proposed based on the model and experiments.

42 ENGINEERING↗

A combined experimental and numerical approach that eliminates the non-uniqueness associated with the Johnson-Cook parameters obtained using inverse methods

Abstract Johnson-Cook constitutive model is a commonly used material model for machining simulations. The model includes five parameters that capture the initial yield stress, strain-hardening, strain-rate hardening, and thermal softening behavior of the material. These parameters are difficult to determine using experiments since the conditions observed during machining (such as high strain-rates of the order of $$10^5$$ 10 5 /sec - $$10^6$$ 10 6 /sec) are challenging to recreate in the laboratory. To address this problem, several researchers have recently proposed inverse approaches where a combination of experiments and analytical models are used to predict the Johnson-Cook parameters. The errors between the measured cutting forces, chip thicknesses and temperatures and those predicted by analytical models are minimized and the parameters are determined. In this work, it is shown that only two of the five Johnson-Cook parameters can be determined uniquely using inverse approaches. Two different algorithms, namely, Adaptive Memory Programming for Global Optimization (AMPGO) and Particle Swarm Optimization (PSO), are used for this purpose. The extended Oxley’s model is used as the analytical tool for optimization. For determining a parameter’s value, a large range for each parameter is provided as an input to the algorithms. The algorithms converge to several different sets of values for the five Johnson-Cook parameters when all the five parameters are considered as unknown in the optimization algorithm. All of these sets, however, yield the same chip shape and cutting forces in FEM simulations. Further analyses show that only the strain-rate and thermal softening parameters can be determined uniquely and the three parameters present in the strain-hardening term of the Johnson-Cook model cannot be determined uniquely using the inverse method. A combined experimental and numerical approach is proposed to eliminate this determine all parameters uniquely.

42 ENGINEERING↗

Development of Thermoacoustic Sensors for Sodium-cooled Fast Reactors

This Final Report refers to the project “Development of Thermoacoustic Sensors for Sodium-cooled Fast Reactor Systems”, which was led by Westinghouse Electric Company (Westinghouse) and carried out in collaboration with Argonne National Laboratory (ANL) and University of Pittsburgh. Thermo-acoustic Power Sensors (TAPS) are self-powered, wireless sensors envisioned for measuring key parameters, such as local temperature and neutron flux, in a nuclear reactor core. This project was intended to specifically investigate their applicability to Sodium-cooled Fast Reactors (SFR). TAPS are non-invasive (wireless) and passive (self-powered) devices. The passivity derives from their ability to use conditions that “naturally” exist in a nuclear reactor, such as gamma and neutron flux, as power sources. They generate oscillating pressure waves (i.e., sound waves) which, with a frequency and amplitude dependent upon these conditions, can travel through the core and associated structures, and reach the outside of the reactor vessel where a properly designed network of receivers can detect and interpret them. These receivers require a very small amount of power which, during loss of power events, can be provided for example by harvesting gamma radiation energy, thus resulting in a monitoring system that can function both during normal operation and during loss of power events. The project aimed at TAPS development through a series of tasks which are listed and briefly discussed as follows. TASK 1 – Sensor hardware design Subtask 1a: Assessment of sensor applications to SFRs Subtask 1b: Development of sensor functional requirements Subtask 1c: Definition of sensor hardware design specifications Task description: TAPS design was informed by considerations on their application (Subtask 1a), both the ultimate one in an SFR and the actual one in the ANL testing facilities that was intended to be used in support of the project. Considerations were made to identify optimum sensor design features that optimize the sensor size, materials, and output signal, for installation inside an SFR core. These considerations led to the development of Functional Requirements (Subtask 1b) and Design Requirements (Subtask 1c). TASK 2 – Sensor Hardware Manufacture Subtask 2a: Sensor hardware construction drawing development Subtask 2b: Sensor manufacture and assembly Task description: TAPS technical drawings were developed (Subtask 2a) using the Design Requirements established under Task 1. Subsequently, in spite of some problems which ultimately caused the program to be delayed, TAPS manufacturing was completed based on drawings (Subtask 2b). TASK 3 – Development of TAPS Signal Measurement System and TAPS Testing in Water Subtask 3a: Design, assembly and testing of signal measurement system, and TAPS testing in water Subtask 3b: Signal prediction-correction methodology development Task description: An assessment was performed on the techniques that can potentially be used to detect the signals emitted by the TAPS, e.g. a fiber-optic based acoustic signal measurement system, a laser vibrometer system, or an accelerometer-based system. The most suited technology, i.e. the accelerometer-based system, was developed further, and tested in water (Subtask 3a). Moreover, efforts were made to develop the methodology required to determine the actual system temperature and neutron flux distribution using differences between the measured and predicted TAPS responses (Subtask 3b). TASK 4 – Sensor System Testing in Sodium Subtask 4a: Test plan development Subtask 4b: Design, assembly and testing in small-scale sodium facility Subtask 4c: Design, assembly and testing in large-scale sodium and structures facility Task description: Upon proper test plan development (Subtask 4a), the fabricated TAPS was planned to be tested in sodium, by using two sodium facilities at ANL having different size and different purpose. The Under Sodium Viewing (USV) small-scale facility was intended to be used to investigate the effect of sodium on the sensor and its performance (Subtask 4b). The Mechanism Engineering Test Loop (METL) large-scale facility was instead intended to be used to assess the additional effect of prototypical SFR structures, such as fuel assembly mockup or parts of the core restrain structure, on sensor performance (Subtask 4c). As discussed in Section 3.2.2.7, unexpected issues during the TAPS manufacturing process resulted in some activities being delayed, with the TAPS and USV facility developed to the point to be ready for testing in sodium, however without the possibility to actually perform such testing (including the testing in METL) due to the end of the program’s performance period. Overall, through the development and testing (in water only) of two TAPS devices (a First-Generation TAPS followed by an optimized Second-Generation TAPS), the project confirmed the capability of this technology to generate acoustic signals proportional to temperature, which can be detected through a network of accelerometers identified as the best-suited type of receivers for acoustic signal detection. Moreover, the project also developed a computational model to predict the characteristics of the acoustic signals being generated, which combines thermal analysis of the TAPS with Finite Element Modeling (FEM)-aided acoustic characterization of the system. This model was benchmarked against experimental data collected during the project and, although general agreement was obtained, some limitations of the modeling methods were identified, which require additional development. Additional testing is needed in order to assess the effect, on TAPS operation and performance, of environmental changes resulting from the transition from water to liquid sodium. Such testing, which is suggested to be performed in the future, should look specifically at 1) both the effect resulting from the different thermoacoustic behavior of sodium (relative to water) and the effects of higher temperature on TAPS performance, and 2) the performance of the sensor-receiver system when multiple TAPS are used simultaneously and prototypical reactor structures are positioned in the testing environment. The latter testing is needed to assess the effects that potential signal attenuation/ distortion phenomena, as well as potential interference between signals emitted simultaneously, have on the performance of the technology for ultimate application in a nuclear reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Measurement of rabbit eardrum vibration through stroboscopic digital holography

In this work, we present a setup for high-power single shot stroboscopic digital holography and demonstrate it in an application on rabbit eardrum vibration measurement. The setup is able to make full-field time-resolved measurements of vibrating surfaces with a precision in the nanometer range in a broad frequency range. The height displacement of the measured object is visualized over the entire surface as a function of time. Vibration magnitude and phase maps can be extracted from these data, the latter proving to be very useful to reveal phase delays across the surface. Such deviations from modal motion indicate energy losses due to internal damping, in contrast to purely elastic mechanics. This is of great interest in middle ear mechanics and finite element modelling. In our setup, short laser pulses are fired at selected instants within the surface vibration period and are recorded by a CCD camera. The timing of the pulses and the exposure of the camera are synchronized to the vibration phase by a microprocessor. The high-power frequency-doubled Nd:YAG laser produces pulses containing up to 5 mJ of energy, which is amply sufficient to record single-shot holograms. As the laser pulse length is 8 ns and the smallest time step of the trigger electronics is 1 μs, vibration measurements of frequencies up to 250 kHz are achievable through this method, provided that the maximum vibration amplitude exceeds a few nanometers. In our application, middle ear mechanics, measuring frequencies extend from 5 Hz to 20 kHz. The experimental setup will be presented, as well as results of measurements on a stretched circular rubber membrane and a rabbit's eardrum. Two of the challenges when measuring biological tissues, such as the eardrum, are low reflectivity and fast dehydration. To increase reflectivity, a coating is applied and to counteract the undesirable effects of tissue dehydration, the measurement setup and software have been optimized for speed without compromising on the quality. Results of a repeatability test will be presented as well. Since the method measures the membrane motion as a function of time in small time steps, we do not only measure vibration amplitude like in time-average holography, but we can also measure non-linear elastic and transient behaviour. In conclusion, the combination of good spatial, depth and time resolution with the fast data acquisition and very wide frequency range make our technique applicable in a number of fields, including biological tissue vibrations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗