Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “oscillatory flows”

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 NTRS1991

High-order ENO schemes applied to two- and three-dimensional compressible flow

High order essentially non-oscillatory (ENO) finite difference schemes are applied to the 2-D and 3-D compressible Euler and Navier-Stokes equations. Practical issues, such as vectorization, efficiency of coding, cost comparison with other numerical methods, and accuracy degeneracy effects, are discussed. Numerical examples are provided which are representative of computational problems of current interest in transition and turbulence physics. These require both nonoscillatory shock capturing and high resolution for detailed structures in the smooth regions and demonstrate the advantage of ENO schemes.

Shu, Chi-Wang↗
DOE OSTI2022

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2004

Summary of the 2004 CFD Validation Workshop on Synthetic Jets and Turbulent Separation Control

A CFD validation workshop for synthetic jets and turbulent separation control (CFDVAL2004) was held in Williamsburg, Virginia in March 2004. Three cases were investigated: synthetic jet into quiescent air, synthetic jet into a turbulent boundary layer crossflow, and flow over a hump model with no-flow-control, steady suction, and oscillatory control. This paper is a summary of the CFD results from the workshop. Although some detailed results are shown, mostly a broad viewpoint is taken, and the CFD state-of-the-art for predicting these types of flows is evaluated from a general point of view. Overall, for synthetic jets, CFD can only qualitatively predict the flow physics, but there is some uncertainty regarding how to best model the unsteady boundary conditions from the experiment consistently. As a result. there is wide variation among CFD results. For the hump flow, CFD as a whole is capable of predicting many of the particulars of this flow provided that tunnel blockage is accounted for, but the length of the separated region compared to experimental results is consistently overpredicted.

Rumsey, C. L.↗
NASA NTRS2006

Summary of the 2004 CFD Validation Workshop on Synthetic Jets and Turbulent Separation Control

A computational fluid dynamics (CFD) validation workshop for synthetic jets and turbulent separation control (CFDVAL2004) was held in Williamsburg, Virginia in March 2004. Three cases were investigated: synthetic jet into quiescent air, synthetic jet into a turbulent boundary layer crossflow, and flow over a hump model with no-flow-control, steady suction, and oscillatory control. This paper is a summary of the CFD results from the workshop. Although some detailed results are shown, mostly a broad viewpoint is taken, and the CFD state-of-the-art for predicting these types of flows is evaluated from a general point of view. Overall, for synthetic jets, CFD can only qualitatively predict the flow physics, but there is some uncertainty regarding how to best model the unsteady boundary conditions from the experiment consistently. As a result, there is wide variation among CFD results. For the hump flow, CFD as a whole is capable of predicting many of the particulars of this flow provided that tunnel blockage is accounted for, but the length of the separated region compared to experimental results is consistently overpredicted.

Rumsey, C. L.↗
NASA NTRS1974

A finite element method for potential aerodynamics around complex configurations

A general formulation for steady and oscillatory, subsonic and supersonic, potential linearized aerodynamic flow around complex configurations is presented. A linear integral equation relating the unknown potential on the surface of the body to the known downwash is used. The formulation is applied to the analysis of the flow field around wings and wing-body combinations. The surface is divided into small quadrilateral elements which are approximated with a hyperboloidal surface. The potential is assumed to be constant within each element. This yields a set of linear algebraic equations. The coefficients are evaluated analytically. Numerical results for steady and oscillatory, subsonic and supersonic flows indicate that the method, intrinsically general and flexible, is also fast, accurate and in excellent agreement with existing results.

Chen, L.-T.↗
NASA NTRS1976

A steady and oscillatory kernel function method for interfering surfaces in subsonic, transonic and supersonic flow

The theory, results and user instructions for an aerodynamic computer program are presented. The theory is based on linear lifting surface theory, and the method is the kernel function. The program is applicable to multiple interfering surfaces which may be coplanar or noncoplanar. Local linearization was used to treat nonuniform flow problems without shocks. For cases with imbedded shocks, the appropriate boundary conditions were added to account for the flow discontinuities. The data describing nonuniform flow fields must be input from some other source such as an experiment or a finite difference solution. The results are in the form of small linear perturbations about nonlinear flow fields. The method was applied to a wide variety of problems for which it is demonstrated to be significantly superior to the uniform flow method. Program user instructions are given for easy access.

Cunningham, A. M., Jr.↗
NASA NTRS1993

Matching Experimental And Analytical Aerodynamic Data

