Search NASA⌕ Search

SEARCH · Search NASA

Results for “solute transport modeling”

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 145 records · Page 8

Continuum Description of the Role of Negative Transference Numbers on Ion Motion in Polymer Electrolytes

New experimental techniques such as electrophoretic NMR (eNMR) are emerging as powerful methods for directly measuring ion velocities in electrolytes under applied electric fields. The aim of this theoretical study is to predict the spatial- and temporal-dependence of these velocities of ions as a function of the magnitude of the applied field and salt concentration. It has recently been shown that mixtures of poly(ethylene oxide)-based (PEO) and lithium bis(trifluoromethanesulfonyl) imide (LiTFSI) electrolytes exhibit negative cation transference numbers in a certain salt concentration range. In this range, the cation motion at early times is directed to the positive electrode at all locations in the cell; ion migration dominates in this regime. As time progresses, the cation velocity in finite zones near both the electrodes changes sign. These zones grow rapidly with time, reflecting the increasing importance of diffusion, and a point in time is reached beyond which the cation velocity in the entire cell is directed toward the negative electrode. Overall, our work reveals the limited time window over which the results of eNMR can be used to determine the transference number. More importantly, it shows how to account for the effect of diffusional flux in such experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗

Tokamak Disruption Simulation (Final Technical Report)

This was a multi-institutional project concerned with simulation of tokamak disruption – rapid termination of the plasma and confining magnetic field in a tokamak fusion reactor – with the aim of quantifying the physics associated with disruption mitigation. The PI’s role is as one of a group of computational scientists concerned with development of effective solution algorithms for the underlying models of magnetohydrodynamics and transport that represent the physical processes, and to develop techniques of uncertainty quantification for treating unknown aspects of the model such as boundary conditions and effects of impurities in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗

A theoretical study of mixing downstream of transverse injection into a supersonic boundary layer

A theoretical and analytical study was made of mixing downstream of transverse hydrogen injection, from single and multiple orifices, into a Mach 4 air boundary layer over a flat plate. Numerical solutions to the governing three-dimensional, elliptic boundary layer equations were obtained using a general purpose computer program. Founded upon a finite element solution algorithm. A prototype three-dimensional turbulent transport model was developed using mixing length theory in the wall region and the mass defect concept in the outer region. Excellent agreement between the computed flow field and experimental data for a jet/freestream dynamic pressure ratio of unity was obtained in the centerplane region of the single-jet configuration. Poorer agreement off centerplane suggests an inadequacy of the extrapolated two-dimensional turbulence model. Considerable improvement in off-centerplane computational agreement occured for a multi-jet configuration, using the same turbulent transport model.

Baker, A. J.↗

Magnetic monopoles and cosmology

The number density of magnetic monopoles is calculated in the framework of big-bang cosmology on the basis of a solution to the Boltzmann transport equation in evolutionary cosmological models. The product of the mass and the transition probability of the annihilation process for magnetic monopoles is obtained relative to those for protons and is shown to be greater than about 10 to the 20th power. Based on this result, two likely conclusions are made: (1) monopoles exist, but possess a pole-antipole bound state with binding energy of the order of the monopole rest mass; or (2) magnetic monopoles do not exist.

Adams, P. J.↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

TranSyT , an innovative framework for identifying transport systems

The importance and rate of development of genome-scale metabolic models have been growing for the last few years, increasing the demand for software solutions that automate several steps of this process. However, since TRIAGE’s release, software development for the automatic integration of transport reactions into models has stalled. Here, in this paper, we present the Transport Systems Tracker (TranSyT). Unlike other transport systems annotation software, TranSyT does not rely on manual curation to expand its internal database, which is derived from highly curated records retrieved from the Transporters Classification Database and complemented with information from other data sources. TranSyT compiles information regarding transporter families and proteins, and derives reactions into its internal database, making it available for rapid annotation of complete genomes. All transport reactions have GPR associations and can be exported with identifiers from four different metabolite databases. TranSyT is currently available as a plugin for merlin v4.0 and an app for KBase.

59 BASIC BIOLOGICAL SCIENCES↗

A radioisotope - enabled reactive transport model for deep vadose zone carbon

