Search NASA⌕ Search

SEARCH · Search NASA

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Surface Instability of Liquid Propellant under Vertical Oscillatory Forcing

Fluid motion in a fuel tank produced during thrust oscillations can circulate sub-cooled hydrogen near the liquid-vapor interface resulting in increased condensation and ullage pressure collapse. The first objective of this study is to validate the capabilities of a Computational Fluid Dynamics (CFD) tool, CFD-ACE+, in modeling the fundamental interface transition physics occurring at the propellant surface. The second objective is to use the tool to assess the effects of thrust oscillations on surface dynamics. Our technical approach is to first verify the CFD code against known theoretical solutions, and then validate against existing experiments for small scale tanks and a range of transition regimes. A 2D axisymmetric, multi-phase model of gases, liquids, and solids is used to verify that CFD-ACE+ is capable of modeling fluid-structure interaction and system resonance in a typical thrust oscillation environment. Then, the 3D mode is studied with an assumed oscillatory body force to simulate the thrust oscillating effect. The study showed that CFD modeling can capture all of the transition physics from solid body motion to standing surface wave and to droplet ejection from liquid-gas interface. Unlike the analytical solutions established during the 1960 s, CFD modeling is not limited to the small amplitude regime. It can extend solutions to the nonlinear regime to determine the amplitude of surface waves after the onset of instability. The present simulation also demonstrated consistent trends from numerical experiments through variation of physical properties from low viscous fluid to high viscous fluids, and through variation of geometry and input forcing functions. A comparison of surface wave patterns under various forcing frequencies and amplitudes showed good agreement with experimental observations. It is concluded that thrust oscillations can cause droplet formation at the interface, which results in increased surface area and enhanced heat transfer between the liquid and gas phases as the ejected droplets travel well into the warmer gas region.

Yang, H. Q.↗
NASA NTRS1993

Phase 2: HGM air flow tests in support of HEX vane investigation

Following the start of SSME certification testing for the Pratt and Whitney Alternate Turbopump Development (ATD) High Pressure Oxidizer Turbopump (HPOTP), cracking of the leading edge of the inner HEX vane was experienced. The HEX vane, at the inlet of the oxidizer bowl in the Hot Gas Manifold (HGM), accepts the HPOTP turbine discharge flow and turns it toward the Gaseous Oxidizer Heat Exchanger (GOX HEX) coil. The cracking consistently initiated over a specific circumferential region of the hex vane, with other circumferential locations appearing with increased run time. Since cracking had not to date been seen with the baseline HPOTP, a fluid-structural interaction involving the ATD HPOTP turbine exit flowfield and the HEX inner vane was suspected. As part of NASA contract NAS8-36801, Pratt and Whitney conducted air flow tests of the ATD HPOTP turbine turnaround duct flowpath in the MSFC Phase 2 HGM air flow model. These tests included HEX vane strain gages and additional fluctuating pressure gages in the turnaround duct and HEX vane flowpath area. Three-dimensional flow probe measurements at two stations downstream of the turbine simulator exit plane were also made. Modifications to the HPOTP turbine simulator investigated the effects on turbine exit flow profile and velocity components, with the objective of reproducing flow conditions calculated for the actual ATD HPOTP hardware. Testing was done at the MSFC SSME Dynamic Fluid Air Flow (Dual-Leg) Facility, at air supply pressures between 50 and 250 psia. Combinations of turbine exit Mach number and pressure level were run to investigate the effect of flow regime. Information presented includes: (1) Descriptions of turbine simulator modifications to produce the desired flow environment; (2) Types and locations for instrumentation added to the flow model for improved diagnostic capability; (3) Evaluation of the effect of changes to the turbine simulator flowpath on the turbine exit flow environment; and (4) Comparison of the experimental turbine exit flow environment to the environment calculated for the ATD HPOTP.

Cox, G. B., Jr.↗
DOE OSTI2020

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗
NASA NTRS2014

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise about Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the airframe and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large flat plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position.The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non-dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A. E.↗
NASA NTRS1997

On spurious behavior of CFD simulations

Spurious behavior in underresolved grids and/or semi-implicit temporal discretizations for four computational fluid dynamics (CFD) simulations are studied. The numerical simulations consist of (a) a 1-D chemically relaxed nonequilibrium model, (b) the direct numerical simulation (DNS) of 2-D incompressible flow over a backward facing step, (c) a loosely-coupled approach for a 2-D fluid-structure interaction, and (d) a 3-D compressible unsteady flow simulation of vortex breakdown in delta wings. Using knowledge from dynamical systems theory, various types of spurious behaviors that are numerical artifacts were systematically identified. These studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by the available computing power. In large scale computations underresolved grids, semi-implicit procedures, loosely-coupled implicit procedures, and insufficiently long time integration in DNS are most often unavoidable. Consequently, care must be taken in both computation and in interpretation of the numerical data. The results presented confirm the important role that dynamical systems theory can play in the understanding of the nonlinear behavior of numerical algorithms and in aiding the identification of the sources of numerical uncertainties in CFD.

