Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid instabilities”

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 847 records · Page 47

Low-dimensional chaos in turbulence

Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Vastano, John A.↗

Gravitational Instability in Suspension Flows

The gravity-driven flow of non-neutrally buoyant suspensions is shown to be unstable to spanwise perturbations when the shearing motion generates a density profile that increases with height. The instability is simply due to having heavier material over light. The wavelength of the perturbation is found to be on the order of the thickness of the suspension layer. The parameters important to the problem are the angle of inclination of the layer relative to gravity, the relative density difference between the particles and fluid, the ratio of the particle size to the suspension layer, and the bulk volume fraction of particles. An example showing the growth rate as a function of wave number is shown.

Carpen, Ileana C.↗

NASA Common Research Model Test Envelope Extension With Active Sting Damping at NTF

The NASA Common Research Model (CRM) high Reynolds number transonic wind tunnel testing program was established to generate an experimental database for applied Computational Fluid Dynamics (CFD) validation studies. During transonic wind tunnel tests, the CRM encounters large sting vibrations when the angle of attack approaches the second pitching moment break, which can sometimes become divergent. CRM transonic test data analysis suggests that sting divergent oscillations are related to negative net sting damping episodes associated with flow separation instability. The National Transonic Facility (NTF) has been addressing remedies to extend polar testing up to and beyond the second pitching moment break point of the test articles using an active piezoceramic damper system for both ambient and cryogenic temperatures. This paper reviews CRM test results to gain understanding of sting dynamics with a simple model describing the mechanics of a sting-model system and presents the performance of the damper under cryogenic conditions.

Rivers, Melissa B.↗

The Onset of Resonance-Controlled Instability in Spherical Bubble Oscillations

Single bubble dynamics are investigated using acoustic techniques for isolation and manipulation. The goal of the investigations is to understand the dynamic origin of the various phenomena that bubbles exhibit: light emission, enhanced mass transport, chaotic and quasiperiodic oscillations and translations. Once understood, acoustically manipulated bubbles can serve as platforms for materials effects on free surfaces, using surfactants to alter surface rheology and observing how that affects both dynamics and also mass transport. The effects of gravity on the problem will be shown to be significant. The first set of observations from 1g experimentation are presented. These observations are of the onset conditions for instability of the spherical shape of the bubble. For the size range 55 - 90 microns in diameter we observe instability governed by resonant mode coupling, which is significantly affected by the buoyant force and its effects.

bubbles↗

A multi-dimensional dataset for two-phase instability in low pressure natural circulation based on direct transient local measurement

Thermal hydraulic instability can challenge the control and the safety of nuclear reactors. This experimental research captures and investigates the instability in a low-pressure (1–10 atmospheres) natural circulation loop with an adiabatic chimney installed above the heated section. Direct transient local measurement of two-phase flow is performed at multiple axial locations along the flow channel with thermocouples and conductivity probes traversed by customized high-precision mechanisms. Based on the periodicity of oscillations, a post-processing method is proposed to align the data captured non-simultaneously to reconstruct the instability behavior in one representative period. The averaging operation over multiple oscillation cycles is performed to reduce the random errors. Cyclic phenomena are clearly demonstrated by the post-processed data, and rapid void fraction change induced by the flashing is observed to be significant in the driving mechanism for low pressure oscillations. Parametric effects of inlet subcooling, heat flux, inlet flow restriction and system pressure are studied within the experiment range, and qualitative results are explained and compared with existing literature. The ratio between the oscillation period and the travelling time of the fluid particle through the test section is compared with the theoretical value for the density wave oscillation. The new dataset contains directly measured transient multi-dimensional two-phase data along the flow channel for natural circulation under density wave oscillations. The generated dataset is valuable for the validation of system analysis codes and CFD codes in the simulation of natural circulation phenomena and dynamic two-phase flow.

42 ENGINEERING↗

Thermal hydraulic analysis of the AFIP-7 irradiation test in the Advanced test Reactor – Model correlation and performance evaluation

