Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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.

1,075 records · Page 60

Heating of the solar middle chromosphere by large-scale electric currents

A global resistive, two-dimensional, time-dependent magnetohydrodynamic (MHD) model is used to introduce and support the hypothesis that the quiet solar middle chromosphere is heated by resistive dissipation of large-scale electric currents which fill most of its volume. The scale height and maximum magnitude of the current density are 400 km and 31.3 m/sq m, respectively. The associated magnetic field is almost horizontal, has the same scale height as the current density, and has a maximum magnitude of 153 G. The current is carried by electrons flowing across magnetic field lines at 1 m/s. The resistivity is the electron contribution to the Pedersen resitivity for a weakly ionized, strongly magnetized, hydrogen gas. The model does not include a driving mechanism. Most of the physical quantities in the model decrease exponentially with time on a resistive timescale of 41.3 minutes. However, the initial values and spatial; dependence of these quantities are expected to be essentially the same as they would be if the correct driving mechanism were included in a more general model. The heating rate per unit mass is found to be 4.5 x 10(exp 9) ergs/g/s, independent of height and latitude. The electron density scale height is found to be 800 km. The model predicts that 90% of the thermal energy required to heat the middle chromosphere is deposited in the height range 300-760 km above the temperature minimum. It is shown to be consistent to assume that the radiation rate per unit volume is proportional to the magnetic energy density, and then it follows that the heating rate per unit volume is also proportional to the energy from the photosphere into the overlying chromosphere are briefly discussed as possible driving mechanisms for establishing and maintaining the current system. The case in which part of or all of the current is carried by protons and metal ions, and the contribution of electron-proton scattering to the current are also considered, with the conclusion that these effects do not change the qualitative prediction of the model, but probably change the quantitative predictions slightly, mainly by increasing the maximum magntiude of the current density and magnetic field to at most approximately 100 mA/m and approximately 484 G, respectively. The heating rate per unit mass, current density scale height, magnetic field scale height, temperatures, and pressures are unchanged or are only slightly changed by including these additional effects due to protons and ions.

Goodman, M. L.↗

Comparison of Classical and Charge Storage Methods for Determining Conductivity of Thin Film Insulators

Conductivity of insulating materials is a key parameter to determine how accumulated charge will distribute across the spacecraft and how rapidly charge imbalance will dissipate. Classical ASTM and IEC methods to measure thin film insulator conductivity apply a constant voltage to two electrodes around the sample and measure the resulting current for tens of minutes. However, conductivity is more appropriately measured for spacecraft charging applications as the "decay" of charge deposited on the surface of an insulator. Charge decay methods expose one side of the insulator in vacuum to sequences of charged particles, light, and plasma, with a metal electrode attached to the other side of the insulator. Data are obtained by capacitive coupling to measure both the resulting voltage on the open surface and emission of electrons from the exposed surface, as well monitoring currents to the electrode. Instrumentation for both classical and charge storage decay methods has been developed and tested at Jet Propulsion Laboratory (JPL) and at Utah State University (USU). Details of the apparatus, test methods and data analysis are given here. The JPL charge storage decay chamber is a first-generation instrument, designed to make detailed measurements on only three to five samples at a time. Because samples must typically be tested for over a month, a second-generation high sample throughput charge storage decay chamber was developed at USU with the capability of testing up to 32 samples simultaneously. Details are provided about the instrumentation to measure surface charge and current; for charge deposition apparatus and control; the sample holders to properly isolate the mounted samples; the sample carousel to rotate samples into place; the control of the sample environment including sample vacuum, ambient gas, and sample temperature; and the computer control and data acquisition systems. Measurements are compared here for a number of thin film insulators using both methods at both facilities. We have found that conductivity determined from charge storage decay methods is 102 to 104 larger than values obtained from classical methods. Another Spacecraft Charging Conference presentation describes more extensive measurements made with these apparatus. This work is supported through funding from the NASA Space Environments and Effects Program and the USU Space Dynamics Laboratory Enabling Technologies Program.

Swaminathan, Prasanna↗

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗

Decision Aid for Conjunction Risk Mitigation by Differential Drag

