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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Geometric multigrid for an implicit-time immersed boundary method

The immersed boundary (IB) method is an approach to fluid-structure interaction that uses Lagrangian variables to describe the deformations and resulting forces of the structure and Eulerian variables to describe the motion and forces of the fluid. Explicit time stepping schemes for the IB method require solvers only for Eulerian equations, for which fast Cartesian grid solution methods are available. Such methods are relatively straightforward to develop and are widely used in practice but often require very small time steps to maintain stability. Implicit-time IB methods permit the stable use of large time steps, but efficient implementations of such methods require significantly more complex solvers that effectively treat both Lagrangian and Eulerian variables simultaneously. Moreover, several different approaches to solving the coupled Lagrangian-Eulerian equations have been proposed, but a complete understanding of this problem is still emerging. This paper presents a geometric multigrid method for an implicit-time discretization of the IB equations. This multigrid scheme uses a generalization of box relaxation that is shown to handle problems in which the physical stiffness of the structure is very large. Numerical examples are provided to illustrate the effectiveness and efficiency of the algorithms described herein. Finally, these tests show that using multigrid as a preconditioner for a Krylov method yields improvements in both robustness and efficiency as compared to using multigrid as a solver. They also demonstrate that with a time step 100–1000 times larger than that permitted by an explicit IB method, the multigrid-preconditioned implicit IB method is approximately 50–200 times more efficient than the explicit method.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

High-order, linearly stable, partitioned solvers for general multiphysics problems based on implicit–explicit Runge–Kutta schemes

This work introduces a general framework for constructing high-order, linearly stable, partitioned solvers for multiphysics problems from a monolithic implicit–explicit Runge–Kutta (IMEX-RK) discretization of the semi-discrete equations. The generic multiphysics problem is modeled as a system of n systems of partial differential equations where the i th subsystem is coupled to the other subsystems through a coupling term that can depend on the state of all the other subsystems. This coupled system of partial differential equations reduces to a coupled system of ordinary differential equations via the method of lines where an appropriate spatial discretization is applied to each subsystem. The coupled system of ordinary differential equations is taken as a monolithic system and discretized using an IMEX-RK discretization with a specific implicit–explicit decomposition that introduces the concept of a predictor for the coupling term. We propose four coupling predictors that enable the monolithic system to be solved in a partitioned manner, i.e., subsystem-by-subsystem, and preserve the IMEX-RK structure and therefore the design order of accuracy of the monolithic scheme. The four partitioned solvers that result from these predictors are high-order accurate, allow for maximum re-use of existing single-physics software, and two of the four solvers allow the subsystems to be solved in parallel at a given stage and time step. We also analyze the stability of a coupled, linear model problem with a specific coupling structure and show that one of the partitioned solvers achieves unconditional linear stability for this problem, while the others are unconditionally stable only for certain values of the coupling strength. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection–diffusion–reaction systems, fluid–structure interaction problems, and particle-laden flows, where we verify the design order of the IMEX schemes and study various stability properties.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Parallel aeroelastic computations for wing and wing-body configurations

The objective of this research is to develop computationally efficient methods for solving fluid-structural interaction problems by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures on parallel computers. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗

DOE OSTI2021

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↗

DOE OSTI2012

Temperature dependent droplet impact dynamics on flat and textured surfaces

Droplet impact dynamics determines the performance of surfaces used in many applications such as anti-icing, condensation, boiling and heat transfer. We study impact dynamics of water droplets on surfaces with chemistry/texture ranging from hydrophilic to superhydrophobic and across a temperature range spanning below freezing to near boiling conditions. Droplet retraction shows very strong temperature dependence especially for hydrophilic surfaces; it is seen that lower substrate temperatures lead to lesser retraction. Physics-based analyses show that the increased viscosity associated with lower temperatures can explain the decreased retraction. The present findings serve to guide further studies of dynamic fluid-structure interaction at various temperatures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1995

Aeroelasticity of wing and wing-body configurations on parallel computers

The objective of this research is to develop computationally efficient methods for solving aeroelasticity problems on parallel computers. Both uncoupled and coupled methods are studied in this research. For the uncoupled approach, the conventional U-g method is used to determine the flutter boundary. The generalized aerodynamic forces required are obtained by the pulse transfer-function analysis method. For the coupled approach, the fluid-structure interaction is obtained by directly coupling finite difference Euler/Navier-Stokes equations for fluids and finite element dynamics equations for structures. This capability will significantly impact many aerospace projects of national importance such as Advanced Subsonic Civil Transport (ASCT), where the structural stability margin becomes very critical at the transonic region. This research effort will have direct impact on the High Performance Computing and Communication (HPCC) Program of NASA in the area of parallel computing.

Byun, Chansup↗

DOE OSTI2017

A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis

A stable partitioned algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and rigid bodies. This added-mass partitioned (AMP) algorithm remains stable, without sub-iterations, for light and even zero mass rigid bodies when added-mass and viscous added-damping effects are large. The scheme is based on a generalized Robin interface condition for the fluid pressure that includes terms involving the linear acceleration and angular acceleration of the rigid body. Added mass effects are handled in the Robin condition by inclusion of a boundary integral term that depends on the pressure. Added-damping effects due to the viscous shear forces on the body are treated by inclusion of added-damping tensors that are derived through a linearization of the integrals defining the force and torque. Added-damping effects may be important at low Reynolds number, or, for example, in the case of a rotating cylinder or rotating sphere when the rotational moments of inertia are small. In this second part of a two-part series, the general formulation of the AMP scheme is presented including the form of the AMP interface conditions and added-damping tensors for general geometries. A fully second-order accurate implementation of the AMP scheme is developed in two dimensions based on a fractional-step method for the incompressible Navier-Stokes equations using finite difference methods and overlapping grids to handle the moving geometry. Here, the numerical scheme is verified on a number of difficult benchmark problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2015

