Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional diffusion models”

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 325 records · Page 18

Theory and observation of a dynamically evolving negative ion plasma

The study described here examines the full range of negative-ion chemistry in the upper ionosphere by using current reaction rate data to investigate the many chemical paths SF6-type injections might take in an F region environment. Special attention is given to the conditions required to create heavy negative ions that persist long enough to affect the dynamical properties of the F region. The ambipolar diffusion characteristics of a three-component plasma are described, and estimates are given of the incoherent scatter spectra obtained from such a plasma. Model calculations using a first-order chemical code are defined and tested in order to investigate the actual types of negative-ion plasma capable of being created under nighttime conditions.

Mendillo, M.↗

Analytic Couple Modeling Introducing Device Design Factor, Fin Factor, Thermal Diffusivity Factor, and Inductance Factor

A set of convenient thermoelectric device solutions have been derived in order to capture a number of factors which are previously only resolved with numerical techniques. The concise conversion efficiency equations derived from governing equations provide intuitive and straight-forward design guidelines. These guidelines allow for better device design without requiring detailed numerical modeling. The analytical modeling accounts for factors such as i) variable temperature boundary conditions, ii) lateral heat transfer, iii) temperature variable material properties, and iv) transient operation. New dimensionless parameters, similar to the figure of merit, are introduced including the device design factor, fin factor, thermal diffusivity factor, and inductance factor. These new device factors allow for the straight-forward description of phenomenon generally only captured with numerical work otherwise. As an example a device design factor of 0.38, which accounts for thermal resistance of the hot and cold shoes, can be used to calculate a conversion efficiency of 2.28 while the ideal conversion efficiency based on figure of merit alone would be 6.15. Likewise an ideal couple with efficiency of 6.15 will be reduced to 5.33 when lateral heat is accounted for with a fin factor of 1.0.

Seebeck effect↗

Triso Analysis Tool For Predictive Source Terms

Source term modeling for TRi-structural ISOtropic (TRISO) fuel has been performed for previous reactor designs, but few are available in the open literature. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that can be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. To meet this need, the TRISO Analysis Tool for Predictive Source terms (TRISO-ATOPS) was developed. This model calculates the release of the key safety-important fission products by diffusion through the kernel, silicon carbide and graphite based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system for removing fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. From this initial condition, the model calculates the fission product release for any transient temperature profile, and the fission product releases can then be used to assess radiological dose to the workers and the public using conventional radiological dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy Advanced Gas Reactor TRISO fuel development program. The example cases in this work demonstrate the flexibility of the model