In the previous five years, the rate of conjunctions that the NASA Conjunction Assessment Risk Analysis (CARA) team processed and analyzed has more than tripled. (NASA CARA, 2024) New missions in the early development phases are now required to plan for dealing with conjunctions under the present space environment, and also projecting forward into a future likely with even further increased utilization of the space environment. Some missions are investigating the possibility of using differential drag to remediate conjunctions without expending limited fuel or for missions without propulsive capabilities. The NASA CARA team studied the historical record of conjunctions to evaluate the circumstances under which differential drag may be successfully applied and have developed a series of tables to use as a decision aid for missions considering differential drag. CARA records all conjunctions of their protected payloads, with historical records starting in 2005 (with significant conjunction events starting to occur on or after 2013). From this record, approximately 7,300 had a probability of collision (Pc) greater than 1 in 10,000, the nominal requirement for a risk mitigation maneuver (RMM) to be made per NASA Procedural Requirements (NPR) 8079.1 (NASA, 2023), at 3 days prior to the time of closest approach, the time analyzed for differential drag efficacy. CARA’s maneuver trade space tool (MTS) was used to prop-agate the primary satellite forward from that decision point with varying degrees of increase to ballistic coefficient (BC), and then recalculate the Pc to evaluate whether or not the conjunction was mitigated (Pc < 3E-6). These results were then binned and sorted along several axes, including altitude, amount of delta-BC, and (pre-maneuver) rate of energy dissipation, to identify underlying patterns. The altitude plot is shown in Figure 1. We found that differential drag was most successful for satellites with perigees below 560 km, and which could adopt an average delta-BC of 2 or greater (that is, increasing their ballistic coefficient by a factor of 3). However, this is a difficult threshold for a mission to clear; very few spacecraft are capable of adopting a high-drag configuration for 72 hours continuously. Planet’s Dove spacecraft use differential drag to remediate conjunctions (Griffith, et al., 2021), and they have a maximum delta-BC factor of 9, but in practice (with mission and charging constraints) they achieve a time-averaged delta-BC that is closer to 2 (Foster, et al., 2017). In contrast, a NASA mission in development reached out to CARA to evaluate the utility of differential drag to remediate potential conjunctions for a spacecraft which is capable of a similar high delta-BC, but with operational constraints that limit their time-averaged delta-BC to less than 1. CARA has developed tables that can be used as decision aids to advise missions-in-development about the best way to utilize their differential drag capabilities. For missions below 560 km with the operational flexibility to devote multiple days to holding a high-drag configuration (or a sufficiently high drag ratio to compensate for limitations on that time), they are able to successfully remediate high-risk conjunctions. Conversely, missions that do not meet these ex-acting criteria – most missions – can instead be advised to use on-board propulsion systems to perform RMMs or to turn their minimum-area face towards the approach vector, thereby reducing Pc at the moment of conjunction. (NASA, 2023)

risk mitigation↗

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

Large Scale Flame Spread Environmental Characterization Testing

Under the Advanced Exploration Systems (AES) Spacecraft Fire Safety Demonstration Project (SFSDP), as a risk mitigation activity in support of the development of a large-scale fire demonstration experiment in microgravity, flame-spread tests were conducted in normal gravity on thin, cellulose-based fuels in a sealed chamber. The primary objective of the tests was to measure pressure rise in a chamber as sample material, burning direction (upward/downward), total heat release, heat release rate, and heat loss mechanisms were varied between tests. A Design of Experiments (DOE) method was imposed to produce an array of tests from a fixed set of constraints and a coupled response model was developed. Supplementary tests were run without experimental design to additionally vary select parameters such as initial chamber pressure. The starting chamber pressure for each test was set below atmospheric to prevent chamber overpressure. Bottom ignition, or upward propagating burns, produced rapid acceleratory turbulent flame spread. Pressure rise in the chamber increases as the amount of fuel burned increases mainly because of the larger amount of heat generation and, to a much smaller extent, due to the increase in gaseous number of moles. Top ignition, or downward propagating burns, produced a steady flame spread with a very small flat flame across the burning edge. Steady-state pressure is achieved during downward flame spread as the pressure rises and plateaus. This indicates that the heat generation by the flame matches the heat loss to surroundings during the longer, slower downward burns. One heat loss mechanism included mounting a heat exchanger directly above the burning sample in the path of the plume to act as a heat sink and more efficiently dissipate the heat due to the combustion event. This proved an effective means for chamber overpressure mitigation for those tests producing the most total heat release and thusly was determined to be a feasible mitigation strategy to incorporate into the microgravity experiment.

Clayman, Lauren K.↗

