Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

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 235 records · Page 13

Numerical Simulations of Flow-Induced Deflections in MITR LEU Fuel Plate Due to Channel Size Disparity

The hydromechanical stability of the fuel plates in parallel coolant channels of a Materials Testing Reactor (MTR) fuel element design is of great importance to the safety of research and test reactors. Previous analytical, experimental, and numerical efforts focused on parallel channels with the same or similar size; also, in the prior numerical simulations, the fuel plate was often assumed to be perfectly flat. This work presents the results of a fluid-structure interaction simulation performed to evaluate the flow-induced deflections of the fuel plates in the low-enriched uranium (LEU, <20 wt% 235 U) fuel element design for the conversion (from highly enriched uranium) of the Massachusetts Institute of Technology Reactor (MITR-II, also referred to as MITR). Various manufacturing and assembly tolerances of the MITR LEU elements are considered in the analysis, and the effects of channel size disparity, nonideal plate shape, and flow rate uncertainty are investigated. Results show that, for all cases analyzed, the deflection occurs toward the larger channel, and the change in any channel stripe remains small (less than 0.021 mm) compared to fabrication tolerances. In addition to simulation work, a hydraulic performance test of the MITR LEU fuel element is currently planned to support conversion to the use of LEU fuel.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Structural modeling and dynamic analysis of condensate storage tanks in nuclear power plants

Condensate storage tanks (CSTs) in nuclear power plants (NPPs) are classified as critical equipment capable of surviving strong shaking in a design basis earthquake to assure the ability to subsequently provide cooling water. Structural modeling and dynamic analysis of CSTs are complicated due to fluid-structure interaction (FSI) and coupling issues, which could have significant effects on seismic response. The capability of a number of approximate models to capture seismic response is investigated using dynamic analysis results from detailed finite element models. Modal analyses and time history analyses are carried out using both 2D and 3D models to predict the dynamic behavior of CSTs as a function of ground motion intensity. Although 2D simplified models can be used to quickly evaluate the dynamic response of CSTs when there is substantial margin to failure, a detailed 3D model is required when it is necessary to examine a limit state associated with the failure mode of the tank. The findings show good agreement between the natural frequencies of convective modes determined from the 2D and 3D CST models. The difference between the frequencies of impulsive modes were found to be equal to or less than 25%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Interaction of fluids and structures for aircraft applications

Strong interactions occur between the flow about an aircraft and its structural components, which result in several important aeroelastic phenomena which can significantly influence performance. For example, aircraft with highly swept wings experience vortex-induced aeroelastic oscillations. The simulation of these complex aeroelastic phenomena requires coupling the fluid and structural analysis. This paper provides a summary of the development of such coupled methods and their applications to aeroelasticity. Results based on the transonic small-perturbation equations and the Euler equations are presented.

Guruswamy, Guru P.↗

DTK C/Fortran Interface Development for NEAMS FSI Simulations

This report documents the development of DataTransferKit (DTK) C and Fortran interfaces for fluid-structure-interaction (FSI) simulations in NEAMS. In these simulations, the codes Nek5000 and Diablo are being coupled within the SHARP framework to study flow-induced vibration (FIV) in reactor steam generators. We will review the current Nek5000/Diablo coupling algorithm in SHARP and the current state of the solution transfer scheme used in this implementation. We will then present existing DTK algorithms which may be used instead to provide an improvement in both flexibility and scalability of the current SHARP implementation. We will show how these can be used within the current FSI scheme using a new set of interfaces to the algorithms developed by this work. These new interfaces currently expose the mesh-free solution transfer algorithms in DTK, a C++ library, and are written in C and Fortran to enable coupling of both Nek5000 and Diablo in their native Fortran language. They have been compiled and tested on Cooley, the test-bed machine for Mira at ALCF.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Simulation and Optimization Framework for Managing Wind-Driven Loading on PV Systems

