Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable density flow”

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 343 records · Page 19

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

CFD↗

An Overset Grid Motion Driver for Resolving Entry Vehicle CFD Simulations with Large Changes in Wake Orientation

Time-accurate computational fluid dynamics (CFD) simulations using overset meshes can enable reduced computational cost compared to a rigid, monolithic mesh for blunt-body atmospheric entry vehicles. If the vehicle has a significant change in attitude, the high-density region of the overset mesh designed for resolving the subsonic wake can be kept in position as the wake orientation changes relative to the body. Current simulations generally keep the wake mesh fixed in attitude while the body rotates, assuming the wake orientation remains relatively constant. However, changes in the relative velocity over a long-duration simulation or effects of the body geometry can invalidate this assumption. This work presents a method for calculating the desired mesh orientation to resolve a blunt body wake using the atmosphere-relative velocity of the vehicle. A motion driver using these calculations is implemented in a CFD-in-the-loop flight model bridging the FUN3D flow solver and POST2 trajectory propagator. Overset meshing is performed using the Yoga extension for FUN3D. Results are presented for simulations with both overset and monolithic meshes under large changes in relative velocity (and therefore wake orientation). The results demonstrate that the motion driver successfully orients an overset wake mesh to capture the subsonic wake with sufficient resolution despite variable body rotation and body-relative freestream velocity.

Computational Fluid Dynamics↗

Why does preferential diffusion strongly affect premixed turbulent combustion?

Combustion of premixed reactants in a turbulent flow is a classical but unresolved problem. The key problem is to explain the following data: the maximal turbulent and laminar burning velocities u(sub t) and u(sub L) occur at different equivalence ratios Phi. It is known that the equivalence ratio varies along a curved flame if molecular diffusivity D(sub fuel) does not equal D(sub oxygen). However, the mean flame radius of curvature is much larger than the laminar flame thickness delta-L. Therefore, significant influence of preferential diffusion should occur only if the flame propagation speed varies with flame curvature. This conclusion agrees with Zel'dovich's long-standing idea about the important role of leading points of a flame. The main objective of this paper is to prove Zel'dovich's hypothesis. An equation for the mean flame surface area density (MFSAD) is employed for this purpose. The second objective of this paper is to suggest a different approach to the derivation of the equation for MFSAD. It is based on the pdf equation for the reaction progress variable C and the relation between the pdf and MFSAD. This treatment suggests an entirely different closure assumption.

Kuznetsov, Vadim R.↗

Development of a High Reliability Compact Air Independent PEMFC Power System