The 3-D description of vertical current sheets with application to solar flares

Following a brief review of the processes which have been suggested for explaining the occurrence of solar flares we suggest a new scenario which builds on the achievements of the previous suggestion that the current sheets, which develop naturally in 3-D cases with gravity from impacting independent magnetic structures (i.e., approaching current systems), do not consist of horizontal currents but are instead predominantly vertical current systems. This suggestion is based on the fact that as the subphotospheric sources of the magnetic field displace the upper photosphere and lower chromosphere regions, where plasma beta is near unity, will experience predominantly horizontal mass motions which will lead to a distorted 3-D configurations of the magnetic field having stored free energy. In our scenario, a vertically flowing current sheet separates the plasma regions associated with either of the subphotospheric sources. This reflects the balanced tension of the two stressed fields which twist around each other. This leads naturally to a metastable or unstable situation as the twisted field emerges into a low beta region where vertical motions are not inhibited by gravity. In our flare scenario the impulsive energy release occurs, initially, not by reconnection but mainly by the rapid change of the magnetic field which has become unstable. During the impulsive phase the field lines contort in such way as to realign the electric current sheet into a minimum energy horizontal flow. This contortion produces very large electric fields which will accelerate particles. As the current evolves to a horizontal configuration the magnetic field expands vertically, which can be accompanied by eruptions of material. The instability of a horizontal current is well known and causes the magnetic field to undergo a rapid outward expansion. In our scenario, fast reconnection is not necessary to trigger the flare, however, slow reconnection would occur continuously in the current layer at the locations of potential flaring. During the initial rearrangement of the field strong plasma turbulence develops. Following the impulsive phase, the final current sheet will experience faster reconnection which we believe responsible for the gradual phase of the flare. The reconnection will dissipate part of the current and will produce sustained and extended heating in the flare region and in the postflare loops.

Fontenla, Juan M.↗

Focal Gray Matter Plasticity as a Function of Long Duration Head-down Tilt Bed Rest

Long duration spaceflight (i.e., > or = 22 days) has been associated with changes in sensorimotor systems, resulting in difficulties that astronauts experience with posture control, locomotion, and manual control. The microgravity environment is an important causal factor for spaceflight induced sensorimotor changes. Whether these sensorimotor changes may be related to structural and functional brain changes is yet unknown. However, experimental studies revealed changes in the gray matter (GM) of the brain after simulated microgravity. Thus, it is possible that spaceflight may affect brain structure and thereby cognitive functioning and motor behavior. Long duration head-down tilt bed rest has been suggested as an exclusionary analog to study microgravity effects on the sensorimotor system. Bed rest mimics microgravity in body unloading and bodily fluid shifts. In consideration of the health and performance of crewmembers both in- and post-flight, we are conducting a prospective longitudinal 70-day bed rest study as an analog to investigate the effects of microgravity on the brain. VBM analysis revealed a progressive decrease from pre- to in- bed rest in GM volume in bilateral areas including the frontal medial cortex, the insular cortex and the caudate. Over the same time period, there was a progressive increase in GM volume in the cerebellum, occipital-, and parietal cortex, including the precuneus. The majority of these changes did not fully recover during the post-bed rest period. Analysis of lobular GM volumes obtained with BRAINS showed significantly increased volume from pre-bed rest to in-bed rest in GM of the parietal lobe and the third ventricle. Temporal GM volume at 70 days in bed rest was smaller than that at the first pre-bed rest measurement. Trend analysis showed significant positive linear and negative quadratic relationships between parietal GM and time, a positive linear relationship between third ventricle volume and time, and a negative linear relationship between cerebellar GM volume and time. FM performance improved from pre-bed rest session 1 to session 2. From the second pre-bed rest measure to the last-day-in-bed rest, there was a significant decrease in performance that only partially recovered post-bed rest. No significant association was observed between changes in brain volume and changes in functional mobility. Extended bed rest, which is an analog for microgravity, can result in local volumetric GM increase and decrease and adversely affect functional mobility. These changes in brain structure and performance were not related in this sample. Whether the effects of bed rest dissipate at longer times post-bed rest, and if they are associated with behavior are important questions that warrant further research including more subjects and longer follow-up times.

Koppelmans, V.↗

The Impact of Drifting Orbits on the Monthly Regional TOA Flux Assuming Constant Meteorology

