Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STAGNATION PRESSURE”

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 631 records · Page 35

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗

In‐operando thermal transport characterization of moving particle bed heat exchanger

The project aimed to develop modulated photothermal radiometry (MPR) for conducting in-operando measurements of thermal transport in moving particle bed heat exchangers (MPBE) located at Sandia National Laboratory (SNL). The plan was to initially conduct lab-scale measurements on a 1-meter-long particle channel at UCSD, followed by in- operando measurements at SNL on a 20 kWth MPBE prototype built for this project, and if time permitting, potential final deployment on SNL’s 1 MWth G3P3 MPBE. The first two goals were successfully achieved but due to changes with commissioning timeline of the G3P3 MPBE, in-operando measurements on the G3P3 MPBE could not be conducted. The measurement campaigns have resulted in some promising results and the MPR tool has been handed over to trained SNL personnel for future use. MPR measurements on UCSD’s 1-meter-long particle channel showed a peculiar length dependence of the particle bed thermal conductivity and near-wall thermal resistance. The bed thermal conductivity increased, and the near-wall thermal resistance decreased, with increasing distance from the channel inlet, both eventually reaching a plateau. This was found to be due to hydrostatic pressure screening by the channel walls, known as the “Janssen effect”, and is the first experimental study on its role in granular heat transfer. A 20 kWth small-scale MPBE prototype, the first of its kind with 3 mm particle channels and 2 mm sCO2 plates, was designed and manufactured by Vacuum Process Engineering (VPE) for this project. This MPBE prototype was installed in SNL’s HEX test stand and the MPR tool was integrated into the test stand for in-operando measurements. Several modifications have been made to the MPR setup for facile installation and easy operation at SNL. Measurements with Carbo HSP40/70 particles were conducted with various particle and sCO2 temperatures and flow rates. Further, apart from regular MPBE operation, fault scenarios such as plate warping, particle stagnation, and choking were also detected and MPR showed excellent sensitivity to the detection of such conditions. Through the course of this project, MPR has evolved from a sensitive laboratory technique to a rugged industrial tool with various applications in high temperature measurements and diagnostics.

14 SOLAR ENERGY↗

Thermal Protection System to Enable Ice Giant Aerocapture Mission for Delivering Both an Orbiter and an In Situ Probe

