Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLUID-STRUCTURE INTERACTIONS”

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 163 records · Page 9

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↗

Highly Accurate Experimental Heave Decay Tests with a Floating Sphere: A Public Benchmark Dataset for Model Validation of Fluid–Structure Interaction

Highly accurate and precise heave decay tests on a sphere with a diameter of 300 mm were completed in a meticulously designed test setup in the wave basin in the Ocean and Coastal Engineering Laboratory at Aalborg University, Denmark. The tests were dedicated to providing a rigorous benchmark dataset for numerical model validation. The sphere was ballasted to half submergence, thereby floating with the waterline at the equator when at rest in calm water. Heave decay tests were conducted, wherein the sphere was held stationary and dropped from three drop heights: a small drop height, which can be considered a linear case, a moderately nonlinear case, and a highly nonlinear case with a drop height from a position where the whole sphere was initially above the water. The precision of the heave decay time series was calculated from random and systematic standard uncertainties. At a 95% confidence level, uncertainties were found to be very low—on average only about 0.3% of the respective drop heights. Physical parameters of the test setup and associated uncertainties were quantified. A test case was formulated that closely represents the physical tests, enabling the reader to do his/her own numerical tests. The paper includes a comparison of the physical test results to the results from several independent numerical models based on linear potential flow, fully nonlinear potential flow, and the Reynolds-averaged Navier–Stokes (RANS) equations. A high correlation between physical and numerical test results is shown. The physical test results are very suitable for numerical model validation and are public as a benchmark dataset.

42 ENGINEERING↗

High-Speed X-Ray Stereo Digital Image Correlation in a Shock Tube

X-ray stereo digital image correlation (DIC) measurements were performed at 10 kHz on the internal surface of a jointed structure in a shock tube at a shock Mach number of 1.42 and compared with optical stereo DIC measurements on the outer, visible surface of the structure. The shock tube environment introduces temperature and density gradients in the gas through which the structure was imaged, resulting in spatial and temporal index of refraction variations. These variations cause bias errors in optical DIC measurements due to beam-steering but have minimal influence on x-ray DIC measurements. These results demonstrate the utility of time-resolved x-ray DIC measurements in complicated environments where optical measurements suffer severe errors and/or are precluded by lack of optical access.

42 ENGINEERING↗

CFD Simulation of Floating Body Motion with Mooring Dynamics: Coupling MoorDyn with OpenFOAM

It is increasingly popular to use Computational Fluid Dynamics (CFD) models to study floating structures subjected to ocean waves, especially when it comes to applications of floating offshore wind turbines and Wave Energy Converters (WECs). Mooring dynamics are currently lacking in most of these applications. This paper presents a coupled simulation study of moored floating body motion by coupling two open-source libraries: a finite volume CFD toolbox, OpenFOAM, and a lumped-mass mooring model, MoorDyn. The instantaneous floating body position and velocity are passed from the body motion solver in the CFD model to the mooring model to calculate the fairlead kinematics. The mooring reaction forces, which are calculated by MoorDyn after updating the mooring system states, are then returned to the body motion solver to update the floating body motion. Both mesh deformation and overset mesh methods are used as the mesh motion solver in the CFD model to account for the floating body motion. The coupled model was validated against experimental measurements for a floating box moored with four catenary lines under the action of regular waves, which came from a preliminary test campaign for WECs. Apart from the lumped-mass mooring model, the present work also coupled a quasi-static mooring model and a finite element model with the floating body motion solver in OpenFOAM. The mooring line tensions predicted by these models were compared. The coupled model equipped with three mooring line codes may be further used to carry out survivability studies of FOWTs and WECs subject to severe sea states.

floating body↗

High fidelity simulation and validation of crossflow through a tube bundle and the onset of vibration