The NASA Clouds and the Earth's Radiant Energy System (CERES) gridded Single Scanner Footprint (SSF1deg) product provides TOA SW and LW monthly 1° regional all-sky fluxes, which are used to monitor the Earth’s energy balance. The CERES long-term climate data record relies on Terra and Aqua satellite sun-synchronous orbits that are maintained at 10:30 and 1:30 local equator crossing times (LECT), respectively. The Terra and Aqua satellites are expected to drift outside of their respective LECT during mid 2022. Both Terra and Aqua will drift over several years towards sunrise and sunset, respectively, and eventually will be deorbited. The CERES SSF1deg product monthly regional fluxes are based on the well calibrated and stable CERES observed fluxes and are temporally interpolated assuming constant meteorology between measurements to resolve the regional diurnal flux cycle to obtain a daily averaged flux. The drifting orbits may impact the monthly regional TOA flux over regions with systematic diurnal cycles, because the observations will shift in local time. The CERES project would like to determine the maximum Terra and Aqua LECT time shift before the monthly regional fluxes are diurnally impacted and become unreliable for long-term climate monitoring. To determine the impact of the drifting orbits on the SSF1deg monthly regional fluxes,15-minute Geostationary Earth Radiation Budget (GERB) broadband observed fluxes over the Meteosat geostationary satellite domain (±60° in longitude and latitude) are used as proxy CERES observations. The drifting orbit sampling pattern is achieved by simply incrementing the observation time by steps of 15 minutes from the CERES footprint time. For each 15-minute time interval, the CERES observed fluxes are replaced by the GERB observed fluxes. The-15 minute incremented monthly regional fluxes based on constant meteorology are compared to the reference 10:30 and 1:30 LECT fluxes. Regions with systematic diurnal cycles, include morning maritime stratus, where the clouds dissipate during the morning, and land afternoon convection, where clouds increase in the afternoon will impact the regional flux differences. Based on January and July 2010 GERB data, even a 15-minute LECT change caused regional monthly flux differences that would impact long term regional trend analysis. Results will be shown at the conference.

David Doelling↗

The Orbital Drift Impact on the Monthly Regional TOA Flux Assuming Constant Meteorology

The NASA Clouds and the Earth's Radiant Energy System (CERES) gridded Single Scanner Footprint (SSF1deg) product provides TOA SW and LW monthly 1° regional all-sky fluxes, which are used to monitor the Earth’s energy balance. The CERES long-term climate data record relies on Terra and Aqua satellite sun-synchronous orbits that are maintained at 10:30 and 1:30 local equator crossing times (LECT), respectively. The Terra and Aqua satellites are expected to drift outside of their respective LECT during mid 2022. Both Terra and Aqua will drift over several years towards sunrise and sunset, respectively, and eventually will be deorbited. The CERES SSF1deg product monthly regional fluxes are based on the well calibrated and stable CERES observed fluxes and are temporally interpolated assuming constant meteorology between measurements to resolve the regional diurnal flux cycle to obtain a daily averaged flux. The drifting orbits may impact the monthly regional TOA flux over regions with systematic diurnal cycles, because the observations will shift in local time. The CERES project would like to determine the maximum Terra and Aqua LECT time shift before the monthly regional fluxes are diurnally impacted and become unreliable for long-term climate monitoring. To determine the impact of the drifting orbits on the SSF1deg monthly regional fluxes,15-minute Geostationary Earth Radiation Budget (GERB) broadband observed fluxes over the Meteosat geostationary satellite domain (±60° in longitude and latitude) are used as proxy CERES observations. The drifting orbit sampling pattern is achieved by simply incrementing the observation time by steps of 15 minutes from the CERES footprint time. For each 15-minute time interval, the CERES observed fluxes are replaced by the GERB observed fluxes. The-15 minute incremented monthly regional fluxes based on constant meteorology are compared to the reference 10:30 and 1:30 LECT fluxes. Regions with systematic diurnal cycles, include morning maritime stratus, where the clouds dissipate during the morning, and land afternoon convection, where clouds increase in the afternoon will impact the regional flux differences. Based on January and July 2010 GERB data, even a 15-minute LECT change caused regional monthly flux differences that would impact long term regional trend analysis. Results will be shown at the conference.

D R Doelling↗

Signatures of bulge triaxiality from kinematics in Baade's window