The Ice Giants have been identified as high priority science destinations in the last Decadal Survey [1] and could benefit from aerocapture as the primary method for orbit insertion [2]. A mass-efficient aerocapture system will enable the delivery of an orbiter along with an atmospheric probe (for in situ measurements to anchor global data collected by the orbiter) and possibly a lander at Triton [3]. Aerocapture could be executed either using low L/D rigid aeroshell with lift modulation (LMA) [4] or using deployable aeroshell using drag modulation (DMA) [3]. Nearly two decades ago, a NASA-funded team performed Neptune-Triton aerocapture studies with a mid-L/D lifting configuration [5] for achieving orbit using LMA. This study showed aerocapture challenges. Due to very high peak entry conditions combined with very high heat-load, a suite of TPS materials was required and this suite was deemed problematic from a qualification perspective, due test facility limitations. In the past 20 years, progress made in GN&C for lift-guided entry missions such as MSL, Orion EFT1, Mars 2020 and the upcoming Artemis missions, and the expertise in blunt body aerodynamics at large scale (~ 5m) has led the EDL community to conclude that aerocapture is a “go do” engineering activity and most technologies are in hand to propose missions with aerocapture [6] [7]. Aerocapture using DMA, currently in development, is an option for Ice Giant Missions. While DMA is simpler in some sense, due to ballistic entry and no need for lift-guided maneuvering, it has challenges and it’s maturity is lower. LMA and DMA both require one or more ablative Thermal Protection System (TPS) materials for the rigid aeroshell element. The ablative TPS needs to be robust and mass efficient due to the high heat loads and size of the rigid aeroshell. Currently, there are capable ablative thermal protection materials, e.g., Heatshield for Extreme Entry Environments Technology (HEEET), 3-D woven Mid-Density Carbon- Phenolic (3MDCP), and PICA (Phenolic-Impregnated Carbon Ablator) that are mature, i.e., at TRL 6 or higher. NASA also invested in Conformal PICA that was matured to TRL 5. Our goal is to evaluate the applicability of high TRL TPS and consider other design options. We first establish bounding aerocapture trajectories for a wide range of arrival conditions and the associated aerothermal environment. Based on the environments, we then determine the predicted TPS mass for the aeroshell [4]. In this presentation, we will outline the process by which we establish bounding aerocapture trajectories for hyperbolic excess velocities ranging from 27 km/s to 35 km/s, which are shown to be a range of velocities that can reduce the trip time from ~14 years to 8 years. The above velocity range translates to ~12 km/s to ~24 km/s arrival velocities at the planetary entry interface [2]. The velocity reduction required to achieve orbit ranges between ~2.5 km/s to 9.5 km/s for both Neptune and Triton. Propulsive insertion alone, due to the amount of fuel required to achieve the required velocity reduction, limits the science returned [2]. We establish the bounding aerocapture trajectories for a low L/D (~ 0.4) configuration for three different ballistic coefficients. The ballistic coefficient range is determined from three different aeroshell diameters of 3m, 4m and 5m and with an entry system mass of 2200 kg. With the above range of design parameters, we then determine conservative/bounding estimates of aerothermal environments by using a combination of CFD simulations and stagnation point heating estimates [7]. This engineering approach allows us to first assess the TPS need vs. TPS capability and determine the applicability of existing TPS. Once an applicable suite of TPS is determined, the TPS thickness and mass are computed. We show that the TPS mass fraction can be as low as 5% to as high as 20%, depending on total trip time reduction and other design parameters for a range of TPS. This is a large range for TPS mass fraction. We show PICA and HEEET can indeed enable aerocapture missions, but the missions incur a mass penalty. TPS mass savings, can be further reduced with the use of conformal PICA. Advancing the development of Conformal PICA to make it robust across the entire aerothermal environment (peak heat-flux, pressure and shear) range will result in TPS mass fractions of < 10% for Ice Giant aerocapture missions such as the Neptune-Triton mission. Aerocapture allows for not only shortening the trip time but enables larger mass to be placed in orbit. Furthermore, probes deployed from orbit will benefit in reduced entry environments allowing for a lower risk TPS implementation as compared to mission designs where the probe is released prior to orbit insertion. One of the challenges for the Ice Giant community is to ensure mission designs that maximize science and allow flexibility in the placement of the entry probe. The traditional approach to release the probe ahead of the orbiter may not optimize returned science. In this presentation, we will make the case for mature TPS such as HEEET and PICA. While these materials can enable aerocapture missions, completing the development of conformal PICA and extending Conformal PICA to be more robust, will have significant impact to TPS mass efficiency and significantly enhance science return for future Gas- and Ice-Giant missions.

E Venkatapathy↗

High Enthalpy Studies of Capsule Heating in an Expansion Tunnel Facility

Measurements were made on an Orion heat shield model to demonstrate the capability of the new LENS-XX expansion tunnel facility to make high quality measurements of heat transfer distributions at flow velocities from 3 km/s (h(sub 0) = 5 MJ/kg) to 8.4 km/s (h(sub 0) = 36 MJ/kg). Thirty-nine heat transfer gauges, including both thin-film and thermocouple instruments, as well as four pressure gauges, and high-speed Schlieren were used to assess the aerothermal environment on the capsule heat shield. Only results from laminar boundary layer runs are reported. A major finding of this test series is that the high enthalpy, low-density flows displayed surface heating behavior that is observed to be consistent with some finite-rate recombination process occurring on the surface of the model. It is too early to speculate on the nature of the mechanism, but the response of the gages on the surface seems generally repeatable and consistent for a range of conditions. This result is an important milestone in developing and proving a capability to make measurements in a ground test environment and extrapolate them to flight for conditions with extreme non-equilibrium effects. Additionally, no significant, isolated stagnation point augmentation ("bump") was observed in the tests in this facility. Cases at higher Reynolds number seemed to show the greatest amount of overall increase in heating on the windward side of the model, which may in part be due to small-scale particulate.

Dufrene, Aaron↗

Photogrammetric and Global Positioning System Measurements of Active Pahoehoe Lava Lobe Emplacement on Kilauea, Hawaii