We report as a system level code, RELAP5 is widely used in the nuclear field for the reactor hydraulic analysis. In the component/experimental level, it is frequently employed as well. This paper demonstrates that the RELAP5 flow simulation of the AFIP-7 irradiation experiment which has a complicated 3D geometry/flow path deviates by ~28% from the flow evaluation through a computational fluid dynamics (CFD) simulation, which has been validated against a flow test performed at Oregon State University. Thermal safety compliance analysis of an experiment planned to be irradiated in the Advanced Test Reactor is usually performed using the finite element analysis code, Abaqus. This research reveals the disadvantages of the Abaqus simulation in the flow instability and departure from nucleate boiling evaluations via a conjugate heat transfer analysis of the AFIP-7 irradiation experiment. As a result, a detailed CFD simulation is suggested for irradiation experiments with complicated flow paths, rather than a simplified RELAP5 simulation or hand calculations. The validated simulation approach will be integrated with the Boehmite correlations to investigate the oxide growth prediction on the fuel cladding in Part II of this research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗

Performance of end-face seals with diametral tilt and coning - Hydrodynamic effects

Hydrodynamic effects in end-face seals with diametral tilt and coning are analyzed. A closed-form solution for the axial separating force and the restoring and transverse moments is presented that covers the whole range from zero to full angular misalignment at various degrees of coning. Both low-pressure seals with cavitating flow and high-pressure seals with full fluid film are considered. The effect of coning is to reduce the axial force and the restoring and transverse moments compared to their magnitude in flat-face seals. Strong coupling between diametral tilt and transverse moment is demonstrated. This transverse moment which is entirely due to hydrodynamic effects can be the source of dynamic instability in the form of seal wobble.

Sharoni, A.↗

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The nature of the sunspot phenomenon. II - Internal overstable modes

It had been pointed out by Parker (1974) that the basic cause of the sunspot phenomenon is the enhanced heat transport in the magnetic field of the sunspot. The enhanced transport occurs through convective overstability which operates as a heat engine generating Alfven waves. The characteristics of the convective forces present are investigated along with questions concerning overstability and convectively driven Alfven waves. Relations regarding instability and convectively driven surface waves are discussed and attention is given to individual overstable Alfven modes. It is found that the form of an Alfven wave in the absence of convective forces is entirely arbitrary, so that waves with any arbitrary profile can be fitted into a vertical column of the field without disturbing the fluid outside. With the introduction of convective forces the situation changes so that the presence of lateral boundaries alters the form of the basic wave modes.

Parker, E. N.↗

An electroviscous damper

A new type of vibration damper for rotor systems was developed and tested. The damper contains electroviscous fluid which solidifies and provides Coulomb damping when an electric voltage is imposed across the fluid. The damping capacity is controlled by the voltage. The damper was incorporated in a flexible rotor system and found to be able to damp out high levels of unbalanced excitation. Other proven advantages include controllability, simplicity, and no requirement for oil supply. Still unconfirmed are the capabilities to eliminate critical speeds and to suppress rotor instabilities.

Nikolajsen, Jorgen L.↗

CFD code evaluation for internal flow modeling

Research on the computational fluid dynamics (CFD) code evaluation with emphasis on supercomputing in reacting flows is discussed. Advantages of unstructured grids, multigrids, adaptive methods, improved flow solvers, vector processing, parallel processing, and reduction of memory requirements are discussed. As examples, researchers include applications of supercomputing to reacting flow Navier-Stokes equations including shock waves and turbulence and combustion instability problems associated with solid and liquid propellants. Evaluation of codes developed by other organizations are not included. Instead, the basic criteria for accuracy and efficiency have been established, and some applications on rocket combustion have been made. Research toward an ultimate goal, the most accurate and efficient CFD code, is in progress and will continue for years to come.

Chung, T. J.↗

Dynamic Control of Aerodynamic Instabilities in Gas Turbine Engines

This lecture discusses the use of closed loop control at the component level to enhance the performance of gas turbine engines. The general theme is the suppression of flow instabilities (rotating stall and surge) through use of feedback, either actively or by means of the aeromechanical coupling provided by tailored structures. The basic concepts that underlie active control of turbomachinery instability, and their experimental demonstration, are first described for a centrifugal compressor. It is shown that the mechanism for stabilization is associated with damping of unsteady perturbations in the compression system, and the steady-state performance can thus remain virtually unaltered. Control of instability using a tailored structure is then discussed, along with experimental results illustrating the flow range extension achievable using this technique. A considerably more complex problem is presented by active control or rotating stall where the multi-dimensional features mean that distributed sensing and actuation are required. In addition, there are basic questions concerning unsteady fluid mechanics; these imply the need to resolve issues connected with identification of suitable signals as well as with definition of appropriate wave launchers for implementing the feedback. These issues are discussed and the results of initial successful demonstrations of active control of rotating stall in a single-stage and a three-stage axial compressor are presented. The lecture concludes with suggestions for future research on dynamic control of gas turbine engines.

