Search NASA⌕ Search

SEARCH · Search NASA

Results for “ZERO VELOCITY CURVE”

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.

53 records · Page 3

Scale Effect on Clark Y Airfoil Characteristics from NACA Full-Scale Wind-Tunnel Tests

This report presents the results of wind tunnel tests conducted to determine the aerodynamic characteristics of the Clark Y airfoil over a large range of Reynolds numbers. Three airfoils of aspect ratio 6 and with 4, 6, and 8 foot chords were tested at velocities between 25 and 118 miles per hour, and the characteristics were obtained for Reynolds numbers (based on the airfoil chord) in the range between 1,000,000 and 9,000,000 at the low angles of attack, and between 1,000,000 and 6,000,000 at maximum lift. With increasing Reynolds number the airfoil characteristics are affected in the following manner: the drag at zero lift decreases, the maximum lift increases, the slope of the lift curve increases, the angle of zero lift occurs at smaller negative angles, and the pitching moment at zero lift does not change appreciably.

Silverstein, Abe↗

Numerical analysis of a NACA0012 airfoil with leading edge ice accretions

Analysis of a NACA0012 airfoil with leading edge ice has been performed using a Navier-Stokes code coupled with a grid generation code. The computed results were compared to experimental information obtained for an airfoil with a well defined artificial ice shape. The computations were performed at angles of attack ranging from zero to ten degrees. This range is sufficient to show the development of the separation bubble aft of the ice shape on both the upper and lower surfaces. Velocity profile plots in the separation bubble are examined in order to determine if recirculation patterns are predicted properly and if separation and reattachment points are found within the resolution of the experimental information. Also, the massive separation near the point of stall is examined in order to more accurately evaluate the lift coefficient curve in that region. Lift, drag, and moment coefficients are computed and compared to experiment.

Potapczuk, Mark G.↗

Magnus induced diode effect for skyrmions in channels with periodic potentials

Using a particle based model, we investigate the skyrmion dynamical behavior in a channel where the upper wall contains divots of one depth and the lower wall contains divots of a different depth. Under an applied driving force, skyrmions in the channels move with a finite skyrmion Hall angle that deflects them toward the upper wall for -x direction driving and the lower wall for +x direction driving. When the upper divots have zero height, the skyrmions are deflected against the flat upper wall for -x direction driving and the skyrmion velocity depends linearly on the drive. For +x direction driving, the skyrmions are pushed against the lower divots and become trapped, giving reduced velocities and a nonlinear velocity-force response. When there are shallow divots on the upper wall and deep divots on the lower wall, skyrmions get trapped for both driving directions; however, due to the divot depth difference, skyrmions move more easily under -x direction driving, and become strongly trapped for +x direction driving. Further, the preferred -x direction motion produces what we call a Magnus diode effect since it vanishes in the limit of zero Magnus force, unlike the diode effects observed for asymmetric sawtooth potentials. We show that the transport curves can exhibit a series of jumps or dips, negative differential conductivity, and reentrant pinning due to collective trapping events. We also discuss how our results relate to recent continuum modeling on a similar skyrmion diode system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Liquid-vapour surface sensors for liquid nitrogen and hydrogen

The present paper identifies devices to serve as liquid-vapor detectors in zero gravity. The testing in LH2 was done in a sealed glass Dewar system to eliminate any chance of mixing H2 and air. Most of the tests were performed with the leads to the sensor horizontal. Some results of rapid cycle testing of LVDG in LH2 are presented. Findings of rapid-cycle testing of LVDG in LH2 are discussed. The sensor crossed the liquid surface when the position sensor registered 1.9 V, which occurred at about 0.4075 s. The delay time was about 1.5 ms. From the estimated slope of the position sensor curve at 1.9 V, the velocity of the sensor through the liquid surface is over 3 m/s. Results of tests of optical sensors are presented as well.

Siegwarth, J. D.↗

Calculation of three-dimensional (3-D) internal flow by means of the velocity-vorticity formulation on a staggered grid