Basalt is the most common rock type on the surface of terrestrial bodies throughout the solar system and -- by total volume and areal coverage -- pahoehoe flows are the most abundant form of basaltic lava in subaerial and submarine environments on Earth. A detailed understanding of pahoehoe emplacement processes is necessary for developing accurate models of flow field development, assessing hazards associated with active lava flows, and interpreting the significance of lava flow morphology on Earth and other planetary bodies. Here, we examine the active emplacement of pahoehoe lobes along the margins of the Hook Flow from Pu'u 'O'o on Kilauea, Hawaii. Topographic data were acquired between 21 and 23 February 2006 using stereo-imaging and differential global positing system (DGPS) measurements. During this time, the average discharge rate for the Hook Flow was 0.01-0.05 cubic m/s. Using stereogrammetric point clouds and interpolated digital terrain models (DTMs), active flow fronts were digitized at 1 minute intervals. These areal spreading maps show that the lava lobe grew by a series of breakouts tha t broadly fit into two categories: narrow (0.2-0.6 m-wide) toes that grew preferentially down-slope, and broad (1.4-3.5 m-wide) breakouts that formed along the sides of the lobe, nearly perpendicular to the down-flow axis. These lobes inflated to half of their final thickness within approx 5 minutes, with a rate of inflation that generally deceased with time. Through a combination of down-slope and cross-slope breakouts, lobes developed a parabolic cross-sectional shape within tens of minutes. We also observed that while the average local discharge rate for the lobe was generally constant at 0.0064 +/- 0.0019 cubic m/s, there was a 2 to 6 fold increase in the areal coverage rate every 4.1 +/- 0.6 minutes. We attribute this periodicity to the time required for the dynamic pressurization of the liquid core of the lava lobe to exceed the cooling-induced strength of the lobe margins. Using DGPS-derived DTMs of the topography before and after pahoehoe lobe emplacement, we observed that the lava typically concentrated within existing topographic lows, with the lobe reaching a maximum thickness of approx 1.2 m above the lowest points of the initial topography and above reverse-facing slopes. Lobe margins were typically controlled by high-standing topography, with the zone directly adjacent to the final flow margin having average relief that is approx 4 cm higher than the lava-inundated region. This suggests that irregularities approx 25% of the height of the smallest breakout elements (i.e., toes) can exert a strong control on the paths of low-discharge pahoehoe lobes, with stagnated toes forming confining margins that allow interior portions of flow to topographically invert the landscape by inflation.

Hamilton, Christopher W.↗

Overview of Additively Manufactured TPS Proposed Flight Test and Earth Re-Entry Capsule Design

A flight mechanics overview is presented of an Earth flight test designed to investigate a novel, 3D printed thermal protection system (TPS) that is currently in development at NASA as part of the Additive Manufacturing of Thermal Protection Systems project. The project is pioneering a method to print a thermal protection system onto an entry vehicle forebody one layer at a time. This method reduces labor and complexity as compared to traditional manufacturing methods while increasing mission-dependent customization of through-depth materials properties. The flight test has three objectives. First, subject the forebody stagnation point of a capsule equipped with additively manufactured TPS (AMTPS) material to peak heat fluxes in excess of 100 $W/cm^2$. Second, capture in-flight data to enable flight reconstruction and AMTPS material thermal response model improvement. Third, recover the capsule with data storage and forebody AMTPS intact to enable post-flight inspection and analysis of AMTPS performance. The flight test trajectory is designed to achieve a peak stagnation point, cold-wall, entry heat flux of 135~$W/cm^2$. Flight mechanics simulations are performed using the Program to Optimize Simulated Trajectories II (POST2) and Monte-Carlo analysis yields statistical percentiles on vehicle performance at key points along the trajectory. Based on the flight mechanics analysis presented in this paper, a prototype capsule was designed, partially fabricated, and underwent preliminary component stress testing in preparation for fabrication of the flight unit capsule. The capsule outer mold line is a modified version of the heritage Mars Microprobe geometry. The capsule has a 0.356~m diameter, a 30~kg mass, and a hypersonic ballistic coefficient of 300~$kg/m^2$. Sensor selection is guided by flight dynamics simulations with the goal of resolving the re-entry heating pulse. On-board instrumentation include forebody and aftbody pressure sensors and thermocouples, a 9-axis IMU, a GPS receiver, and an Iridium satellite modem, all of which collect and store data throughout flight via on-board avionics systems. A two-stage parachute system is designed to decelerate the capsule to touchdown velocities that will not result in significant fracture or deformation of the charred AMTPS material at ground impact.

Flight Mechanics↗