Autonomous Underwater Vehicles (AUV's) have received increasing attention in recent years as military and commercial users look for means to maintain a mobile and persistent presence in the undersea world. Compact, neutrally buoyant power systems are needed for both small and large vehicles. Historically, batteries have been employed in these applications, but the energy density and therefore mission duration are limited with current battery technologies. Vehicles with stored energy requirements greater than approximately 10 kWh have an alternate means to get long duration power. High efficiency Proton Exchange Membrane (PEM) fuel cell systems utilizing pure hydrogen and oxygen reactants show the potential for an order of magnitude energy density improvement over batteries as long as the subsystems are compact. One key aspect to achieving a compact and energy dense system is the design of the fuel cell balance of plant (BOP). Recent fuel cell work, initially focused on NASA applications requiring high reliability, has developed systems that can meet target power and energy densities. Passive flow through systems using ejector driven reactant (EDR) circulation have been developed to provide high reactant flow and water management within the stack, with minimal parasitic losses compared to blowers. The ejectors and recirculation loops, along with valves and other BOP instrumentation, have been incorporated within the stack end plate. In addition, components for water management and reactant conditioning have been incorporated within the stack to further minimize the BOP. These BOP systems are thermally and functionally integrated into the stack hardware and fit into the small volumes required for AUV and future NASA applications to maximize the volume available for reactants. These integrated systems provide a compact solution for the fuel cell BOP and maximize the efficiency and reliability of the system. Designs have been developed for multiple applications ranging from less than 1 kWe to 70 kWe. These systems occupy a very small portion of the overall energy system, allowing most of the system volume to be used for reactants. The fuel cell systems have been optimized to use reactants efficiently with high stack efficiency and low parasitic losses. The resulting compact, highly efficient fuel cell system provides exceptional reactant utilization and energy density. Key design variables and supporting test data are presented. Future development activities are described.

Vasquez, Arturo↗

Winds Measured by the Rover Environmental Monitoring Station (REMS) During the Mars Science Laboratory (MSL) Rover's Bagnold Dunes Campaign and Comparison with Numerical Modeling Using MarsWRF

A high density of REMS wind measurements were collected in three science investigations during MSL's Bagnold Dunes Campaign, which took place over approx. 80 sols around southern winter solstice (Ls approx. 90deg) and constituted the first in situ analysis of the environmental conditions, morphology, structure, and composition of an active dune field on Mars. The Wind Characterization Investigation was designed to fully characterize the near-surface wind field just outside the dunes and confirmed the primarily upslope/downslope flow expected from theory and modeling of the circulation on the slopes of Aeolis Mons in this season. The basic pattern of winds is 'upslope' (from the northwest, heading up Aeolis Mons) during the daytime (approx. 09:00-17:00 or 18:00) and 'downslope' (from the southeast, heading down Aeolis Mons) at night (approx. 20:00 to some time before 08:00). Between these times the wind rotates largely clockwise, giving generally westerly winds mid-morning and easterly winds in the early evening. The timings of these direction changes are relatively consistent from sol to sol; however, the wind direction and speed at any given time shows considerable intersol variability. This pattern and timing is similar to predictions from the MarsWRF numerical model, run at a resolution of approx. 490 m in this region, although the model predicts the upslope winds to have a stronger component from the E than the W, misses a wind speed peak at approx. 09:00, and under-predicts the strength of daytime wind speeds by approx. 2-4 m/s. The Namib Dune Lee Investigation reveals 'blocking' of northerly winds by the dune, leaving primarily a westerly component to the daytime winds, and also shows a broadening of the 1 Hz wind speed distribution likely associated with lee turbulence. The Namib Dune Side Investigation measured primarily daytime winds at the side of the same dune, in support of aeolian change detection experiments designed to put limits on the saltation threshold, and also appears to show the influence of the dune body on the local flow, though less clearly than in the lee. Using a vertical grid with lower resolution near the surface reduces the relative strength of nighttime winds predicted by MarsWRF and produces a peak in wind speed at approx. 09:00, improving the match to the observed diurnal variation of wind speed, albeit with an offset in magnitude. The annual wind field predicted using this grid also provides a far better match to observations of aeolian dune morphology and motion in the Bagnold Dunes. However, the lower overall wind speeds than observed and disagreement with the observed wind direction at approx. 09:00 suggest that the problem has not been solved and that alternative boundary layer mixing schemes should be explored which may result in more mixing of momentum down to the near-surface from higher layers. These results demonstrate a strong need for in situ wind data to constrain the setup and assumptions used in numerical models, so that they may be used with more confidence to predict the circulation at other times and locations on Mars.

Newman, Claire E.↗

Airborne Observations of Aerosol Properties in Southeast Asia: Overview of Emissions, Cloud Processing, and Long Range Transport during NASA CAMP2Ex

The Southwest Monsoon represents an important synoptic scale meteorological feature that modulates precipitation throughout Southeast Asia. In addition, southwesterly flow can funnel large amounts of biomass burning aerosol into heavily populated areas throughout the Philippines, mixing with urban emissions and representing a significant air quality issue. From a climate perspective, understanding interactions between these particles and clouds has important societal ramifications through frequent extreme periods of flood and drought. Here we present airborne observations made in support of the NASA CAMP2Ex (Cloud and Aerosol Monsoonal Processes-Philippines Experiment) mission aboard the Wallops Flight Facility P-3 research aircraft, based at Clark International Airport in the Luzon region of Northern Philippines. Nineteen research flights were conducted in August-October of 2019. An extensive suite of aerosol microphysical, optical, and chemical instrumentation was deployed, along with trace gas and remote sensing measurements, to characterize aerosol emissions and their vertical/spatial variability associated with moderately convective marine cloud systems in the region. Analysis will focus on characterizing disparate aerosol sources with respect to derived chemical and optical properties relevant to global circulation modeling. For example, higher than expected mass scattering efficiencies for aged biomass burning emissions, derived from a combination of nephelometer and aerosol mass spectrometer measurements, are attributed to increased particle size and decreased density. Results will also illustrate the evolution of aerosols from Indonesian agricultural biomass burning activities into the Philippines region, in comparison with urban emissions from the megacity Manila and transported air masses from mainland Asia.

Luke David Ziemba↗

Improved Prediction of Losses with Large Eddy Simulation in a Low-Pressure Turbine

There is a need to improve predictions of losses resulting from large eddy simulations (LES) of low-pressure turbines (LPT) in gas turbines. This may be done by assessing the accuracy of predictions against validation data and understanding the source of any inaccuracies. LES is a promising approach for capturing the laminar/turbulent transition process in a LPT. In previous studies, the authors utilized LES to model the flow field over a Variable Speed Power Turbine (VSPT) blade and successfully captured characteristic features of separation/reattachment and transition on the suction side at both the cruise (positive incidence) and take-off conditions (negative incidence) and as well, simulated the effect of free-stream turbulence (FST) on those phenomena. The predicted pressure loading profiles agreed well with the experimental data for both a high and a low FST case at a Reynolds number of Re(ex)=220,000. In this paper, we present wake profiles resulting from computations for a range of FST values. Although the predicted wake profiles for the lowest FST case (Tu=0.5%) matched the experimental data, at higher FST (Tu=10-15%,) the wake was wider than the experimentally measured wake and for both cases were displaced laterally when compared to the experimental measurements. In our investigation of the causes of the said discrepancies we have identified important effects which could strongly influence the predicted wake profile. Predicted losses were improved by assuring the validity of the flow solution. This was done by utilizing spectral analysis to scrutinize the dynamic behavior of the wake and determine solution accuracy resulting from low mesh density and low accuracy of convective modeling.

LES-LDKM↗

Efficient Gradient-Based Shape Optimization Methodology Using Inviscid/Viscous CFD

The formerly developed preconditioned-biconjugate-gradient (PBCG) solvers for the analysis and the sensitivity equations had resulted in very large error reductions per iteration; quadratic convergence was achieved whenever the solution entered the domain of attraction to the root. Its memory requirement was also lower as compared to a direct inversion solver. However, this memory requirement was high enough to preclude the realistic, high grid-density design of a practical 3D geometry. This limitation served as the impetus to the first-year activity (March 9, 1995 to March 8, 1996). Therefore, the major activity for this period was the development of the low-memory methodology for the discrete-sensitivity-based shape optimization. This was accomplished by solving all the resulting sets of equations using an alternating-direction-implicit (ADI) approach. The results indicated that shape optimization problems which required large numbers of grid points could be resolved with a gradient-based approach. Therefore, to better utilize the computational resources, it was recommended that a number of coarse grid cases, using the PBCG method, should initially be conducted to better define the optimization problem and the design space, and obtain an improved initial shape. Subsequently, a fine grid shape optimization, which necessitates using the ADI method, should be conducted to accurately obtain the final optimized shape. The other activity during this period was the interaction with the members of the Aerodynamic and Aeroacoustic Methods Branch of Langley Research Center during one stage of their investigation to develop an adjoint-variable sensitivity method using the viscous flow equations. This method had algorithmic similarities to the variational sensitivity methods and the control-theory approach. However, unlike the prior studies, it was considered for the three-dimensional, viscous flow equations. The major accomplishment in the second period of this project (March 9, 1996 to March 8, 1997) was the extension of the shape optimization methodology for the Thin-Layer Navier-Stokes equations. Both the Euler-based and the TLNS-based analyses compared with the analyses obtained using the CFL3D code. The sensitivities, again from both levels of the flow equations, also compared very well with the finite-differenced sensitivities. A fairly large set of shape optimization cases were conducted to study a number of issues previously not well understood. The testbed for these cases was the shaping of an arrow wing in Mach 2.4 flow. All the final shapes, obtained either from a coarse-grid-based or a fine-grid-based optimization, using either a Euler-based or a TLNS-based analysis, were all re-analyzed using a fine-grid, TLNS solution for their function evaluations. This allowed for a more fair comparison of their relative merits. From the aerodynamic performance standpoint, the fine-grid TLNS-based optimization produced the best shape, and the fine-grid Euler-based optimization produced the lowest cruise efficiency.

Baysal, Oktay↗

Impact of U-10Mo HALEU fuel element tolerances on the Massachusetts Institute of Technology reactor safety and operational performance – Thermal hydraulics

The U.S. is coordinating efforts for the conversion of six U.S. High Performance Research Reactors (USHPRR) including one critical facility from highly enriched uranium (HEU) to low-enriched uranium (LEU). In order to continue the mission of these reactors, including the Massachusetts of Institute of Technology Reactor (MITR), and achieve similar performance, high-assay low-enriched uranium (HALEU) with a high-density metallic alloy of uranium with 10 wt% molybdenum (U-10Mo) is being evaluated. The impact of the fabrication specification and tolerances was assessed following the preliminary design of the MITR LEU fuel elements using the U-10Mo monolithic alloy. This research focuses on the analysis of fabrication specification impact on thermal hydraulics (TH) characteristic of the MITR LEU core as a function of the variation of the relevant fuel specification parameters (e.g., coolant channel gap thickness, fuel plate thickness, etc.). The analyses are performed based on an all-fresh LEU fuel conversion plan identified in a preliminary safety analysis report submitted to the Nuclear Regulatory Commission. The reactor power margin to the onset of nucleate boiling (ONB) is assessed under the limiting safety system settings (LSSS), where a scram occurs, to ensure there is sufficient margin to the reactor safety limit, which is defined by the onset of flow instability that occurs after the ONB. The best estimate plus uncertainty approach is employed to analyze this TH characteristic, which yields realistic results while maintaining adequate conservatism, utilizing a statistical uncertainty propagation method with the STAT7 code. The TH characteristic is analyzed as a function of the variability of the specification parameters resulting from the fabrication process. The main findings of this study show that the MITR core can meet the TH safety and operational requirements at the all-LEU initial core startup (cycle 1), selected transition cycles (most reactive cycle and most limiting cycle: cycle 3 and 5, respectively) and equilibrium (cycle 14) cores under all limiting fabrication parameter combinations considered. In addition, the analyses show that the dependency of the core power margin to ONB on those specification parameters that have the most direct impact on TH performance is non-linear but monotonically decreasing within the specification tolerances. The third order polynomial fit curves are reported in detail for selected limiting cases and can serve as a powerful tool for future MITR fuel management in cases such as when HALEU supply is established that may allow additional cycle length or other operational benefits.

Conversion↗

Paper or Plastic? Multiscale Material Handling Properties of Two Model Municipal Solid Waste Streams

Purpose: Municipal Solid Waste (MSW) is a potentially valuable sustainable feedstock for fuel and chemical production due to its carbon-rich content and low cost. This study aims to assess the material handling properties of paperand plastic-rich MSW feedstocks to mitigate equipment failure and processing downtime. Methods: The material handling properties of crumbled MSW feedstocks were measured using apowder rheometer with mass flow hopper calculations to assess handling performance. Inverse gas chromatography was use to measure the surface energy differences between feedstocks. Electron microscopy and Raman spectroscopy was used to evaluate microscale features that may contribute to material handling differences. Results: Plastic and paper rich feedstocks crumbled to a nominal 2 mm particle size were observed to have similar flow and handling characteristics with reasonable hopper outlets. 2 mm plastic rich crumbles, with their higher bulk density, exhibited superior flow performance. By contrast, 4 mm material required significantly larger hopper outlets, indicating poor flowability. Paper rich and 4 mm plastic rich samples displayed broad particle size distributions, which contributed to particle interlocking, jamming, and other flow issues. Electron microscopy revealed that plastic rich samples were significantly smoother, enhancing their flowability compared to the rougher, paper rich materials. Conclusions: This study establishes critical material handling baselines for processing MSW as a viable feedstock for fuel and chemical production. The findings highlight the importance of optimizing particle size and feedstock composition to improve flowability and handling performance.

09 BIOMASS FUELS↗

Growth and Morphology of Supercritical Fluids, a Fluid Physics Experiment Conducted on Mir, Complete

The Growth and Morphology of Supercritical Fluids (GMSF) is an international experiment facilitated by the NASA Glenn Research Center and under the guidance of U.S. principal investor Professor Hegseth of the University of New Orleans and three French coinvestigators: Daniel Beysens, Yves Garrabos, and Carole Chabot. The GMSF experiments were concluded in early 1999 on the Russian space station Mir. The experiments spanned the three science themes of near-critical phase separation rates, interface dynamics in near-critical boiling, and measurement of the spectrum of density fluctuation length scales very close to the critical point. The fluids used were pure CO2 or SF6. Three of the five thermostats used could adjust the sample volume with the scheduled crew time. Such a volume adjustment enabled variable sample densities around the critical density as well as pressure steps (as distinct from the usual temperature steps) applied to the sample. The French-built ALICE II facility was used for these experiments. It allows tightly thermostated (left photograph) samples (right photograph) to be controlled and viewed/measured. Its diagnostics include interferometry, shadowgraph, high-speed pressure measurements, and microscopy. Data were logged on DAT tapes, and PCMCIA cards and were returned to Earth only after the mission was over. The ground-breaking near critical boiling experiment has yielded the most results with a paper published in Physical Review Letters (ref. 1). The boiling work also received press in Science Magazine (ref. 2). This work showed that, in very compressible near-critical two-phase pure fluids, a vapor bubble was induced to temporarily overheat during a rapid heating of the sample wall. The temperature rise in the vapor was 23-percent higher than the rise in the driving container wall. The effect is due to adiabatic compression of the vapor bubble by the rapid expansion of fluid near the boundary during heatup. Thermal diffusivity is low near the critical point, so getting heat out of the compressed bubble is observably slow. This gives the appearance of a backward heat flow, or heat flow from a cold surface to a warm fluid.

Wilkinson, R. Allen↗

Estimating Resolution Lengths of Hybrid Turbulence Models

A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.

Abdol-Hamid, Khaled S.↗

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22, 2015 was one of the largest storms in the current solar cycle, resulting from an active region on the Sun that produced numerous coronal mass ejections (CMEs) and associated interplanetary shock waves. On June 22 at 18:36 UT the magnetosphere was impacted by the leading-edge shock wave and a sheath carrying a large and highly variable interplanetary magnetic field (IMF) Bz with values ranging from +25 to -40 nT. During the subsequent interval from 0000 to 0800 UT, there was a second intensification of the geomagnetic storm resulting from the impact of the CME. We present dramatic responses of simultaneous particle measurements from the high-altitude Magnetospheric Multiscale Mission (MMS) at high altitudes in the magnetosphere (approx. 9-12 Re) and from the low-altitude (F-region) Floating Potential Measurement Unit (FPMU) on board the International Space Station (ISS). We analyze potential causes of these dramatic particle flux dropouts by putting them in the context of storm-time electrodynamics, and support our results with numerical simulations of the global magnetosphere and ionosphere. During the sheath phase of the storm, the MMS spacecraft in the near-earth equatorial plane observed a rapid reconfiguration of the magnetic field near 1923 UT. Initially in the warm plasmasheet, particle flux dropouts were observed as they tracked the plasma-sheet to lobe transitions with the stretching and thinning of the plasmasheet. Anti-sunward flowing O+ ions of ionospheric origin were also measured during this period, confirming that the MMS spacecraft temporarily was in a lobe.

Coffey, Victoria↗

Coupled ice-ocean variability in the Greenland Sea

A major surface feature of the Greenland Sea during winter is the frequent eastward extension of sea ice south of 75 deg N and an associated embayment to the north. These features are nominally connected with the East Greenland Current, and both the promontory) and the embayment are readily apparent on climatic ice charts. However, there are significant changes in these features on time-scales as short as a few days. Using a combination of satellite microwave images (SSM/I) of ice cover, meteorological data and in situ velocity, temperature and salinity records, we relate the ice distribution and its changes to the developing structure and circulation of the upper ocean during winter 1988-1989. Our measurements illustrate the preconditioning that leads to convective overturn, which in turn brings warmer water to the surface and results in the rapid disappearance of ice. In particular, the surface was cooled to the freezing point by early December and the salinity then increased through ice formation (about 0.016 m/d) and brine rejection. Once the vertical density gradient was sufficiently eroded, a period of high heat flux (greater than 300 W/sq m) in late January provided enough buoyancy loss to convectively mix the upper water column to at least 200 m. We estimate vertical velocities at about 3 cm/s downward during the initial sinking. The deepening of the thermocline raised surface temperatures by over 1 C resulting in nearly 1.5 x 105 km of ice-melt within two days. Average rates of ice retreat are about 2 km/d southwestward, generally consistent with a wind-driven flow. Comparison of hydrographic surveys from before and after the overturning indicate the fresh water was advected out of the area, possibly to the south and east of our moorings.