In mountainous regions, which constitute the principle source of recharge to major rivers and regional aquifers, infiltration occurs through fractured, partially saturated, weathered bedrock that acts as a boundary layer between saturated aquifers and surface soil. Commonly this deep vadose zone (DVZ) is many meters thick, and yet its role in regulating the generation, retention and mobility of reactive solutes, including nutrients, contaminants and weathering products, is largely unknown. In particular, many of the key reactions that drive the formation of the weathered DVZ and the quality of water moving through it are redox processes, regulated by the availability of organic carbon and oxygen below the soil layer. The role of the DVZ is thus also poorly constrained in the context of carbon stocks and mobility, particularly in lithologies that are naturally high in organic carbon, such as shales. The overarching hypothesis of this study is that upland regions developed in geologic settings with abundant petrogenic carbon store and actively cycle carbon in the weathered DVZ below the soil and above the water table at rates that are significant and currently unconstrained. In order to quantify this cycling, the current study combines novel instrumentation techniques allowing new direct sampling of DVZ systems with advanced numerical reactive transport simulations of carbon transport and transformation. Critically, these simulations will explicitly treat the three isotopes of carbon (the abundant 12C, the stable rare 13C and the radioactive 14C) in a unified framework, thus clearly parsing between the contributions of modern surface derived carbon and lithologic carbon sources in integrated measurements of fluid and gas phase fluxes. This novel model capability will be applied to test the role of DVZ carbon cycling as a regulator of water quality and geological weathering in two complementary field sites both located in organic carbon rich shale lithologies. The first is the Eel River Critical Zone Observatory (ERCZO) in Mendocino County, California, and the second is the Lawrence Berkeley National Laboratory Watershed Function Scientific Focus Area (SFA) in the East River watershed, near Crested Butte, Colorado. At the ERCZO, a novel vadose zone monitoring system has been installed in a 20 m thick, partially saturated, weathered shale hillslope, and preliminary data already indicate substantial CO2 flux generated many meters below the soil surface. At the SFA field site, an instrumented hillslope transect indicates a more complex multi-dimensional fluid and solute transport regime, which will serve as a key test of the calibrated models. Collectively, this project will advance understanding of the cycling of carbon belowground and in relation to transport pathways across the poorly constrained DVZ characteristic of primary water recharge areas. The key product of this work will be enhanced isotope simulation capabilities that are robust and publicly available for application across a broad diversity of systems.

58 GEOSCIENCES↗

Analytical solutions to the problem of transient heat transfer in living tissue.

An analytical model of transient heat transfer in living biological tissue is considered. The model includes storage, generation, conduction, and convective transport of heat in the tissue. Solutions for rectangular and cylindrical coordinates are presented and discussed. Transient times for reaching the ?locally fully developed' temperature profile were found to be of the order of 5 to 25 min. These transients are dominated by a geometrical parameters and, to a lesser extent, by a parameter representing the ratio of heat supplied by blood flow to heat conducted in the tissue.

Shitzer, A.↗

Oceanic Fluxes of Mass, Heat and Freshwater: A Global Estimate and Perspective

Data from fifteen globally distributed, modern, high resolution, hydrographic oceanic transects are combined in an inverse calculation using large scale box models. The models provide estimates of the global meridional heat and freshwater budgets and are used to examine the sensitivity of the global circulation, both inter and intra-basin exchange rates, to a variety of external constraints provided by estimates of Ekman, boundary current and throughflow transports. A solution is found which is consistent with both the model physics and the global data set, despite a twenty five year time span and a lack of seasonal consistency among the data. The overall pattern of the global circulation suggested by the models is similar to that proposed in previously published local studies and regional reviews. However, significant qualitative and quantitative differences exist. These differences are due both to the model definition and to the global nature of the data set.

MacDonald, Alison Marguerite↗

Static Aeroelastic and Longitudinal Trim Model of Flexible Wing Aircraft Using Finite-Element Vortex-Lattice Coupled Solution

This paper presents a static aeroelastic model and longitudinal trim model for the analysis of a flexible wing transport aircraft. The static aeroelastic model is built using a structural model based on finite-element modeling and coupled to an aerodynamic model that uses vortex-lattice solution. An automatic geometry generation tool is used to close the loop between the structural and aerodynamic models. The aeroelastic model is extended for the development of a three degree-of-freedom longitudinal trim model for an aircraft with flexible wings. The resulting flexible aircraft longitudinal trim model is used to simultaneously compute the static aeroelastic shape for the aircraft model and the longitudinal state inputs to maintain an aircraft trim state. The framework is applied to an aircraft model based on the NASA Generic Transport Model (GTM) with wing structures allowed to flexibly deformed referred to as the Elastically Shaped Aircraft Concept (ESAC). The ESAC wing mass and stiffness properties are based on a baseline "stiff" values representative of current generation transport aircraft.

Aeroelasticity↗

Fully consistent Eulerian Monte Carlo fields method for solving probability density function transport equations in turbulence modeling

