SEARCH · Engineering Papers

Results for “equation of fluid dynamics”

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 1,297 records · Page 72

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Thermal Fault Tolerance Analysis of Carbon Fiber Rope Barrier Systems for Use in the Reusable Solid Rocket Motor ( RSRM) Nozzle Joints

Carbon Fiber Rope (CFR) thermal barrier systems are being considered for use in several RSRM (Reusable Solid Rocket Motor) nozzle joints as a replacement for the current assembly gap close-out process/design. This study provides for development and test verification of analysis methods used for flow-thermal modeling of a CFR thermal barrier subject to fault conditions such as rope combustion gas blow-by and CFR splice failure. Global model development is based on a 1-D (one dimensional) transient volume filling approach where the flow conditions are calculated as a function of internal 'pipe' and porous media 'Darcy' flow correlations. Combustion gas flow rates are calculated for the CFR on a per-linear inch basis and solved simultaneously with a detailed thermal-gas dynamic model of a local region of gas blow by (or splice fault). Effects of gas compressibility, friction and heat transfer are accounted for the model. Computational Fluid Dynamic (CFD) solutions of the fault regions are used to characterize the local flow field, quantify the amount of free jet spreading and assist in the determination of impingement film coefficients on the nozzle housings. Gas to wall heat transfer is simulated by a large thermal finite element grid of the local structure. The employed numerical technique loosely couples the FE (Finite Element) solution with the gas dynamics solution of the faulted region. All free constants that appear in the governing equations are calibrated by hot fire sub-scale test. The calibrated model is used to make flight predictions using motor aft end environments and timelines. Model results indicate that CFR barrier systems provide a near 'vented joint' style of pressurization. Hypothetical fault conditions considered in this study (blow by, splice defect) are relatively benign in terms of overall heating to nozzle metal housing structures.

Clayton, J. Louie↗

NASA NTRS2010

Parameterized Linear Longitudinal Airship Model

A parameterized linear mathematical model of the longitudinal dynamics of an airship is undergoing development. This model is intended to be used in designing control systems for future airships that would operate in the atmospheres of Earth and remote planets. Heretofore, the development of linearized models of the longitudinal dynamics of airships has been costly in that it has been necessary to perform extensive flight testing and to use system-identification techniques to construct models that fit the flight-test data. The present model is a generic one that can be relatively easily specialized to approximate the dynamics of specific airships at specific operating points, without need for further system identification, and with significantly less flight testing. The approach taken in the present development is to merge the linearized dynamical equations of an airship with techniques for estimation of aircraft stability derivatives, and to thereby make it possible to construct a linearized dynamical model of the longitudinal dynamics of a specific airship from geometric and aerodynamic data pertaining to that airship. (It is also planned to develop a model of the lateral dynamics by use of the same methods.) All of the aerodynamic data needed to construct the model of a specific airship can be obtained from wind-tunnel testing and computational fluid dynamics

Kulczycki, Eric↗

DOE OSTI2023

Nonlinear dynamics of small-scale Alfvén waves

Herein we study the nonlinear evolution of very oblique small-scale Alfvén waves with k ⊥ d i ≳1. At these scales, the waves become significantly compressive, unlike in magnetohydrodynamics, due to the Hall term in the equations. We demonstrate that when frequencies are small compared to the ion gyrofrequency and amplitudes are small compared to unity, no new nonlinear interaction appears due to the Hall term alone at the lowest non-trivial order, even when k ⊥ d i ~1. However, at the second non-trivial order, we discover that the Hall physics leads to a slow but resonant nonlinear interaction between co-propagating Alfvén waves, an inherently three-dimensional effect. Including the effects of finite temperature, finite frequency, and electron inertia, the two-fluid Alfvén wave also becomes dispersive once one or more of k ⊥ ρ s , k ⊥ d e , or k ∥ d i becomes significant: for oblique waves at low β as studied here, this can be at a much smaller scale than d i . We show that the timescale for one-dimensional steepening of two-fluid Alfvén waves is only significant at these smaller dispersive scales, and also derive an expression for the amplitude of driven harmonics of a primary wave. Importantly, both new effects are absent in gyrokinetics and other commonly used reduced two-fluid models. Our calculations have relevance for the interpretation of laboratory Alfvén wave experiments, as well as shedding light on the physics of turbulence in the solar corona and inner solar wind, where the dominant nonlinear interaction between counter-propagating waves is suppressed, allowing these new effects to become important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗

