Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flux shape deformation”

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.

8 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Inferring small-sample-reactivity worth from oscillation experiment with auto-corrected local flux

To estimate reactivity changes as a function of time from the signal of detectors count rates, one usually uses point kinetic equations. However, when detectors are located inside the core, it is often necessary to compute a correction factor to account for the flux-shape deformation (i.e. the change of detector efficiency that may have been induced by neutron absorbers movements). Various simulation methods have been developed depending on transients observed. They range from the resolution of static equations up to full-time-dependent neutron transport in 3D. However, those methods require high computational-cost calculations of local reaction rates that rely on nuclear data evaluations. Therefore, a 'measured' reactivity is never determined through pure measurements. Taking advantage of the known oscillation movement of a small sample at the center of the MINERVE reactor, the innovative algorithm presented in this paper derives the correction factor needed to get no time drift on reactivity steps. Derived correction factors are in good agreement with those obtained independently with TRIPOLI-4 through sensitivity analysis. Low sensitivity to the sample composition and moderate sensitivity to the detectors position reinforce the reliability of those results. This application paves the way to almost pure nuclear data reactivity measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

(Doublon) Benchmarking of Different Inverse Point Kinetics Implementations for an Autocorrected Reactimeter Algorithm

In November 2017, the Transient Reactor Test Facility returned to operation. Since that time, many transient test series have been completed, such as the Transient Heatsink Overpower Response capsule (THOR), the Transient Water Irradiation System for TREAT (TWIST), and Sirius. Each has provided valuable data for materials performance and reactor safety that can be applied in future designs. During each experimental series, detector count rates provided important information on the core behavior during transients. However, a limitation of these data is that variations in the neutron distribution during experiments can cause errors when attempting to infer reactivity evolution from detector signals. Neutron physics codes can be used to compute the flux shape variations. However, this is a poor solution when the experimental data is used for code verification, validation and uncertainty quantification. Indeed, if the output of the code is used both as a reference and to correct what the reference is compared to, the circular dependency limits the quality of the verification, validation and uncertainty quantification approach. To overcome this problem, the autocorrected reactimeter algorithm (ACRA) has been developed. This approach infers a time-dependent reactivity evolution by testing different spatial corrections and selecting the one that minimizes reactivity variations when the core is in a frozen configuration (i.e., when there is no variation in parameters affecting reactivity). However, the scope of this method was limited to transients where there were negligible thermal feedback. Indeed, the core is never in a frozen configuration when the fuel temperature varies during the whole transient. This is our motivation for developing an improved version of the ACRA that does not require frozen configurations. To develop this new algorithm, we need a precise and unbiased implementation of the inverse point kinetic equations (IPKEs) as any error in the reactivity evaluation will be propagated into the choice of the optimal spatial correction. Indeed, the previous reactimeter algorithm would use approximations, such as a negligible flux amplitude derivative, to focus on rapidity. For the numerical validation of ACRA, we aim at absolute error under for reactivity derived from signals similar to the one of this study. In this summary, we test eight different IPKE implementations. Each will process a mockup signal built for this study, similar to those that the future ACRA will process. Each reactivity output will be compared to the reference reactivity that has been used to generate the mockup signal. The implementation minimizing the difference with the reference reactivity will be used in the development of a new ACRA formulation.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Combined Dimensional and Topology Optimization of Interior Permanent Magnet Synchronous Machine Rotors Using a Permanent Magnet Function Interpolation Method

This paper presents a magneto-structural combined dimensional and topology optimization technique for interior permanent magnet synchronous machine (IPMSM) rotors. Dimensional changes to the permanent magnet (PM) location or size are accomplished by interpolating or projecting a smoothed Heaviside rectangular function representing the presence of PM material onto the IPMSM rotor design domain mesh. A density based Solid Isotropic with Material Penalization (SIMP) topology optimization approach is then used to vary the presence of electrical steel in mesh elements to form flux barriers around the PM. The proposed method enforces a defined shape for the PM without requiring the mesh in the design domain to be deformed. Three examples are presented to demonstrate the technique: two flat bar IPMSM and one V-shaped IPMSM.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Combined Dimensional and Topology Optimization of Synchronous Machine Rotors Using a Material Density Interpolation Method