The probability density function (PDF) transport equation method is a sophisticated model for the closure of turbulent mixing and turbulent reactive flows. An efficient solution approach for solving the PDF transport equation has been vital for the method to be widely used in applications. The Eulerian Monte Carlo fields (EMCF) method has been developed to solve the PDF transport equation efficiently for decades. A recent work by Wang et al. revealed a serious issue of the EMCF method for not being fully consistent with the PDF transport equation for which the method is designed to solve. This work advances the state of the art by introducing fully consistent EMCF methods for solving the PDF transport equation. The fully consistent EMCF formulations are derived for two different PDF equation forms. The consistency of the EMCF formulations is mathematically confirmed by examining the derived moment transport equations from the EMCF formulations and from the PDF transport equation. The method of manufactured solutions is employed to further verify the consistency and convergence of the different EMCF formulations numerically. The newly introduced EMCF formulations bring the EMCF method to full consistency with the PDF transport equations for the first time.

42 ENGINEERING↗

Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution

Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.

36 MATERIALS SCIENCE↗

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of large-scale effects on mass transfer and concentration polarization in Reverse Osmosis membrane systems

Here, we present well-resolved computational fluid dynamics simulations of a large-scale reverse osmosis membrane-spacer configuration (1 m). Our computational model solves the flow and transport equations with variable solute-dependent properties. We utilize a high resolution computational mesh to resolve all relevant length scales associated with spacer-induced mixing and thin concentration boundary layers. An important contribution of this work is the development of a modified mass-transfer correlation that accounts for the development of the concentration boundary layer along the channel. A set of 2D axisymmetric simulations were performed for a spiral wound module layer with varying cross-flow conditions and spacer diameters which indicate a significant entrance length effect for concentration profile development at lower flow rates while mixing effects dominate at higher flow rates. The mass-transfer correlations at higher flow rates compare well with published correlations while a surrogate model for Sherwood number was obtained that depends on an additional similarity variable that accounted for entrance length effects at lower flow rates. Finally, a large-scale membrane-spacer design relevant to high-pressure reverse osmosis is studied with a non-uniform arrangement of spacers, which indicate a substantial saving in pressure drop (~40%) compared to traditional uniformly spaced pattern with minor variations (~2%) in concentration polarization, product water quality (~1%) and water recovery (~7%) compared to a uniform spacer pattern.

97 MATHEMATICS AND COMPUTING↗

Photoelectrons in the Quiet Polar Wind

This study presents a newly coupled model capable of treating the superthermal electron population in the global polar wind solution. The model combines the hydrodynamic Polar Wind Outflow Model (PWOM) with the kinetic SuperThermal Electron Transport (STET) code. The resulting PWOM-STET coupled model is described and then used to investigate the role of photoelectrons in the polar wind. We present polar wind results along single stationary field lines under dayside and nightside conditions, as well as the global solution reconstructed from nearly 1000 moving field lines. The model results show significant day-night asymmetries in the polar wind solution owing to the higher ionization and photoelectron fluxes on the dayside compared to the nightside. Field line motion is found to modify this dependence and create global structure by transporting field lines through different conditions of illumination and through the localized effects of Joule heating.

wind↗

Dendritic Alloy Solidification Experiment (DASE)

A space experiment, and supporting ground-based research, is proposed to study the microstructural evolution in free dendritic growth from a supercooled melt of the transparent model alloy succinonitrile-acetone (SCN-ACE). The research is relevant to equiaxed solidification of metal alloy castings. The microgravity experiment will establish a benchmark for testing of equiaxed dendritic growth theories, scaling laws, and models in the presence of purely diffusive, coupled heat and solute transport, without the complicating influences of melt convection. The specific objectives are to: determine the selection of the dendrite tip operating state, i.e. the growth velocity and tip radius, for free dendritic growth of succinonitrile-acetone alloys; determine the growth morphology and sidebranching behavior for freely grown alloy dendrites; determine the effects of the thermal/solutal interactions in the growth of an assemblage of equiaxed alloy crystals; determine the effects of melt convection on the free growth of alloy dendrites; measure the surface tension anisotropy strength of succinon itrile -acetone alloys establish a theoretical and modeling framework for the experiments. Microgravity experiments on equiaxed dendritic growth of alloy dendrites have not been performed in the past. The proposed experiment builds on the Isothermal Dendritic Growth Experiment (IDGE) of Glicksman and coworkers, which focused on the steady growth of a single crystal from pure supercooled melts (succinonitrile and pivalic acid). It also extends the Equiaxed Dendritic Solidification Experiment (EDSE) of the present investigators, which is concerned with the interactions and transients arising in the growth of an assemblage of equiaxed crystals (succinonitrile). However, these experiments with pure substances are not able to address the issues related to coupled heat and solute transport in growth of alloy dendrites.

Beckermann, C.↗