A method has been developed to accurately compute the viscous flow in three-dimensional (3-D) enclosures. This method is the 3-D extension of a two-dimensional (2-D) method developed for the calculation of flow over airfoils. The 2-D method has been tested extensively and has been shown to accurately reproduce experimental results. As in the 2-D method, the 3-D method provides for the non-iterative solution of the incompressible Navier-Stokes equations by means of a fully coupled implicit technique. The solution is calculated on a body fitted computational mesh incorporating a staggered grid methodology. In the staggered grid method, the three components of vorticity are defined at the centers of the computational cell sides, while the velocity components are defined as normal vectors at the centers of the computational cell faces. The staggered grid orientation provides for the accurate definition of the vorticity components at the vorticity locations, the divergence of vorticity at the mesh cell nodes and the conservation of mass at the mesh cell centers. The solution is obtained by utilizing a fractional step solution technique in the three coordinate directions. The boundary conditions for the vorticity and velocity are calculated implicitly as part of the solution. The method provides for the non-iterative solution of the flow field and satisfies the conservation of mass and divergence of vorticity to machine zero at each time step. To test the method, the calculation of simple driven cavity flows have been computed. The driven cavity flow is defined as the flow in an enclosure driven by a moving upper plate at the top of the enclosure. To demonstrate the ability of the method to predict the flow in arbitrary cavities, results will he shown for both cubic and curved cavities.

Stremel, Paul M.↗

LAMOST Time-Domain survey: first results of four K2 plates

From Oct. 2019 to Apr. 2020, LAMOST performed a time-domain (TD) spectroscopic survey of four K2 plates with both low- and medium-resolution observations. The low-resolution spectroscopic survey acquired 282 exposures ( ≈ 46.6 h) over 25 nights, yielding a total of about 767 000 spectra, and the medium-resolution survey took 177 exposures ( ≈ 49.1 h) over 27 nights, collecting about 478 000 spectra. More than 70%/50% of low-resolution/medium-resolution spectra have signal-to-noise ratio higher than 10. We determine stellar parameters (e.g., T eff , log g, [Fe/H]) and radial velocity (RV) with different methods, including LASP, DD-Payne and SLAM. In general, these parameter estimations from different methods show good agreement, and the stellar parameter values are consistent with those of APOGEE. We use the Gaia DR2 RV values to calculate a median RV zero point (RVZP) for each spectrograph exposure by exposure, and the RVZP-corrected RVs agree well with the APOGEE data. The stellar evolutionary and spectroscopic masses are estimated based on the stellar parameters, multi-band magnitudes, distances and extinction values. Finally, we construct a binary catalog including about 2700 candidates by analyzing their light curves, fitting the RV data, calculating the binarity parameters from medium-resolution spectra and cross-matching the spatially resolved binary catalog from Gaia EDR3. In conclusion, the LAMOST TD survey is expected to represent a breakthrough in various scientific topics, such as binary systems, stellar activity, stellar pulsation, etc.

79 ASTRONOMY AND ASTROPHYSICS↗

Synthetic Jet Interaction With A Turbulent Boundary Layer Flow

Perhaps one of the more notable advances to have occurred in flow control technology in the last fifteen years is the application of surface-issuing jets for separation control on aerodynamic surfaces. The concept was introduced by Johnston and Night (1990) who proposed using circular jets, skewed and inclined to the wall, to generate streamwise vortices for the purpose of mitigating boundary layer separation. The skew and inclination angles have subsequently been shown to affect the strength and sign of the ensuing vortices. With a non-circular orifice, in addition to skew and inclination, the yaw angle of the major axis of the orifice can influence the flow control effectiveness of the jet. In particular, a study by Chang arid Collins (1997) revealed that a non-circular orifice, yawed relative to the freestream, can be used to control the size and strength of the vortices produced by the control jet. This early work used jets with only a steady injection of mass. Seifert et al. revealed that an unsteady blowing jet, could be as effective at separation control as a steady jet but with less mass flow. Seifert et al. showed that small amplitude blowing oscillations superimposed on a low momentum steady jet Was the most effective approach to delaying separation on a NACA 0015 airfoil at post-stall angles of attack. More recent work suggests that perhaps the most efficient jet control effect comes from a synthetic (oscillatory) jet where the time-averaged mass flux through the orifice is zero, but the net wall normal momentum is non-zero. The control effectiveness of synthetic jets has been demonstrated for several internal and external flow fields used synthetic jet control on a thick, blunt-nosed airfoil to delay stall well beyond the stall angles for the uncontrolled airfoil and with a dramatic increase in the lift-to-drag performance. Amitay et al. used an array of synthetic jets to mitigate flow separation in curved and diffusing ducts. While the control effectiveness of synthetic jet actuators in an application setting has been confirmed through a number of studies, the fluid dynamic mechanism by which control is effected is not well understood. Issues like yaw angle arid velocity ratio that have been studied for steady control jets have not been investigated for synthetic jets. Moreover, the role played by the inherently unsteady nature of a synthetic jet in the interaction with the controlled flow is not known. Some recent work by Rinehart and Glazer and Smith suggests that away from tile immediate vicinity of the jet orifice the flow field is steady in the mean, and that the interaction of a synthetic jet and a boundary layer creates a secondary flow in the boundary layer consisting primarily of streamwise vortices.

