Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

Flux-vector splitting for unsteady calculations on dynamic meshes

The method of flux vector splitting used is that of Van Leer. The fluxes split in this manner have the advantage of being continuously differentiable at eigenvalue sign changes and this allows normal shocks to be captured with at most two interior zones, although in practice only one zone is usually observed. The fluxes as originally derived, however did not include the necessary terms appropriate for calculations on a dynamic mesh. The extension of the splitting to include these terms while retaining the advantages of the original splitting is the main purpose of this investigation. In addition, the use of multiple grids to reduce the computer time is investigated. A subiterative procedure to eliminate factorization and linearization error so that larger time steps can be used is also investigated.

Anderson, W. Kyle↗

Aerothermal Analysis of the Dragonfly Titan Entry

The Dragonfly mission will send a rotorcraft lander to the surface of Saturn’s moon Titan as part of the New-Frontiers program. This will be the first spacecraft to land on Titan since the Huygens probe’s descent in January of 2005 and only the second spacecraft to enter Titan’s atmosphere. The Dragonfly entry capsule is significantly larger than the Huygens probe and will experience higher aerothermal environments. This poster will provide an overview of the aerothermodynamic models used in the Dragonfly aeroshell design process. Titan provides a unique entry environment that has several fundamental differences from the environments at the more common entry destinations of Earth and Mars. The lack of atmospheric oxygen significantly reduces heatshield recession (expected to be negligible), which simplifies material response analysis but also reduces the efficiency of the thermal protection system. The atmos-phere is composed of mostly nitrogen and trace amounts of methane. During hypersonic entry this methane dissociates and leads to the formation of the molecule CN, which is known to radiate strongly in the shock layer. Current Dragonfly heating predictions estimate that radiation contributes about 50% of the total heat flux along the forebody and up to 90% on the aftbody, making radiation a key element of aerothermal analysis. To properly address the importance of radiation, a unique radiative-heating correlation was developed that incorporates normal-shock equilibrium chemistry predictions based on work from the Mars 2020 program. These new heating correlations agree well with Computational Fluid Dynamics (CFD) predictions and enable large-scale trajectory analyses by providing rapid heating predictions along candidate trajectories. Dragonfly aerothermal environments are generated with the Data Parallel Line Relaxation (DPLR) finite-volume Navier-Stokes solver alongside the Nonequilibrium Radiative Transport and Spectra (NEQAIR) radiation transport code. DPLR produces shock-aligned flowfields that incorporate the effects of both chemical and thermal nonequilibrium while NEQAIR solves the radiation transport equations along lines-of-sight through the flowfield solutions to predict the radiative heating on the aeroshell. Titan’s atmosphere produces complex thermal and chemical nonequilibria and so classical radiation approximations that extrapolate from a single line-of-sight, such as the Tangent-Slab method, were found to overpredict peak radiative heat fluxes by 10-20%. Therefore, all radiative heating is directly calculated using 3-D transport for both the forebody and aftbody to reduce unnecessary conservatism. A 21-species finite-rate chemistry model, including electrons and ionized species, is used in the CFD to ensure accurate vehicle heating predictions (the level of ionization is particularly for accurate radiative flux predictions). This poster will cover: 1) The models and assumptions used to simulate entry into Titan’s atmosphere, 2) Dragonfly aerothermal engineering correlations alongside statistical heating values, and 3) Driving aerothermal considerations for a Titan entry such as CN formation and radiation-flowfield coupling.

EDL↗

Numerical simulation by TVD schemes of complex shock reflections from airfoils at high angle of attack

The shock-capturing capability of total variation diminishing (TVD) schemes is demonstrated for a more realistic complex shock-diffraction problem for which the experimental data are available. Second-order explicit upwind and symmetric TVD schemes are used to solve the time-dependent Euler equations of gas dynamics for the interaction of a blast wave with an airfoil at high angle-of-attack. The test cases considered are a time-dependent moving curved-shock wave and a contant moving planar-shock wave impinging at an angle-of-attack 30 deg on a NACA 0018 airfoil. Good agreement is obtained between isopycnic contours computed by the TVD schemes and those from experimental interferograms. No drastic difference in flow-field structure is found between the curved- and planar-shock wave cases, except for a difference in density level near the lower surface of the airfoil. Computation for cases with higher shock Mach numbers is also possible. Numerical experiments show that the symmetric TVD scheme is less sensitive to the boundary conditions treatment than the upwind scheme.