FY24 Report on Water NSTF Testing: Lower Tank Inlet Piping Configuration

The following report serves as a summary of the accomplishments and testing results by the Natural convection Shutdown heat removal Test Facility (NSTF) experimental program over the past 12-month period. A major activity included reconfiguration of the chimney piping geometry which altered the discharge position into the tank from the original 50% tank height to a new lower position at a 10.9% tank height. This modification increases the volume of available coolant thus extending long-term operating capacity, however also results in a decreased liquid driving head which has the potential to reduce the natural circulation efficiency and decrease overall performance. Examination of the tradeoffs for these two configurations is an area of interest for RCCS designers and drove planned test activities this year. Ten matrix tests were performed in FY24, totaling 195 hours of heated operations and 12.9 MWh of electrical heating throughout the year, with eight classified as Accepted per NQA-1 and two classified as Trending. Testing prior to the facility reconfiguration examined the facility response to throttling at the tank inlet, demonstrating increased sensitivity of flow instabilities to throttling within the two-phase region when compared to throttling at the lower sensitivity single-phase inlet region. Testing after the chimney reconfiguration began with a baseline test at conditions of 80% inventory fill and prototypic thermal power input of 2.1 MW t . In addition to establishing a reference for nominal system behavior and performance, repeat testing with multiple subsequent runs demonstrated strong repeatability between tests performed at the same conditions in this new configuration. These tests also examined if influences would occur to system behavior with installation of new higher-resolution instrumentation within the upper chimney. In this critical region where boiling and flashing phenomena dominate, the test results provided confidence that the new instrumentation does not uncharacteristically influence the observed behaviors. An inventory parametric series was then initiated to examine system behavior with the new chimney configuration at six varying initial inventory levels, ranging from high 80% to low 20% fill. Two-phase oscillations with similar peak and mean flow rates were observed when comparing to the previous mid tank configuration. Generally, similar system response trends with inventory were also observed, such as two phase oscillations suppressing as inventories were lowered. However, in one absolute inventory comparison at the highest fill of 80%, flow oscillations saw a gradual growth over the 4 hours of two-phase operation in the new lower tank inlet configuration, opposite of the gradual dampening observed in the mid tank inlet configuration. This can be attributed to, in part, a greater hydrostatic head pressure above the two-phase discharge region where the boiling front is developed. Furthermore, the change resulted in greater sensitivity to liquid degassing phenomena during single-phase heating, causing loop instabilities to form which trigged moderate flow degradation during the period approaching saturation and boiling conditions. This behavior had been observed previously under some conditions but was a common occurrence in recent testing with the newer lower tank inlet configuration. Initial observations from these first data sets suggest an overall larger window of stability for the mid tank inlet configuration compared to the lower tank inlet. Lastly, the lowest inventory fill test was repeated over an extended testing window to examine depletion behavior. Natural circulation flow and effective heat removal performance were observed during most of the testing period; only after the tank became fully drained (0% fill) did flow stagnate and violent geysering events occur. This early observation confirms one relative advantage over the mid tank inlet configuration, which stagnated under comparable conditions at ~20% inventory remaining in the tank.

42 ENGINEERING↗

CFD Simulations of the IHF Arc-Jet Flow: 9-Inch Nozzle, Flow Surveys, LEAF Wedge Calibration Data

This paper reports computational analyses and flow characterization studies in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using a wedge model placed in a free jet downstream of new 9-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. Both the nozzle and wedge model were specifically designed for testing in the new Laser-Enhanced Arc-jet Facility. Data were obtained using stagnation calorimeters and wedge models placed downstream of the nozzle exit. Two instrumented wedge calibration plates were used: one water-cooled and the other RCG-coated tile plate. Experimental surveys of arc-jet test flow with pitot and heat flux probes were also performed at three arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. By taking into account nonuniform total enthalpy and mass flux profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields, these simulations approximately reproduce the probe survey data and predict the wedge model surface pressure and heat flux measurements.

Arc-jets↗

Exit plane H2O concentration measurements correlated with OH PLIF near-injector mixing measurements for scramjet flows