As PV modules continue to trend toward larger, thinner, and more flexible forms they grow more susceptible to damage from dynamic wind loading. As a result, understanding the impact of wind on PV systems, particularly when mounted on solar-tracking hardware, and identifying robust, stable array layouts and stow strategies is becoming increasingly important for the PV community. We are developing an open-source software package, PVade (PV aerodynamic design engineering), to simulate the cascading fluid-structure interaction that occurs within solar-tracking arrays to enable researchers to test hardware, layout, and tracker control changes, leading to enhanced stability and a reduction in wind-driven damage. We will give an overview of the PVade software and present the latest outcomes from our ongoing validation campaign in which we compare statistical structural responses with field data. From there, we will present simulated results from a larger, multi-row array and highlight the effect of varying tracker angles on stability as measured by both acceleration and deformation, focusing on the stability differences between positive and negative tilt angles.

fluid structure interaction↗

Deflection predictions of involute-shaped fuel plates using a fully-coupled numerical approach

This paper describes the modeling and simulation of fluid structure interactions (FSI) of involute-shaped fuel plates used in nuclear research reactors. We believe this to be the first time that this type of application is described in the literature using a fully-coupled, and monolithic, finite element approach. The simulations are validated against plate deflection data for the conceptual design of the Advanced Neutron Source Reactor (ANSR), which was envisioned to be the world’s most powerful nuclear research reactor for neutron scattering and other applications, but was ultimately never completed. The high performance of the ANSR creates a bounding envelope for involute-shaped research reactors such as that used in the High Flux Isotope Reactor (HFIR) at the Oak Ridge National Laboratory (ORNL) which is undergoing research for the conversion from highly-enriched uranium (HEU) to low-enriched uranium (LEU) fuel. As such, the findings from the present FSI analyses carried out herein for the ANSR plates provide good guidelines and inform designers what should be expected for the next generation of plates in the HFIR. Furthermore it is shown herein that the current approach can accurately capture the leading-edge deflections of the involute-shaped plates and simulations can predict the ‘S-shaped’ deflection of the first mode instilling confidence in the methodology.

42 ENGINEERING↗

Failure investigation and mitigation after experimental research reactor fuel plate deformation in an irradiation device

Experimental research reactor fuel testing is conducted in the Belgian Reactor 2 (BR2) of the Belgian Nuclear Research Centre (SCK CEN) in dedicated irradiation vehicles or rigs. One such vehicle allows flat full-size fuel plates to be irradiated by inserting them into slotted baskets that captures a narrow portion of the longitudinal edges of the plates. The motion of the fuel plates within the baskets is possible within the narrow slots and thus, the plate is considered to be unattached. The design intentionally omits fixing mechanisms of the fuel plates to the baskets to facilitate the inspection and repositioning of the plates between the irradiation cycles and the accommodation of thermal expansion of the plates in the lateral direction. However, loosely inserted fuel plates have weak structural boundary conditions allowing for larger out-of-plane deflections caused by hydrodynamic loads exerted by the flowing coolant, as compared to those of fixed plates. Unexpected large deformations of plates occurred in several irradiation cycles that further resulted in a loss of cladding integrity. These deformations could not be attributed to a single source. This triggered a series of thermal hydraulic, structural, and fluid-structure interaction analyses aiming at understanding the observed phenomenon. The analyses revealed that, for a certain combination of unfavorable manufacturing and assembly tolerances, fuel plate edges could escape out of the slots in the irradiation basket due to the hydrodynamic load. Subsequently, the plate could become wedged inside the basket coolant channel opening. This resulted in reduced coolant flow and accelerated temperature increase and thermal expansion of the plate while under irradiation. This unfavorable feedback loop could then lead to excessive plate surface temperatures, deformed plates and cladding failure, as was observed in the experiments. These analyses not only provided a probable cause of the fuel plate failures, but also resulted in a new and improved design of the irradiation basket to avoid these issues in the future. In conclusion, a series of recent successful irradiations confirm that the sources of failures were identified correctly, and the implemented mitigations were adequate.

BR2↗

Fully Coupled Aeroelastic Stability Analysis of Adaptive Shape Memory Alloy Structural Technologies for Airframe Noise Reduction