Smith, Douglas R.↗

EP 250108a/SN 2025kg: Observations of the Most Nearby Broad-line Type Ic Supernova Following an Einstein Probe Fast X-Ray Transient

With a small sample of fast X-ray transients (FXTs) with multiwavelength counterparts discovered to date, their progenitors and connections to γ-ray bursts (GRBs) and supernovae (SNe) remain ambiguous. Here, we present photometric and spectroscopic observations of SN 2025kg, the SN counterpart to the FXT EP 250108a. At z = 0.17641, this is the closest known SN discovered following an Einstein Probe (EP) FXT. We show that SN 2025kg’s optical spectra reveal the hallmark features of a broad-lined Type Ic SN. Its light-curve evolution and expansion velocities are comparable to those of GRB-SNe, including SN 1998bw, and two past FXT-SNe. We present JWST/NIRSpec spectroscopy taken around SN 2025kg’s maximum light, and find weak absorption due to He I 1.0830 μm and 2.0581 μm and a broad, unidentified emission feature at ∼4–4.5 μm. Further, we observe broadened Hα in optical data at 42.5 days that is not detected at other epochs, indicating interaction with H-rich material. From its light curve, we derive a 56 Ni mass of 0.2–0.6 M ⊙ . Together with our companion Letter, our broadband data are consistent with a trapped or low-energy (≲10 51 erg) jet-driven explosion from a collapsar with a zero-age main-sequence mass of 15–30 M ⊙ . Finally, we show that the sample of EP FXT-SNe supports past estimates that low-luminosity jets seen through FXTs are more common than successful (GRB) jets, and that similar FXT-like signatures are likely present in at least a few percent of the brightest Type Ic-BL SNe.

79 ASTRONOMY AND ASTROPHYSICS↗

Experimental data and model for the turbulent boundary layer on a convex, curved surface

Experiments were performed to determine how boundary layer turbulence is affected by strong convex curvature. The data gathered on the behavior of the Reynolds stress suggested the formulation of a simple turbulence model. Data were taken on two separate facilities. Both rigs had flow from a flat surface, over a convex surface with 90 deg of turning and then onto a flat recovery surface. The geometry was adjusted so that, for both rigs, the pressure gradient along the test surface was zero. Two experiments were performed at delta/R approximately 0.10, and one at weaker curvature with delta/R approximately 0.05. Results show that after a sudden introduction of curvature the shear stress in the outer part of the boundary layer is sharply diminished and is even slightly negative near the edge. The wall shear also drops off quickly downstream. When the surface suddenly becomes flat again, the wall shear and shear stress profiles recover very slowly towards flat wall conditions. A simple turbulence model, which was based on the theory that the Prandtl mixing length in the outer layer should scale on the velocity gradient layer, was shown to account for the slow recovery.

Gillis, J. C.↗

Flight Investigation to Determine the Aerodynamic Characteristics of Rocket-Powered Models Representative of a Fighter-Type Airplane Configuration Incorporating an Inverse-Taper Wing and a Vee Tail

Two rocket-powered models representative of a fighter-type airplane were investigated in flight at Mach numbers up to 1.01 and 1.07 by the Langley Pilotless Aircraft Research Division at its testing station at Wallops Island, Va. These models incorporated an inverse-taper wing and a vee tail and were flown with controls undeflected and wing and stabilizer set at 0 deg incidence. Values of lateral acceleration, normal acceleration velocity, and drag were obtained by use of telemeters and a Doppler velocimeter radar unit. The results of this investigation indicated no unusual variation in the lateral acceleration characteristics. After the cessation of powered flight, the lateral oscillation quickly damped to zero. The data indicated that the airplane, at low lift coefficients, should not experience any abrupt trim changes until it attains a Mach number of 0.97. The change in normal-force coefficient associated with this trim change will amount to about 0.03 with the center of gravity located at 4.48% of the mean aerodynamic chord. At higher lift coefficients, on the basis of other data, the Mach number at which this trim change occurs would be expected to be decreased. The neutral point of the model at Mach numbers near 1.05 was estimated to fall at 45% of the mean aerodynamic chord, assuming a lift-curve slope of 0.05. A value of the static-directional-stability parameter dCn/d(psi) of approximately -0.002 was estimated for a Mach number of 0.93. The values of drag coefficient obtained from both model flights were in a good comparative agreement. The highest drag coefficient occurred at a Mach number of 1.01 and was equal to 0.044.

