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 199 records · Page 11

Modeling Concept and Numerical Simulation of Ultrasonic Wave Propagation in a Moving Fluid-Structure Domain based on a Monolithic Approach

Here in the present study, we propose a novel multiphysics model that merges two time-dependent problems – the Fluid-Structure Interaction (FSI) and the ultrasonic wave propagation in a fluid-structure domain with a one directional coupling from the FSI problem to the ultrasonic wave propagation problem. This model is referred to as the “eXtended fluid-structure interaction (eXFSI)” problem. This model comprises isothermal, incompressible Navier-Stokes equations with nonlinear elastodynamics using the Saint-Venant Kirchhoff solid model. The ultrasonic wave propagation problem comprises monolithically coupled acoustic and elastic wave equations. To ensure that the fluid and structure domains are conforming, we use the ALE technique. The solution principle for the coupled problem is to first solve the FSI problem and then to solve the wave propagation problem. Accordingly, the boundary conditions for the wave propagation problem are automatically adopted from the FSI problem at each time step. The overall problem is highly nonlinear, which is tackled via a Newton-like method. The model is verified using several alternative domain configurations. To ensure the credibility of the modeling approach, the numerical solution is contrasted against experimental data.

97 MATHEMATICS AND COMPUTING↗

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗

An operator splitting approach for the interaction between a fluid and a multilayered poroelastic structure

We develop a loosely coupled fluid‐structure interaction finite element solver based on the Lie operator splitting scheme. The scheme is applied to the interaction between an incompressible, viscous, Newtonian fluid, and a multilayered structure, which consists of a thin elastic layer and a thick poroelastic material. The thin layer is modeled using the linearly elastic Koiter membrane model, while the thick poroelastic layer is modeled as a Biot system. We prove a conditional stability of the scheme and derive error estimates. Theoretical results are supported with numerical examples. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1054–1100, 2015

Bukač, Martina↗

Fluid-thermal-structural interaction of aerodynamically heated leading edges

A two-dimensional finite element approach is presented for the integrated fluid-thermal-structural analysis of aerodynamically heated leading edges. The approach is combined with an adaptive unstructured remeshing technique to solve the Navier-Stokes equations for high speed compressible flow, the energy equation for the structure thermal response, and the quasi-static equilibrium equations for the structural response. Coupling and interaction between the three disciplines are demonstrated using two applications for high speed flow over a cylinder and a simulated engine leading edge verification test.

Dechaumphai, Pramote↗

Method advancements for fuel assessments during seismic and LOCA events

The evaluation of fuel assembly structural responses of light-water reactors to seismic and LOCA excitations has received increased attention in the last five years. Real seismic events such as those that affected the Fukushima Daiichi and North Anna plants in 2011 have caused the industry to more closely examine our preparedness for seismic events. Furthermore, the review of design certifications for new plants has raised questions about the methods used to evaluate these events and is leading to revised guidance from regulators (e.g. US NRC Information Notice 2012-09). As a result, the industry is seeing fuel structural design margins based on current methods begin to erode for seismic and LOCA events. This paper will review points for consideration as the industry responds to these changes. To begin, it is worthwhile for the industry to understand the history and the critical thought behind the current methods such that improved criteria and methods can be formulated. This review identifies three known areas of weakness in the current framework for this methodology. First, a broader understanding of spacer grid behavior is needed. Spacer grid behavior in response to impact loads is of paramount importance, but the industry needs to re-frame our characterization of this behavior in order to focus on the actual critical safety function of this component. Secondly, despite guidance from the Standard Review Plan, the effects of irradiation, on spacer grid characteristics, in particular, must be considered. Thirdly, the area of fluid-structure-interaction may be the greatest frontier to conquer in this methodology and may provide our best possibility to recover eroded margins. However, fluid-structure-interactions must be cautiously evaluated and applied with conservatism. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Dynamic responses of submerged components in advanced reactors: experimental and numerical studies

The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.

Engineering↗

Elastic response of water-filled fiber composite tubes under shock wave loading