We study a sample of 62 Baade's Window, (l,b) = (1, -4)deg, K giants that have published proper motions, radial velocity, and metallicity. Using R(sub 0) = 8 kpc, we construct the velocity ellipsoids, namely the 3x3 velocity dispersion tensors, for the metal rich stars ((Fe/H) greater than or equal to 0) and metal poor stars ((Fe/H) less than or equal to -0.2). After diagonalizing the tensor, we find a vertex deviation characteristic of a nonaxisymmetric system. Eigenvalues for the two velocity ellipsoids (sigma(sub 1), sigma(sub 2), sigma(sub 3)) are (126, 89, 65) +/- 13 km/s for the metal rich sample and (154, 77, 83) +/- 25 km/s for the metal poor sample with their long axes pointing to two nearly perpendicular directions (l(sub v), b(sub v)) = (-65 +/- 9 deg, +14 +/- 9 deg) and (l(sub v), b(sub v)) = (25 +/- 14 deg, -11 +/- 14 deg), respectively. The vertex deviations of the velocity ellipsoids cannot be consistently explained by any oblate model. We are able to reject the hypothesis that the metal poor and metal rich populations are drawn from the same distribution at better than the 97% confidence level. We populate orbits in a realistic bar potential with a Gaussian velocity distribution, allowing us to simulate and interpret observations. We conclude that the data are consistent with a triaxial bulge pointing towards (l,b) with l less than 0 deg and b = 0 deg as suggested by earlier work on gas dynamics and the observed light distribution. We also predict that low latitude (absolute value of b less than or equal to 4 deg) bulge fields should show the vertex deviation more strongly and would therefore be the best locations for future proper motion studies. In the classification scheme of Athanassoula et al. (1983) the metal rich stars appear to occupy the B-family orbits which rotate in the prograde sense in the rest frame and have boxy shapes that are aligned with and supporting the bar. The metal poor stars in the sample lag behind the metal rich bulge and appear to occupy R-family orbits which rotate in the retrograde sense in the rest frame. They have nearly round loop shapes and are aligned perpendicularly to the bar, hence limit the triaxiality of the bar potential. The correlations between the metallicity and the orbit families can develop if the bulge forms dissipatively on a sufficiently long time scale. However, it is difficult to explain such correlations if most stars in the inner Galaxy form during the violent relaxation phase.

Zhao, Hongsheng↗

Preliminary Assessment of the Impact on the V-Band Oxygen Channels From Satellite Communication Uplinks