Stoyer, Benjamin [Idaho National Laboratory (INL),↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

Combined Experimental and Modeling Study of the Interactions of Acid Gas with Common Spacecraft Surfaces for Fire Safety Applications

A fire in a spacecraft poses detrimental consequences and risks mission success in addition to crew safety. This is compounded during long-duration missions when the crew has limited options to recover from a fire. A common spacecraft fire concern is the smoldering of wire insulation, typically made from Polyvinyl chloride (PVC) or Polytetrafluoroethylene (PTFE). This creates acid gases such as Hydrogen Chloride (HCl), Hydrogen Fluoride (HF) and Hydrogen Cyanide (HCN). These poisonous gases are hazardous to the crew. They also interact with common surfaces within the spacecraft more than dominant combustion products such as CO2 and H2O. This makes them more difficult to track for potential fire detection techniques, or for postfire clean-up. It is imperative to be able to understand and predict the fate of these poisonous species in a microgravity environment in order to design a safe vehicle. HCl interacts with a number of materials inside a spacecraft. Primary among these materials is aluminum, which is abundantly used due to its strong and light weight nature. Aluminum has a natural oxide layer that protects it from corrosion but is typically treated to enhance this oxide layer. Among these treatments is a chromate conversion coating (CCC), which provides a thin enough protective oxide layer to still conduct electricity, and a traditional anodized material that has a thicker oxide layer that does not conduct electricity. Nomex is another common material found inside a spacecraft. It is a flame-resistant woven polymer that is related to nylon. This commercially available material is used for cargo storage bags and as a fire barrier. Physics-based models were developed to predict the uptake of HCl by these materials. The ultimate objective of these models is to predict the fate of HCl within the spacecraft so that sensors can be placed in meaningful locations in future missions based on the model predictions. To support these modeling efforts, experiments were performed in a cast acrylic test cell that measured the difference between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. Different uptake capacities were realized for each type of sample tested. A computational fluid dynamics model (CFD) model of the reactor was then constructed that used a one-step global reaction rate with calibratable reaction (or kinetic) constants. These constants were calibrated to match the HCl uptake on the CCC aluminum samples, and the same kinetic constants were then tested for the stock and anodized aluminum samples. Model predictions matched the experimental data for the stock aluminum, and to a much lesser extent, the anodized aluminum. The model was additionally validated at different flow rates, sample surface areas, and inlet concentrations, and showed good agreement for all stock and CCC samples. The model did not accurately predict the HCl uptake in the anodized samples compared to the other two types of aluminum. Adjusting the kinetic constants and transport properties did little to improve the prediction. X-Ray Photoelectron Spectroscopy (XPS) was used to determine that the oxide layer thickness of anodized aluminum is approximately 5,000 nm, compared to 250 nm for CCC and 50 nm for stock. XPS also revealed presence of chlorine further down in the aluminum oxide layer in anodized samples than CCC and stock samples after the samples were saturated with HCl, indicating that accounting for diffusion of HCl into the oxide layer is important for accurate prediction of HCl uptake onto anodized aluminum. Consequently, a multi-scale model was developed and tested. First, a single pore inside the anodized aluminum oxide layer was modeled and is referred to as the pore-scale model. In this model, HCl diffused through the pore and reacted with the aluminum oxide pore wall to create aluminum chloride. The sample was then saturated when the mass transfer resistance through the growing aluminum chloride layer became too large for the HCl to reach the aluminum oxide wall and continue the reaction. This pore-scale model was coupled to the reactor-scale model using a concentration-dependent diffusion coefficient, resulting in much more accurate predictions (approximately half the sum square error of the aforementioned reactor-scale model that produced good agreement for stock and CCC) for a variety of operating conditions. The amount of water vapor or relative humidity (RH) in the flow during a reactor experiment was determined to influence HCl uptake. Experiments were performed to understand the interaction of gaseous HCl with aluminum surfaces in the presence of water vapor. The results show that increasing levels of RH increased the capacity of aluminum to adsorb HCl but decreased the capacity of Nomex to uptake HCl. A series of tests were performed on individual aluminum samples after they had been saturated with a fixed concentration of HCl in dry air conditions with the goal of determining how their HCl uptake capacity changes after various treatments with water relative to the original saturation tests. HCl-saturated aluminum samples subjected to a second dry air flow at the same HCl concentration as the original test had an uptake of 23.5% of the original sample with no treatment in between. Saturated aluminum samples subjected to an air flow with a RH of 90% in between tests had an uptake of 35.6% of the original. Saturated aluminum samples submerged in distilled water for 12 hours in between tests had an uptake of 82.2% of the original sample. Previously saturated aluminum tested with HCl and a 50% RH air flow resulted in similar uptake characteristics in multiple repeated tests. The experiments show the profound effect water vapor has on HCl uptake onto aluminum surfaces. In the samples subjected to water vapor or liquid water, capillary condensation and capillary diffusion alters the transport of HCl significantly. A model was proposed that developed a relationship between RH and the coefficient of HCl diffusion in aluminum chloride. This produced an “S-shaped” curve with diffusion coefficient as a function of RH, with 45% RH represented as the point where the diffusion coefficient is halfway between no water saturation and 100% water saturation in the aluminum chloride product layer. No difference in uptake characteristics for the experiment or model were realized between 50% and 62% RH. The results from the large-scale microgravity experiment, Saffire, are discussed as they pertain to the fate of HCl throughout a spacecraft. HCl was released, both as a standalone event, and in concurrence with the burning of a structured cloth. These events only produced a small response in the far field HCl sensor, while a PMMA burn that did not produce HCl had a significantly greater response. A ground-based large-scale facility was constructed to flow acid gas at the scale and configuration realized in the Saffire experiments. A CFD model of this duct was constructed to test kinetic parameters developed in this work at a larger scale and different geometric configuration and to predict the results of the large-scale facility. The models developed in this work were used to interpret the results of the microgravity tests and lead the discussion on what further experiments and models are needed in order to predict the fate of acid gas in a spacecraft environment. To summarize, the major contributions of this work are as follows: the capacity to uptake HCl, with and without the presence of water vapor, was measured for a variety of real spacecraft surfaces. Several different models (single reactor-scale, multiscale, spacecraft-scale) were developed and with the aid of modeling, the rate of uptake for those surfaces was also predicted and validated. The kinetic parameters determined from the small-scale reactor experiments and models were used to predict large-scale and microgravity tests. Conclusions from this research will be used in the design of spacecraft vehicles and large-scale microgravity fire safety experiments. The models built by this work will aid designers in sensor placement and could be used to predict acid gas transport from fires in partial gravity, as would be seen in Lunar and Martian habitats.

fire safety↗

Cosmic-ray gradients in the outer heliosphere

We report on a study of the expected spatial variation of the galactic cosmic-ray intensity in the vicinity of the solar-wind termination shock. Model simulations predict that the radial gradients change abruptly at the shock, and that the nature of the effect varies significantly with particle energy. At low energies, the radial gradient changes abruptly from a lower value inside the shock to a higher value outside, whereas at high energies, the higher value of the gradient is inside the shock. This effect, which is a consequence of the matching conditions at the shock and is closely related to diffusive shock acceleration, is qualitatively the same for both heliospheric magnetic polarity states and remains much the same in one-dimensional, two-dimensional and three-dimensional models. Hence drifts do not change the nature of this phenomenon, although they change it quantitatively. The effect may prove to be an important diagnostic tool for the study of the termination of the solar wind.

Kota, J.↗

Parametric studies with an atmospheric diffusion model that assesses toxic fuel hazards due to the ground clouds generated by rocket launches

Parametric studies were made with a multilayer atmospheric diffusion model to place quantitative limits on the uncertainty of predicting ground-level toxic rocket-fuel concentrations. Exhaust distributions in the ground cloud, cloud stabilized geometry, atmospheric coefficients, the effects of exhaust plume afterburning of carbon monoxide CO, assumed surface mixing-layer division in the model, and model sensitivity to different meteorological regimes were studied. Large-scale differences in ground-level predictions are quantitatively described. Cloud alongwind growth for several meteorological conditions is shown to be in error because of incorrect application of previous diffusion theory. In addition, rocket-plume calculations indicate that almost all of the rocket-motor carbon monoxide is afterburned to carbon dioxide CO2, thus reducing toxic hazards due to CO. The afterburning is also shown to have a significant effect on cloud stabilization height and on ground-level concentrations of exhaust products.

Stewart, R. B.↗

Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Quenching

For about a half century, the stabilization of a steady planar deflagration on a heat-sink-type flat-flame burner has been of extraordinary service for the theoretical modeling and diagnostic probing of combusting gaseous mixtures. However, most engineering devices and most unwanted fire involve the burning of initially unmixed reactants. The most vigorous burning of initially separated gaseous fuel and oxidizer is the diffusion flame. In this useful idealization (limiting case), the reactants are converted to product at a mathematically thin interface, so no interpenetration of fuel and oxidizer occurs. This limit is of practical importance because it often characterizes the condition of optimal performance (and sometimes environmentally objectionable operation) of a combustor. A steady planar diffusion flame is most closely approached in the laboratory in the counterflow apparatus. The utility of this simple-strain-rate flow for the modeling and probing of diffusion flames was noted by Pandya and Weinberg 35 years ago, though only in the last decade or so has its use become internationally common place. However, typically, as the strain rate a is reduced below about 20 cm(exp -1), and the diffusion-flame limit (reaction rate much faster than the flow rate) is approached, the burning is observed to become unstable in earth gravity. The advantageous steady planar flow is not available in the diffusion-flame limit in earth gravity. This is unfortunate because the typical spatial scale in a counterflow is (k/a)(sup 1/2), where k denotes a characteristic diffusion coefficient; thus, the length scale becomes large, and the reacting flow is particularly amenable to diagnostic probing, as the diffusion-flame limit is approached. The disruption of planar symmetry is owing the fact that, as the strain rate a decreases, the residence time (l/a) of the throughput in the counterflow burner increases. Observationally, when the residence time exceeds about 50 msec, the inevitably present convective (Rayleigh-Benard) instabilities, associated with hot-under-cold (flame-under-fresh-reactant) stratification of fluid in a gravitational field, have time to grow to finite amplitude during transit of the burner.

Fendell, Francis↗

Elucidating the role of Cr migration in Ni-Cr exposed to molten FLiNaK via multiscale characterization

Structural materials used in nuclear reactor environments are subjected to coupled extremes such as high temperature, irradiation, and corrosion which act in concert to degrade their functional performance. Connecting alloy microstructure such as grain boundaries, and accumulating point defects with corrosion attack and pore morphology is critical to understanding underlying mechanisms. We uncover the compositional variations and morphology at multiple length scales in corrosion-damaged Ni-Cr alloy after exposure to oxidants in molten fluoride salts. A complex network of dense corrosion pores is detected by surface-level SEM observations. The corrosion pores take on a 1-dimensional morphology and are enriched with Ni and depleted of Cr 1–5 µm from the pore surface. STEM-EDS and 4D-STEM strain maps acquired simultaneously highlight the local variations in composition and structure of a ≤ 200 nm Cr-rich layer identified from a cross-section taken at the bottom of an isolated corrosion pore between the Ni-Cr alloy matrix and the embedded salt. However, the absence of an observed interface between the Ni-Cr alloy matrix and the FCC-structured Cr-rich layer suggests that Cr plating from the salt did not transpire. These findings support a proposed Cr lattice diffusion mechanism rather than Cr-precipitation from the salt to accommodate temperature transient conditions during sample cooling. Through the development of a 1D phase field model, these results are rationalized by formation energies for the Ni- and Cr-oxidation into the molten salt. This study reveals the locally altered microstructure caused by high temperature corrosion in non-steady-state molten salt nuclear reactor environments.

36 MATERIALS SCIENCE↗

Non-Invasive Electrode Impedance Estimation for Optimized Charge Profile Parameterization of Lithium-Ion Batteries

This work presents a non-invasive method for parameterizing a physically motivated equivalent circuit model of lithium-ion batteries using operando electrochemical impedance spectroscopy and time-domain data. The proposed model consists exclusively of linear circuit elements, enabling computationally efficient simulation and real-time implementation on battery management system chips. By integrating frequency- and time-domain measurements, the model accurately estimates internal states such as the negative electrode potential, achieving a root mean square error of 12.3 mV during fast charging validation. Parameterization requires only rate tests with sinusoidal perturbations at three different ambient temperatures, making the approach experimentally accessible. The model reveals key insights into battery behavior, including rate-dependent overpotentials primarily governed by charge transfer kinetics at the positive electrode, and temperature-dependent impedance contributions from both charge transfer and solid-state diffusion processes. Validation using reference electrodes confirms the model’s ability to detect lithium plating onset and reproduce impedance behavior across a wide range of operating conditions. The approach enables in situ optimization of fast charging profiles and lays the foundation for future extensions incorporating aging effects and plating dynamics.

25 ENERGY STORAGE↗

Data-Driven Method for Groundwater-Level Mapping and Monitoring-Well Network Optimization at Hanford

This report summarizes the initial results and outcomes of a physics-informed, data-driven groundwater level (GWL) mapping capability for the Hanford Site. GWL mapping at Hanford is typically conducted annually and requires a significant amount of computational and expert resources, and it does not allow assessment of the informational value of specific monitoring wells. The proposed method produces spatially and temporally resolved fields consistent with sparse, irregularly sampled, and nonuniformly distributed well measurements. Implemented successfully, this capability will allow rapid mapping of groundwater levels and provide an opportunity to optimize monitoring activities (both location and sampling frequency) based on data information value evaluation. The approach integrates a diffusion-based generative model – trained on MODFLOW simulation data from the Plateau-to-River (P2R) model – with score-based data assimilation (SDA), allowing observation-conditioned mapping without retraining for each monitoring-network layout.

54 ENVIRONMENTAL SCIENCES↗

Experimental measurements of the ground cloud effluents and cloud growth for the May 20, 1975, Titan 3C launch at Air Force Eastern Test Range, Florida

The experiment included surface level and airborne in situ cloud measurements of the exhaust effluents from the Titan IIIC solid rocket boosters. Simultaneous visible spectrum photographic pictures of the ground cloud as well as infrared imaging of the cloud were obtained to study the cloud rise, growth, and direction of travel within the earth's surface mixing layer. The NASA multilayer diffusion model predictions of cloud growth, direction of travel, and expected surface level effluent concentrations were made prior to launch and after launch using measured meteorological conditions. Prelaunch predictions were used to position the effluent monitoring instruments, and the postlaunch predictions were compared with the measured data. Measurement results showed that surface level effluent values were low, often below the detection limits of the instrumentation. The maximum surface level hydrogen chloride concentration measured 50 parts per billion at about 8 km from the launch pad. The maximum observed in-cloud (airborne measurement) hydrogen chloride concentration was 7 per million.

Gregory, G. L.↗

Ground based experiments on the growth and characterization of L-Arginine Phosphate (LAP) crystals

L-Arginine Phosphate (LAP) is a new nonlinear optical material with higher efficiency for harmonic generation compared to KDP. Crystals of LAP were grown in the laboratory from supersaturated solutions by temperature lowering technique. Investigations revealed the presence of large dislocation densities inside the crystals which are observed to produce refractive index changes causing damage at high laser powers. This is a result of the convection during crystal growth from supersaturated solutions. It is proposed to grow these crystals in a diffusion controlled growth condition under microgravity environment and compare the crystals grown in space with those grown on ground. Physical properties of the solutions needed for modelling of crystal growth are also presented.

Rao, S. M.↗

Turbulent Jet Flames Into a Vitiated Coflow

Examined is the vitiated coflow flame, an experimental condition that decouples the combustion processes of flows found in practical combustors from the associated recirculating fluid mechanics. The configuration consists of a 4.57 mm diameter fuel jet into a coaxial flow of hot combustion products from a lean premixed flame. The 210 mm diameter coflow isolates the jet flame from the cool ambient, providing a hot environment similar to the operating conditions of advanced combustors; this important high temperature element is lacking in the traditional laboratory experiments of jet flames into cool (room) air. A family of flows of increasing complexity is presented: 1) nonreacting flow, 2) all hydrogen flame (fuel jet and premixed coflow), and 3) set of methane flames. This sequence of experiments provides a convenient ordering of validation data for combustion models. Laser Raman-Rayleigh-LIF diagnostics at the Turbulent Diffusion Flame laboratory of Sandia National Laboratories produced instantaneous multiscalar point measurements. These results attest to the attractive features of the vitiated coflow burner and the well-defined boundary conditions provided by the coflow. The coflow is uniform and steady, isolating the jet flame from the laboratory air for a downstream distance ranging from z/d = 50-70. The statistical results show that differential diffusion effects in this highly turbulent flow are negligible. Complementing the comprehensive set of multiscalar measurements is a parametric study of lifted methane flames that was conducted to analyze flame sensitivity to jet and coflow velocity, as well as coflow temperature. The linear relationship found between the lift-off height and the jet velocity is consistent with previous experiments. New linear sensitivities were found correlating the lift-off height to coflow velocity and temperature. A blow-off study revealed that the methane flame blows off at a common coflow temperature (1260 K), regardless of coflow or jet velocity. An explanation for this phenomenon is that entrainment of ambient air at the high lift-off heights prevents autoignition. Analysis of the results suggests that flame stabilization occurs through a combination of flame propagation, autoignition, and localized extinction processes. Proposed is an expanded view of distributed reaction combustion based on analysis of the distributions of probe volume conditions at the stabilization region of the lifted hydrogen and methane flames. Turbulent eddies the size of the flame thickness mix fuel and hot coflow across the flame front, thereby enhancing the reaction zone with autoignition of reactants at elevated temperatures; this is the reverse effect of turbulent flames in ambient air, where intense turbulence in cool mixtures result in localized extinction. Each of the three processes (i.e., flame propagation, autoignition and localized extinction) contributes to flame stabilization in varying degrees, depending on flow conditions.