Roach, A. T.↗

Observations of Deep Ionospheric F-Region Density Depletions with FPMU Instrumentation and Their Relationship with the Global Dynamics of the June 22-23, 2015 Geomagnetic Storm

The magnetic storm that commenced on June 22, 2015 was one of the largest storms in the current solar cycle, resulting from an active region on the Sun that produced numerous coronal mass ejections (CMEs) and associated interplanetary shock waves. On June 22 at 18:36 UT the magnetosphere was impacted by the leading-edge shock wave and a sheath carrying a large and highly variable interplanetary magnetic field (IMF) Bz with values ranging from +25 to -40 nT. During the subsequent interval from 0000 to 0800 UT, there was a second intensification of the geomagnetic storm resulting from the impact of the CME. We present dramatic responses of simultaneous particle measurements from the high-altitude Magnetospheric Multiscale Mission (MMS) at high altitudes in the magnetosphere (approx. 9-12 Re) and from the low-altitude (F-region) Floating Potential Measurement Unit (FPMU) on board the International Space Station (ISS). We analyze potential causes of these dramatic particle flux dropouts by putting them in the context of storm-time electrodynamics, and support our results with numerical simulations of the global magnetosphere and ionosphere. During the sheath phase of the storm, the MMS spacecraft in the near-earth equatorial plane observed a rapid reconfiguration of the magnetic field near 1923 UT. Initially in the warm plasma sheet, particle flux dropouts were observed as they tracked the plasma-sheet to lobe transitions with the stretching and thinning of the plasma sheet. Anti-sunward flowing O+ ions of ionospheric origin were also measured during this period, confirming that the MMS spacecraft temporarily was in a lobe.

