Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Rheological diagnostic tools for state transitions

Thermodynamic or kinetic state transitions of six polymeric systems were probed using rheology in isothermal and nonisothermal isochronal experiments. The polymeric systems were representative of three different types of transitions: Gelation, macrophase separation, and order-disorder transitions. A phenomenological dimensionless parameter termed the relative elasticity, R, was found to be particularly sensitive to these transitions and was finitely bound by the viscous limit and the elastic limit with values of zero and one, respectively. In addition, the (time or temperature) derivative of R and tan δ (with δ being the phase angle) proved to be more sensitive to transitions as compared to the respective pure values. When compared to other methods found in literature, these new phenomenological protocols are found to be robust while providing a convenient and precise means of detecting transitions. Hence, they may complement existing methods, whereas they are beneficial to industrial or other systems exhibiting multiple transitions due to their high sensitivity.

42 ENGINEERING↗
DOE OSTI2021

TRUST Nonlinear Dynamics Testbed Assessment

The following assessment evaluates the efficacy of the control script for carrying out a linear signal to generate a linear mechanical response of the system: the TRUST nonlinear dynamics (TRUSTND) testbed. The hardware has three main components: the controller (NI PXIe 8861) embedded in the National Instruments chassis (NI PXIe 1092), the signal amplifier (The Modal Shop Linear Power Amplifier 2050E09) to amplify the output of a custom LabVIEW script, and a shaker (The Modal Shop Electrodynamic Exciter 2075E) where wave spring specimens live and are tested within the attached aluminum testbed (center mass). Different kinds of signals (single-tone, swept frequency and white noise) were sent through this hardware in the E-1 lab space at TA-53. First, the linearity of the setup was tested by specifying a sine wave in the control script and checking the quality of the oscillations in a wave spring (McMaster-Carr 9714K19 [1]). When passing a sine wave through the amplifier was succeeded, the frequency response of the wave spring was tested with white noise for comparison to finite-element model predictions. The testbed outputs a repeatable, linear response to the input excitation when used with the wave spring. These measurements reduce the uncertainty associated with the testbed and test procedure to better characterize uncertainty associated with material behavior (SX358 foam samples) in future tests

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A Lagrangian discontinuous Galerkin hydrodynamic method

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for solving the two-dimensional gas dynamic equations on unstructured hybrid meshes. The physical conservation laws for the momentum and total energy are discretized using a DG method based on linear Taylor expansions. Three different approaches are investigated for calculating the density variation over the element. The first approach evolves a Taylor expansion of the specific volume field. The second approach follows certain finite element methods and uses the strong mass conservation to calculate the density field at a location inside the element or on the element surface. The third approach evolves a Taylor expansion of the density field. The nodal velocity, and the corresponding forces, are explicitly calculated by solving a multidirectional approximate Riemann problem. An effective limiting strategy is presented that ensures monotonicity of the primitive variables. This new Lagrangian DG hydrodynamic method conserves mass, momentum, and total energy. Results from a suite of test problems are presented to demonstrate the robustness and expected second-order accuracy of this new method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Insights on finite size effects in ab initio study of CO adsorption and dissociation on Fe 110 surface

Adsorption and dissociation of hydrocarbons on metallic surfaces represent crucial steps on the path to carburization, eventually leading to dusting corrosion. While adsorption of CO molecules on Fe surface is a barrier-less exothermic process, this is not the case for the dissociation of CO into C and O adatoms and the diffusion of C beneath the surface that are found to be associated with large energy barriers. In practice, these barriers can be affected by numerous factors that combine to favour the CO-Fe reaction such as the abundance of CO and other hydrocarbons as well as the presence of structural defects. From a numerical point of view, studying these factors is challenging and a step-by-step approach is necessary to assess, in particular, the influence of the finite box size on the reaction parameters for adsorption and dissociation of CO on metal surfaces. Here, we use density functional theory (DFT) total energy calculations with the climbing-image nudged elastic band method to estimate the adsorption energies and dissociation barriers for different CO coverages with surface supercells of different sizes. We further compute the effect of periodic boundary condition for DFT calculations and find that the contribution from van der Waals interaction in the computation of adsorption parameters is important as they contribute to correcting the finite-size error in small systems. The dissociation process involves carbon insertion into the Fe surface causing a lattice deformation that requires a larger surface system for unrestricted relaxation. We show that, in the larger surface systems associated with dilute CO-coverages, C-insertion is energetically more favourable, leading to a significant decrease in the dissociation barrier. This observation suggests that a large surface system with dilute coverage is necessary for all similar metal-hydrocarbon reactions in order to study their fundamental electronic mechanisms, as an isolated phenomenon, free from finite-size effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Quantization and anomalous structures in the conductance of Si/SiGe quantum point contacts