A simplified model for dynamics of cell rolling and cell-surface adhesion

We propose a three dimensional model for the adhesion and rolling of biological cells on surfaces. We study cells moving in shear flow above a wall to which they can adhere via specific receptor-ligand bonds based on receptors from selectin as well as integrin family. The computational fluid dynamics are governed by the lattice-Boltzmann method. The movement and the deformation of the cells is described by the immersed boundary method. Both methods are fully coupled by implementing a two-way fluid-structure interaction. The adhesion mechanism is modelled by adhesive bonds including stochastic rules for their creation and rupture. We explore a simplified model with dissociation rate independent of the length of the bonds. We demonstrate that this model is able to resemble the mesoscopic properties, such as velocity of rolling cells.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

An Evaluation of An Asynchronous Task Based Dataflow Approach For Uintah

The challenge of running complex physics code on the largest computers available has led to dataflow paradigms being explored. While such approaches are often applied at smaller scales, the challenge of extreme-scale data flow computing remains. The Uintah dataflow framework has consistently used dataflow computing at the largest scales on complex physics applications. At present Uintah contains two main dataflow models. Both are based upon asynchronous communication. One uses a static graph-based approach with asynchronous communication and the other uses a more dynamic approach that was introduced almost a decade ago. Subsequent changes within the Uintah runtime system combined with many more large scale experiments, has necessitated a reevaluation of these two approaches, comparing them in the context of large scale problems. While the static approach has worked well for some large-scale simulations, the dynamic approach is seen to offer performance improvements over the static case for a challenging fluid-structure interaction problem at large scale that involves fluid flow and a moving solid represented using particle method on an adaptive mesh.

Humphrey, Alan↗

DOE OSTI2012

IMPOSING A LAGRANGIAN PARTICLE FRAMEWORK ON AN EULERIAN HYDRODYNAMICS INFRASTRUCTURE IN FLASH

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid-structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

New investigations of the Phenix negative reactivity events

In 1989 and 1990, four very fast and high amplitude decreases of the measured neutron flux, known as 'AURN' events, occurred in the Phenix fast neutron reactor and were not well explained despite a significant investigation program carried out in the 90 s. The perspective of the Phenix end of life led the French Nuclear Energy Div. to reinitiate AURN studies in particular with dedicated Phenix tests. It comes out from simulations with improved fluid-structure interaction model of the core movement that an outward movement of the sub-assemblies can generate a power variation consistent with the AURN signals. 'AURN' Phenix tests confirm the abnormal thermal behavior of an experimental sub-assembly, a new A URN scenario has been proposed on this basis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

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↗

DOE OSTI2019

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. 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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Predicting Instability and the Effect of Wind Loading on Single-Axis Trackers

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. In our ongoing DuraMAT project, 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, highlighting recent user-interface and algorithm developments, before presenting the latest outcomes from our ongoing validation campaign in which we analyze and compare with experimental data obtained from a DuraMAT 1 project. From there, we will present simulated results from a larger, multi-row array and highlight relationships between varying tracker angles and stability as measured by different experimentally validated metrics.

fluid structure interaction↗

DOE OSTI2016

Moving Beyond Traditional Tornado Missile Impact Analysis

On March 11, 2011 a 9.0 magnitude earthquake occurred off the coast of Japan creating a tsunami, which together severely damaged northeastern Japan including the Fukushima Dai-ichi nuclear power plant site. In response to this event, the United States (U.S.) Nuclear Regulatory Commission (NRC) issued Order EA-12-049, 'Issuance of Order to Modify Licenses with Regard to Requirements for Mitigation Strategies for Beyond-Design-Basis External Events'. In order to meet the guidance of the NRC Order EA-12- 049, the a central U.S. nuclear generating station (Plant) was required to demonstrate that the currently installed Borated Water Storage Tank (BWST), Refueling Water Tank (RWT), and Qualified Condensate Storage Tank (QCST) could provide the necessary capacity following a Beyond Design Basis External Event (BDBEE). These tanks were not originally designed to meet current site design basis requirements for tornado wind and missile impacts. To save the Plant from having to build costly new storage tanks off site, the existing tanks were credited as robust for tornado missile impact through advanced finite element impact analyses using the ANSYS Autodyn solver and fluid-structure interaction (FSI). This procedure is the first known NRC acceptance of advanced analysis and FSI for to satisfy the FLEX tornado missile impact criteria. Extensive verification and validation (V and V) was performed to verify perforation of steel plates which included FSI. These V and V analyses were confirmed against large databases of experimental and theoretical results to ensure accurate modeling and results. The use of this advanced analytical method benefitted the Plant in several ways. A new storage tank was not required to be built, saving not only initial construction costs but also ongoing maintenance costs and monitoring activities. By determining specific vulnerable components of each tank a set of modifications and protection schemes were developed. These modifications were much quicker to implement saving the Plant on downtime costs. (authors)

42 ENGINEERING↗

DOE OSTI2015

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 ↗

DOE OSTI2023

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↗

DOE OSTI2020

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↗