The effect of wall friction on the strength of shock waves in tubes and hydraulic jumps in channels
Explore the source record for details and available documents.
HYDRODYNAMICS↗
SEARCH · Engineering Papers
Results for “fluid wave”
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.
Computational Fluid Dynamics Study on the Effects of RATO Timing on the Scale Model Acoustic Test
The Scale Model Acoustic Test (SMAT) is a 5% scale test of the Space Launch System (SLS), which is currently being designed at Marshall Space Flight Center (MSFC). The purpose of this test is to characterize and understand a variety of acoustic phenomena that occur during the early portions of lift off, one being the overpressure environment that develops shortly after booster ignition. The SLS lift off configuration consists of four RS-25 liquid thrusters on the core stage, with two solid boosters connected to each side. Past experience with scale model testing at MSFC (in ER42), has shown that there is a delay in the ignition of the Rocket Assisted Take Off (RATO) motor, which is used as the 5% scale analog of the solid boosters, after the signal to ignite is given. This delay can range from 0 to 16.5ms. While this small of a delay maybe insignificant in the case of the full scale SLS, it can significantly alter the data obtained during the SMAT due to the much smaller geometry. The speed of sound of the air and combustion gas constituents is not scaled, and therefore the SMAT pressure waves propagate at approximately the same speed as occurs during full scale. However, the SMAT geometry is much smaller allowing the pressure waves to move down the exhaust duct, through the trench, and impact the vehicle model much faster than occurs at full scale. To better understand the effect of the RATO timing simultaneity on the SMAT IOP test data, a computational fluid dynamics (CFD) analysis was performed using the Loci/CHEM CFD software program. Five different timing offsets, based on RATO ignition delay statistics, were simulated. A variety of results and comparisons will be given, assessing the overall effect of RATO timing simultaneity on the SMAT overpressure environment.
Nielsen, Tanner↗
Computational fluid mechanics
Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.
Hassan, H. A.↗
Real-Time Nonlinear Model Predictive Controller for Multiple Degrees of Freedom Wave Energy Converters with Non-Ideal Power Take-Off
An increase in wave energy converter (WEC) efficiency requires not only consideration of the nonlinear effects in the WEC dynamics and the power take-off (PTO) mechanisms, but also more integrated treatment of the whole system, i.e., the buoy dynamics, the PTO system, and the control strategy. It results in an optimization formulation that has a nonquadratic and nonstandard cost functional. This article presents the application of real-time nonlinear model predictive controller (NMPC) to two degrees of freedom point absorber type WEC with highly nonlinear PTO characteristics. The nonlinear effects, such as the fluid viscous drag, are also included in the plant dynamics. The controller is implemented on a real-time target machine, and the WEC device is emulated in real-time using the WECSIM toolbox. The results for the successful performance of the design are presented for irregular waves under linear and nonlinear hydrodynamic conditions.
Haider, Ali S. (ORCID:0000000246513157)↗
Laboratory and theoretical models of planetary-scale instabilities and waves
Meteorologists and planetary astronomers interested in large-scale planetary and solar circulations recognize the importance of rotation and stratification in determining the character of these flows. In the past it has been impossible to accurately model the effects of sphericity on these motions in the laboratory because of the invariant relationship between the uni-directional terrestrial gravity and the rotation axis of an experiment. Researchers studied motions of rotating convecting liquids in spherical shells using electrohydrodynamic polarization forces to generate radial gravity, and hence centrally directed buoyancy forces, in the laboratory. The Geophysical Fluid Flow Cell (GFFC) experiments performed on Spacelab 3 in 1985 were analyzed. Recent efforts at interpretation led to numerical models of rotating convection with an aim to understand the possible generation of zonal banding on Jupiter and the fate of banana cells in rapidly rotating convection as the heating is made strongly supercritical. In addition, efforts to pose baroclinic wave experiments for future space missions using a modified version of the 1985 instrument led to theoretical and numerical models of baroclinic instability. Rather surprising properties were discovered, which may be useful in generating rational (rather than artificially truncated) models for nonlinear baroclinic instability and baroclinic chaos.
Hart, John E.↗
An improved ten-moment closure for reconnection and instabilities
The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Room temperature depinning of the charge-density waves in quasi-two-dimensional 1T-TaS 2 devices
Here we report on the depinning of nearly commensurate charge-density waves in 1T-TaS 2 thin films at room temperature. A combination of the differential current–voltage measurements with the low-frequency noise spectroscopy provides unambiguous means for detecting the depinning threshold field in quasi-2D materials. The depinning process in 1T-TaS 2 is not accompanied by an observable abrupt increase in electric current—in striking contrast to depinning in the conventional charge-density-wave materials with quasi-1D crystal structure. We explained it by the fact that the current density from the charge-density waves in the 1T-TaS 2 devices is orders of magnitude smaller than the current density of the free carriers available in the discommensuration network surrounding the commensurate charge-density wave islands. The depinning fields in 1T-TaS 2 thin-film devices are several orders of magnitude larger than those in quasi-1D van der Waals materials. Obtained results are important for the proposed applications of the charge-density wave devices in electronics.
2D materials↗
DPFEHM: a differentiable subsurface physics simulator
The Earth’s subsurface is a key resource that provides energy via fossil fuels and geothermal sources, stores drinking water, and is used in the fight against climate change via carbon sequestration. Simulating the physical processes that occur in the Earth’s subsurface with computers enables better use of this resource. DPFEHM is a Julia package that includes computer models with a focus on the Earth’s subsurface, especially fluid flow, which is critical for the aforementioned applications. DPFEHM is able to solve the groundwater flow equations (single phase flow), Richards equation (air/water), the advection-dispersion equation, and the 2d wave equation. One of the key features of DPFEHM is that it supports automatic differentiation, so it can be integrated into machine learning workflows using frameworks such as Flux or PyTorch. The automatic differentiation capabilities give it the same performance as adjoint methods.
54 ENVIRONMENTAL SCIENCES↗
Nonthermal turbulent heating in the solar envelope.
It is shown that MHD pulses, in the form of fast magnetosonic waves or solitons, can produce a strong electron-ion coupling capable of maintaining electron-proton temperature equilibrium in the solar envelope. The mechanism producing the nonthermal heating is the fluid-like modified two-stream instability, which, since it is essentially independent of the electron-proton temperature ratio and the value of beta, becomes a prime candidate for the anomalous collisions required by the fluid models inside a distance less than 30 solar radii, in order to explain the dominant features of the solar-wind flow.
Papadopoulos, K.↗
Preferential concentration of particles by turbulence
Direct numerical simulation of isotropic turbulence was used to investigate the effect of turbulence on the concentration fields of heavy particles. The hydrodynamic field was computed using 64-cubed points, and a statistically stationary flow was obtained by forcing the low-wave-number components of the velocity field. The particles used in the simulations were time advanced according to Stokes drag law and were also assumed to be much more dense than the fluid. Properties of the particle cloud were obtained by following the trajectories of 1,000,000 particles through the simulated flow fields. Three values of the ratio of the particle time constant to large-scale turbulence time scale were used in the simulations: 0.075, 0.15, and 0.52. The simulations show that the particles collect preferentially in regions of low vorticity and high strain rate. This preferential collection was most pronounced for particles with the instantaneous number density as much as 25 times the mean value. The fact that dense particles collect in regions of low vorticity and high strain in turn implies that turbulence may actually inhibit rather than enhance mixing of particles.
Squires, Kyle D.↗
NAS (Numerical Aerodynamic Simulation Program) technical summaries, March 1989 - February 1990
Given here are selected scientific results from the Numerical Aerodynamic Simulation (NAS) Program's third year of operation. During this year, the scientific community was given access to a Cray-2 and a Cray Y-MP supercomputer. Topics covered include flow field analysis of fighter wing configurations, large-scale ocean modeling, the Space Shuttle flow field, advanced computational fluid dynamics (CFD) codes for rotary-wing airloads and performance prediction, turbulence modeling of separated flows, airloads and acoustics of rotorcraft, vortex-induced nonlinearities on submarines, and standing oblique detonation waves.
Source record↗
Contributions from geomagnetic inverse theory to the study of hydromagnetic conditions near the core-mantle boundary
The Final Report on contributions from geomagnetic inverse theory to the study of hydromagnetic conditions near the core-mantle boundary (CMB) is presented. The original proposal was to study five questions concerning what the surface and satellite magnetic data imply about hydromagnetic and electromagnetic conditions near the CMB. The five questions are: (1) what do the surface and satellite data imply about the geomagnetic field B near the surface of the earth; (2) how does one extrapolate B down through the conducting mantle to the CMB; (3) if B on the CMB is visible, how accurately does it satisfy the frozen-flux approximation; (4) if frozen flux is a good approximation on the CMB, what can be inferred about the fluid velocity v in the upper core; and (5) if v at the CMB is visible, does it suggest any dynamical properties of the core, such as vertical advection, Alfven-inertial waves, link instabilities, or mantle effects. A summary of the research is provided.
Backus, George E.↗
New techniques for diffusing-wave spectroscopy
We present two new types of measurements that can be made with diffusing-wave spectroscopy (DWS), a form of dynamic light scattering that applies in limit of strong multiple scattering. The first application is to measure the frequency-dependent linear viscoelastic moduli of complex fluids using light scattering. This is accomplished by measuring the mean square displacement of probe particles using DWS. Their response to thermal fluctuations is determined by the fluctuation-dissipation relation, and is controlled by the response of the surrounding complex fluid. This response can be described in terms of a memory function, which is directly related to the complex elastic modulus of the system. Thus by measuring the mean square displacement, we are able to determine the frequency dependent modulus. The second application is the measurement of shape fluctuations of scattering particles. This is achieved by generalizing the theory for DWS to incorporate the effects if amplitude fluctuations in the scattering intensity of the particles. We apply this new method to study the thermally induced fluctuations in the shape of spherical emulsion droplets whose geometry is controlled by surface tension.
Mason, T. G.↗
Increasing Accuracy in Computed Inviscid Boundary Conditions
A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.
Dyson, Roger↗
Design and Analysis of A Jet Stretcher for High-Altitude Supersonic Free Jet Testing
A jet stretcher can be used to extend the region of clean flow for a test article in a supersonic free jet. The device acts to shield the article from shock waves reflected from the free jet shear layer. It can also allow for overexpanded nozzle operation through blocking shocks emanating from the nozzle lip, enlarging the potential test envelope. The FUN3D computational fluid dynamics solver was used to test a design process for developing jet stretchers. The inner surface of the axisymmetric jet stretcher was obtained from a streamline extracted from a free flight simulation around a notional test article at the intended test condition of Mach 3.2, 53,000 ft. The test article and jet stretcher were then simulated in a free jet nozzle exhaust, and the flowfield within a notional inlet aperture was compared between the installed and free flight calculations. Ambient pressure was increased until separation within the free jet nozzle finally induced disturbances at the inlet. The jet stretcher was shown to provide clean flow to the inlet at cell pressure ratios up to 2.8. A process for trimming the axisymmetric jet stretcher was also tested. This was intended to provide blockage relief for complex test articles, aiding with starting, as well as reducing weight and fabrication cost. The trimmed jet stretcher performed identically to the axisymmetric model at slightly-overexpanded conditions, but at cell pressure ratios above 1.3 allowed disturbances from the trailing edge to reach the inlet.
free jet↗
Design and Analysis of A Jet Stretcher for High-Altitude Supersonic Free Jet Testing
A jet stretcher can be used to extend the region of clean flow for a test article in a supersonic free jet. The device acts to shield the article from shock waves reflected from the free jet shear layer. It can also allow for overexpanded nozzle operation through blocking shocks emanating from the nozzle lip, enlarging the potential test envelope. The FUN3D computational fluid dynamics solver was used to test a design process for developing jet stretchers. The inner surface of the axisymmetric jet stretcher was obtained from a streamline extracted from a free flight simulation around a notional test article at the intended test condition of Mach 3.2, 53,000 ft. The test article and jet stretcher were then simulated in a free jet nozzle exhaust, and the flowfield within a notional inlet aperture was compared between the installed and free flight calculations. Ambient pressure was increased until separation within the free jet nozzle finally induced disturbances at the inlet. The jet stretcher was shown to provide clean flow to the inlet at cell pressure ratios up to 2.8. A process for trimming the axisymmetric jet stretcher was also tested. This was intended to provide blockage relief for complex test articles, aiding with starting, as well as reducing weight and fabrication cost. The trimmed jet stretcher performed identically to the axisymmetric model at slightly overexpanded conditions, but at cell pressure ratios above 1.3 allowed disturbances from the trailing edge to reach the inlet.
free jet↗
Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry
As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.
Eve Papajak↗
Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry
As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.
Eve Papajak↗