This final report documents work performed by ATA Engineering, Inc., (ATA) to develop computational models and analyze the coupled fluid-structure response of two types of noise treatments applied to the leading-edge-slat component of a high-lift system typical of modern transport aircraft. The first treatment is a slat-gap filler (SGF), which closes the gap between the suction surfaces of a deployed slat and an aircraft main wing, and the second treatment is a slat-cove filler (SCF), which replaces the recirculating flow on the slat cove with a surface that promotes flow attachment. The representative airframe chosen for this work was NASA’s High-Lift Common Research Model (CRM-HL) in a baseline high-lift configuration. Superelastic shape memory alloys (SMAs) have been identified as enabling materials for these structural treatments. Since the technology elements rely upon having a highly reconfigurable structure, designs were assessed for their static aeroelastic deflection as well as their dynamic aeroelastic stability using coupled computational fluid dynamics (CFD) and nonlinear computational structural dynamics (NL-CSD) tools. Specifically, fluid-structure interaction (FSI) problems were solved computationally using the CFD solver Loci/CHEM and the NL-CSD solver Abaqus. As a part of the overall project, a similar capability was implemented using the CFD solver FUN3D coupled to Abaqus, although that work is documented in a separate report and that FSI framework was not used to analyze any of NASA’s SGF and SCF configurations. The technical approach consisted of solving for the flow field around the entire vehicle using a global CFD model, followed by extraction of relevant local subdomain data for CFD and NL-CSD co-simulations. The SGF design was analyzed using both 2D and 3D co-simulations to predict quasi-static aeroelastic deformations and to assess dynamic aeroelastic stability, whereas the SCF was analyzed in 2D only. SGF static aeroelastic response predictions focused on characterizing the deformed shape, with maximum displacements predicted to be on the order of magnitude of the technology element panel thickness. SGF dynamic aeroelastic response predictions used Partial Floquet analysis of the temporal evolution of selected nodal displacements to quantify the sign and magnitude of aeroelastic damping. Results suggest that the CRM-HL operating conditions would result in a dynamically stable response. The simulated dynamic pressure was also increased up to a factor of about four, and resulting responses suggest that predicted dynamic stability would be achieved with some margin.

Fluid Structure Interaction↗

Fluid and gas dynamics

A clearer understanding of the interactions of the hot gas flow with the structure in the duct system, the flow passages of the rotating machinery, and the thrust chamber nozzle is sought for the purpose of finding ways and means to increase the life and performance of the systems. The objective of the transient nozzle test is to understand the shock-boundary layer interactions during transient flow and the resulting large side forces acting on the nozzle skirt. The Fluctuating Pressures in Ducts study deals mainly with the fluid-structural interactions of the lox post tube banks in the injector. The Ablative Nozzle Insert Study (H4) has as its goal the performance improvement of the space shuttle main engine.

Povinelli, L.↗

Bayesian inference of nonlinear unsteady aerodynamics from aeroelastic limit cycle oscillations

A Bayesian model selection and parameter estimation algorithm is applied to investigate the influence of nonlinear and unsteady aerodynamic loads on the limit cycle oscillation (LCO) of a pitching airfoil in the transitional Reynolds number regime. At small angles of attack, laminar boundary layer trailing edge separation causes negative aerodynamic damping leading to the LCO. The fluid–structure interaction of the rigid, but elastically mounted, airfoil and nonlinear unsteady aerodynamics is represented by two coupled nonlinear stochastic ordinary differential equations containing uncertain parameters and model approximation errors. Several plausible aerodynamic models with increasing complexity are proposed to describe the aeroelastic system leading to LCO. The likelihood in the posterior parameter probability density function (pdf) is available semi-analytically using the extended Kalman filter for the state estimation of the coupled nonlinear structural and unsteady aerodynamic model. The posterior parameter pdf is sampled using a parallel and adaptive Markov Chain Monte Carlo (MCMC) algorithm. The posterior probability of each model is estimated using the Chib–Jeliazkov method that directly uses the posterior MCMC samples for evidence (marginal likelihood) computation. The Bayesian algorithm is validated through a numerical study and then applied to model the nonlinear unsteady aerodynamic loads using wind-tunnel test data at various Reynolds numbers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Evaluation of flow-induced plate deflection for University of Missouri research reactor low-enriched uranium fuel element