Three methods of calculating correction factors developed to match steady lifting distributions of pressure, forces and moments on airfoil sections, and total forces and moments. Data for rectangular supercritical wing tested in NASA Langley Research Center Transonic Dynamics Tunnel used to determine correction factors to match distributions of surface pressure for mach numbers ranging from 0.266 to 0.8. These correction factors also applied to calculations of unsteady aerodynamic flow, and comparisons made with oscillatory experimental data for range of reduced frequencies at several mach numbers. Methods of combining analytical and experimental data should lead to more-accurate predictions of performance of airfoils. Also applied to such airfoils as those in gas-turbine engines, fans, and wind turbines.

Wieseman, Carol D.↗
NASA NTRS1996

Deformation of Fluid Column by Action of Axial Vibration and Some Aspects of High-Rate Thermocapillary Convection

The deformation of the fluid column by an action of a low-frequency vibration is considered. It is shown that behavior of the free fluid surface depends on the frequency of applied vibration and its amplitude. In the area of very low frequencies when fluid has time to comment on travel of bounding solid walls limiting column, the harmonical oscillations of free surface with given frequency are observed. With increase of vibration frequency the steady-state relief on free fluid surface is formed. If the amplitude of vibration is very small and the frequency corresponding to the first peak in the vibration spectrum on the Mir orbital station, the deformation of free surface tends to zero. Fluid flow induced thermocapillary effect on deformed free surface is more unstable as in the case of smooth cylindrical surface. It was shown that width of heating zone affects very essentially the flow pattern and transition to oscillatory regime of thermocapillary convection.

Feonychev, Alexander I.↗
NASA NTRS1976

Unsteady aerodynamic flow field analysis of the space shuttle configuration. Part 3: Unsteady aerodynamics of bodies with concave nose geometries

An analysis of the unsteady aerodynamics of bodies with concave nose geometries was performed. The results show that the experimentally observed pulsating flow on spiked bodies and in forward facing cavities can be described by the developed simple mathematical model of the phenomenon. Static experimental data is used as a basis for determination of the oscillatory frequency of spike-induced flow pulsations. The agreement between predicted and measured reduced frequencies is generally very good. The spiked-body mathematical model is extended to describe the pulsations observed in forward facing cavities and it is shown that not only the frequency but also the pressure time history can be described with the accuracy needed to predict the experimentally observed time average effects. This implies that it should be possible to determine analytically the impact of the flow pulsation on the structural integrity of the nozzles for the jettisoned empty SRM-shells.

Ericsson, L. E.↗
NASA NTRS1992

An experimental study of oscillatory thermocapillary convection in cylindrical containers

An experimental study of oscillatory thermocapillary in small cylindrical containers with a heating wire placed along the center axis is performed by investigating the flow structures and temperature distributions under various conditions. To supplement the flow visualization the surface is scanned using an infrared imager. Here, 2 cS viscosity (Pr = 27) silicone oil is used as the test fluid. It is observed that beyond a certain temperature difference between the container wall and the heating wire, a distinctive unsteady flow pattern appears. This unsteady phenomenon is identified as oscillatory thermocapillary. After the onset of oscillations the flow structure becomes nonaxisymmetric and wave motion is observed at the free surface. It is shown that the critical temperature difference is independent of container dimensions if the aspect ratio is fixed.

Kamotani, Y.↗
NASA NTRS1975

Fully unsteady subsonic and supersonic potential aerodynamics for complex aircraft configurations with applications to flutter

A general formulation is presented for the analysis of steady and unsteady, subsonic and supersonic aerodynamics for complex aircraft configurations. The theoretical formulation, the numerical procedure, the description of the program SOUSSA (steady, oscillatory and unsteady, subsonic and supersonic aerodynamics) and numerical results are included. In particular, generalized forces for fully unsteady (complex frequency) aerodynamics for a wing-body configuration, AGARD wing-tail interference in both subsonic and supersonic flows as well as flutter analysis results are included. The theoretical formulation is based upon an integral equation, which includes completely arbitrary motion. Steady and oscillatory aerodynamic flows are considered. Here small-amplitude, fully transient response in the time domain is considered. This yields the aerodynamic transfer function (Laplace transform of the fully unsteady operator) for frequency domain analysis. This is particularly convenient for the linear systems analysis of the whole aircraft.

Tseng, K.↗
NASA NTRS2013

Propellant-Flow-Actuated Rocket Engine Igniter