Flow induced vibrations can be detrimental to all manner of engineering applications with fluid flow. Heat exchangers with bundles of thin tubes experiencing a crossflow are especially susceptible to the unexpected onset of vibration. Heat-exchanger design necessarily involves extensive modeling and testing to ensure that significant vibrations cannot develop for any expected flow conditions. A properly validated numerical simulation can work in conjunction with physical experiments to identify problematic vibrations and allow for rapid iteration of design at relatively low expense. A high-fidelity fluid–structure interaction code has been developed by fully coupling CFD LES/DNS code Nek5000 and CSM code Diablo. The coupled code is used to simulate crossflow through a tube bundle in a geometry recreated after a physical experiment. Validation against this experiment involves comparing the amplitude and frequency spectrum for three different flow velocities. The velocities compared straddle the onset of large magnitude vibrations. The simulation accurately captures the onset of the vibrations with a marginally lower predicted frequency of vibration. When the onset velocity is compensated for the difference in natural frequency, the simulation results closely match the experiment. Finally, the continuous nature of the simulation measurements helps illustrate the fluid mechanics behind the pin motion and reveals the front to back propagation of the vibrations as the fluid velocity increases.

42 ENGINEERING↗

An added-mass partition algorithm for fluid–structure interactions of compressible fluids and nonlinear solids

We describe an added-mass partitioned (AMP) algorithm for solving fluid–structure interaction (FSI) problems involving inviscid compressible fluids interacting with nonlinear solids that undergo large rotations and displacements. The computational approach is a mixed Eulerian–Lagrangian scheme that makes use of deforming composite grids (DCG) to treat large changes in the geometry in an accurate, flexible, and robust manner. The current work extends the AMP algorithm developed in Banks et al. [1] for linearly elasticity to the case of nonlinear solids. To ensure stability for the case of light solids, the new AMP algorithm embeds an approximate solution of a nonlinear fluid–solid Riemann (FSR) problem into the interface treatment. The solution to the FSR problem is derived and shown to be of a similar form to that derived for linear solids: the state on the interface being fundamentally an impedance-weighted average of the fluid and solid states. Numerical simulations demonstrate that the AMP algorithm is stable even for light solids when added-mass effects are large. The accuracy and stability of the AMP scheme is verified by comparison to an exact solution using the method of analytical solutions and to a semi-analytical solution that is obtained for a rotating solid disk immersed in a fluid. The scheme is applied to the simulation of a planar shock impacting a light elliptical-shaped solid, and comparisons are made between solutions of the FSI problem for a neo-Hookean solid, a linearly elastic solid, and a rigid solid. The ability of the approach to handle large deformations is demonstrated for a problem of a high-speed flow past a light, thin, and flexible solid beam.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A cut-cell finite volume – finite element coupling approach for fluid–structure interaction in compressible flow

We present a loosely coupled approach for the solution of fluid–structure interaction problems between a compressible flow and a deformable structure. The method is based on staggered Dirichlet–Neumann partitioning. The interface motion in the Eulerian frame is accounted for by a conservative cut-cell Immersed Boundary method. The present approach enables sub-cell resolution by considering individual cut-elements within a single fluid cell, which guarantees an accurate representation of the time-varying solid interface. The cut-cell procedure inevitably leads to non-matching interfaces, demanding for a special treatment. A Mortar method is chosen in order to obtain a conservative and consistent load transfer. We validate our method by investigating two-dimensional test cases comprising a shock-loaded rigid cylinder and a deformable panel. Moreover, the aeroelastic instability of a thin plate structure is studied with a focus on the prediction of flutter onset. Finally, we propose a three-dimensional fluid–structure interaction test case of a flexible inflated thin shell interacting with a shock wave involving large and complex structural deformations.

42 ENGINEERING↗

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↗

A fast platform for simulating semi-flexible fiber suspensions applied to cell mechanics