Quantum point contacts (QPCs) are fabricated on modulation-doped Si/SiGe heterostructures and ballistic transport is studied at low temperatures. We observe quantized conductance with subband separations up to 4 meV and anomalies in the first conductance plateau at 4e{sup 2}/h. At a temperature of T = 22 mK in the linear transport regime, a weak anomalous kink structure arises close to 0.5(4e{sup 2}/h), which develops into a distinct plateau-like structure as temperature is raised up to T = 4 K. Under magnetic field parallel to the wire up to B = 14 T, the anomaly evolves into the Zeeman spin-split level at 0.5(4e{sup 2}/h), resembling the '0.7 anomaly' in GaAs/AlGaAs QPCs. Additionally, a zero-bias anomaly (ZBA) is observed in nonlinear transport spectroscopy. At T = 22 mK, a parallel magnetic field splits the ZBA peak up into two peaks. At B = 0, elevated temperatures lead to similar splitting, which differs from the behavior of ZBAs in GaAs/AlGaAs QPCs. Under finite dc bias, the differential resistance exhibits additional plateaus approximately at 0.8(4e{sup 2}/h) and 0.2(4e{sup 2}/h) known as '0.85 anomaly' and '0.25 anomaly' in GaAs/AlGaAs QPCs. Unlike the first regular plateau at 4e{sup 2}/h, the 0.2(4e{sup 2}/h) plateau is insensitive to dc bias voltage up to at least V{sub DS} = 80 mV, in-plane magnetic fields up to B = 15 T, and to elevated temperatures up to T = 25 K. We interpret this effect as due to pinching off one of the reservoirs close to the QPC. We do not see any indication of lifting of the valley degeneracy in our samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Equivalent equilibrium ion distributions in ion sources and beams

Ion sources and ion beams often contain ions of more than one charge-to-mass ratio, due to the presence of multiple ionization states and/or multiple atomic and molecular ion species in the source. We show that, under certain conditions, analysis of these systems can be simplified using the notion of “equivalent equilibrium ion distributions” (EEIDs). These distributions have identical ion trajectories, space-charge distributions, and electric-field distributions, but can be composed of ion species with arbitrary charge-to-mass ratios. We use this result to obtain multi-species generalizations for the perveance of an ion diode, which relates the diode voltage to the extracted beam current, and for the perveance of a drifting ion beam, which determines the axial length-scale for beam expansion due to the self-electric field. Applied to a quasi-neutral plasma ion source containing cold ions and hot electrons, we find that the location of the beam-plasma interface (meniscus), where ions are extracted from the plasma, depends on the electron density and ion energy per unit charge at the interface, independent of the ion composition. For certain types of numerical simulations, EEIDs allow one to calculate the beam optics of an intense multi-species ion source by doing a simulation for a single ion species with an artificial charge-to-mass ratio chosen to minimize computation time. We illustrate these results using the 3D particle-in-cell, finite-element code Aleph. The simulations show how ions in a plasma source self-consistently accelerate to the supersonic speed needed for steady, radially expanding flow. Multi-species equilibria are subject to streaming instabilities because, while the different species have the same kinetic energy per unit charge, their velocities scale with the square root of the charge-to-mass ratio. The simulations show that a strong ion streaming instability can occur in a multi-species plasma ion source.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2003

Applications of Space-Filling-Curves to Cartesian Methods for CFD