Numerical Solution↗
NASA NTRS2014

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise About Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the air-frame and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large at plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position. The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non- dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A.↗
DOE OSTI2020

Assessment of Survey Results from Advanced Reactor Industry Domain Experts on Nonlinear Soil-Structure Interaction Analysis Software Verification and Validation

The seismic load case has a significant impact on the design and cost of nuclear power plants. Given the need to substantially reduce cost, advanced reactor designers are looking to leverage numerical tools that enable seismic analysis of an integrated vessel/support/structure/soil system. Modern nonlinear analysis tools provide a solution, capturing dynamic coupling between components (including soil-structure interaction (SSI)) concurrent with nonlinear behavior in one or more parts of the system. Although these methods have a rich history of technical development and implementation, software quality assurance (SQA) following nuclear industry standards remains a significant burden for those wishing to adopt such methods for advanced reactor design and licensing. To reduce the SQA burden, this project will develop the guidance for SQA verification and validation (V&V) of coupled nonlinear SSI analysis tools, to include a test matrix of software features and test problems to support commercial grade dedication (CGD). The guidance is intended to be technology-neutral: supporting designers of all advanced reactors. To maximize the value of the project to reactor designers, the project team performed focused surveys of and interviews with advanced reactor designers. Additionally, the project surveyed members of the broader industry involved in reactor design and licensing. The aggregated industry feedback consists of written survey responses, live polling responses, focused interviews, and informal feedback provided after outreach presentations, collectively referred to as “survey results”. Herein, the survey results are reported, reviewed, and assessed to inform follow-on project activities. The survey results confirmed the familiarity of the industry with the coupled nonlinear SSI analysis. They also affirm the project need based upon the expressed intent to (1) incorporate nonlinear features and (2) pursue the coupled nonlinear SSI analysis to reduce seismic demands and construction costs. All reactor designers indicated their intent to incorporate two or more nonlinear features, and all reactor designer respondents opined that the coupled nonlinear SSI analysis would allow for the optimization of the reactor design and construction. However, the programmatic challenges, whether real or perceived, present a significant barrier to reactor designers. The barrier most commonly identified by the reactor designer population was regulatory risk, with 70% citing this as a reason not to pursue the approach. The perceived regulatory risk identified by the reactor designers underscores the importance of regulator engagement and dialogue in this project. Additionally, half of the reactor designers identified cost and lack of guidance as a deterrent. The survey results also provide insights to tailor specific aspects of the guidance document and test matrix. The project plans to prepare both a formal referenceable guidance document and a collaborative, web-based test matrix and problem set. The project will focus on more complete test problem definitions for the prioritized nonlinear features over shallower problem descriptions for a larger set of features. The nonlinearities prioritized as high based upon survey feedback include fluid-structure interaction and seismic isolation and energy dissipation devices. The nonlinearities prioritized as intermediate include nonlinear geomaterials, nonlinear concrete, nonlinear steel, and interface/contact nonlinearity. Deep embedment and the associated nonlinear phenomena are assigned the lowest priority based upon survey feedback.

42 ENGINEERING↗
DOE OSTI2025

High fidelity blade-resolved and actuator line data from a 16 turbine wind farm simulation using ExaWind