We experimentally and numerically investigate the response of fluid-filled filament-wound composite tubes subjected to axial shock wave loading in water. Our study focuses on the fluid structure interaction occurring when the shock wave in the fluid propagates parallel to the axis of the tube, creating pressure waves in the fluid coupled to flexural waves in the shell. The in-house-developed computational scheme couples an Eulerian fluid solver with a Lagrangian shell solver, which includes a new and simple material model to capture the response of fiber composites in finite kinematics. In the experiments and simulations we examine tubes with fiber winding angles equal to 45 and 60 , and we measure the precursor and primary wave speeds, hoop and longitudinal strains, and pressure. The experimental and computational results are in agreement, showing the validity of the computational scheme in complex fluid structure interaction problems involving fiber composite materials subjected to shock waves. The analyses of the measured quantities show the strong coupling of axial and hoop deformations and the significant effect of fiber winding angle on the composite tube response, which differs substantially from that of a metal tube in the same configuration.

Perotti, Luigi E.↗

Fluid-Thermal-Structural Interactions in Ramp-Induced Shock-Wave Boundary-Layer Interactions at Mach6

Shock-wave boundary-layer interactions (SWBLI) pose unique challenges for the slender airframes and sharp control surfaces characteristic of high-speed airbreathing aircraft. Spikes in pressure and thermal loading have been shown to occur in SWBLI regions; thus, predicting their location and behavior proves critical for assuring aircraft structural integrity and fatigue resilience. SWBLIs frequently lead to severe turbulent boundary-layer separation and corresponding unsteady motions characterized by a wide range of frequencies. Low-frequency motion has been shown to couple with the structure of a compliant panel, and turbulent boundary layers thus present a particular danger for thin control surfaces. Previous studies have shown that the thermal loading on a compliant panel in a ramp-induced SWBLI may have a signi cant impact on the panel vibrations. A multitude of studies has investigated the heating on rigid structures, but a charactrization of the heating effects on compliant structures is necessary to fully understand the fluid-thermal-structural interactions (FTSI) occurring in ramp-induced SWBLIs. This studuses infrared thermography and high-frequency Kulite pressure transducers to examine the heating loads on rigid and compliant ramps at various angles and flow conditions at Mach 6.

Hypersonic↗

Fluid-Thermal-Structural Interactions in Ramp-Induced Shock-Wave Boundary-Layer Interactions at Mach 6

Shock-wave boundary-layer interactions (SWBLI) pose unique challenges for the slender airframes and sharp control surfaces characteristic of high-speed airbreathing aircraft. Spikes in pressure and thermal loading have been shown to occur in SWBLI regions; thus, predicting their location and behavior proves critical for assuring aircraft structural integrity and fatigue resilience. SWBLIs frequently lead to severe turbulent boundary-layer separation and corresponding unsteady motions characterized by a wide range of frequencies. Low-frequency motion has been shown to couple with the structure of a compliant panel, and turbulent boundary layers thus present a particular danger for thin control surfaces. Previous studies have shown that the thermal loading on a compliant panel in a ramp-induced SWBLI may have a signi cant impact on the panel vibrations. A multitude of studies has investigated the heating on rigid structures, but a charactrization of the heating effects on compliant structures is necessary to fully understand the fluid-thermal-structural interactions (FTSI) occurring in ramp-induced SWBLIs. This studuses infrared thermography and high-frequency Kulite pressure transducers to examine the heating loads on rigid and compliant ramps at various angles and flow conditions at Mach 6.

Hypersonic↗

KIVA-Update June 2012