Alexander, Sidney R.↗

The Atacama Cosmology Telescope: Large-scale velocity reconstruction with the kinematic Sunyaev-Zel'dovich effect and DESI LRGs

The kinematic Sunyaev-Zel'dovich (kSZ) effect induces a non-zero density-density-temperature bispectrum, which we can use to reconstruct the large-scale velocity field from a combination of cosmic microwave background (CMB) and galaxy density measurements, in a procedure known as “kSZ velocity reconstruction”. This method has been forecast to constrain large-scale modes with future galaxy and CMB surveys, improving their measurement beyond what is possible with the galaxy surveys alone. Such measurements will enable tighter constraints on large-scale signals such as primordial non-Gaussianity, deviations from homogeneity, and modified gravity. In this work, we demonstrate a statistically significant measurement of kSZ velocity reconstruction for the first time, by applying quadratic estimators to the combination of the ACT DR6 CMB+kSZ map and the DESI LRG galaxies (with photometric redshifts) in order to reconstruct the velocity field. We do so using a formalism appropriate for the 2-dimensional projected galaxy fields that we use, which naturally incorporates the curved-sky effects important on the largest scales. We find evidence for the signal by cross-correlating with an external estimate of the velocity field from the spectroscopic BOSS survey and rejecting the null (no-kSZ) hypothesis at 3.8σ. Our work presents a first step towards the use of this observable for cosmological analyses.

Sunyaev-Zeldovich effect↗

An extended theory of thin airfoils and its application to the biplane problem

The report presents a new treatment, due essentially to von Karman, of the problem of the thin airfoil. The standard formulae for the angle of zero lift and zero moment are first developed and the analysis is then extended to give the effect of disturbing or interference velocities, corresponding to an arbitrary potential flow, which are superimposed on a normal rectilinear flow over the airfoil. An approximate method is presented for obtaining the velocities induced by a 2-dimensional airfoil at a point some distance away. In certain cases this method has considerable advantage over the simple "lifting line" procedure usually adopted. The interference effects for a 2-dimensional biplane are considered in the light of the previous analysis. The results of the earlier sections are then applied to the general problem of the interference effects for a 3-dimensional biplane, and formulae and charts are given which permit the characteristics of the individual wings of an arbitrary biplane without sweepback or dihedral to be calculated. In the final section the conclusions drawn from the application of the theory to a considerable number of special cases are discussed, and curves are given illustrating certain of these conclusions and serving as examples to indicate the nature of the agreement between the theory and experiment.

Millikan, Clark B↗

Isolating Curvature Effects in Computing Wall-Bounded Turbulent Flows

The flow over the zero-pressure-gradient So-Mellor convex curved wall is simulated using the Navier-Stokes equations. An inviscid effective outer wall shape, undocumented in the experiment, is obtained by using an adjoint optimization method with the desired pressure distribution on the inner wall as the cost function. Using this wall shape with a Navier-Stokes method, the abilities of various turbulence models to simulate the effects of curvature without the complicating factor of streamwise pressure gradient can be evaluated. The one-equation Spalart-Allmaras turbulence model overpredicts eddy viscosity, and its boundary layer profiles are too full. A curvature-corrected version of this model improves results, which are sensitive to the choice of a particular constant. An explicit algebraic stress model does a reasonable job predicting this flow field. However, results can be slightly improved by modifying the assumption on anisotropy equilibrium in the model's derivation. The resulting curvature-corrected explicit algebraic stress model possesses no heuristic functions or additional constants. It lowers slightly the computed skin friction coefficient and the turbulent stress levels for this case (in better agreement with experiment), but the effect on computed velocity profiles is very small.

Rumsey, Christopher L.↗

Behavior of Langmuir Probes in Non-Equilibrium Plasmas