Coffey, Victoria↗

Direct Numerical Simulations of Turbulent Convection: Variable Gravity and Uniform Rotation - Part 1

Turbulent convection may have played a major role in determining the structure and evolution of the primordial solar nebula, but current, incomplete models of convection and turbulence give very different results and remain largely untested in the absence of detailed astronomical observations. Numerical simulations provide an 'experimental' database for comparison with these models, and, to this end, direct numerical simulations of turbulent convection were performed with modifications intended to mimic some of the unique physical features of thin accretion disks, such as the primordial solar nebula: internal heating, a gravitational acceleration that is linearly proportional to the distance from midplane of the nebula, and rapid rotation. Peclet numbers in the simulations are comparable to those in solar nebula models; Rossby numbers in the simulations are an order of magnitude larger than those in solar nebula models because of the unrealistically high Prandtl and low, Reynolds numbers required to resolve all scales of the convective flow. We find that, despite the loss of buoyancy at midplane, turbulent motions easily penetrate the midplane region with little loss of intensity, providing efficient transport of heat and turbulent kinetic energy throughout the interior. A simple mixing length model modified to include rotation is found to give convective heat fluxes for the interior flow in rough agreement with the numerical simulations. We discuss these preliminary results with regard to assumptions about heating distributions and convective heat fluxes made in standard solar nebula models. More definitive comparisons with solar nebula modelling will become possible when more realistic effects of shear, density, stratification, and compression are included.