Development of a fractional step, a Predictor-Corrector Split (PCS), or what is often known as a projection method combining hp-adaptive system in a Finite Element Method (FEM) for combustion modeling has been achieved. This model will advance the accuracy and range of applicability of the KIVA combustion model and software used typically for internal combustion engine modeling. This abstract describes a PCS hp-adaptive FEM model for turbulent reactive flow spanning all velocity regimes and fluids that is being developed for the new KIVA combustion algorithm, particularly for internal combustion engines. The method and general solver is applicable to Newtonian and non- Newtonian fluids and also for incompressible solids, porous media, solidification modeling, and fluid structure interaction problems. The fuel injection and injector modeling could easily benefit from the capability of solving the fluid structure interaction problem in an injector, helping to understand cycle to cycle variation and cavitation. This is just one example where the new algorithm differs from the old, in addition to handling Conjugate Heat Transfer (CHT), although there a numerous features that makes the new system more robust and accurate. In these ways, the PCS hp-adaptive algorithm does not compete with commercial software packages, those often used in conjunction with the currently distributed KIVA codes for engine combustion modeling. In addition, choosing a local ALE method on immersed moving parts represented by overset grid that is 2nd order spatially accurate, allows for easy grid generation from CAD to fluid grid while also provide for robustness in handling any possible moving parts configuration without any code modifications. The combined methods employed produce a minimal amount of computational effort as compared to fully resolved grids at the same accuracy. We demonstrate the solver on benchmark problems for the all flow regimes as follows: (1) 2-D backward-facing step using h-adaption, (2) 2-D driven cavity, (3) 2-D natural convection in a differentially heat cavity with h-adaptation, (4) NACA 0012 airfoil in 2-D, (5) supersonic flows over compression ramps, (6) 2-D natural convection in a differentially heat cavity with hp-adaptation, (7) 3-D natural convection in a differentially heat sphere with hp-adaptation. In addition, we show the new moving parts algorithm for working for a 2-D piston; the immersed moving parts method also for valves and pistons, vanes, etc... The movement is performed using an overset grid method and is 2nd order accurate in space, and never produces a tangle grid, that is, robust system at any resolution and any parts configuration. We also show CHT for the currently distributed KIVA-4mpi software and some fairly automatic grid generation using Sandia's Cubit unstructured grid generator. A new electronic web-based manual for KIVA-4 has been developed as well.

02 PETROLEUM↗

Structural Loads Analysis for Wave Energy Converters: Preprint

This study explores and verifies the generalized body-modes method for evaluating the structural loads on a wave energy converter (WEC). Historically, WEC design methodologies have focused primarily on accurately evaluating hydrodynamic loads, while methodologies for evaluating structural loads have yet to be fully considered and incorporated into the WEC design process. As wave energy technologies continue to advance, however, it has become increasingly evident that an accurate evaluation of the structural loads will enable an optimized structural design, as well as the potential utilization of composites and flexible materials, and hence reduce WEC costs. Although there are many computational fluid dynamics, structural analyses and fluid-structure-interaction (FSI) codes available, the application of these codes is typically too computationally intensive to be practical in the early stages of the WEC design process. The generalized body-modes method, however, is a reduced order, linearized, frequency-domain FSI approach, performed in conjunction with the linear hydrodynamic analysis, with computation times that could realistically be incorporated into the WEC design process. The objective of this study is to verify the generalized body-modes approach in comparison to high-fidelity FSI simulations to accurately predict structural deflections and stress loads in a WEC. Two verification cases are considered, a free-floating barge and a fixed-bottom column. Details for both the generalized body-modes models and FSI models are first provided. Results for each of the models are then compared and discussed. Finally, based on the verification results obtained, future plans for incorporating the generalized body-modes method into the WEC simulation tool, WEC-Sim, and the overall WEC design process are discussed.

17 WIND ENERGY↗

On the continuity of mean total normal stress in geometrical multiscale cardiovascular problems

In this work an iterative strategy to implicitly couple dimensionally-heterogeneous blood flow models accounting for the continuity of mean total normal stress at interface boundaries is developed. Conservation of mean total normal stress in the coupling of heterogeneous models is mandatory to satisfy energetic consistency between them. Nevertheless, existing methodologies are based on modifications of the Navier–Stokes variational formulation, which are undesired when dealing with fluid–structure interaction or black box codes. The proposed methodology makes possible to couple one-dimensional and three-dimensional fluid–structure interaction models, enforcing the continuity of mean total normal stress while just imposing flow rate data or even the classical Neumann boundary data to the models. This is accomplished by modifying an existing iterative algorithm, which is also able to account for the continuity of the vessel area, when required. Comparisons are performed to assess differences in the convergence properties of the algorithms when considering the continuity of mean normal stress and the continuity of mean total normal stress for a wide range of flow regimes. Finally, examples in the physiological regime are shown to evaluate the importance, or not, of considering the continuity of mean total normal stress in hemodynamics simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Steady-State Solution of a Flexible Wing