Langmuir probes are diagnostic tools used to determine electron temperature, number density, and plasma potential. Irving Langmuir first used an electrostatic probe in the 1920s to find these characteristics in ionized gases. Single, double, and triple Langmuir probes are commonly used in plasma diagnostics because of their relative simplicity. In the single probe, a swept voltage is applied between the probe tip and circuit common to acquire a waveform showing the collected current as a function of applied voltage. A double Langmuir probe consists of two tips, both inserted into the plasma, with a voltage applied between them. As this voltage is swept, a current-voltage characteristic is measured. In a triple probe three probe tips are electrically coupled to each other with constant non-swept voltages applied between each of the tips. The voltages are selected to represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Triple Langmuir probe measurements have been widely employed for various types of plasmas, including pulsed and time-varying plasmas such as those seen in pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. The typical Langmuir probe analysis for determining electron temperature and number density of the plasma (for a single, double, or triple Langmuir probe) includes an assumption that the plasma is in thermal equilibrium. While the this assumption may be justified for some applications, it is unlikely that it is fully justifiable for pulsed and time-varying plasmas or for the entire time a plasma device is in use. In the present work, we model the responses of Langmuir probes as they are inserted into a range of simple equilibrium and non-equilibrium plasmas. We return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. Comparing the distribution function that was assumed a priori with the plasma density and temperature one would infer when applying standard probe theory to analyze the collected currents serves to illustrate the effect a non- Maxwellian plasma would have on results interpreted using the equilibrium probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

The triple Langmuir probe is an electrostatic probe in which three probe tips collect current when inserted into a plasma. The triple probe differs from a simple single Langmuir probe in the nature of the voltage applied to the probe tips. In the single probe, a swept voltage is applied to the probe tip to acquire a waveform showing the collected current as a function of applied voltage (I-V curve). In a triple probe three probe tips are electrically coupled to each other with constant voltages applied between each of the tips. The voltages are selected such that they would represent three points on the single Langmuir probe I-V curve. Elimination of the voltage sweep makes it possible to measure time-varying plasma properties in transient plasmas. Under the assumption of a Maxwellian plasma, one can determine the time-varying plasma temperature T(sub e)(t) and number density n(sub e)(t) from the applied voltage levels and the time-histories of the collected currents. In the present paper we examine the theory of triple probe operation, specifically focusing on the assumption of a Maxwellian plasma. Triple probe measurements have been widely employed for a number of pulsed and timevarying plasmas, including pulsed plasma thrusters (PPTs), dense plasma focus devices, plasma flows, and fusion experiments. While the equilibrium assumption may be justified for some applications, it is unlikely that it is fully justifiable for all pulsed and time-varying plasmas or for all times during the pulse of a plasma device. To examine a simple non-equilibrium plasma case, we return to basic governing equations of probe current collection and compute the current to the probes for a distribution function consisting of two Maxwellian distributions with different temperatures (the two-temperature Maxwellian). A variation of this method is also employed, where one of the Maxwellians is offset from zero (in velocity space) to add a suprathermal beam of electrons to the tail of the main Maxwellian distribution (the bump-on-the-tail distribution function). For a range of parameters in these non-Maxwellian distributions, we compute the current collection to the probes. We compare the distribution function that was assumed a priori with the distribution function one would infer when applying standard triple probe theory to analyze the collected currents. For the assumed class of non-Maxwellian distribution functions this serves to illustrate the effect a non-Maxwellian plasma would have on results interpreted using the equilibrium triple probe current collection theory, allowing us to state the magnitudes of these deviations as a function of the assumed distribution function properties.

Polzin, Kurt A.↗

Turbulence Scales, Rise Times, Caustics, and the Simulation of Sonic Boom Propagation

The general topic of atmospheric turbulence effects on sonic boom propagation is addressed with especial emphasis on taking proper and efficient account of the contributions of the portion oi the turbulence that is associated with extremely high wavenumber components. The recent work reported by Bart Lipkens in his doctoral thesis is reexamined to determine whether the good agreement between his measured rise times with the 1971 theory of the author is fortuitous. It is argued that Lipken's estimate of the distance to the first caustic was a gross overestimate because of the use of a sound speed correlation function shaped like a gaussian curve. In particular, it is argued that the expected distance to the first caustic varies with the kinematic viscosity nu and the energy epsilon dissipated per unit mass per unit time, and the sound speed c as : d(sub first caustic) = nu(exp 7/12) c(exp 2/3)/ epsilon(exp 5/12)(nu x epsilon/c(exp 4))(exp a), where the exponent a is greater than -7/12 and can be argued to be either O or 1/24. In any event, the surprising aspect of the relationship is that it actually goes to zero as the viscosity goes to zero with s held constant. It is argued that the apparent overabundance of caustics can be grossly reduced by a general computational and analytical perspective that partitions the turbulence into two parts, divided by a wavenumber k(sub c). Wavenumbers higher than kc correspond to small-scale turbulence, and the associated turbulence can be taken into account by a renormalization of the ambient sound speed so that the result has a small frequency dependence that results from a spatial averaging over of the smaller-scale turbulent fluctuations. Selection of k(sub c). can be made so large that only a very small number of caustics are encountered if one adopts the premise that the frequency dispersion of pulses is caused by that part of the turbulence spectrum which lies in the inertial range originally predicted by Kolmogoroff. The acoustic propagating wave's dispersion relation has the acoustic wavenumber being of the form k = (omega/c) + F(omega), where c is a spatially averaged sound speed and where, for mechanical turbulence, the extra term F(omega) must depend on only the angular frequency omega, the sound speed c, and the turbulent energy dissipation epsilon per unit fluid mass and per unit time. If the turbulence is weak, then the quantity F(omega) has to be of second order in the portions of the turbulent fluid velocity in the inertial range, so, following Kolmogoroff's reasoning, it must vary with epsilon as epsilon(exp 2/3). Simple dimensional analysis then reveals that F(omega) is K epsilon(exp 2/3) c(exp -7/3) omega(exp l/3), K being a universal dimensionless complex constant.