E M Greitzer↗

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Explaining the lack of power degradation of energy confinement in wide pedestal quiescent H-modes via transport modeling

Wide pedestal quiescent H (WPQH)-mode is an attractive scenario for future burning plasmas as they operate without ELMs. WPQH is characterized by formation of a wider and higher pedestal (than quiescent H-mode), and broadband fluctuations in the pedestal. Unlike conventional H-modes, where the energy confinement time reduces with increasing heating power, the WPQH plasmas reported in this paper do not show power degradation of the energy confinement. As the injected neutral beam power was increased, reduced core (ρ ≤ 0.45) transport calculated by transp, as well as increased core temperatures, pressure gradient and diamagnetic E × B shear rate were observed. The reduction in the heat transport and rapid decrease in the ion temperature gradient scale length suggest the formation of an ion internal transport barrier (ITB) that was accompanied by increased stored energy in the core. Quasilinear turbulent transport modeling using the trapped gyro Landau fluid (tglf) code was used to predict the ITB and its turbulence stability properties. By using profiles and equilibria produced by matching the transp transport fluxes with the tglf transport model within the tgyro transport solver, the energy confinement time captures the experimentally observed insensitivity to the increased PNBI. Linear stability analysis reveals that drift-wave instabilities in the core are stabilized by E × B shear, T i /T e ratio and Shafranov shift; the latter was found to have the strongest effect on the turbulence suppression at the highest heating level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Combustion Instabilities Modeled

NASA Lewis Research Center's Advanced Controls and Dynamics Technology Branch is investigating active control strategies to mitigate or eliminate the combustion instabilities prevalent in lean-burning, low-emission combustors. These instabilities result from coupling between the heat-release mechanisms of the burning process and the acoustic flow field of the combustor. Control design and implementation require a simulation capability that is both fast and accurate. It must capture the essential physics of the system, yet be as simple as possible. A quasi-one-dimensional, computational fluid dynamics (CFD) based simulation has been developed which may meet these requirements. The Euler equations of mass, momentum, and energy have been used, along with a single reactive species transport equation to simulate coupled thermoacoustic oscillations. A very simple numerical integration scheme was chosen to reduce computing time. Robust boundary condition procedures were incorporated to simulate various flow conditions (e.g., valves, open ends, and choked inflow) as well as to accommodate flow reversals that may arise during large flow-field oscillations. The accompanying figure shows a sample simulation result. A combustor with an open inlet, a choked outlet, and a large constriction approximately two thirds of the way down the length is shown. The middle plot shows normalized, time-averaged distributions of the relevant flow quantities, and the bottom plot illustrates the acoustic mode shape of the resulting thermoacoustic oscillation. For this simulation, the limit cycle peak-to-peak pressure fluctuations were 13 percent of the mean. The simulation used 100 numerical cells. The total normalized simulation time was 50 units (approximately 15 oscillations), which took 26 sec on a Sun Ultra2.

Paxson, Daniel E.↗

Assessment of Artemis-1 Pogo Flight Instrumentation System

As the space industry continues to strive for more efficient launch vehicles, it must relyon increasingly accurate predictive models. Verification of models typically requires physical testing. Flight data measurements offer the most real and therefore the most accurate data for model correlation. As NASA prepares for the inaugural launch of Space Launch System(SLS), Artemis-1, they must rely heavily on predictive system models to ensure flight safety.NASA has implemented a Development Flight Instrumentation (DFI) system in hopes of recovering useful flight data to aid in model correlation. Historically, some of the most important flight measurements are those that monitor the potentially destructive dynamic interaction of the structural and propellant modes – a phenomenon known as pogo. Pogo is a dynamic instability that can occur on a launch vehicle during any phase of ascent. During this investigation an end-to-end assessment of the Artemis-I Pogo-related DFI was performed to identify any obstacles inherent in the current instrumentation system which may prevent successful measurement of the data necessary to validate the current predictive fluid-structural models. Input drive signals were approximated and applied to a state-space system model to derive predicted pressure and acceleration responses. These predicted responses were then fed through a simulation of the data acquisition process in order to recover predicted measurements. Finally, a mock Flight Data Analysis (FDA) was performed to assess the ability of these measurements to meet the Flight Test Objectives(FTO).

Adam Johnson↗