This article presents a magneto-structural combined dimensional and topology optimization technique for interior permanent magnet synchronous machine (IPMSM) and wound field synchronous machine (WFSM) rotors. Dimensional changes to the permanent magnet (PM) and rotor winding location or size are accomplished by interpolating or projecting a smoothed Heaviside rectangular function representing the presence of PM material onto the IPMSM rotor design domain mesh, whereas the presence of copper material for WFSM rotor is interpolated in the same manner. A density based Solid Isotropic Material with Penalization (SIMP) topology optimization approach is then used to vary the presence of electrical steel in mesh elements to form flux barriers around the PM. The proposed method enforces a defined shape for the PM without requiring the mesh in the design domain to be deformed. Four examples are presented to demonstrate the technique: Two flat bar IPMSM, one V-shaped IPMSM, and a WFSM. Here, a comparative study is performed on one of the flat bar IPMSMs using Metamodel of Optimal Prognosis (MOP) based method.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Examining astrophysical gas cloud collapse using an optical depth-scaled, x-ray-irradiated, carbon-foam sphere

When stellar radiation interacts with a molecular cloud, the cloud's fate depends on the strength of the incident radiation and the radiation's mean-free-path within the cloud [F. Bertoldi, Astrophys. J. 346, 735–755 (1989)]. Under the right conditions, the radiation compresses the cloud and a star formation may occur. Where and when the stellar formation occurs in the cloud's collapse are open questions. Direct observation of the complete star–cloud lifecycle is nearly impossible due to the immense timescales and distances over which the interaction occurs. Laboratory astrophysics offers a way to investigate such a system by scaling the important astrophysical parameters to the laboratory. This work describes laboratory experiments to study the radiation-driven implosion of clouds, using x rays from a laser-irradiated, thin, gold foil as a surrogate star and a carbon-foam sphere as a surrogate cloud. An optically thick system, theoretically corresponding to a star-forming regime, was selected by choice of the foam density. Gold foil and sphere motions were imaged by x-ray radiography. Radiographic images show the formation of an interface between rarefied gold and carbon plasmas, a shock moving into the sphere, and a blunting of the initial sphere's shape. Measurements show that the shock moved linearly around 64 μm/ns into the sphere, and the gold–carbon interface formed by 2 ns at the sphere edge remained stationary. The deformation of the sphere was driven by the incident radiation and not by mechanical pressures applied by gold plasma. The blunting of the sphere was likely due to the geometric reduction of flux near the sphere's poles. Higher x-ray flux near the sphere's equator caused high compression and a faster shock, which flattened the sphere. We will discuss the results and implications of our observations.

VanDervort, R. W. [University of Michigan 1 , Ann ↗
DOE OSTI2023

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗
DOE OSTI2020

Primal interface debonding formulation for finite strain isotropic plasticity

In this work, a framework is developed for modeling ductile damage of nonlinear materials whose plastic deformation is characterized using rate independent classical plasticity. This method relies on the assumption that the free energy can be decomposed into elastic, plastic and damage parts. A thermodynamically consistent method is derived which satisfies the second law of thermodynamics in the Clausius–Duhem inequality form. The dissipation associated with plasticity takes place in the domain only, while damage dissipation is localized to the interface. The method is developed using Variational Multiscale ideas to obtain definitions of the interface fluxes within a primal formulation analogous to the Discontinuous Galerkin method, which ensures weakly vanishing interface gap prior to reaching a damage initiation criterion. The local nonlinear problem to calculate both plastic deformation gradient and damage variable follows an incremental approach similar to classical plasticity return mapping algorithm. This elastoplastic damage formulation is developed for material undergoing finite strain, and it naturally accommodates a trapezoidal traction separation law (TSL) whose shape can be varied to model either ductile interface behavior or brittle interface behavior. The formulation's performance is assessed through modeling a patch test and a compact tension specimen.

42 ENGINEERING↗
DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