Mixing and combusting high enthalpy flows, similar to those encountered in scramjet engines, were investigated using a shock tunnel to produce the flow in conjunction with non-intrusive optical diagnostics which monitored the performance of two injector configurations. The shock tunnel is configured to produce Mach 3 flow and stagnation enthalpies corresponding to flight equivalent Mach numbers between 7 and 11. A pulsed hydrogen injection capability and interchangeable injector blocks provide a means of examining high speed, high enthalpy reacting flows. Planar laser induced fluorescence (PLIF) of OH molecules in the near injector region produced images which show the combusting and mixing zones for the reacting flow. Line-of-sight exit plane measurement of water concentration and temperature were used to provide a unique method of monitoring exit plane products. These results demonstrated that a velocity matched axial injection system produced a fuel jet that lifted off the floor of the duct. Mixing was observed to increase for this system as a velocity mismatch was introduced. Comparison of exit plane water concentrations for a wall jet injection system and a velocity matched injection system indicated similar mixing performance but an accurate pressure measurement is necessary to further validate the result. In addition, exit plane measurements indicated an approximate steady-state condition was achieved during the 1 to 2 ms test times.

Parker, T. E.↗

Prediction of Unsteady Transitional Layers in Turbomachinery Using Navier Stokes Equations

The objective of the research reported in this presentation is to develop computational techniques for the prediction of unsteady transitional flows associated with the rotor stator interaction in turbomachinery. Three low-Reynolds number turbulence models are incorporated in two unsteady Navier-Stokes codes (one is pressure based and the other is time marching with Runge-Kutta time stepping) and evaluated for accuracy in predicting the onset and the end of unsteady transitional patches due to wake passing. The best model is then used for modification and improvement for the leading edge effect. An existing steady Navier-Stokes code was modified to include pseudo-time stepping, which provided acceleration from 5 to 25 times that of the original code. A systematic validation procedure was implemented to assess the effects of the grid, artificial dissipation, physical, and the pseudo-time step for an accurate prediction of transitional flows resulting from the rotor-stator interaction. The ability of the Navier-Stokes code to predict the unsteady transitional flow on a turbomachinery blade is demonstrated. The unsteady pressure and velocity fields are in good agreement with the experimental data and the prediction from the Euler/boundary layer approach. The numerical solver was able to capture all zones (wake induced transitional strip, wake induced turbulent strip, calmed region, etc.) associated with wake induced transition in a compressor cascade. Another significant step is the assessment of k-epsilon turbulence models, including the leading edge modifications. Best results were obtained from the FLB model. The LB model predicted earlier inception of the transition and shorter transition length. Modification of the k-epsilon model was found to be essential for an accurate prediction of the unsteady transitional flow in a compressor cascade. The CH model failed to predict the unsteady transitional flow. Predicted boundary layer was turbulent from the leading edge, even with the modification of the k-epsilon model near the stagnation point. A comparison of the instantaneous shape factor, the skin friction coefficient, and the momentum thickness indicates that the Navier-Stokes predictions are reasonably good. Interaction between the upstream wake and stator wake results in shedding of unsteady vortices from the trailing edge and increased dissipation in the stator wake and, as a consequence, increased rate of decay of the stator wake. The procedure developed at Penn State should aid the designers of turbomachinery in allowing for the unsteady transitional flows due to rotor-stator interaction.

Lakshminarayana, B.↗

Heat Transfer Measurements and Predictions on a Power Generation Gas Turbine Blade

Detailed heat transfer measurements and predictions are given for a power generation turbine rotor with 129 deg of nominal turning and an axial chord of 137 mm. Data were obtained for a set of four exit Reynolds numbers comprised of the design point of 628,000, -20%, +20%, and +40%. Three ideal exit pressure ratios were examined including the design point of 1.378, -10%, and +10%. Inlet incidence angles of 0 deg and +/-2 deg were also examined. Measurements were made in a linear cascade with highly three-dimensional blade passage flows that resulted from the high flow turning and thick inlet boundary layers. Inlet turbulence was generated with a blown square bar grid. The purpose of the work is the extension of three-dimensional predictive modeling capability for airfoil external heat transfer to engine specific conditions including blade shape, Reynolds numbers, and Mach numbers. Data were obtained by a steady-state technique using a thin-foil heater wrapped around a low thermal conductivity blade. Surface temperatures were measured using calibrated liquid crystals. The results show the effects of strong secondary vortical flows, laminar-to-turbulent transition, and also show good detail in the stagnation region.

Giel, Paul W.↗

Viscous shock-layer solutions for the low-density hypersonic flow past long slender bodies