Moon, Young J.↗

Bennett ion mass spectrometers on the Pioneer Venus Bus and Orbiter

Identical Bennett radio-frequency ion mass spectrometer instruments on the Pioneer Venus Bus and Orbiter have provided the first in-situ measurements of the detailed composition of the planet's ionosphere. The sensitivity, resolution, and dynamic range are sufficient to provide measurements of the solar-wind-induced bow-shock, the ionopause, and highly structured distributions of up to 16 thermal ion species within the ionosphere. The use of adaptive scan and detection circuits and servo-controlled logic for ion mass and energy analysis permits detection of ion concentrations as low as 5 ions/cu cm and ion flow velocities as large as 9 km/sec for O(+). A variety of commandable modes provides ion sampling rates ranging from 0.1 to 1.6 sec between measurements of a single constituent. A lightweight sensor and electronics housing are features of a compact instrument package.

Taylor, H. A., Jr.↗

Role of computational fluid dynamics in unsteady aerodynamics for aeroelasticity

In the last two decades there have been extensive developments in computational unsteady transonic aerodynamics. Such developments are essential since the transonic regime plays an important role in the design of modern aircraft. Therefore, there has been a large effort to develop computational tools with which to accurately perform flutter analysis at transonic speeds. In the area of Computational Fluid Dynamics (CFD), unsteady transonic aerodynamics are characterized by the feature of modeling the motion of shock waves over aerodynamic bodies, such as wings. This modeling requires the solution of nonlinear partial differential equations. Most advanced codes such as XTRAN3S use the transonic small perturbation equation. Currently, XTRAN3S is being used for generic research in unsteady aerodynamics and aeroelasticity of almost full aircraft configurations. Use of Euler/Navier Stokes equations for simple typical sections has just begun. A brief history of the development of CFD for aeroelastic applications is summarized. The development of unsteady transonic aerodynamics and aeroelasticity are also summarized.

Guruswamy, Guru P.↗

The benchmark aeroelastic models program: Description and highlights of initial results

An experimental effort was implemented in aeroelasticity called the Benchmark Models Program. The primary purpose of this program is to provide the necessary data to evaluate computational fluid dynamic codes for aeroelastic analysis. It also focuses on increasing the understanding of the physics of unsteady flows and providing data for empirical design. An overview is given of this program and some results obtained in the initial tests are highlighted. The tests that were completed include measurement of unsteady pressures during flutter of rigid wing with a NACA 0012 airfoil section and dynamic response measurements of a flexible rectangular wing with a thick circular arc airfoil undergoing shock boundary layer oscillations.

Bennett, Robert M.↗

The benchmark aeroelastic models program: Description and highlights of initial results

An experimental effort was implemented in aeroelasticity called the Benchmark Models Program. The primary purpose of this program is to provide the necessary data to evaluate computational fluid dynamic codes for aeroelastic analysis. It also focuses on increasing the understanding of the physics of unsteady flows and providing data for empirical design. An overview is given of this program and some results obtained in the initial tests are highlighted. The tests that were completed include measurement of unsteady pressures during flutter of a rigid wing with an NACA 0012 airfoil section and dynamic response measurements of a flexible rectangular wing with a thick circular arc airfoil undergoing shock boundary layer oscillations.

Bennett, Robert M.↗

Crossing shock wave turbulent boundary layer interactions - Variable angle and shock generator length geometry effects at Mach 3

By comparing the detailed wall static pressure distributions for 9 inch and 11 inch long fins generating a crossing shock configuration at M = 2.93, the high resolution results of the 9 inch fins are shown to be free of exit effects. Analysis of the static pressure profiles have delineated the limited regions where the single fin results are valid. The characteristics of the complex interaction, with varying shock wave strength, have been described. The data provide a critical test for computational fluid dynamics which, in its initial phase, has performed poorly in predicting the measured wall static pressure distributions.

Bogdonoff, S. M.↗

On the spatial range of validity of the gas dynamic model in the magnetosheath of Venus