Holdeman, James D.↗

The role of eddy turbulence in the development of self-consistent models of the lower and upper thermosphere

Numerical solutions of mutually coupled time dependent equations of continuity, momentum and energy balance are presented to illustrate the effect of eddy turbulence on the neutral composition and temperature of the lower and upper atmosphere. The procedure adopted eliminates the necessity of making ad-hoc assumptions about the conditions at the turbopause level and allows the study of the two regions in a self-consistent manner. From the illustrative examples comprising parametric changes in the eddy diffusion coefficient, the specific roles of eddy turbulence in the development of theoretical models of the thermosphere are discussed.

Chandra, S.↗

The gradient of galactic cosmic rays at the solar-wind termination shock

We report on a study of the expected spatial variation of the galactic cosmic-ray intensity in the outer heliosphere, in the vicinity of the solar-wind termination shock. Model simulations which contain all of the transport effects, including drifts, predict that the radial gradients change abruptly at the shock, and that the nature of the effect varies significantly with particle energy. At low energies, the radial gradient changes abruptly from a lower value inside the shock to a higher value outside, whereas at high energies, the higher value of the gradient is inside the shock. This effect, which is a consequence of the matching conditions at the shock and is closely related to diffusive shock acceleration, is qualitatively the same for both helisopheric magnetic polarity states and remains much the same in one-dimensional, two-dimensional, and three-dimensional models. Hence drifts do not change the nature of this phenomenon, although they change it quantitatively. The effect is, of course, not present in the absence of a terminal shock and may prove to be an important diagnostic tool for the study of the termination of the solar wind.