Results are obtained for the surface pressure, drag, heat-transfer, and skin-friction coefficients for hyperboloids and sphere cones. Body half angles from 5 to 22.5 degrees are considered for various low-density flow conditions. Recently obtained surface-slip and shock-slip equations are employed to account for the low-density effects. The method of solution employed for the viscous shock-layer (VSL) equations is a partially coupled spatial-marching implicit finite-difference technique. The flow cases analyzed include highly cooled long slender bodies in high Mach number flows. The present perfect-gas VSL calculations compare quite well with available experimental data. Results have also been obtained from the steady-state Navier-Stokes (NS) equations by successive approximations. Comparison between the NS and VSL results indicates that VSL equations even with body and shock-slip boundary conditions may not be adequate in the stagnation region at altitudes greater than about 75 km for the cases analyzed here.

Gupta, R. N.↗

Ablation of ATJ Graphite at High Temperatures

Results of an extensive experimental program on the ablation of ATJ graphite in air at surface pressures of 0.3 to 4.4 atm and surface temperatures of 2570 to 4030 K are presented. The measured mass-loss rates are correlated with surface temperature, pressure, and effective nose radius. The results are compared with several equilibrium thermochemical ablation theories which differ basically in the chemical species that are considered and in the thermodynamic properties assigned to these species. All the theories predict about the same mass-loss rate in the diffusion-controlled oxidation regime and are in good agreement with the experimental results. At higher temperatures, however, the experimental and theoretical results do not agree. At temperatures above 3700 K, the experimental mass-loss rate becomes independent of pressure and an exponential function of temperature. A comparison of the high-temperature results with the theories indicates that all the theories, except the one based on the thermodynamic properties from the JANAF tables, are invalid because they overpredict the mass-loss rate. On the other hand, the theory based on the JANAF properties underpredicts the measured mass-loss rate by a factor of three at 4000 K. At least a portion of this difference between experiment and theory is attributed to particulate mass loss which is visually observed.

Thermochemical ablation↗

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗

SEPAC data analysis in support of the environmental interaction program

Data analyses of the Space Experiments with Particle Accelerators (SEPAC) data and computer modeling were conducted to investigate spacecraft environmental effects associated with injection of electron beams, plasma clouds, and neutral gas clouds from the Shuttle orbiter. The data analysis indicates that Extremely Low Frequency oscillations from 150 to 200 Hz were seen in the Langmuir probe current when the beam was fired in a continuous mode. The strongest oscillations occurred when the ambient pressure was augmented by neutral gas releases from the SEPAC plasma accelerator magnetoplasma-dynamic (MPD) arcjet. To understand the dependence of spacecraft charging potential on beam density and other plasma parameters, a two-dimensional electrostatic particle code was used to simulate the injection of electron beams from an infinite conductor into a plasma. The simulations show that the conductor charging potential depends critically on the reflection coefficient of the conductor surface, which is defined as the percentage of incident particles reflected by the conductor. The ionization effects on spacecraft charging were examined by including interactions of electrons with neutral gas. The simulations show that the conductor charging potential decreases with increasing neutral background density due to the production of secondary electrons near the conductor surface. The simulations also indicate that the beam radius is generally proportional to the beam electron gyroradius when the conductor is charged to a large potential. It appears that the charge buildup at the beam stagnation point causes the beam radial expansion. A survey of the simulation results suggests that the ratio of the beam radius to the beam electron gyroradius increases with the square root of beam density and decreases inversely with beam injection velocity. These results are useful for explaining the spacecraft charging phenomena observed during SEPAC experiments from Spacelab 1.

Lin, Chin S.↗

Computational and Physics-Based Modeling for the Development of in-Space Welding Technology