NASA NTRS1993

Dynamic ocean-tide effects on Earth's rotation

This article develops 'broad-band' Liouville equations which are capable of determining the effects on the rotation of the Earth of a periodic excitation even at frequencies as high as semi-diurnal; these equations are then used to predict the rotational effects of altimetric, numerical and 32-constituent spherical harmonic ocean-tide models. The rotational model includes a frequency-dependent decoupled core, the effects of which are especially marked near retrograde diurnal frequencies; and a fully dynamic oceanic response, whose effects appear to be minor despite significant frequency dependence. The model also includes solid-earth effects which are frequency dependent as the result of both anelasticity at long periods and the fluid-core resonance at nearly diurnal periods. The effects of both tidal inertia and relative angular momentum on Earth rotation (polar motion, length of day, 'nutation' and Universal Time) are presented for 32 long- and short-period ocean tides determined as solutions to the author's spherical harmonic tide theory. The lengthening of the Chandler wobble period by the pole tide is also re-computed using the author's full theory. Additionally, using the spherical harmonic theory, tidal currents and their effects on rotation are determined for available numerical and altimetric tide height models. For all models, we find that the effects of tidal currents are at least as important as those of tide height for diurnal and semi-diurnal constituents.

Dickman, S. R.↗

NASA NTRS1994

Design of a variable area diffuser for a 15-inch Mach 6 open-jet tunnel

The Langley 15-inch Mach 6 High Temperature Tunnel was recently converted from a Mach 10 Hypersonic Flow Apparatus. This conversion was effected to improve the capability of testing in Mach 6 air at relatively high reservoir temperatures not previously possible at Langley. Elevated temperatures allow the matching of the Mach numbers, Reynolds numbers, and ratio of wall-to-adiabatic-wall temperatures (TW/Taw) between this and the Langley 20-inch Mach 6 CF4 Tunnel. This ratio is also matched for Langley's 31-inch Mach 10 Tunnel and is an important parameter useful in the simulation of slender bodies such as National Aerospace Plane (NASP) configurations currently being studied. Having established the nozzle's operating characteristics, the decision was made to install another test section to provide model injection capability. This test section is an open-jet type, with an injection system capable of injecting a model from retracted position to nozzle centerline between 0.5 and 2 seconds. Preliminary calibrations with the new test section resulted in Tunnel blockage. This blockage phenomenon was eliminated when the conical center body in the diffuser was replaced. The issue then, is to provide a new and more efficient variable area diffuser configuration with the capability to withstand testing of larger models without sending the Tunnel into an unstart condition. Use of the 1-dimensional steady flow equation with due regard to friction and heat transfer was employed to estimate the required area ratios (exit area / throat area) in a variable area diffuser. Correlations between diffuser exit Mach number and area ratios, relative to the stagnation pressure ratios and diffuser inlet Mach number were derived. From these correlations, one can set upper and lower operating pressures and temperatures for a given diffuser throat area. In addition, they will provide appropriate input conditions for the full 3-dimensional computational fluid dynamics (CFD) code for further simulation studies.

Loney, Norman W.↗

NASA NTRS1996

A Burst Mode, Ultrahigh Temperature UF4 Vapor Core Reactor Rankine Cycle Space Power System Concept