In the past, the global solar wind interaction with Venus has been treated with gas dynamic models which, while successful in modeling some of the global characteristics of the interaction, do not include the magnetic barrier in a self-consistent manner. This magnetic barrier is formed in the inner magnetosheath where it transfers solar wind momentum flux to the obstacle via magnetic pressure. In this study, we examine the extent to which the gas dynamic fluid approximation describes the magnetic field in the dayside Venus magnetosheath by comparing with two gas dynamic models, one which matches the observed ionopause location and one which matches the bow shock location. We find that each model predicts the field profile reasonably well in the vicinity of the matched bow shock or ionopause, but neither model provides an adequate model over the entire range from the ionopause to the bow shock.

Zhang, T. L.↗

The shock-absorbed system of the airplane landing gear

A discussion is given of the behavior of the shock-absorbing system, consisting of elastic struts and tires, under landing, take-off, and taxying conditions, and a general formula derived for obtaining the minimum stroke required to satisfy the conditions imposed on the landing gear. Finally, the operation of some typical shock-absorbing systems are examined and the necessity brought out for taking into account, in dynamic landing-gear tests, the effect of the wing lift at the instant of contact with the ground.

Callerio, Pietro↗

Structure of relativistic shocks in pulsar winds: A model of the wisps in the Crab Nebula

We propose a model of a optical 'wisps' of the Crab Nebula, features observed in the nebular synchrotron surface brightness near the central pulsar, as manifestations of the internal structure of the shock terminating the pulsar wind. We assume that this wind is composed of ions and a much denser plasma of electrons and positrons, frozen together to a toroidal magnetic field and flowing relativistically. We construct a form of solitary wave model of the shock structure in which we self-consistently solve for the ion orbits and the dynamics of the relativistically hot, magnetized e(+/-) background flow. We ignore dispersion in the ion energies, and we treat the pairs as an adiabatic fluid. The synchrotron emission enhancements, observed as the wisps, are then explained as the regions where reflection of the ions in the self-consistent magnetic field causes compressions of the e(+/-).

Gallant, Yves A.↗

Pressure Distribution Measurements at High Speed and Oblique Incidence of Flow

The present report contains the results of a series of observations obtained for a wing of symmetrical profile for different angles of yaw. The shock tunnel with 0.4m x 0.4m cross section, of the Institute for Gas Dynamics at L.F.A. - Braunchweig was available for this work. The profile used was a 9 percent, thick hyperbolic profile with maximum thickness 40 percent aft, which was calculated by the method of F. Ringleb (2). Since the conformal transformation of this profile is known, the theoretical pressure distribution could be determined exactly for the case of incompressible flow. Then by the use of the Prandtl rule one may extend the comparison of theory and experiment to the case of higher velocity of incident flow.

Lippisch, A.↗

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↗

Updates to Implicit Edge-Based Gradient Methods

In this paper, we report updates to the implicit edge-based gradient methods originally introduced in [H. Nishikawa, AIAA Paper 2020-3048, 2020]. First, we clarify the relationship between gradient accuracy and truncation error and show that the quadratic method involves a free parameter. Then, we provide a complete description and a simplified matrix form of the implicit gradient systems including a consistent boundary treatment, and derive a set of parameters for achieving fourth-order gradient accuracy on regular tetrahedral grids. A stability analysis is performed for a relaxation scheme used to solve the implicit gradient systems, and the result serves as a guide for choosing parameters. Numerical results are shown for accuracy verification and also for realistic inviscid flow problems in three dimensions, including flows with shock waves.

Weighted Least-Squares↗

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗

A new procedure for dynamic adaption of three-dimensional unstructured grids

A new procedure is presented for the simultaneous coarsening and refinement of three-dimensional unstructured tetrahedral meshes. This algorithm allows for localized grid adaption that is used to capture aerodynamic flow features such as vortices and shock waves in helicopter flowfield simulations. The mesh-adaption algorithm is implemented in the C programming language and uses a data structure consisting of a series of dynamically-allocated linked lists. These lists allow the mesh connectivity to be rapidly reconstructed when individual mesh points are added and/or deleted. The algorithm allows the mesh to change in an anisotropic manner in order to efficiently resolve directional flow features. The procedure has been successfully implemented on a single processor of a Cray Y-MP computer. Two sample cases are presented involving three-dimensional transonic flow. Computed results show good agreement with conventional structured-grid solutions for the Euler equations.

Biswas, Rupak↗