A fluid-structure interaction code, ENSAERO, has been used to compute the aerodynamic loads on a swept-tapered wing. The code has the capability of using Euler or Navier-Stokes equations. Both options have been used and compared in the present paper. In the calculation of the steady-state solution, we are interested in knowing how the flexibility of the wing influences the lift coefficients. If the results of a flexible wing are not affected by the flexibility of the wing significantly, one could consider the wing to be rigid and reduce the problem from fluid-structure interaction to a fluid problem.

Karkehabadi, Reza↗

Reduction of vortex induced forces and motion through surface roughness control

Roughness is added to the surface of a bluff body in a relative motion with respect to a fluid. The amount, size, and distribution of roughness on the body surface is controlled passively or actively to modify the flow around the body and subsequently the Vortex Induced Forces and Motion (VIFM). The added roughness, when designed and implemented appropriately, affects in a predetermined way the boundary layer, the separation of the boundary layer, the level of turbulence, the wake, the drag and lift forces, and consequently the Vortex Induced Motion (VIM), and the fluid-structure interaction. The goal of surface roughness control is to decrease/suppress Vortex Induced Forces and Motion. Suppression is required when fluid-structure interaction becomes destructive as in VIM of flexible cylinders or rigid cylinders on elastic support, such as underwater pipelines, marine risers, tubes in heat exchangers, nuclear fuel rods, cooling towers, SPAR offshore platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Finite-Volume method for compressible non-equilibrium two-phase flows in networks of elastic pipelines using the Baer–Nunziato model

A novel Finite-Volume scheme for the numerical computations of compressible two-phase flows in pipelines is proposed for the fully non-equilibrium Baer–Nunziato model. The present FV approach is the extension of the method proposed in Daude and Galon (2018) in the context of the Euler equations to the Baer–Nunziato model. In addition, proper approximations of the non-conservative terms are proposed to consider jumps of volume fraction as well as jumps of cross-section in order to respect uniform pressure and velocity profiles preservation. In particular, focus is given to the numerical treatment of abrupt changes in area and to networks wherein several pipelines are connected at junctions. The proposed method makes it possible to avoid the use of an iterative procedure for the solution of the junction problem. The present approach can also deal with general Equations Of State. In addition, the fluid–structure interaction of compressible fluid flowing in flexible pipes is also considered. Lastly, the proposed scheme is then assessed on a variety of shock-tubes and other transient flow problems and experiments demonstrating its capability to resolve such problems efficiently, accurately and robustly.

42 ENGINEERING↗

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗

Miniature Flow and Viscosity Sensors

Integration of flow sensing into microfluidic systems such as automated sample preparation and processing has seen limited applications due to the large cost and size of commercially available flow sensors. Nevertheless, flow sensing is a critical missing components for many high precision microfluidic applications. Los Alamos National Laboratory researchers have developed a miniature flow sensor technology based on a novel microfluidic method developed at the Laboratory. The method takes advantage of the fluid-structure interaction between the fluid and an elastic membrane to provide precision measurement. Depending on how it is integrated, it is also possible to use the technology to measure viscosity of the fluid. The flow sensor can work with a wide range of fluids.

47 OTHER INSTRUMENTATION↗

New developments in adaptive methods for computational fluid dynamics

New developments in a posteriori error estimates, smart algorithms, and h- and h-p adaptive finite element methods are discussed in the context of two- and three-dimensional compressible and incompressible flow simulations. Applications to rotor-stator interaction, rotorcraft aerodynamics, shock and viscous boundary layer interaction and fluid-structure interaction problems are discussed.

Oden, J. T.↗