Static and dynamic neutronic analyses have been performed on an innovative burst mode (100's of MW output for a few thousand seconds) Ulvahigh Temperature Vapor Core Reactor (UTVR) space nuclear power system. The NVTR employs multiple, neutronically-coupled fissioning cores and operates on a direct, closed Rankine cycle using a disk Magnetohydrodynamic (MHD) generater for energy conversion. The UTVR includes two types of fissioning core regions: (1) the central Ultrahigh Temperature Vapor Core (UTVC) which contains a vapor mixture of highly enriched UF4 fuel and a metal fluoride working fluid and (2) the UF4 boiler column cores located in the BeO moderator/reflector region. The gaseous nature of the fuel the fact that the fuel is circulating, the multiple coupled fissioning cores, and the use of a two phase fissioning fuel lead to unique static and dynamic neutronic characteristics. Static neutronic analysis was conducted using two-dimensional S sub n, transport theory calculations and three-dimensional Monte Carlo transport theory calculations. Circulating-fuel, coupled-core point reactor kinetics equations were used for analyzing the dynamic behavior of the UTVR. In addition to including reactivity feedback phenomena associated with the individual fissioning cores, the effects of core-to-core neutronic and mass flow coupling between the UTVC and the surrounding boiler cores were also included in the dynamic model The dynamic analysis of the UTVR reveals the existence of some very effectlve inherent reactivity feedback effects that are capable of quickly stabilizing this system, within a few seconds, even when large positive reactivity insertions are imposed. If the UTVC vapor fuel density feedback is suppressed, the UTVR is still inherently stable because of the boiler core liquid-fuel volume feedback; in contrast, suppression of the vapor fuel density feedback in 'conventional" gas core cavity reactors causes them to become inherently unstable. Due to the strength of the negative reactivity feedback in the UTVR, it is found that external reactivity insertions alone are inadequate for bringing about significant power level changes during normal reactor operations. Additional methods of reactivity control such as variations in the gaseous fuel mass flow rate, are needed to achieve the desired power level oontrol.

Dugan, E. T.↗

NASA NTRS2004

A Hybrid Approach To Tandem Cylinder Noise

Aeolian tone generation from tandem cylinders is predicted using a hybrid approach. A standard computational fluid dynamics (CFD) code is used to compute the unsteady flow around the cylinders, and the acoustics are calculated using the acoustic analogy. The CFD code is nominally second order in space and time and includes several turbulence models, but the SST k - omega model is used for most of the calculations. Significant variation is observed between laminar and turbulent cases, and with changes in the turbulence model. A two-dimensional implementation of the Ffowcs Williams-Hawkings (FW-H) equation is used to predict the far-field noise.

Lockard, David P.↗

DOE OSTI2020

Local Truncation Error-Informed Code Verification

The method of manufactured solutions (MMS) has become increasingly popular in conducting code verification studies on predictive codes, such as nuclear power system codes and computational fluid dynamic codes. The reason for the popularity of this approach is that it can be used when an analytical solution is not available. Using MMS, code developers are able to verify that their code is free of coding errors that impact the observed order of accuracy. While MMS is still an excellent tool for code verification, it does not identify coding errors that are of the same order as the numerical method. This paper presents a method that combines MMS with modified equation analysis (MEA), which calculates the local truncation error (LTE) to identify coding error up to and including the order of the numerical method. This method is referred to as modified equation analysis methd of manufactured solutions (MEAMMS). MEAMMS is then applied to a custom-built code, which solves the shallow water equations, to test the performance of the code verification method. MEAMMS is able to detect all coding errors that impact the implementation of the numerical scheme. To show how MEAMMS is different than MMS, they are both applied to the same first-order numerical method test problem with a first-order coding error. When there are first-order coding errors, only MEAMMS is able to identify them. Finally, this shows that MEAMMS is able to identify a larger set of coding errors while still being able to identify the coding errors MMS is able to identify.

42 ENGINEERING↗

DOE OSTI2023

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗

DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗

NASA NTRS1990

A mathematical model of the structure and evolution of small scale discrete auroral arcs

A three dimensional fluid model which includes the dispersive effect of electron inertia is used to study the nonlinear macroscopic plasma dynamics of small scale discrete auroral arcs within the auroral acceleration zone and ionosphere. The motion of the Alfven wave source relative to the magnetospheric and ionospheric plasma forms an oblique Alfven wave which is reflected from the topside ionosphere by the negative density gradient. The superposition of the incident and reflected wave can be described by a steady state analytical solution of the model equations with the appropriate boundary conditions. This two dimensional discrete auroral arc equilibrium provides a simple explanation of auroral acceleration associated with the parallel electric field. Three dimensional fully nonlinear numerical simulations indicate that the equilibrium arc configuration evolves three dimensionally through collisionless tearing and reconnection of the current layer. The interaction of the perturbed flow and the transverse magnetic field produces complex transverse structure that may be the origin of the folds and curls observed to be associated with small scale discrete arcs.

Seyler, C. E.↗

NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗

NASA NTRS1996

Application of the kappa-omega Turbulence Model to Quasi-Three-Dimensional Turbomachinery Flows

Many computational fluid dynamics codes for turbomachinery use the Baldwin-Lomax (B-L) turbulence model. It is easy to implement in two dimensions and works well for predicting overall turbomachinery performance. However, it is awkward to implement in three dimensions, often has difficulty finding the length scale, has a crude transition model, and neglects freestream turbulence, surface roughness, and mass injection. The kappa-omega model developed by Wilcox is an appealing alternative for several reasons. First, it is the only two-equation model that can be integrated to the wall without requiring damping functions or the distance to the wall, and hence, should behave well numerically. Second, the effects of freestream turbulence, surface roughness, and mass injection are easily included in the model. Finally, transition can be simulated using the low Reynolds number version of the model. Menter applied the kappa-model to external flows and showed very good results for flows with adverse pressure gradients. Liu and Zheng described their implementation of the kappa-model in a cascade code that included an area change term to account for endwall convergence. They validated the model for a flat plate, and compared the B-L and kappa-models to measured surface pressures for a low-pressure turbine cascade. Since they did not use the low Reynolds number version of the model, their results showed problems resulting from early transition. In this Note the low Reynolds number kappa-model was incorporated in the author's quasi-three-dimensional turbomachinery analysis code. The code includes the effects of rotation, radius change, and stream-surface thickness variation, and also includes the B-L turbulence model. The kappa-omega model was implemented using many of Menter's recommendations and an implicit approximate-factorization scheme described by Baldwin and Barth. The model was tested for a transonic compressor with rotation and variable stream-surface radius and height, and for a transonic turbine vane with transition and heat transfer. Results were compared to the B-L model and to experimental data.

Chima, Rodrick V.↗

NASA NTRS2002

Opportunities for Breakthroughs in Large-Scale Computational Simulation and Design

Opportunities for breakthroughs in the large-scale computational simulation and design of aerospace vehicles are presented. Computational fluid dynamics tools to be used within multidisciplinary analysis and design methods are emphasized. The opportunities stem from speedups and robustness improvements in the underlying unit operations associated with simulation (geometry modeling, grid generation, physical modeling, analysis, etc.). Further, an improved programming environment can synergistically integrate these unit operations to leverage the gains. The speedups result from reducing the problem setup time through geometry modeling and grid generation operations, and reducing the solution time through the operation counts associated with solving the discretized equations to a sufficient accuracy. The opportunities are addressed only at a general level here, but an extensive list of references containing further details is included. The opportunities discussed are being addressed through the Fast Adaptive Aerospace Tools (FAAST) element of the Advanced Systems Concept to Test (ASCoT) and the third Generation Reusable Launch Vehicles (RLV) projects at NASA Langley Research Center. The overall goal is to enable greater inroads into the design process with large-scale simulations.

Alexandrov, Natalia↗

NASA NTRS2004

Linearized Aeroelastic Solver Applied to the Flutter Prediction of Real Configurations

A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code.

Reddy, Tondapu S.↗

NASA NTRS2014

A Coupled Fluid-Structure Interaction Analysis of Solid Rocket Motor with Flexible Inhibitors

A capability to couple NASA production CFD code, Loci/CHEM, with CFDRC's structural finite element code, CoBi, has been developed. This paper summarizes the efforts in applying the installed coupling software to demonstrate/investigate fluid-structure interaction (FSI) between pressure wave and flexible inhibitor inside reusable solid rocket motor (RSRM). First a unified governing equation for both fluid and structure is presented, then an Eulerian-Lagrangian framework is described to satisfy the interfacial continuity requirements. The features of fluid solver, Loci/CHEM and structural solver, CoBi, are discussed before the coupling methodology of the solvers is described. The simulation uses production level CFD LES turbulence model with a grid resolution of 80 million cells. The flexible inhibitor is modeled with full 3D shell elements. Verifications against analytical solutions of structural model under steady uniform pressure condition and under dynamic condition of modal analysis show excellent agreements in terms of displacement distribution and eigen modal frequencies. The preliminary coupled result shows that due to acoustic coupling, the dynamics of one of the more flexible inhibitors shift from its first modal frequency to the first acoustic frequency of the solid rocket motor.

Yang, H. Q.↗

NASA NTRS1978

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