We calculate the percentage of time that an ATMS-like instrument [1] will be illuminated by the uplink beam of one of the proposed V-band communication system and estimate the damage resulting from such exposure. Using a combination of openly available information and educated guesses about the location and characteristics of the up/down link terminals, we constructed the ground segment of a hypothetical high-speed communication network. The space segment of the network was constructed from the orbital data of the existing Starlink constellation [2] of 6223 communication satellites (comsats) which is used as strawman to represent any other possible constellation of communication satellites. It is shown that without a very delicate balance of frequency allocations (science vs telecommunications), coupled with extremely steep and deep bandpass-defining filters, and strict adherence to the agreed limits (i.e. no out-of-band transmissions) the deployment of the telecommunication network leads to almost-complete loss of some important geophysical data. For the analysis we use the spectral characteristics of the ATMS instrument with the ephemeris for the NOAA-21 satellite [3]. The analysis is conducted for the USA and the simulation covers 8 consecutive days in July 2024. Effective and accurate microwave remote sensing of the atmosphere depends on the availability of interference-free spectrum windows at frequencies which are prescribed by physical processes [e.g. 4]. The family of resonant lines of the oxygen molecule near 60 GHz provides a unique opportunity to sample the vertical distribution of temperature and density from space, and it has been exploited for weather and climate studies from polar-orbiting satellites since 1978 (MSU on TIROS-N [5]). It remains a staple in the payloads operated by Russia, China, USA, Japan, France, India, UK, Ukraine [6] which are built around a common blueprint: a few wide-band (hundreds of MHz) channels around 50 GHz to sample the atmosphere and the surface while several more channels with high spectral resolution (few MHz) sample the individual resonant lines. Accurate retrieval of the environmental parameters depends upon the data provided by both sets of channels, and the their location in frequency space is not arbitrary and cannot be altered at will [7, 8]. The introduction of 5G technology in 2019 has driven telecommunication companies to request more bandwidth to be dedicated to their devices. This additional bandwidth is only available in spectral regions traditionally reserved for environmental and astrophysical research, such as the V-band between 50 and 60 GHz for up/downlink between satellites in low-earth orbits and terminals connected to fiberoptics network for distribution to high-speed local internet services. The power broadcast by the uplink communication leg is many orders of magnitude greater than the natural thermal signal emitted from the Earth scene. If the ground antenna were to perfectly align with the passive instrument’s antenna, the spaceborne receiver would suffer permanent, irreparable damage. While a direct boresight-to-boresight conjunction is extremely unlikely (even with a large constellation of satellites the fraction of the celestial sphere occupied by the satellites remains minuscule) the finite size of the ground station’s antenna beam in the sky suggests that the ATMS will be in the near background (as seen from the ground station) of one of the communication satellites and will be illuminated by either the main lobe or the near sidelobes of the uplink antenna more often than it is desirable. For our analysis we first calculate the position of the ATMS with respect to each of the ground stations at a resolution of 0.2 sec, then calculate the position of each of the comsats which are at least 25 deg above the station’s local horizon; finally we calculate the angle between the line-of-sight of the ATMS and the line-of-sight of the comsat. We assume that the gain pattern of the ground station is circularly symmetric; the angle-off-station-boresight then provides an attenuation of the uplink power which we use to assess the likely effect upon the passive instrument’s operations. We assume that each ground station can communicate with all the comsats in its field of view; this implies that, on average, a ground station can engage with 46 comsats simultaneously. The analysis is repeated for the case when the uplink broadcast within the ATMS passive channels (in-band scenario) and for the case when the uplink is limited to frequencies adjacent to the ATMS channels (out-of-band scenario). The antenna of the ground station is modelled as having a HPBW (Half-Power Beam Width) of 0.16 deg and EIRP (Equivalent Isotropic Radiated Power) of 70 dBW. We account for the geometric dissipation of the signal caused by the satellite orbital altitude, the attenuation induced by atmospheric gasses at 51 GHz and the mismatch between the circular polarization of the ground-based transmitting antenna and the linear polarization of the satellite-borne receiving antenna. The damages on ATMS are estimated from bench-level measurement conducted at the ATMS’ manufacturer facilities [unpublished].

passive microwave↗

Development and Testing of Pulse Guns for Combustion Instability Testing