The proposed paper presents a variety novel uses of Space-Filling-Curves (SFCs) for Cartesian mesh methods in 0. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, most are applicable on general body-fitted meshes -both structured and unstructured. We demonstrate the use of single O(N log N) SFC-based reordering to produce single-pass (O(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations. Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 512 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 10% of ideal even with only around 50,000 cells in each subdomain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with O(max(M,N)) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for finite-difference-based gradient design methods.

Aftosmis, Michael J.↗
NASA NTRS1996

Development of Methodologies for the Estimation of Thermal Properties Associated with Aerospace Vehicles

A thermal stress analysis is an important aspect in the design of aerospace structures and vehicles such as the High Speed Civil Transport (HSCT) at the National Aeronautics and Space Administration Langley Research Center (NASA-LaRC). These structures are complex and are often composed of numerous components fabricated from a variety of different materials. The thermal loads on these structures induce temperature variations within the structure, which in turn result in the development of thermal stresses. Therefore, a thermal stress analysis requires knowledge of the temperature distributions within the structures which consequently necessitates the need for accurate knowledge of the thermal properties, boundary conditions and thermal interface conditions associated with the structural materials. The goal of this proposed multi-year research effort was to develop estimation methodologies for the determination of the thermal properties and interface conditions associated with aerospace vehicles. Specific objectives focused on the development and implementation of optimal experimental design strategies and methodologies for the estimation of thermal properties associated with simple composite and honeycomb structures. The strategy used in this multi-year research effort was to first develop methodologies for relatively simple systems and then systematically modify these methodologies to analyze complex structures. This can be thought of as a building block approach. This strategy was intended to promote maximum usability of the resulting estimation procedure by NASA-LARC researchers through the design of in-house experimentation procedures and through the use of an existing general purpose finite element software.

Scott, Elaine P.↗
DOE OSTI2018

Magnonic Floquet Hofstadter butterfly

We introduce the magnonic Floquet Hofstadter butterfly in the two-dimensional (2D) insulating honeycomb ferromagnet using a combination of spin wave theory and quantum field theory. We show that hopping magnetic dipole moment (i.e. magnon quasiparticles) in the 2D insulating honeycomb ferromagnet accumulates the Aharonov–Casher phase when irradiated by an oscillating space- and time-dependent electric field. We study two different cases. In the first case, we consider the effects of only space-dependent electric field. In this case, we realize the magnonic Hofstadter spectrum with similar fractal structure as graphene subject to a perpendicular magnetic field, but with no spin degeneracy due to broken time-reversal symmetry by the ferromagnetic order. In addition, the magnonic Dirac points and Landau levels occur at finite energy as expected in a bosonic system. Remarkably, this discrepancy does not affect the topological invariant of the system. Consequently, the magnonic Chern number assumes odd values and the magnon Hall conductance gets quantized by odd integers. In the second case, we study the effects of both space- and time-dependent electric field. In this case, the theoretical framework is studied by the Floquet formalism. We show that the magnonic Floquet Hofstadter spectrum emerges entirely from the oscillating space- and time-dependent electric field, which is in stark contrast to electronic Floquet Hofstadter spectrum, where irradiation by circularly polarized light and a perpendicular magnetic field are applied independently. We study the deformation of the fractal structure at different laser frequencies and amplitudes, and analyze the topological phase transitions associated with gap openings in the magnonic Floquet Hofstadter butterfly.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

CODE-TO-CODE BENCHMARK STUDY FOR THERMAL STRESS MODELING AND PRELIMINARY ANALYSIS OF THE HIGH-TEMPERATURE SINGLE HEAT PIPE EXPERIMENT

Microreactors are very small nuclear reactors with typical thermal-energy output of up to 20 MWth. The microreactor concept is gaining more and more attention for safe, robust, and reliable supply of electricity for remote locations and to meet industrial process heat needs. Heat-pipe-cooled microreactor is one of the microreactors being investigated at Idaho National Laboratory. In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios for the heat-pipe-cooled microreactor designs. This paper discusses finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory, using two commercial software packages, Abaqus and STAR-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the locations where the largest temperature gradients appeared between heat pipe and electric heater. The temperature fields have good agreements between Abaqus and STAR-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Assessing the structure of classical molecular optimal control landscapes

Here we consider laser-driven optimal control landscape of a molecule from a classical mechanical perspective. The goal of optimal control in the present work is to steer the molecule from an initial state to a target state, denoted by two distinct points in phase space. Thus, a particular control objective is given as the difference between the final achieved phase space point and the target. The corresponding control landscape is defined as the latter control objective as a functional of the control field. While previous examination of the landscape critical points (i.e., a suboptimal point on the landscape where there is a zero gradient) has shown that the landscape topology is generally trap-free, the structure of the landscape away from these critical points is not well understood. We explore the landscape structure by examining an underlying metric defined as the ratio R of the gradient-based optimization path length of the control field evolution to the Euclidean distance between a given initial control field and the resultant optimal control field, where the latter field corresponds to a point at the top of the landscape. We analyze the path length-to-distance ratio R analytically for a linear forced harmonic oscillator and numerically for a nonlinear forced Morse oscillator. For the linear forced harmonic oscillator, we find that R ≤ $\sqrt{2}$ and reaches its minimum value of 1 (i.e., corresponding to “a straight shot” through control space) in the large target time limit, as well as at special finite target times. The ratio R is similarly small for Morse oscillator simulations when following a steepest-ascent path to the top of the landscape, implying that the landscape is quite smooth and devoid of gnarled features. This conclusion is exemplified for a path discovered with R≃ 1.0 where simply following the initial gradient direction takes the climb very close to the top of the landscape. These findings are consistent with a variety of previous like simulations examining R in quantum control scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Iterative stability analysis of spatial domain decomposition based on block Jacobi algorithm for the diamond-difference scheme

Here, we study convergence of the integral transport matrix method (ITMM) based on a parallel block Jacobi (PBJ) iterative strategy for solving particle transport problems. The ITMM is a spatial domain decomposition method proposed for massively parallel computations. AFourier analysis of the PBJ-based iterations applied to SN diamond-difference equations in 1D slab and 2D Cartesian geometries is performed. It is carried out for infinite-medium problems with homogeneous material properties. To analyze the performance of the ITMM with the PBJ algorithm and evaluate its potential in scalability we consider a limiting case of one spatial cell per subdomain. The analysis shows that in such limit the spectral radius of the iteration method is one without regard to values of the scattering ratio and optical thickness of the spatial cells. This implies lack of convergence in infinite medium. Numerical results of finite-medium problems are presented. They demonstrate effects of finite size of spatial domain on the performance of the iteration algorithm as well as its asymptotic behavior when the extent of the spatial domain increases. Finally, these numerical experiments also show that for finite domains iterative convergence to a finite criterion is achievable in a multiple of the sum of number of cells in each dimension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Efficient nanorod-based amorphous silicon solar cells with advanced light trapping

We present a simple, low-cost, and scalable approach for the fabrication of efficient nanorod-based solar cells. Templates with arrays of self-assembled ZnO nanorods with tunable morphology are synthesized by chemical bath deposition using a low process temperature at 80 °C. The nanorod templates are conformally coated with hydrogenated amorphous silicon light absorber layers of 100 nm and 200 nm thickness. An initial efficiency of up to 9.0% is achieved for the optimized design. External quantum efficiency measurements on the nanorod cells show a substantial photocurrent enhancement both in the red and the blue parts of the solar spectrum. Key insights in the light trapping mechanisms in these arrays are obtained via a combination of three-dimensional finite-difference time-domain simulations, optical absorption, and external quantum efficiency measurements. Front surface patterns enhance the light incoupling in the blue, while rear side patterns lead to enhanced light trapping in the red. The red response in the nanorod cells is limited by absorption in the patterned Ag back contact. With these findings, we develop and experimentally realize a further advanced design with patterned front and back sides while keeping the Ag reflector flat, showing significantly enhanced scattering from the back reflector with reduced parasitic absorption in the Ag and thus higher photocurrent generation. Many of the findings in this work can serve to provide insights for further optimization of nanostructures for thin-film solar cells in a broad range of materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Modeling Nonlinear Energy Dissipation of the Ministack Assembly

An assessment of two methodologies used at Sandia National Laboratories to model mechanical interfaces is performed on the Ministack finite element model. One method uses solid mechanics models to model contacting surfaces with Coulomb frictional contact to capture the physics. The other, termed the structural dynamics reduced order model, models the interface with a simplified whole joint model using four-parameter Iwan elements. The solid mechanics model resolves local kinematics at the interface while the simplified structural dynamics model is significantly faster to simulate. One of the current challenges to using the whole joint model is that it requires calibration to data. A novel approach is developed to calibrate the reduced structural dynamics model using data from the solid mechanics model to match the global dynamics of the system. This is achieved by calibrating to amplitude dependent frequency and damping of the system modes, which are estimated using three different approaches.

42 ENGINEERING↗
NASA NTRS2014

A Technique for Mapping Characteristic Lengths to Preserve Energy Dissipated via Strain Softening in a Multiscale Analysis

It is often advantageous to account for the microstructure of the material directly using multiscale modeling. For computational tractability, an idealized repeating unit cell (RUC) is used to capture all of the pertinent features of the microstructure. Typically, the RUC is dimensionless and depends only on the relative volume fractions of the different phases in the material. This works well for non-linear and inelastic behavior exhibiting a positive-definite constitutive response. Although, once the material exhibits strain softening, or localization, a mesh objective failure theories, such as smeared fracture theories, nodal and element enrichment theories (XFEM), cohesive elements or virtual crack closure technique (VCCT), can be utilized at the microscale, but the dimensions of the RUC must then be defined. One major challenge in multiscale progressive damage modeling is relating the characteristic lengths across the scales in order to preserve the energy that is dissipated via localization at the microscale. If there is no effort to relate the size of the macroscale element to the microscale RUC, then the energy that is dissipated will remain mesh dependent at the macroscale, even if it is regularized at the microscale. Here, a technique for mapping characteristic lengths across the scales is proposed. The RUC will be modeled using the generalized method of cells (GMC) micromechanics theory, and local failure in the matrix constituent subcells will be modeled using the crack band theory. The subcell characteristic lengths used in the crack band calculations will be mapped to the macroscale finite element in order to regularize the local energy in a manner consistent with the global length scale. Examples will be provided with and without the regularization, and they will be compared to a baseline case where the size and shape of the element and RUC are coincident (ensuring energy is preserved across the scales).

softening↗
NASA NTRS2003

Applications of Space-Filling-Curves to Cartesian Methods for CFD

This paper presents a variety of novel uses of space-filling-curves (SFCs) for Cartesian mesh methods in CFD. While these techniques will be demonstrated using non-body-fitted Cartesian meshes, many are applicable on general body-fitted meshes-both structured and unstructured. We demonstrate the use of single theta(N log N) SFC-based reordering to produce single-pass (theta(N)) algorithms for mesh partitioning, multigrid coarsening, and inter-mesh interpolation. The intermesh interpolation operator has many practical applications including warm starts on modified geometry, or as an inter-grid transfer operator on remeshed regions in moving-body simulations Exploiting the compact construction of these operators, we further show that these algorithms are highly amenable to parallelization. Examples using the SFC-based mesh partitioner show nearly linear speedup to 640 CPUs even when using multigrid as a smoother. Partition statistics are presented showing that the SFC partitions are, on-average, within 15% of ideal even with only around 50,000 cells in each sub-domain. The inter-mesh interpolation operator also has linear asymptotic complexity and can be used to map a solution with N unknowns to another mesh with M unknowns with theta(M + N) operations. This capability is demonstrated both on moving-body simulations and in mapping solutions to perturbed meshes for control surface deflection or finite-difference-based gradient design methods.

Aftosmis, M. J.↗
DOE OSTI2024

Environmentally Assisted Fatigue in Light Water Reactor Environment

This report summarizes the Environmentally Assisted Fatigue (EAF) research conducted at ANL under the US DOE Light Water Reactor Sustainability (LWRS) program. Starting from a rich background in theoretical and experimental EAF, ANL previously developed an approach to evaluate fatigue performance of reactor materials in light water reactor environments with the correction factor F en . The approach was based on a large body of experimental work performed at ANL and elsewhere, and was consistent with American Society of Mechanical Engineers (ASME)’s methodology governing the design and construction of reactor components. In recent years, the program was focused on component fatigue prediction and made several major and fundamental contributions in this area. These accomplishments help meet the needs identified by the industry concerning component level fatigue predictions in complex, transient conditions. The main contribution of the ANL program involved the development of a system-level model for estimating residual strain and life of nuclear reactor coolant system components under connected-system-thermal-mechanical boundary conditions. The goal was to predict the stress hotspots, strain residuals, strain amplitudes and the resulting fatigue lives. Thermal-mechanical stress analysis was performed considering thermal stratification and a design-basis reactor loading cycle. Based on the finite element (FE) model results, the strain residuals, strain amplitudes and resulting fatigue lives of reactor coolant system (RCS) components were predicted. The results show that some of the RCS components can have significantly different strain amplitudes, residual strain, and fatigue lives, despite having similar geometry and material. In addition, the simulated component-level strain profile can guide the selection of appropriate test inputs for conducting laboratory-scale EAF tests. Building upon the system-level model, ANL developed a digital twin (DT) framework to predict the structural states and associated fatigue life of components in real-time. This framework is a comprehensive system designed to predict the structural states and fatigue lives of reactor components. It includes multiple models and integrates artificial intelligence (AI), machine learning (ML), and FE based modeling tools to evaluate the structural states and fatigue lives.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