We present a novel platform for the large-scale simulation of three-dimensional fibrous structures immersed in a Stokesian fluid and evolving under confinement or in free-space in three dimensions. One of the main motivations for this work is to study the dynamics of fiber assemblies within biological cells. For this, we also incorporate the key biophysical elements that determine the dynamics of these assemblies, which include the polymerization and depolymerization kinetics of fibers, their interactions with molecular motors and other objects, their flexibility, and hydrodynamic coupling. This work, to our knowledge, is the first technique to include many-body hydrodynamic interactions (HIs), and the resulting fluid flows, in cellular assemblies of flexible fibers. We use non-local slender body theory to compute the fluid–structure interactions of the fibers and a second-kind boundary integral formulation for other rigid bodies and the confining boundary. A kernel-independent implementation of the fast multipole method is utilized for efficient evaluation of HIs. The deformation of the fibers is described by nonlinear Euler–Bernoulli beam theory and their polymerization is modeled by the reparametrization of the dynamic equations in the appropriate non-Lagrangian frame. We use a pseudo-spectral representation of fiber positions and implicit time-stepping to resolve large fiber deformations, and to allow time-steps not excessively constrained by temporal stiffness or fiber–fiber interactions. The entire computational scheme is parallelized, which enables simulating assemblies of thousands of fibers. We use our method to investigate two important questions in the mechanics of cell division: (i) the effect of confinement on the hydrodynamic mobility of microtubule asters; and (ii) the dynamics of the positioning of mitotic spindle in complex cell geometries. Finally to demonstrate the general applicability of the method, we simulate the sedimentation of a cloud of semi-flexible fibers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An immersed boundary method for fluid–structure interaction with compressible multiphase flows

This paper presents a two-dimensional immersed boundary method for fluid–structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid–structure interaction coupling. In the flow solver, the compressible multiphase Navier–Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge–Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A computational risk assessment approach to the integration of seismic and flooding hazards with internal hazards

Probabilistic risk assessment (PRA) of nuclear power plants historically focuses on internal events at the plant, rather than external hazards. Although the importance of external hazards risk analysis is now well recognized, the methods for assessing the risk associated with low probability external hazards rely heavily on subjective judgment of specialists, often using conservative elements in the analysis, as is the case with the Conservative Deterministic Failure Margin approach commonly used for seismic PRA. The U.S. Department of Energy Light Water Reactor Sustainability program has been investigating the use of first principles analyses to provide a more realistic assessment of the risk of external hazards, in a process called computational risk assessment (CRA). Efforts undertaken are described in the development of CRA methods to integrate the risk associated with seismic and flooding events into traditional PRA. Further, the results of four case studies are presented.

99 GENERAL AND MISCELLANEOUS↗

Validation of numerical models for seismic fluid-structure-interaction analysis of nuclear, safety-related equipment

Seismic design, qualification, and risk assessment of nuclear safety-related vessels filled with liquid will have to consider the interaction between the vessel (tank), the contained liquid, and submerged components, if any. Seismic fluid-structure-interaction (FSI) analysis of nuclear vessels will rely on numerical models, which are required to be verified and validated. This study validates previously-verified numerical models using test data generated from earthquake-simulator experiments involving a cylindrical tank. Two solvers in LS-DYNA for FSI simulations are used for the numerical analysis: Arbitrary-Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD). Numerical and test results are compared for responses critical to seismic design of advanced reactor vessels: hydrodynamic pressures on the tank wall, reactions at the support, and wave heights of the contained liquid. Analysis is performed for one-, two-, and three-directional seismic inputs with a range of intensity, and rocking motions. The accuracy of the numerical results and different methods for outputting wave heights are discussed. Recommendations for validation of seismic FSI numerical models of advanced reactors are provided. As a result, the validation exercise presented in this study is broadly applicable to cylindrical tanks, regardless of industry sector.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A beam-membrane structure micromachined differential pressure flow sensor

A beam-membrane structure micromachined flow sensor is designed, depending on the principle of differential pressure caused by the mass flow, which is directly proportional to the square flow rate. The FSI (fluid structure interaction) characteristics of the differential pressure flow sensor are investigated via numerical analysis and analog simulation. The working mechanism of the flow sensor is analyzed depending on the FSI results. Then, the flow sensor is fabricated and calibrated. The calibration results show that the beam-membrane structure differential pressure flow sensor achieves ideal static characteristics and works well in the practical applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Size modulated transition in the fluid–structure interaction losses in nano mechanical beam resonators