This data was generated with the ExaWind code suite (https://github.com/Exawind) as a demonstration of a large, 16 turbine wind farm simulation, calculated using two different levels of fidelity. The lower level of fidelity approach uses an actuator line approach to represent the turbines, and was simulated with AMR-Wind (https://github.com/Exawind/amr-wind/) as the background flow solver, coupled to OpenFAST (https://github.com/OpenFAST/openfast). The higher level of fidelity simulation uses a blade-resolved approach, and is done using AMR-Wind, Nalu-Wind (https://github.com/Exawind/nalu-wind), OpenFAST, and TIOGA (https://github.com/Exawind/tioga). In the blade-resolved simulation, ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. OpenFAST handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. In the actuator line simulation, a mesh of 295M elements was used for a 5km x 5km domain, and it was simulated using 256 nodes (2048 GPU's) on the Oak Ridge Leadership Computing Facility Frontier supercomputer. For the blade-resolved simulation, 1.5B element mesh was used in the AMR-Wind background 5km x 5km domain, and 16M elements were used for each turbine in the Nalu-Wind domains, for a total of 1.7B elements. This was simulated using 384 nodes on Frontier, with each node using 56 cores for Nalu-Wind and 8 GPU cores. The data in this archive includes the turbine outputs from OpenFAST, 2D sampling planes from AMR-Wind, and full-field solution files from AMR-Wind and Nalu-Wind.

17 WIND ENERGY↗
DOE OSTI2025

High temporal frequency data from a four turbine, blade-resolved wind farm simulation with ExaWind

The data was generated with ExaWind (https://github.com/Exawind) which couples AMR-Wind (https://github.com/Exawind/amr-wind/), Nalu-Wind (https://github.com/Exawind/nalu-wind), TIOGA (https://github.com/Exawind/tioga), and OpenFAST (https://github.com/OpenFAST/openfast). This is a large-scale simulation of a blade-resolved wind farm using the ExaWind software stack. ExaWind couples together a background flow solver, AMR-Wind, and a near-body solver, Nalu-Wind, through an overset technique from the TIOGA application. Another application, OpenFAST, handles the structural dynamics of the turbine blades and towers, which informs the fluid-structure interaction of the wind turbines with the flow solvers. This particular simulation includes four blade-resolved wind turbines operating in a turbulent atmospheric boundary layer. The AMR-Wind solver uses 500 million cells and is being solved on 256 AMD GPUs of the Oakridge Leadership Computing Facility Frontier supercomputer. Each turbine is assigned its own Nalu-Wind solver with over 13 million elements per turbine and solved using 448 CPU cores, for a total of 1792 CPU cores. For each node, 56 cores contain Nalu-Wind, while 8 cores correspond to AMR-Wind operations on the GPUs. Consequently, ExaWind is entirely utilizing the CPUs and the GPUs of the nodes concurrently. The data used in the visualization is full flow field data output from the simulation. It is lossy-compressed to a specific accuracy using ZFP and written to disk every 16 time-steps to enable real-time flow visualization. The flow fields are sampled at a high temporal frequency to enable real-time, 24fps visualization. The flow fields are sampled every 12 simulation time steps (every 0.04132s).

17 WIND ENERGY↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2020

High-order partitioned spectral deferred correction solvers for multiphysics problems

We present an arbitrarily high-order, conditionally stable, partitioned spectral deferred correction (SDC) method for solving multiphysics problems using a sequence of pre-existing single-physics solvers. This method extends the work in [1], [2], which used implicit-explicit Runge-Kutta methods (IMEX) to build high-order, partitioned multiphysics solvers. We consider a generic multiphysics problem modeled as a system of coupled ordinary differential equations (ODEs), coupled through coupling terms that can depend on the state of each subsystem; therefore the method applies to both a semi-discretized system of partial differential equations (PDEs) or problems naturally modeled as coupled systems of ODEs. The sufficient conditions to build arbitrarily high-order partitioned SDC schemes are derived. Based on these conditions, various of partitioned SDC schemes are designed. The stability of the first-order partitioned SDC scheme is analyzed in detail on a coupled, linear model problem. We show that the scheme is conditionally stable, and under conditions on the coupling strength, the scheme can be unconditionally stable. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems with moderate coupling strength. They include a stiff linear system of ODEs, advection-diffusion-reaction systems, and fluid-structure interaction problems with both incompressible and compressible flows, where we verify the design order of the SDC schemes and study various stability properties. We also directly compare the accuracy, stability, and cost of the proposed partitioned SDC solver with the partitioned IMEX method in [1], [2] on this suite of test problems. The results suggest that the high-order partitioned SDC solvers are more robust than the partitioned IMEX solvers for the numerical examples considered in this work, while the IMEX methods require fewer implicit solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Gas-driven short disconnection mitigates thermal runaway in Li-ion batteries under mechanical abuse

Mechanical abuse poses a critical safety risk to Li-ion batteries by inducing internal short circuits that initiate thermal runaway. Remarkably, voltage recovery frequently emerges during thermal runaway initiation, a phenomenon that conventional theories fail to explain. Here, our study develops a new multiphysics mechanism to explain voltage recovery, termed gas-driven short disconnection, whereby internal gas pressure mechanically disengages short-circuit contacts and causes the cell voltage to rebound. This mechanism incorporates gas generation and its structural impact on the short circuit. Real-time optical and thermal imaging and X-ray computed tomography reveal fluid-structure interaction between internal gas flow and adjacent shorting contacts. We establish a mechanistic framework linking gas-driven short disconnection to cell-level voltage and temperature responses, elucidating the extension–truncation pattern of voltage recovery. Furthermore, thermal regime maps show that a voltage recovery duration exceeding 5 s correlates with limited temperature rise below 150 °C, indicating that sustained short-circuit disconnection suppresses Joule heating. Additionally, a dimensionless criterion is deduced from scaling analysis for physical plausibility of gas-driven short disconnection in mechanically abused cells. This finding inspires smart venting control, which regulates gas release to maintain the internal pressure while dissipating gas enthalpy, thereby providing a device-level strategy for thermal runaway mitigation.

Gas-driven short disconnection↗
DOE OSTI2021

Combined selection of the dynamic model and modeling error in nonlinear aeroelastic systems using Bayesian Inference

Here, we report a Bayesian framework for concurrent selection of physics-based models and (modeling) error models. We investigate the use of colored noise to capture the mismatch between the predictions of calibrated models and observational data that cannot be explained by measurement error alone within the context of Bayesian estimation for stochastic ordinary differential equations. Proposed models are characterized by the average data-fit, a measure of how well a model fits the measurements, and the model complexity measured using the Kullback–Leibler divergence. The use of a more complex error models increases the average data-fit but also increases the complexity of the combined model, possibly over-fitting the data. Bayesian model selection is used to find the optimal physical model as well as the optimal error model. The optimal model is defined using the evidence, where the average data-fit is balanced by the complexity of the model. The effect of colored noise process is illustrated using a nonlinear aeroelastic oscillator representing a rigid NACA0012 airfoil undergoing limit cycle oscillations due to complex fluid–structure interactions. Several quasi-steady and unsteady aerodynamic models are proposed with colored noise or white noise for the model error. The use of colored noise improves the predictive capabilities of simpler models.

42 ENGINEERING↗
DOE OSTI2026

Margin to onset of nucleate boiling and flow instability studies for preliminary MITR design-demonstration element thermal-hydraulics

This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗
DOE OSTI2022

Bio-inspired bistable piezoelectric energy harvester for powering animal telemetry tags: Conceptual design and preliminary experimental validation

Here, this paper presents the conceptual design, preliminary experimental validation, and performance evaluation of a novel bio-inspired bi-stable piezoelectric energy harvester for self-powered animal telemetry tags. The overall conceptual design, which includes a bio-inspired attachment and a bi-stable piezoelectric energy harvester, is introduced firstly with a specific application example of marine fish tracking. The self-powered telemetry tag can be externally deployed on fish (dorsal fin) to monitor fish habitats, population, and underwater environment. Inspired by the Venus flytrap's rapid shape transition, a bi-stable piezoelectric energy harvester is developed to scavenge energy from fish maneuvering and the surrounding fluid flow for a sustainable power supply. The bistability of the harvester is characterized by the measured force-displacement curve and double potential wells. A bluff body is integrated to the free end of the bistable piezoelectric energy harvester to enhance the structure-fluid interaction for the large-amplitude snap-through vibrations and higher voltage output. Controlled laboratory experiments are conducted in a water tank on the bio-inspired bi-stable piezoelectric energy harvester using a servo motor system to simulate fish swing motion at various conditions to evaluate the power generation performance. The preliminary underwater experimental results demonstrated that the proposed bio-inspired bi-stable piezoelectric energy harvester could effectively convert fish swing motions into electricity. The device collected 17.25 mJ of energy over 130 s under a peak-to-peak swing angle of 30° at 1.5 Hz in the capacitor charging experiments.

Animal tracking↗
DOE OSTI2020

CU-BENs: A structural modeling finite element library

The present work discusses capabilities within the finite element library CU-BENs. CU-BENs focuses on applying the finite element method to structural mechanics problems encountered within the context of inverse problems and partitioned fluid–structure interaction; thus element formulations are primarily of a structural type — truss, frame, and triangular discrete Kirchhoff theory shells. CU-BENs defaults to the skyline sparse storage scheme for the system matrix, but also supports other storage schemes when using external libraries such as LAPACK and UMFPACK. CU-BENs includes built-in nonlinear solution strategies, such as the Newton Raphson method and modified spherical arc length method, that are available within static analyses as well within the context of transient dynamic analyses involving a generalized-α implementation of the Newmark implicit time integration scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Modeling of the spray-induced wall stress acting on the ignition assistance device

This research introduces a novel wall-stress model called the Spray-Induced Wall Stress (SIWS) model, which considers the effects of spray-wall impingement and the resulting formation of wall stress within the Lagrangian spray modeling framework. The primary objective of this paper is to provide a mathematical description of the fundamental physics underlying the model. Subsequently, the proposed model is validated using existing experimental data. Here, the remainder of the study focuses on the practical application of the model to an ignition assistance device. Specifically, this device is installed in a compression ignition engine and designed to enhance ignition in aviation-fueled high-altitude aircraft propulsion systems. The research sheds light on the mechanical impulse caused by the high-speed impact of the spray jet, leading to the accumulation of mechanical stress on the rigid body of the ignition assistance device. Previous studies on fluid-structure interaction have only considered the interaction between the gas phase and the solid wall. However, the SIWS model incorporates the additional impact of the impinging liquid spray jet. Consequently, the simulated stress distribution on the ignition assistance device can be estimated by considering both the gas-phase-induced term and the spray-induced term simultaneously.

33 ADVANCED PROPULSION SYSTEMS↗