Cabot, W.↗

Recovery of postconsumer mechanically recycled polymers

Mechanical recycling plays a key role in reducing landfill bound plastics that pollute our environment. This process converts plastic waste into marketable pellets by sorting, cleaning, grinding into flakes, compounding in the molten state, and ultimately pelletizing. A primary restriction for the widescale usage of mechanical recycling is the highly variable quality and mechanical properties of the plastic waste feedstock. Degradation can occur during the plastic life cycle with the consumer, during the mechanical processing itself, or during the complex sorting process required to produce the feedstock. This study explores how rheological characterization can mitigate the batch-to-batch variability and identify a potential application for each batch. Shear and extensional rheology of “application-specific” virgin high-density polyethylene (HDPE) and virgin polypropylene (PP) was used as the control for this categorization process. Recycled HDPE and PP from three different streams were then measured and compared to the results from the control study. Rheological measurements proved to be very effective at providing sufficient differentiation to categorize the recycled polymer as suitable for different applications such as injection molding, blow molding, or thermoforming. Finally, the usage of an additional step to sort the recycled polymers by their initial use application was found to achieve a remarkably consistent recovery of application-specific material properties. Furthermore, this secondary sorting could provide significant added value for mechanically recycled polymers.

Differential scanning calorimetry↗

Modeling, Simulation and Control of a Propellant Mixer

This paper addresses some modeling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer must combine high pressure liquid hydrogen (LH2) and gaseous hydrogen (GH2) to produce and output flow that meets certain thermodynamic properties before it is fed into a test article. More precisely, this paper considers that the quantities to be tracked and/or regulated are mixer internal pressure, exit mass flow, and exit temperature. The available control inputs are given by three value positions, namely those of the GH2, LH2 and exit valves. The mixer is modelled by a system of two nonlinear ordinary differential equations having density and internal energy as states. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as controllability and uniqueness of the equilibrium point are shown to hold. Validation of the model against real data is also provided. As a first control approach, a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions. A feedback linearization controller is designed and used to achieve tracking and regulation of the outputs over an extended range of the variables of interest.

Richter, Hanz↗