A rocket engine igniter has been created that uses a pneumatically driven hammer that, by specialized geometry, is induced into an oscillatory state that can be used to either repeatedly impact a piezoelectric crystal with sufficient force to generate a spark capable of initiating combustion, or can be used with any other system capable of generating a spark from direct oscillatory motion. This innovation uses the energy of flowing gaseous propellant, which by means of pressure differentials and kinetic motion, causes a hammer object to oscillate. The concept works by mass flows being induced through orifices on both sides of a cylindrical tube with one or more vent paths. As the mass flow enters the chamber, the pressure differential is caused because the hammer object is supplied with flow on one side and the other side is opened with access to the vent path. The object then crosses the vent opening and begins to slow because the pressure differential across the ball reverses due to the geometry in the tube. Eventually, the object stops because of the increasing pressure differential on the object until all of the kinetic energy has been transferred to the gas via compression. This is the point where the object reverses direction because of the pressure differential. This behavior excites a piezoelectric crystal via direct impact from the hammer object. The hammer strikes a piezoelectric crystal, then reverses direction, and the resultant high voltage created from the crystal is transferred via an electrode to a spark gap in the ignition zone, thereby providing a spark to ignite the engine. Magnets, or other retention methods, might be employed to favorably position the hammer object prior to start, but are not necessary to maintain the oscillatory behavior. Various manifestations of the igniter have been developed and tested to improve device efficiency, and some improved designs are capable of operation at gas flow rates of a fraction of a gram per second (0.001 lb/s) and pressure drops on the order of 30 to 50 kilopascal (a few psi). An analytical model has been created and tested in conjunction with a precisely calibrated reference model. The analytical model accurately captures the overall behavior of this innovation. The model is a simple "volume-orifice" concept, with each chamber considered a single temperature and pressure "node" connected to adjacent nodes, or to vent paths through flow control orifices. Mass and energy balances are applied to each node, with gas flow predicted using simple compressible flow equations.

Wollen, Mark↗
NASA NTRS1999

Fluid Physics in a Fluctuating Acceleration Environment

Our program of research aims at developing a stochastic description of the residual acceleration field onboard spacecraft (g-jitter) to describe in quantitative detail its effect on fluid motion. Our main premise is that such a statistical description is necessary in those cases in which the characteristic time scales of the process under investigation are long compared with the correlation time of g-jitter. Although a clear separation between time scales makes this approach feasible, there remain several difficulties of practical nature: (i), g-jitter time series are not statistically stationary but rather show definite dependences on factors such as active or rest crew periods; (ii), it is very difficult to extract reliably the low frequency range of the power spectrum of the acceleration field. This range controls the magnitude of diffusive processes; and (iii), models used to date are Gaussian, but there is evidence that large amplitude disturbances occur much more frequently than a Gaussian distribution would predict. The lack of stationarity does not constitute a severe limitation in practice, since the intensity of the stochastic components changes very slowly during space missions (perhaps over times of the order of hours). A separate analysis of large amplitude disturbances has not been undertaken yet, but it does not seem difficult a priori to devise models that may describe this range better than a Gaussian distribution. The effect of low frequency components, on the other hand, is more difficult to ascertain, partly due to the difficulty associated with measuring them, and partly because they may be indistinguishable from slowly changing averages. This latter effect is further complicated by the lack of statistical stationarity of the time series. Recent work has focused on the effect of stochastic modulation on the onset of oscillatory instabilities as an example of resonant interaction between the driving acceleration and normal modes of the system, and on cavity flow as an example of how an oscillatory response under periodic driving becomes diffusive if the forcing is random instead. This paper describes three different topics that illustrate behavior that is peculiar to a stochastic acceleration field. In the first case, we show that g-jitter can induce effective attractive or repulsive forces between a pair of spherical particles that are suspended in an incompressible fluid of different density provided that the momentum diffusion length is larger than the interparticle separation (as in the case in most colloidal suspensions). Second, a stochastic modulation of the control parameter in the vicinity of a pitchfork or supercritical bifurcation is known not to affect the location of the threshold. We show, however, that resonance between the modulation and linearly stable modes close to onset can lead to a shift in threshold. Finally, we discuss the classical problem of vorticity diffusion away from a plane boundary that is being vibrated along its own plane. Periodic motion with zero average vorticity production results in an exponential decay of the vorticity away from the boundary. Random vibration, on the other hand, results in power law decay away from the boundary even if vorticity production averages to zero.

Drolet, Francois↗
NASA NTRS1996

Unsteady, Transonic Flow Around Delta Wings Undergoing Coupled and Natural Modes Response: A Multidisciplinary Problem