Jokipii, J. R.↗

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗

Physiological Response of Plants Grown on Porous Ceramic Tubes

This research involves the manipulation of the root-zone water potential for the purposes of discriminating the rate limiting step in the inorganic nutrient uptake mechanism utilized by higher plants. This reaction sequence includes the pathways controlled by the root-zone conditions such as water tension and gradient concentrations. Furthermore, plant based control mechanisms dictated by various protein productions are differentiated as well. For the nutrients limited by the environmental availability, the kinetics were modeled using convection and diffusion equations. Alternatively, for the nutrients dependent upon enzyme manipulations, the uptakes are modeled using Michaelis-Menten kinetics. In order to differentiate between these various mechanistic steps, an experimental apparatus known as the Porous Ceramic Tube - Nutrient Delivery System (PCT-NDS) was used. Manipulation of the applied suction pressure circulating a nutrient solution through this system imposes a change in the matric component of the water potential. This compensates for the different osmotic components of water potential dictated by nutrient concentration. By maintaining this control over the root-zone conditions, the rate limiting steps in the uptake of the essential nutrients into tomato plants (Lycopersicon esculentum cv. Cherry Elite) were differentiated. Results showed that the uptake of some nutrients were mass transfer limited while others were limited by the enzyme kinetics. Each of these were adequately modeled with calculations and discussions of the parameter estimations provided.

Tsao, David↗