Pierce, Allan D.↗

Hall Magnetohydrodynamic Power Generation and Drag Augmentation Using A Coaxial-Electrode Configuration During Hypersonic Entry

For interplanetary missions, a spacecraft must reduce its relative velocity from orbital speeds of multiple kilometers per second to zero in order to safely land on the surface of a planetary body. During this phase known as planetary reentry, the spacecraft undergoes high speed above Mach 5 up to Mach 30. This hyper-sonic flight regime imposes a strong bow shock in front of the vehicle that imparts extreme aerodynamic and thermal effects. In addition, during portions of this regime, a plasma is formed in the post-shock flow-field around the vehicle due to these high temperatures. Interaction between an applied magnetic field with this electrically conductive fluid, or magnetohydrodynamic (MHD) interaction, can be utilized to convert kinetic energy from the flow into storable electrical energy. Onboard MHD energy generation can benefit the spacecraft by increasing the power of active thermal control components and vehicle control systems during reentry. MHD interaction also exerts a Lorentz body force on the plasma flow that can be manipulated to act as an additional drag force through the vehicle body. MHD drag augmentation can benefit the spacecraft with an additional control mechanism that operates without moving parts. This can be utilized at higher altitudes and subsequently reduce convective heat transfer to the vehicle. MHD interaction with plasma flow suggests that both power generation and drag augmentation can be integrated together in one design. Modern numerical studies have only demonstrated the viability and potential of MHD energy generation and drag modulation separately for implementation into future spacecraft designs [1] [2] [3]. Steeves et. al. presented one possible generator design specifically for reentry vehicles that utilized a modular panel with two embedded electromagnets and two extruding electrodes for air plasma to flow between. Fujino et. al. presented results that indicated increased total drag with a permanent magnet embedded in the blunt body of the Orbital Reentry Experiment (OREX) trajectory. Furthermore, experimental investigations of these two MHD applications have been limited and separate due to the difficult technical nature of experimental de-sign and analysis. For MHD energy generation, there have been two designs which used a cylindrical fore-body with two embedded permanent magnets and two outward facing curved electrodes for tangency to an artificially ionized plasma that flowed along the body length [4] [5]. The first design was studied using an air microwave (MW) ionized supersonic plasma wind tunnel and the second design used radio frequency (RF) ionization of Argon, but both physically demonstrated the feasibility of MHD energy generation in reentry plasma conditions. For MHD drag modulation, one study involved a small-scale cylindrical blunt body embedded with permanent magnets immersed in plasma flow in an arcjet tunnel [6], and another utilized a spherical permanent magnet enclosed with a spherical forebody that was evaluated in a shock expansion tunnel [7]. Both experimental studies provided physical evidence of drag augmentation due to MHD interactions. Therefore, because of the multi-faceted effects of MHD along with both applications having been deemed feasible and beneficial during reentry, this indicates that a coupled power generation and drag modulation design could be achievable. The axi-symmetric nature of blunt bodies for reentry suggests that a spherical set of two ring electrodes near the nose along the body axis. Based on first principles, flow along the walls of the spherical body can then produce MHD interactions between the electrodes. As a result, two currents are generated: a Hall current that produces power and a primary inductive current that augments the drag force. The goal of this work is to provide a simulated design of a coaxial-electrode con-figuration for power generation and drag modulation for experimental testing.

E. Leong↗