To test liquid rocket engines (LREs) for combustion instabilities, devices such as bombs are often used to create pressure wave disturbances. Bombs, while effective, are inherently dangerous, expensive, and difficult to procure. Over the years, pulse guns have been used as a safer and more cost-effective way to generate controlled pressure disturbances in engine chambers. In anticipation of the need for stability testing of prototype LREs at NASA Marshall Space Flight Center (MSFC), a set of pulse guns have been designed, fabricated, tested, and characterized. The pulse gun program is funded by the RS-25 Engine Program managed out of MSFC and funded by NASA’s Space Launch System (SLS) through the MSFC Liquid Engines Office (LEO). A pulse gun is a simple device – like a traditional gun, it has a breech and barrel section. However, unlike a traditional gun, there is no bullet, as the purpose of the pulse gun is strictly to fire a high pressure pulse. Instead of a firing pin and primer that would normally be used in a traditional gun, an initiator is used to activate the gun powder. The initiator is a highly reliable pyrotechnic initiator. For this study, clone versions of the NASA Standard Initiator (NSI) were used. The initiator is used to ignite a pre-measured amount of gun powder loaded into the breech. When the pressure of the burnt gun powder exceeds the set point of a downstream burst disk, the disk ruptures, allowing the high pressure pulse to travel downstream through the barrel section. A ballistic pressure transducer located in the breech section is used to measure the short duration, high pressure pulse. Some configurations of the pulse gun have barrel sections that include one or more additional ballistic pressure transducers. These additional pressure measurements help track the degradation and damping of the pulse as it travels out of the barrel section. The measurements may also be used to determine the velocity of pulse. The objective of this paper is to present the different variants of this newly-developed pulse gun and characterize performance over a range of parameters. The parameters include breech diameter, barrel length, amount of gun powder used, the way the gun powder is packed, and the pressure setting of the burst disk. There are a total of six pulse guns: three with the 0.25 inch breech and three different length barrels, and three with the 0.40 inch breech and three different length barrels. For both breech sizes, the associated barrels are referred to as standard, one-port, and two-port barrels. The standard barrel has no instrumentation and is likely the barrel that will be used for engine stability testing. The one- and two-port barrels were designed specifically for pulse gun component testing to allow measurements of the magnitude and timing of the high pressure pulse as it makes its way through the pulse gun. The burst disks tested were commercially-available burst disks designed to rupture at 8,000, 16,000, and 24,000 psid. Testing was accomplished by firing the pulse gun into a test chamber pressurized with nitrogen at about 2300 psig. A total of four Model 113B23 High frequency ICP® pressure sensors (10k psi) were mounted in the test chamber, in the same plane as the pulse gun. Two of the sensors had “trimmed” adapters, and two did not. Given the symmetrical configuration of the sensors within the test chamber, different amplitudes of pressure measurements are attributed to the use of trimmed versus untrimmed adapters. The untrimmed adapters, with their narrower passages, tended to amplify the pressure amplitudes by as much as 50%. In total, 41 pulse gun tests have been conducted.Data are still being analyzed, but some trends are apparent. For example, measurements taken within the pulse gun are shown in Figures 3 and 4 for Hot-fires (HFs) # 11 and 39, respectively. Both tests were identical in that the 0.40 inch ID breech, two-port barrel, 8,000 psid burst disk, and same amount of gunpowder (6.639 grains for HF# 11, and 6.576 grains for HF# 39) were used, and the back pressure in the nitrogen test chamber was ~2300 psig. Despite keeping all these variables constant, the results from these two tests look quite different. For both tests, the first pressure peak shown in the P2111 trace within the breech is the firing of the initiator. The second peak and any subsequent peaks are from the combustion of the gunpowder. It should be noted that because the sensors are dynamic, the ~2300 psig baseline static pressure is not shown in the figures. For HF# 11, the burst disk clearly ruptured in a little over 0.1 ms from the time the initiator was fired. This is evident from the abrupt rise in pressure for P2112 followed by another abrupt rise in pressure for P2113. The sensors track the movement of the high pressure pulse through the pulse gun. For HF# 39, however, there was nearly 1 ms delay between the initiator firing and the burst disk rupturing. The only difference between these two tests was the way in which the gunpowder was packed within the breech. For HF# 11, the gunpowder was poured directly into the breech, atop the initiator, and held in place with a vegetable fiber wad appropriately sized for the barrel diameter. Once in place, the wad firmly held the gunpowder in place. For HF# 39, the gunpowder was rolled within cigarette paper with the ends of the paper twisted. The rolled gunpowder was gently pushed down into the breech, toward the initiator, and no wad was used. This method of loading the gunpowder consistently produced longer delays in the rupture of the burst disk, presumably because there was a steady, fuller burning of the gunpowder. The burst disk was rated for 8,000 psid, which means that the burst disk should not have opened until the pressure in the breech reached ~10,300 since there was ~2300 psig back pressure in the test chamber. In the case of HF# 11, the breech pressure only reached roughly 4200 psi before there was leakage either through or past the burst disk, but for HF#39, the breech pressure reached 10,600 psi prior to the burst disk opening. The temperature of the burst disk is certainly a factor in determining when the burst disk will rupture. A sharp rise in the burst disk temperature could weaken it and cause it to stray from the designed set point burst pressure. The burst disks are manufactured of Inconel to try to reduce their sensitivity to temperature, but how the burst disks perform above 900°F is not documented. All burst disks were engineered and tested by the manufacturer for ambient temperature conditions. One theory for why rolling the gunpowder in cigarette paper may produce better results is related to a more controlled temperature environment. The cigarette paper may act as a sheath that protects the burst disk from excessive temperatures during the rapid buildup of pressure in the breech. The ultimate benefit of not breeching the burst disk prematurely is shown in Figures 5 and 6, which graphs the resultant pressure pulse produced in the test chamber. For HF #11, the peak pressure on the P2123 trimmed adapter is about 270 psi, while for HF# 39 it is 343 psi, which is about 1.27 times greater in magnitude. In either case, the pressure pulse is ~20 μs in width once it reaches the test chamber. Another pressure rise about 100 μs later can be seen on the sensors located across from the pulse gun port, P2124 and P2121. These pressures are lower due to attenuation of the pressure wave as it moves across the test chamber. After the wave encounters the wall, it reflects back and forth within the chamber until it completely dissipates.

Pulse Gun↗