An understanding of the dominant dissipative mechanisms is crucial for the design of a high-Q doubly clamped nanobeam resonator to be operated in air. We focus on quantifying analytically the viscous losses—the squeeze film damping and drag force damping—that limit the net quality factor of a beam resonator, vibrating in its flexural fundamental mode with the surrounding fluid as air at atmospheric pressure. Specifically, drag force damping dominates at smaller beam widths and squeeze film losses dominate at larger beam widths, with no significant contribution from structural losses and acoustic radiation losses. The combined viscous losses agree well with the experimentally measured Q of the resonator over a large range of beam widths, within the limits of thin beam theory. We propose an empirical relation between the maximum quality factor and the ratio of maximum beam width to the squeeze film air gap thickness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

Significance Fluid–structure interactions are ubiquitous in many natural and man-made environments. They are difficult to study analytically, and therefore accurate and flexible computational methods are an indispensable tool in the field. Typically, fluids are simulated with a fixed background computational mesh, whereas a solid is simulated with a mesh that moves with it, making it challenging to couple the two. Here we develop a three-dimensional computational method where both fluid and solid can be represented on a fixed computational grid, which simplifies the coupling between the two phases considerably. Our method can simulate scenarios that are difficult to do with other methods, such as complex suspensions containing soft particles that are both heavier and lighter than the fluid.

3D fluid-structure interaction↗

Reducing Uncertainty in Hydrodynamic Modeling of ATR Experiments Via Flow Testing, Validation, and Optimization

The development, characterization, and qualification testing of nuclear fuel at Idaho National Laboratory’s Advanced Test Reactor (ATR) requires extensive design and analysis activities prior to the insertion of an irradiation experiment in-pile. Significant effort is made in the design and development phase of all in-pile experiments to ensure that the maximum feasible impacts of all necessary experimental requirements are satisfied. The advancement of fuel, cladding, and in-reactor materials technology in recent years has introduced complexities associated with the design and construct of in-pile experiments necessitating deeper understanding of boundary conditions and increasingly comprehensive observations resulting from the experiment. Each unique experiment must be assessed for neutronics response, thermal/hydraulic/hydrodynamic performance, and structural integrity. This is accomplished either analytically, computationally, or experimentally, or some combination thereof, prior to insertion into the ATR. The various effects are interrelated to various degrees, such as the case with the experiment temperature affecting the thermal cross section of the fuel or the increased temperature of the experiment’s materials reducing the mechanical strength of the assemblies. Additionally, the feedback between the experiment’s response to a reactor transient could alter the neutron flux profile of the reactor during the transient. Each experiment must therefore undergo a barrage of analyses to assure the ATR operational safety review committee that the insertion and irradiation of the experiment will not detrimentally affect the safe operational envelope of the reactor. In many cases, the nuclear fuel being tested can be double-encapsulated to ensure safety margins are adequately addressed, whereas failed fuel would be encased in a protective capsule. In other cases, the experiments can be inserted in a self-contained loop that passes through the reactor core, remaining isolated from the primary coolant. In the case of research reactor fuel, however, the fuel plates must be tested in direct contact with the reactor coolant, and being fuel designed for high neutron fluxes, they are inherently power-dense plates. The combination of plate geometry, high-power density, and direct contact with primary coolant creates a scenario where the neutronic/thermomechanic/hydrodynamic characteristics of the fuel plates are tightly coupled, necessitating as complete characterization as possible to support the safety and programmatic assessments, thus enabling a successful experiment. This paper explores the efforts of the U.S. High-Performance Research Reactor program to thermomechanically/hydromechanically characterize the program’s wide variety of experiments, which range from stacks of miniplate capsules to full-sized, geometrically representative curved plates. Special attention is given to instances where the combination of experimental characterization and analytical assessment has reduced uncertainties of the safety margins, allowing experiments to be irradiated that would otherwise not have passed the rigorous qualification process for irradiation in the ATR. In some cases, the combined processes have exposed flow and heat transfer characteristics that would have been missed using historical methods, which allows for more accurate and representative postirradiation assessments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