A Monte Carlo approach to transition and free molecular flow problems.
Monte Carlo approach to transition and free molecular flow problems having mass as well as thermal motion of molecules
SEARCH · Engineering Papers
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.
Monte Carlo approach to transition and free molecular flow problems having mass as well as thermal motion of molecules
This work concerns both the theoretical prediction and measurement of structural parameters in open-cell highly porous polyurethane foams. Of particular interest are the dynamic flow resistance, thermal time constant, and mass structure factor and their dependence on frequency and geometry of the cellular structure. The predictions of cell size parameters, static flow resistance, and heat transfer as accounted for by a Nusselt number are compared with measurement. Since the static flow resistance and inverse thermal time constant are interrelated via the 'mean' pore size parameter of Biot, only two independent measurements such as volume porosity and mean filament diameter are required to make the predictions for a given fluid condition. The agreements between this theory and nonacoustical experiments are excellent.
The 8- by 6-Foot Supersonic Wind Tunnel (SWT) at Lewis Research Center is available for use by qualified researchers. This manual contains tunnel performance maps which show the range of total temperature, total pressure, static pressure, dynamic pressure, altitude, Reynolds number, and mass flow as a function of test section Mach number. These maps are applicable for both the aerodynamic and propulsion cycle. The 8- by 6-Foot Supersonic Wind Tunnel is an atmospheric facility with a test section Mach number range from 0.36 to 2.0. General support systems (air systems, hydraulic system, hydrogen system, infrared system, laser system, laser sheet system, and schlieren system are also described as are instrumentation and data processing and acquisition systems. Pretest meeting formats are outlined. Tunnel user responsibility and personal safety requirements are also stated.
Neutral interstellar atoms penetrate deeply into the inner heliosphere before they become ionized by various processes. As ions they are picked-up by the frozen-in magnetic fields and are convected outwards with the solar wind plasma. Thereby the primary plasma flow is mass, momentum, and energy-loaded. The dynamics of the distant multi-constituent solar wind is, however, not solely determined by these loading processes, but is also affected by the wave-mediated pick-up ion pressure gradients derivable from the pick-up ion distribution function. The action of the radial components of these pressures essentially counter balances the decelerating effect of the solar wind momentum loading, diminishing strongly the deceleration of the distant solar wind. Furthermore the latitudinal components of the pick-up ion pressures induce latitudinal forces acting on the multiconstituent solar plasma outflow and inducing nonradial bulk flow components. The enforced nonradial outflow geometry on the upwind hemisphere may partly be responsible for the magnetic flux deficit which was claimed since several years in the PIONEER-10 magnetic flux data.
We use a variety of visualization techniques to display the interior and surface flows in a double white dwarf binary undergoing direct-impact mass transfer and evolving dynamically to a merger. The structure of the flow can be interpreted in terms of standard dynamical, cyclostrophic, and geostrophic arguments. We describe and showcase some visualization and analysis techniques of potential interest for astrophysical hydrodynamics. In the context of R Coronae Borealis stars, we find that mixing of accretor material with donor material at the shear layer between the fast accretion belt and the slower rotating accretor body will always result in some dredge-up. We also discuss briefly some potential applications to other types of binaries.
A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.
An Integrated Refrigeration and Storage (IRAS) experimental system called the Ground Operations and Demonstration Unit for Liquid Hydrogen (GODU-LH2) demonstrated the ability store cryogens in a zero boiloff (ZBO) process, and to densify liquid hydrogen by reducing the temperature and pressure down to the triple point resulting in solidification. The incompressible Navier-Stokes was utilized to simulate ZBO, densification, and solidification of liquid hydrogen (LH2) in the IRAS tank. The simulations were performed using a commercially available Computational Fluid Dynamics (CFD) pressure-based mass and momentum flow model and an enthalpy-porosity energy model. Results demonstrated the simulation’s ability to predict time-dependent flow and temperature fields and solid-liquid phase locations for hydrogen during ZBO, densification, and solidification. The simulations showed good agreement with experimental data, with errors within acceptable ranges for temperature and pressure predictions, provided detailed insights into natural convection and solid-liquid phase change dynamics. These findings are critical for the design of future cryogenic fluid management systems.
The advent of novel microreactor technology has driven a focused effort to explore safety and efficiency improvements that can be achieved through the use of automated system control. Development of control strategies, especially for initial demonstration, requires an adequate surrogate environment to safely research failure modes and control integration with realistic hardware delay. However, efficiency gains from control strategies are improved when the scope of controller action is expanded to include system-level dynamics such as downstream heat extraction and mass flow. For this reason, a balance-of-plant (BOP) model of a representative microreactor system has been developed using the TRANsient Simulation Framework of Reconfigurable Models library in Modelica. This model captures a reactor and primary NaK coolant loop that represent corresponding system components of the Microreactor Applications Research Validation and EvaLuation (MARVEL) design as well as a secondary coolant loop and heat extraction representative of the Microreactor Agile Non-Nuclear Experimental Test Bed (MAGNET). This model configuration allows for hardware-in-the-loop (HIL) integration with microreactor automated control system (MACS) hardware in real time through a Python-based gRPC client. Real-time simulation of model performance with emulated hardware and communication delay suggests that under independent proportional-integral-derivative control of BOP model drum dynamics and downstream heat extraction, stable power load following is achievable. A slight delay in load following, filtering of high-frequency dynamics, and localized temperature fluctation suggest room for improvement through the development of higher-level control strategies. The simulated coupling of the MAGNET facility lays the groundwork for future digital twin analysis with a coupled MACS-MAGNET HIL demonstration.
A brief review is given of non-flare investigations using data obtained by the Ultraviolet Spectrometer and Polarimeter on the Solar Maximum Mission. The major topics described are sunspot research including magnetic field measurements, oscillations, and models; mass motions in quiet and active regions including steady flows and acoustic waves; and prominence research including physical conditions, dynamics, and mass motions around prominences. Also discussed are studies of UV bursts, the formation of the Cl I line at 1351 A, ozone in the terrestrial atmosphere, and active regions using correlated observations from other instruments on the spacecraft or on the ground.
Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.
Sprays occur in a wide variety of industrial and power applications and in materials processing. A liquid spray is a two-phase flow with a gas as the continuous phase and a liquid as the dispersed phase in the form of droplets or ligaments. The interactions between the two phases--which are coupled through exchanges of mass, momentum, and energy--can occur in different ways at disparate time and length scales involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the ratecontrolling processes associated with turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates of the spray, among many other factors. With the aim of developing an efficient solution procedure for use in multidimensional combustor modeling, researchers at the NASA Glenn Research Center have advanced the state-of-the-art in spray computations in several important ways.
The oxidation of pyrolytic carbon (PyC) deposited via fluidized bed chemical vapor deposition was characterized and compared with that of standard nuclear-grade graphite. The materials were heated at 700 to 1000 °C in a thermogravimetric analysis system under 20% v/v O 2 flow, allowing for direct comparison of dynamic oxidative mass change in each material. Further, three different PyC samples fabricated under different conditions exhibited variation in total mass loss and mass loss rate, varying by as much as 709 mg/cm 2 in total mass loss and 14.2 (mg/cm 2 )/min in mass loss rate at a single temperature. These variations highlight the correlation between PyC microstructure/defect density and oxidation susceptibility. Additionally, changes in the microstructure and composition between PyC and graphite were characterized via scanning electron microscopy and correlated to the mass loss results. The results of this work have implications toward the safety of tristructural isotropic (TRISO) and other coated particle fuels, especially under off-normal conditions, given the limited information that exists about the oxidation behavior of PyC.
An analytical model for the orifice pulse tube is developed, in which the system is described in terms of average enthalpy flow with such simplifying assumptions as an ideal gas and sinusoidal pressure variation. Phasor analysis is used to represent temperature, pressure, and mass-flow rate waves in vector form. The model predictions, namely, that the refrigeration power of a pulse-tube refrigerator is proportional to the average pressure, the pulse frequency, the mass-flow ratio, and the square of the dynamic pressure ratio, were verified by experimental measurements. It was found that, as a result of the simplifying assumptions, the magnitudes of the refrigeration power predicted by the model were between three and five times greater than experimental values.
A numerical, and experimental investigation to study the effects of flow distortion on a Mass Flow Plug (MFP) used to control and measure mass-flow during an inlet test has been conducted. The MFP was first calibrated using the WIND-US flow solver for uniform (undistorted) inflow conditions. These results are shown to compare favorably with an experimental calibration under similar conditions. The effects of distortion were investigated by imposing distorted flow conditions taken from an actual inlet test to the inflow plane of the numerical simulation. The computational fluid dynamic (CFD) based distortion study only showed the general trend in mass flow rate. The study used only total pressure as the upstream boundary condition, which was not enough to define the flow. A better simulation requires knowledge of the turbulence structure and a specific distortion pattern over a range of plug positions. It is recommended that future distortion studies utilize a rake with at least the same amount of pitot tubes as the AIP rake.
A test program was conducted on a two-dimensional supersonic inlet. Internal disturbances in diffuser exit mass flow were produced by oscillating overboard bypass doors. Open-loop dynamic responses of shock position, throat exit and diffuser exit static pressures are presented. The steady-state and dynamic coupling between ducts were also obtained. The experimental results from the two-dimensional inlet are compared to results from a similar size axisymmetric inlet and also to a transfer function synthesis program.
A critical barrier to the deployment of MSRs is the absence of a well-defined nuclear material control and accounting (MC&A) approach, a vital prerequisite to meet NRC licensing requirements as well as facilitating future international exports. Liquid-fueled MSR variants (especially those incorporating online refueling or reprocessing) present a unique set of challenges to traditional nuclear material control and accountancy. Unlike solid-fueled light-water reactors, traditional item-counting methods cannot be applied as an accountancy strategy. Rather, MC&A approaches to MSR variants (including both uranium and thorium fueled designs) are more analogous to bulk material handling facilities (e.g., enrichment and reprocessing); yet further complicating matters is the fact that fuel medium is also highly radioactive. Meanwhile, the space of MSRs covers a broad range of design parameters, including thermal and fast spectra designs; operation in actinide breeder or burner modes, choice of actinide fuel used (e.g., 235 U, 232 Th / 233 U, denatured 233 U, etc.), pool or loop-type configuration, and even different salt chemistry. Each of these design choices introduce significant challenges to MC&A approaches within MSR facilities. We propose to bridge this gap for liquid-fueled MSRs by developing a modular, component-based test framework for evaluating viable process monitoring and MC&A techniques specifically suited to liquid-fueled MSR system variants employing online refueling or reprocessing. This test platform will consist of a toolbox of independent process modules representing discrete physical units (such as the reactor core, off-gas processing, decay tanks, and actinide separation units), each with its own self-contained physics responsive to the input mass flow, along with appropriate measurement models that can be coupled to key flow points. These dynamic physical signatures thus afford the ability to test the viability and efficacy of potential accountancy techniques under the full range of reactor operating conditions. As process modules are connected via mass flows, the result is a reconfigurable, generic MSR mass flow model capable of serving as an MC&A test platform for a broad spectrum of possible MSR configurations. The resulting MSR MC&A toolbox will enable robust assessment of accountancy strategies for this unique facility type, including analysis of physical feedbacks arising both from depletion of the fuel over time as well as from potential off-normal events, including those introduced by equipment failures (e.g., a pump failure) as well as by malicious action (i.e., attempts to divert material). The proposed toolbox both addresses a critical needs area for the MPACT analysis toolkit while leveraging existing MPACT-sponsored tools, especially with respect to simulation of measurement and accountancy techniques for advanced fuel cycle facilities. Beyond enabling new analysis capabilities for MSR systems, the design of this toolbox will be to such to maximize compatibility with existing MPACT tools, such to enhance existing facility MC&A analysis capabilities.
Abstract This paper analyzes the dynamics of unstable azimuthal thermoacoustic modes in a lean premixed combustor. Azimuthal modes can be decomposed into two counter rotating waves where they can either compete and potentially suppress one of them (spinning) or coexist (standing), depending on the operating conditions. This paper describes experimental results of the dynamical behaviors of these two waves. The experimental data were taken at different mass flow rates as well as different azimuthal fuel staging in a multi-nozzle can combustor. It is shown that at a low flow rate with uniform fuel distribution, the two waves have similar amplitudes, giving rise to a standing wave. However, the two amplitudes are slowly oscillating out of phase to each other, and the phase difference between the two waves also shows oscillatory behavior. For an intermediate flow rate, the dynamics show intermittency between standing and spinning waves, indicating that the system is bistable. In addition, the phase difference dramatically shifts when the mode switches between standing and spinning waves. For a high flow rate, the system stabilizes at a spinning wave most of the time. These experimental observations demonstrate that not only the amplitudes of two waves but also the phase difference plays an important role in the dynamics of azimuthal mode. For non-uniform azimuthal fuel staging, the modal dynamics exhibit only an oscillatory standing wave behavior regardless of the mass flow rate. Compared to the uniform fuel staging, however, the pressure magnitude is considerably reduced, which provides a potential strategy to mitigate and/or suppress the instabilities.
The microbe, Zymomonas mobilis, can efficiently convert sugars to 2,3 butanediol (BDO), an important intermediate for downstream chemical products, only in a well-controlled microaerated environment. However, controlling oxygen distribution in industrial-scale bioreactors is challenging, and at-scale BDO production is hard to optimize using traditional engineering methods. This study takes a step towards addressing this problem through a computational model for reacting multiphase flows in large-scale bioreactors. A phenomenological metabolic model was first developed and validated against experiments, and then it was coupled to a multiphase computational fluid dynamics (CFD) solver using a subcycling algorithm for simulating long conversion times (~30 hours). Large-scale (~500 m 3 ) bubble column simulations using this coupled approach demonstrated a 25% improvement in BDO yield with low sparging rates (0.002 m/s) and low height-to-diameter ratio (0.875) compared to a baseline bubble column (35 m height, 5 m diameter) with a superficial gas velocity of 0.06 m/s.