The unsteady, three-dimensional Navier-Stokes equations coupled with the Euler equations of rigid-body dynamics are sequentially solved to simulate and analyze the aerodynamic response of a high angle of attack delta wing undergoing oscillatory motion. The governing equations of fluid flow and dynamics of the multidisciplinary problem are solved using a time-accurate solution of the laminar, unsteady, compressible, full Navier- Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The primary model under consideration consists of a 65 deg swept, sharp-edged, cropped delta wing of zero thickness at 20 deg angle of attack. In a freestream of Mach 0.85 and Reynolds number of 3.23 x 10(exp 6), the flow over the upper surface of the wing develops a complex shock system which interacts with the leading-edge primary vortices producing vortex breakdown. The effect of the oscillatory motion of the wing on the vortex breakdown and overall aerodynamic response is detailed to provide insight to the complicated physics associated with unsteady flows and the phenomenon of wing rock. Forced sinusoidal single and coupled mode rolling and pitching motion is presented for the wing in a transonic freestream. The Reynolds number, frequency of oscillation, and the phase angle are varied. Comparison between the single and coupled mode forced rolling and pitching oscillation cases illustrate the effects of coupling the motion. This investigation shows that even when coupled, forced rolling oscillation at a reduced frequency of 2(pi) eliminates the vortex breakdown which results in an increase in lift. The coupling effect for in phase forced oscillations show that the lift coefficient of the pitching-alone case and the rolling-moment coefficient of the rolling-alone case dominate the resulting response. However, with a phase lead in the pitching motion, the coupled motion results in a non-periodic response of the rolling moment. The second class of problems involve releasing the wing in roll to respond to the flowfield. Two models of sharp-edged delta wings, the previous 65 deg swept model and an 80 deg swept, sharp-edged delta wing, are used to observe the aerodynamic response of a wing free to roll in a transonic and subsonic freestream, respectively. These cases demonstrate damped oscillations, self-sustained limit cycle oscillations, and divergent rolling oscillations. Ultimately, an active control model using a mass injection system was applied on the surface of the wing to suppress the self-sustained limit cycle oscillation known as wing rock. Comparisons with experimental investigations complete this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Menzies, Margaret Anne↗
NASA NTRS2003

Flow Control Analysis on the Hump Model with RANS Tools

A concerted effort is underway at NASA Langley Research Center to create a benchmark for Computational Fluid Dynamic (CFD) codes. both unstructured and structured, against a data set for the hump model with actuation. The hump model was tested in the NASA Langley 0.3-m Transonic Cryogenic Tunnel. The CFD codes used for the analyses are the FUN2D (Full Unstructured Navier-Stokes 2-Dimensional) code, the structured TLNS3D (Thin-Layer Navier-Stokes 3-Dimensional) code, and the structured CFL3D code, all developed at NASA Langley. The current investigation uses the time-accurate Reynolds-Averaged Navier-Stokes (RANS) approach to predict aerodynamic performance of the active flow control experimental database for the hump model. Two-dimensional computational results verified that steady blowing and suction and oscillatory suction/blowing can be used to significantly reduce the separated flow region on the model. Discrepancies do exist between the CFD results and experimental data in the region downstream of the slot with the largest differences in the oscillatory cases. Overall, the structured CFD codes exhibited similar behavior with each other for a wide range of control conditions, with the unstructured FUN2D code showing moderately different results in the separated flow region for the suction and oscillatory cases.

Viken, Sally A.↗
NASA NTRS1994

Convection in the Physical Vapor Transport Process: Thermosolutal - Part 2

We consider the effect of an inert gas on the diffusive-convective physical vapor transport process. We investigate the case when the temperature gradient is stabilizing and the concentration gradient is destabilizing for a wide parametric range. When an inert gas is present, the thermal and solutal convection oppose each other. The solutal field is destabilizing while the thermal field and the advective-diffusive flux stabilize the flow field. When the pressure of the inert component is increased, the stabilizing effect of the advective-diffusive flux is decreased. The intensity of convection as well as the oscillatory transient time increases. Below, the critical Rayleigh number, the nonlinear dynamics of the flow field show an oscillatory approach to steady state. For parametric values in the neighborhood of the critical Rayleigh number, the flow field undergoes a chaotic transient which settles to a periodic state. The asymptotic state of the flow field shows that growth and amalgamation of cells yields an overturning motion which results in an asymmetric cellular structure. The low gravity environment yields the stabilizing advective-diffusive flow which results in uniform temperature and concentration gradients near the crystal interface.

Duval, Walter M. B.↗