The University of Missouri Research Reactor (MURR), located on the campus of the University of Missouri in Columbia, Missouri, is one of the six United States (U.S.) High Performance Research Reactors (USHPRR), including one critical facility, that are actively collaborating with the U.S. Department of Energy (DOE) National Nuclear Security Administration (NNSA) Office of Material Management and Minimization (M3) Reactor Conversion Program to convert from highly enriched uranium (HEU, ≥20 wt% U-235) fuel to low-enriched uranium (LEU, <20 wt% U-235) fuel. A new type of very high-density LEU fuel based on a monolithic alloy of uranium and 10 wt% molybdenum (U-10Mo) is expected to allow conversion of some USHPRR, including MURR. In the design of its fuel elements, MURR is using thin parallel curved fuel plates separated by coolant channels. In this work, fluid-structure interaction (FSI) analysis of the MURR LEU fuel element is performed at the element level (as compared to the plate level analysis), which models all components of the LEU fuel element, including fuel plates and the supporting structures. Therefore, the effect of supporting structures on the flow distribution within the element and the fuel plate deflection are evaluated. In addition to the element nominal flow rate and dimensions, the tolerances in the geometry of the coolant channel and plate thickness, the effect of a comb on plate deflection, and the uncertainty of the flow rate per element are evaluated. For the LEU fuel plates, which are thinner than the current HEU plates, the predicted plate deflection is found to be small compared to the fabrication and assembly tolerances. Thus, the FSI-induced deflections are not expected to noticeably reduce the coolant flow rate or predicted safety margins in the limiting channels for the MURR LEU fuel element. In addition to the simulation work, a hydraulic performance test of the MURR LEU fuel element is currently being planned to support conversion to the use of LEU fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Immersed Boundary Method Halo Exchange in a Hemodynamics Application

In recent years, highly parallelized simulations of blood flow resolving individual blood cells have been demonstrated. Simulating such dense suspensions of deformable particles in flow often involves a partitioned fluid-structure interaction (FSI) algorithm, with separate solvers for Eulerian fluid and Lagrangian cell grids, plus a solver - e.g., immersed boundary method - for their interaction. Managing data motion in parallel FSI implementations is increasingly important, particularly for inhomogeneous systems like vascular geometries. In this study, we evaluate the influence of Eulerian and Lagrangian halo exchanges on efficiency and scalability of a partitioned FSI algorithm for blood flow. We describe an MPI+OpenMP implementation of the immersed boundary method coupled with lattice Boltzmann and finite element methods. We consider how communication and recomputation costs influence the optimization of halo exchanges with respect to three factors: immersed boundary interaction distance, cell suspension density, and relative fluid/cell solver costs.

Gounley, John↗

ExaWind: A multifidelity modeling and simulation environment for wind energy

We introduce the open-source ExaWind modeling and simulation environment for wind energy. The primary physics codes of ExaWind are Nalu-Wind and OpenFAST. Nalu-Wind is a wind-focused computational fluid dynamics (CFD) code that is coupled to the whole-turbine simulation code OpenFAST. The ExaWind environment was created under U.S. Department of Energy funding to achieve the highest-fidelity simulations of wind turbines and wind farms to date, with the goal of enabling disruptive changes to turbine and plant design and operation. Innovation will be gleaned through better understanding of the complex flow dynamics in wind farms, including wake evolution and the impact of wakes on downstream turbines and turbulent flow from complex terrain. High-fidelity predictive simulations employ hybrid turbulence models, geometry/boundary-layer-resolving CFD meshes, atmospheric turbulence, nonlinear structural dynamics, and fluid-structure interaction. While there is an emphasis on very high-fidelity simulations (e.g., blade resolved with full fluid-structure coupling), the ExaWind environment supports lower-fidelity modeling capabilities including actuator-line and -disk methods. Important in the development of ExaWind codes is that the codes scale well on today's largest petascale supercomputers and on the next-generation platforms that will enable exascale computing.

17 WIND ENERGY↗

Analysis and control of low-speed forced unsteady flow

A capability for numerically simulating 2-D flows in temporally deforming geometries is described, with emphasis on flow with forced unsteadiness, particularly on the simulation and analysis of these flows. The simulation of forced unsteady flows makes the examination of fundamental unsteady flow mechanisms, such as dynamic stall and unsteady separation, possible. A turbulence model is being incorporated into the analysis so as to obtain solutions for the higher Reynolds numbers used in the experiments. The analysis is also of utility in studying fluid-structure interactions, free surfaces, metal-forming, and bio-fluid mechanics involving flow through passages with flexible walls.

Ghia, U.↗

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗

Infrastructure for Multiphysics Software Integration in High Performance Computing-Aided Science and Engineering

The project described in this report constructed and exercised an innovative multiphysics coupling toolkit called the Illinois Rocstar MultiPhysics Application Coupling Toolkit (IMPACT). IMPACT is an open source, flexible, natively parallel infrastructure for coupling multiple uniphysics simulation codes into multiphysics computational systems. IMPACT works with codes written in several high-performance-computing (HPC) programming languages, and is designed from the beginning for HPC multiphysics code development. It is designed to be minimally invasive to the individual physics codes being integrated, and has few requirements on those physics codes for integration. The goal of IMPACT is to provide the support needed to enable coupling existing tools together in unique and innovative ways to produce powerful new multiphysics technologies without extensive modification and rewrite of the physics packages being integrated. There are three major outcomes from this project: 1) construction, testing, application, and open-source release of the IMPACT infrastructure, 2) production of example open-source multiphysics tools using IMPACT, and 3) identification and engagement of interested organizations in the tools and applications resulting from the project. This last outcome represents the incipient development of a user community and application echosystem being built using IMPACT. Multiphysics coupling standardization can only come from organizations working together to define needs and processes that span the space of necessary multiphysics outcomes, which Illinois Rocstar plans to continue driving toward. The IMPACT system, including source code, documentation, and test problems are all now available through the public gitHUB.org system to anyone interested in multiphysics code coupling. Many of the basic documents explaining use and architecture of IMPACT are also attached as appendices to this document. Online HTML documentation is available through the gitHUB site. There are over 100 unit tests provided that run through the Illinois Rocstar Application Development (IRAD) lightweight testing infrastructure that is also supplied along with IMPACT. The package as a whole provides an excellent base for developing high-quality multiphysics applications using modern software development practices. To facilitate understanding how to utilize IMPACT effectively, two multiphysics systems have been developed and are available open-source through gitHUB. The simpler of the two systems, named ElmerFoamFSI in the repository, is a multiphysics, fluid-structure-interaction (FSI) coupling of the solid mechanics package Elmer with a fluid dynamics module from OpenFOAM. This coupling illustrates how to combine software packages that are unrelated by either author or architecture and combine them into a robust, parallel multiphysics system. A more complex multiphysics tool is the Illinois Rocstar Rocstar Multiphysics code that was rebuilt during the project around IMPACT. Rocstar Multiphysics was already an HPC multiphysics tool, but now that it has been rearchitected around IMPACT, it can be readily expanded to capture new and different physics in the future. In fact, during this project, the Elmer and OpenFOAM tools were also coupled into Rocstar Multiphysics and demonstrated. The full Rocstar Multiphysics codebase is also available on gitHUB, and licensed for any organization to use as they wish. Finally, the new IMPACT product is already being used in several multiphysics code coupling projects for the Air Force, NASA and the Missile Defense Agency, and initial work on expansion of the IMPACT-enabled Rocstar Multiphysics has begun in support of a commercial company. These initiatives promise to expand the interest and reach of IMPACT and Rocstar Multiphysics, ultimately leading to the envisioned standardization and consortium of users that was one of the goals of this project.

42 ENGINEERING↗

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Wind Turbine Rotor Design Using High-Fidelity Aerostructural Optimization

Large wind turbines yield more energy but demand careful aeroelastic blade design. Coupled multiphysics design strategies can reduce wind energy costs by exploiting fluid-structure interactions. This work presents the first high-fidelity aerostructural optimization study of a large wind turbine rotor. We use blade-resolved fluid dynamics and structural solvers in a monolithic gradient-based optimization framework to explore steady-state torque and blade mass tradeoffs. The coupled-adjoint approach computes gradients efficiently, enabling the optimization of over 100 structural and geometric parameters simultaneously. Our optimization study modifies a DTU 10 MW benchmark with a simplified structure and isotropic material properties. The tightly coupled optimizations increase torque by 14% while reducing rotor mass by 9% or reduce blade mass by 27% while maintaining torque. Blade-resolved models provide greater design freedom, enabling 5% higher mass reductions than conventional parameterizations at equal torque. This framework paves the way for more detailed high-fidelity optimization studies to complement conventional design approaches.

17 WIND ENERGY↗