SEARCH · Search NASA
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.
The dynamics of a one-parameter family of mesh equidistribution schemes coupled with finite difference discretisations of linear and nonlinear convection-diffusion model equations is studied numerically. It is shown that, when time marched to steady state, the grid adaption not only influences the stability and convergence rate of the overall scheme, but can also introduce spurious dynamics to the numerical solution procedure.
A completely coupled one-dimensional radiative-convective (RC) and photochemical-diffusion (PC) model has been developed recently and used to study the climate-chemical interactions. The importance of radiative-chemical interactions within the troposphere and stratosphere has been examined in some detail. We find that increases of radiatively and/or chemically active trace gases such as CO2, CH4 and N2O have both the direct effects and the indirect effects on climate change by changing the atmospheric O3 profile through their interaction with chemical processes in the atmosphere. It is also found that the climatic effect of ozone depends strongly on its vertical distribution throughout the troposphere and stratosphere, as well on its column amount in the atmosphere.
The primary objective is to understand how buoyancy affects the structure of the shear layer, the development of fluid dynamic instabilities, and formation of the coherent structures in the near-nozzle regions of gas jets. The secondary objectives are to study the role of buoyancy in lifting and reattachment process of diffusion flames, to evaluate the scaling behavior of diffusion flames, and to aid development and/or validation of theoretical models by providing quantitative data in the absence of buoyancy. Fast reacting hydrogen or hydrogen-inert fuels are used to isolate the effects of buoyancy on fluid dynamics without masking the flame behavior by soot and radiative heat transfer. This choice of fuel also permits an evaluation of simulating low gravity in low pressure ground experiments because the similarity constraints are relaxed for the fast reacting, nonsooting diffusion flames. The diagnostics consists primarily of a color schlieren system coupled with computer generated rainbow filters, video recording, and image analysis. The project involves (1) drop tower experiments, (2) ground experiments, and (3) theoretical analysis.
Hadronization models based on local string-breaking dynamics are typically Markovian by construction, yet the physical ensemble of final states is shaped by global constraints that couple the entire fragmentation trajectory. Recasting hadronization as a conditioned stochastic diffusion process provides a precise mathematical resolution to this tension. In particular, this language reveals explicitly that constraints stemming from conservation laws induce non-Markovian correlations between otherwise independent fragmentation steps, and that these correlations can be absorbed exactly into a renormalization of the local dynamics through a Doob $h$-transform. We develop this formalism for a $q\bar{q}$ string in the chiral limit, where the longitudinal-transverse factorization of the Lund kernel becomes exact, enabling systematic power counting and clean ultraviolet (UV)/infrared (IR) separation. The dynamics organize naturally into a tower of effective theories distinguished by the remaining string mass, spanning a UV fixed point with scale-invariant transport coefficients, an intermediate regime where transverse phase space induces controlled running, and an IR boundary layer where non-local effects enter at leading order. The tower exhibits genuine Wilsonian structure, including $β$-functions, anomalous dimensions, and systematic matching conditions. The resulting framework achieves a clean factorization of universal microscopic fragmentation dynamics from infrared constraint effects, and opens new directions for both the theoretical analysis and practical simulation of hadronization.
An experimental evaluation of the effects of free-stream turbulence on the performance of a subsonic two-dimensional diffuser has been made. Increases of the diffuser's static pressure recovery coefficient of 11.3 and 23.9 percent at total included divergence angles of 12 and 20 degrees respectively were obtained when the value of the inlet integral free-stream scale of turbulence in the flow direction was at least 7.2 times larger than the inlet boundary layer displacement thickness, when the inlet total free-stream turbulence intensity was at least 3.5 percent, and when the axes of upstream rods used to generate turbulence were perpendicular to the flow and parallel to the diverging walls of the diffuser. It is hypothesized that a larger scale of turbulence with the specified eddy axis orientation transmits the free-stream energy to the walls more effectively and, when coupled with large turbulence intensities, are mechanisms which act to decrease the distortion and delay separation within the diffuser.
Unraveling collective modes arising from coupled degrees of freedom is crucial for understanding complex interactions in solids and developing new functionalities. Unique collective behaviors emerge when two degrees of freedom, ordered on distinct length scales, interact. Polar skyrmions, three-dimensional electric polarization textures in ferroelectric superlattices, disrupt the lattice continuity at the nanometer scale with nontrivial topology, leading to previously unexplored collective modes. Here, using terahertz-field excitation and femtosecond x-ray diffraction, we discover subterahertz collective modes, dubbed “skyrons”, which appear as swirling patterns of atomic displacements functioning as atomic-scale gearsets. The key to activating skyrons is the use of the THz field that couples primarily to skyrmion domain walls. Momentum-resolved time-domain measurements of diffuse scattering reveal an avoided crossing in the dispersion relation of skyrons. Atomistic simulations and dynamical phase-field modeling provide microscopic insights into the three-dimensional crystallographic and polarization dynamics. The amplitude and dispersion of skyrons are demonstrated to be controlled by sample temperature and electric-field bias. The discovery of skyrons and their coupling with terahertz fields opens avenues for ultrafast control of topological polar structures.
Understanding the physical phenomena controlling the ignition and spread of flames in microgravity has importance for space safety as well as for characterizing dynamical and chemical combustion processes which are normally masked by buoyancy and other gravity-related effects. Unfortunately, combustion is highly complicated by fluid mechanical and chemical kinetic processes, requiring the use of numerical modeling to compare with carefully designed experiments. More sophisticated diagnostic methods are needed to provide the kind of quantitative data necessary to characterize the properties of microgravity combustion as well as provide accurate feedback to improve the predictive capabilities of the models. Diode lasers are a natural choice for use under the severe conditions of low gravity experiments. Reliable, simple solid state operation at low power satisfies the operational restrictions imposed by drop towers, aircraft and space-based studies. Modulation wavelength absorption spectroscopy (WMS) provides a means to make highly sensitive and quantitative measurements of local gas concentration and, in certain cases, temperature. With near-infrared diode lasers, detection of virtually all major combustion species with extremely rapid response time is possible in an inexpensive package. Advancements in near-infrared diode laser fabrication technology and concurrent development of optical fibers for these lasers led to their use in drop towers. Since near-infrared absorption line strengths for overtone and combination vibrational transitions are weaker than the mid-infrared fundamental bands, WMS techniques are applied to increase detection sensitivity and allow measurement of the major combustion gases. In the first microgravity species measurement, Silver et al. mounted a fiber-coupled laser at the top of the NASA 2.2-sec drop tower and piped the light through a single-mode fiber to the drop rig. A fiber splitter divided the light into eight channels that directed the laser beam across a methane or propane diffusion jet flame. The light beams were recaptured by a set of gradient index lenses, coupled back into separate fiber optic lines, and transmitted back to detectors and electronics in the instrument package. In these experiments a 6-mm od fiber cable (containing the nine optical fibers) fell with the drop rig. Using separate detection and demodulation channels, spatial and temporal (up to 20 Hz) maps of water vapor and methane concentrations were obtained at differing heights in the flames. While this apparatus was useful from a demonstration standpoint, several drawbacks needed attention before useful scientific measurements could be obtained. First, eight lines of sight are somewhat insufficient for detailing the spatial profiles of the gas. Second, multiple detection channels operating in parallel are both expensive and present a challenge for accurate calibration. As a result, a newer scanning system was developed in our first contract under this program. The primary characteristic of this system is that it contains a single detection channel and achieves "continuous" spatial resolution by scanning the laser beam across the flame region, then directing this beam onto a single detector. Thus spatial measurements are converted to a temporal series of data. The true spatial resolution is limited only by the beam diameter and width of the sweep. In these experiments the beam is focused to about 1-mm diameter and scans across a region up to 4-cm wide.
Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.
Human cognitive responses, behavioral responses, and disease dynamics co-evolve over the course of any disease outbreak, and can result in complex feedbacks. We present a dynamic agent-based model that explicitly couples the spread of disease with the spread of fear surrounding the disease, implemented within the EpiCast simulation framework. EpiCast models transmission within a realistic synthetic population, capturing individual-level interactions. In our model, fear propagates through both in-person contact and broadcast media, prompting individuals to adopt protective behaviors that reduce disease spread. In order to better understand these coupled dynamics, we create and compare a range of compartmental models to ensure that introducing additional disease states does not prevent the emergence of multiple waves in these simpler models. Additionally, we compare a range of behavioral scenarios within EpiCast, varying the level and intensity of fear and behavior change. Our results show that the addition of asymptomatic, exposed, and pre-symptomatic disease states can impact both the rate at which an outbreak progresses and its overall trajectory in compartmental models. In EpiCast, the combination of non-local fear spread via broadcasters and strong behavioral responses by fearful individuals generally leads to multiple epidemic waves, an outcome that occurs only within a narrow parameter range when fear spreads purely through local contact. Accounting for the coupled spread of fear and disease is critical for understanding disease dynamics and designing timely, targeted responses to emerging infectious threats.
In this work, we present validation test results of fully coupled neutronics and thermal-hydraulics models of the Molten Salt Reactor Experiment (MSRE) against experimental data of the zero power pump transients and the natural circulation tests at low power. To capture the strong coupling between neutronics and thermal-hydraulics due to fuel circulation, and to account for the delayed neutron precursor (DNP) distribution, the porous media thermal-hydraulics solver Pronghorn was fully coupled to the spatial neutron dynamics code Griffin, which solves the neutron diffusion equation, and to the 0-D point kinetics solver Squirrel, using a 2-D homogenized representation of the MSRE. The validation test results show very good agreement with experimental data for both point kinetics and spatial dynamics simulations, capturing the strong feedback effect and DNP losses in the MSRE. The 0-D code Squirrel accurately predicted the time-dependent behavior in the MSRE given the steady-state spatial dynamics solution of Griffin.
Directly pressurizing seawater for desalination with reverse osmosis membranes via wave motion is a promising and sustainable method for producing freshwater in coastal regions. However, such a system could result in significant pressure fluctuations and a departure from conventional steady-state desalination operations. This study sought to assess if membrane transport parameters (apparent water and salt permeability) should be modeled as transient or constant in solution–diffusion-based modeling efforts of dynamically operated desalination systems, such as those coupled to wave power. Two approaches were used to model membrane transport parameters: one considered each parameter to be a function of the net driving pressure of the system, and the other assumed they were constant across all conditions. A pilot-scale system was used to conduct steady-state and controlled ramping experiments. Data from steady-state experiments were used to calculate transient and constant transport parameters. Parameter combinations were used in a simulation model to predict water flux and effective permeate salinity, and simulation outcomes were compared against experimental ramping results. The transient relationships for both water and salt permeability produced the most accurate results for water flux and comparable results for effective permeate salinity. Development of such relationships would be unique to a specific system but could be valuable in modeling wave-driven desalination systems across the wide range of operating conditions they experience.
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 mathematical model of isoelectric focusing at the steady state has been developed for an M-component system of electrochemically defined ampholytes. The model is formulated from fundamental principles describing the components' chemical equilibria, mass transfer resulting from diffusion and electromigration, and electroneutrality. The model consists of ordinary differential equations coupled with a system of algebraic equations. The model is implemented on a digital computer using FORTRAN-based simulation software. Computer simulation data are presented for several two-component systems showing the effects of varying the isoelectric points and dissociation constants of the constituents.
The efficiency with which heat may be converted into resonance radiation in a cesium thermionic diode is investigated theoretically. An analytical model of a thermionic converter is used which combines the coupled effects of line radiation transport, excited-state kinetics, and plasma diffusion. Operating regimes are established for various degrees of optical density in the plasma. The results indicate that monochromatic radiation can be produced with efficiencies on the order of 30 percent, provided that there is an adequate voltage drop across the plasma. A drop of 1 V was used since it can be maintained without any electrical power input to the device. It is found that high efficiencies are ude to the higher interelectrode distances which the solutions accommodate, and that radiation can be generated efficiently, even with optically dense gases.
Teflon has been identified as suitable material for use in oxygen-enriched atmospheres because of its low specific enthalpy of combustion that is, less than 1500 cal/gram. However, once ignited, Teflon burns in a heterogeneous reaction until total consumption or depletion of oxygen occurs. A model is developed for the burning of Teflon in pure, high-pressure oxygen, 3.4 to 68.9 MPa. The Teflon polymer chain dissociates via monomer units due to pyrolysis. These monomer units diffuse to the surface due to free convection. The model consists of coupled mass and energy balances for the polymer and an energy balance for the free convection of product gases. The model is used to obtain appropriate kinetic parameters for the dissociation and surface reactions. The model is validated against experimental measurements of Teflon discs supported on aluminum rods (2216) at oxygen pressures of 3.4, 6.8, 34.4, and 68.9 MPa. The model simulates the temperature distribution in the product gases above burning Teflon.
One of the main challenges in modeling living systems is to distinguish a random walk of physical origin (for instance, Brownian motions) from those of biological origin and that will constitute the starting point of the proposed approach. As conjectured, the biological random walk must be nonlinear. Indeed, any stochastic Markov process can be described by linear Fokker-Planck equation (or its discretized version), only that type of process has been observed in the inanimate world. However, all such processes always converge to a stable (ergodic or periodic) state, i.e., to the states of a lower complexity and high entropy. At the same time, the evolution of living systems directed toward a higher level of complexity if complexity is associated with a number of structural variations. The simplest way to mimic such a tendency is to incorporate a nonlinearity into the random walk; then the probability evolution will attain the features of diffusion equation: the formation and dissipation of shock waves initiated by small shallow wave disturbances. As a result, the evolution never "dies:" it produces new different configurations which are accompanied by an increase or decrease of entropy (the decrease takes place during formation of shock waves, the increase-during their dissipation). In other words, the evolution can be directed "against the second law of thermodynamics" by forming patterns outside of equilibrium in the probability space. Due to that, a specie is not locked up in a certain pattern of behavior: it still can perform a variety of motions, and only the statistics of these motions is constrained by this pattern. It should be emphasized that such a "twist" is based upon the concept of reflection, i.e., the existence of the self-image (adopted from psychology). The model consists of a generator of stochastic processes which represents the motor dynamics in the form of nonlinear random walks, and a simulator of the nonlinear version of the diffusion equation which represents the mental dynamics. It has been demonstrated that coupled mental-motor dynamics can simulate emerging self-organization, prey-predator games, collaboration and competition, "collective brain," etc.
The specific anisotropic extinction coefficient e* couples the effective Radiative properties of a fibrous insulation into the radiation diffusion equation. This coefficient can be calculated using Mie scattering theory if fiber diameters, refractive indices and fiber orientation distributions are known. In general, fiber orientation distributions are not readily accessible and past calculations have considered fibers as either randomly distributed or normal to the heat flow direction. In certain rigid fibrous ceramic insulations neither of these cases apply well and a simple procedure is described for approximating e* from values calculated for the random and normal orientation cases. The intrinsic error associated with this scaling procedure is investigated. Numerical computations for several test structures and fiber materials show the average error to be less than 5% for net heat flux and radiation conductivity calculations.
When targeting physical understanding of space plasmas, our focus is gradually shifting away from discovery-type investigations to missions and studies that address our basic understanding of processes we know to be important. For these studies, theory and models provide physical predictions that need to be verified or falsified by empirical evidence. Within this paradigm, a tight integration between theory, modeling, and space flight mission design and execution is essential. NASA's Magnetospheric MultiScale (MMS) mission is a pathfinder in this new era of space research. The prime objective of MMS is to understand magnetic reconnection, arguably the most fundamental of plasma processes. In particular, MMS targets the microphysical processes, which permit magnetic reconnection to operate in the collisionless plasmas that permeate space and astrophysical systems. More specifically, MMS will provide closure to such elemental questions as how particles become demagnetized in the reconnection diffusion region, which effects determine the reconnection rate, and how reconnection is coupled to environmental conditions such as magnetic shear angles. Solutions to these problems have remained elusive in past and present spacecraft missions primarily due to instrumental limitations - yet they are fundamental to the large-scale dynamics of collisionless plasmas. Owing to the lack of measurements, most of our present knowledge of these processes is based on results from modern theory and modeling studies of the reconnection process. Proper design and execution of a mission targeting magnetic reconnection should include this knowledge and have to ensure that all relevant scales and effects can be resolved by mission measurements. The SMART mission has responded to this need through a tight integration between instrument and theory and modeling teams. Input from theory and modeling is fed into all aspects of science mission design, and theory and modeling activities are tailored to SMART needs during mission development and science analysis. In this presentation, we will present an overview of SMART theory and modeling team activities. In particular, we will provide examples of science objectives derived from state-of-the art models, and of recent research results that continue to be utilized in SMART mission development.