Welding will be an enabling technology for in-space assembly and manufacturing (ISAM) and in-space repair operations in support of the growing Space Economy. Prospective use cases of welding include metal repair on long-duration lunar surface missions and assembly of large structures that exceed dimensions of the currently available launch vehicle envelopes such as space arrays, tall lunar transmission towers, or even large transit vehicles powered by nuclear propulsion. There is a brief history of beam welding for in-space use. Soviet and NASA astronauts performed electron beam welds (EBW) in space in the late-1960s through mid-1980s and several parabolic flights of laser beam welding (LBW) were undertaken in the USA and Canada in the early-1990s. Work has since stagnated until the recent boom of the in-space economy has demonstrated need for continued development of EBW and LBW processes, which are desirable for their compactness, low mass, low power consumption, and high wall plug efficiency. No space-based manufacturing laboratory capable of demonstrating, developing, and applying welding processes exists, and costs are high to develop and fly individual space-based trials. This drives the current approach to develop welding technology with computational and physics-based models in tandem with earth-based experiments. Several approaches are taken to complement and drive the selection of instrumentation for planned parabolic and suborbital flights. First, a thermal model of Skylab EBW experiments is established and compared to empirical data. Next, the model is used to show variation between Earth-based and in-space weld structure and properties through exploration of key variables that define the Low earth orbit and lunar surface space environments, namely temperature (40 – 400 K), pressure (10-10 to 104 Pa), and gravity (10-6 g to 1 g). Finally, physical processes associated with heat and mass transport during welding are evaluated through several case studies across the above variable ranges demonstrating need for further research.

in-space welding↗

Mars Science Laboratory Heatshield Aerothermodynamics: Design and Reconstruction

The Mars Science Laboratory heatshield was designed to withstand a fully turbulent heat pulse based on test results and computational analysis on a pre-flight design trajectory. Instrumentation on the flight heatshield measured in-depth temperatures in the thermal protection system. The data indicate that boundary layer transition occurred at 5 of 7 thermocouple locations prior to peak heating. Data oscillations at 3 pressure measurement locations may also indicate transition. This paper presents the heatshield temperature and pressure data, possible explanations for the timing of boundary layer transition, and a qualitative comparison of reconstructed and computational heating on the as-flown trajectory. Boundary layer Reynolds numbers that are typically used to predict transition are compared to observed transition at various heatshield locations. A uniform smooth-wall transition Reynolds number does not explain the timing of boundary layer transition observed during flight. A roughness-based Reynolds number supports the possibility of transition due to discrete or distributed roughness elements on the heatshield. However, the distributed roughness height would have needed to be larger than the pre-flight assumption. The instrumentation confirmed the predicted location of maximum turbulent heat flux near the leeside shoulder. The reconstructed heat flux at that location is bounded by smooth-wall turbulent calculations on the reconstructed trajectory, indicating that augmentation due to surface roughness probably did not occur. Turbulent heating on the downstream side of the heatshield nose exceeded smooth-wall computations, indicating that roughness may have augmented heating. The stagnation region also experienced heating that exceeded computational levels, but shock layer radiation does not fully explain the differences.

Edquist, Karl T.↗

Welding in Space: Past, Present, and Future

It has been over fifty years since the first welds were made in space by Soviet cosmonauts on Soyuz-6 in October of 1969. The United States performed bead-on-plate welding, brazing, and metal melting experiments onboard the Skylab orbital space station several years later in 1973. Finally, Soviet cosmonauts departed their Salyut-7 capsule and made the first (and last) welds in open space in 1984. Progress on further demonstrations of welding in space stagnated, and subsequent microgravity welding research work shifted to lower-cost earth-based experiments that include drop towers and parabolic flights. With the advent of the International Space Station, relevant microgravity research was undertaken in the field of metal solidification science, and limited studies of brazing and soldering were undertaken; however, no welds have ever been on the Space Station. Therefore, it has been almost 40 years since a weld was made in space. Much like terrestrial construction and manufacturing industries, welding, joining, and allied processes will be enabling technologies for In-space Assembly and Manufacturing (ISAM) in the nascent “Space Economy”, a sector that is expected to approach a value of 10^12 USD within the next decade. It is critical that the welding and joining research community, along with the welding industry, engages the space industry to advance the understanding of those critical manufacturing processes which must be evaluated and matured in the extreme environments of space. Such environments include variable gravity (microgravity in low earth orbit, 0.17 gravity on the lunar surface, and 0.38 gravity on the martian surface), reduced pressure (extreme vacuum in space and the lunar surface to a predominantly CO2 atmosphere on Mars), and extreme temperatures (between 40 and 400 K). Past experiments and analyses are reviewed to suggest requirements that the welding and joining community should target to make strides on closing the current space welding gaps. Efforts to continue welding in space are currently underway at NASA in conjunction with academia and industry. Those are explored along with future agency goals which are seen as opportunities to engage the welding community on this historic effort.

in-space assembly and manufacturing↗