Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion model”

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 415 records · Page 23

On the nature of the cosmic ray positron spectrum

A new calculation has been made of the flux of secondary positrons above 100 MeV expected for various propagation models. The models investigated are the leaky box or homogeneous model, a disk-halo diffusion model, a dynamical halo model, and the closed galaxy model. The parameters of these models have, in each case, been adjusted for agreement with the observed secondary/primary ratios and Be-10 abundance. The positron flux predicted for these models is compared with the available data. The possibility of a primary positron component is considered.

Protheroe, R. J.↗

The upper ionospheres of Jupiter and Saturn

The electron-density profiles of Jupiter's and Saturn's ionospheres are modeled with a 1D chemical diffusive model incorporating measured parameters of the neutral atmospheric structure. The Voyager RSS ion-density data are modeled by accounting for ion chemistry, the H3(+) recombination-rate coefficient, and H2O chemistry. Electron-density peaks for both ionospheres are 900-1000 km lower in the model than the measured values, and the role of vibrational excitation of H2 is discussed in converting H(+) to H2(+) and H3(+). Vertical ion flow is also considered which can maintain the plasma peaks under the conditions of electrical fields or horizontal neutral winds. An assumed influx of H2O molecules of a specific quantity is theorized to reproduce the measured values on Saturn when combined with a vertical plasma drift.

Majeed, Tariq↗

Spiral waves are stable in discrete element models of two-dimensional homogeneous excitable media

The spontaneous breakup of a single spiral wave of excitation into a turbulent wave pattern has been observed in both discrete element models and continuous reaction-diffusion models of spatially homogeneous 2D excitable media. These results have attracted considerable interest, since spiral breakup is thought to be an important mechanism of transition from the heart rhythm disturbance ventricular tachycardia to the fatal arrhythmia ventricular fibrillation. It is not known whether this process can occur in the absence of disease-induced spatial heterogeneity of the electrical properties of the ventricular tissue. Candidate mechanisms for spiral breakup in uniform 2D media have emerged, but the physical validity of the mechanisms and their applicability to myocardium require further scrutiny. In this letter, we examine the computer simulation results obtained in two discrete element models and show that the instability of each spiral is an artifact resulting from an unphysical dependence of wave speed on wave front curvature in the medium. We conclude that spiral breakup does not occur in these two models at the specified parameter values and that great care must be exercised in the representation of a continuous excitable medium via discrete elements.

NASA Discipline Cardiopulmonary↗

Effect of Rain-Adjusted Satellite Sea Surface Salinity on ENSO Predictions from the GMAO S2S Forecast System

ENSO has a significant impact on climate variability throughout the world and so has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of altimetry and in situ data lead to improved short-term forecasts of the coupled system. However, few studies have focused on improving the near-surface density/mixing through assimilation of satellite sea surface salinity (SSS). For expediency, most projects that do assimilate SSS do so as if these data were observed at the top model layer instead of at the surface. In rainy regions where buoyant water sits as a fresh lens at the surface, this assumption is likely invalid. Therefore, we adjust SSS so that it more accurately represents the salinity at the first model layer. The Rain Impact Model (RIM) uses a simple diffusion model to determine the near surface salinity gradient (i.e., 1 cm to 5 m). Satellite SSS data are modified using this near-surface salinity gradient, so the salinity values are now valid at the first model layer (we call this SSS@5m). We assess the impact of satellite SSS observations for near-surface dynamics within ocean reanalyses and how these impact dynamical ENSO forecasts using the NASA GMAO S2S coupled system. For all reanalysis experiments, all available along-track sea level and in situ observations are assimilated for 2011 to 2020 using the LETKF scheme (Penny et al., 2013). One reanalysis assimilates Aquarius /SMAP SSS as before. An additional reanalysis is performed assimilating the SSS@5m data. Validation statistics are compared for experiments that assimilate SSS (sub-optimally as before) versus the SSS@5m. We also compare results of coupled forecasts that are initialized from these reanalyses in spring. For all but the big La Niña in 2017, all NINO3.4 forecasts were improved by using the RIM. We will show that improved SSS estimates upgrades density and near-surface mixing leading to more accurate coupled air/sea interaction and better forecasts.

ENSO↗

Taking the CCDs to the ultimate performance for low threshold experiments

Scientific grade CCDs show atractive capabilities for the detection of particles with small energy deposition in matter. Their very low threshold of approximately 40 eV and their good spatial reconstruction of the event are key properties for currently running experiments: CONNIE and DAMIC. Both experiments can benefit from any increase of the detection efficiency of nuclear recoils at low energy. In this work we present two different approaches to increase this efficiency by increasing the SNR of events. The first one is based on the reduction of the readout noise of the device, which is the main contribution of uncertainty to the signal measurement. New studies on the electronic noise from the integrated output amplifier and the readout electronics will be presented together with result of a new configuration showing a lower limit on the readout noise which can be implemented on the current setup of the CCD based experiments. A second approach to increase the SNR of events at low energy that will be presented is the studies of the spatial conformation of nuclear recoil events at different depth in the active volume by studies of new effects that differ from expected models based on not interacting diffusion model of electrons in the semiconductor.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A New Photocatalyst Bismuth Oxo Citrate: Synthesis, Characterization, and Photocatalytic Performance

A new photocatalyst bismuth oxo citrate was synthesized by facile precipitation process with calcination at 200 8C. The photocatalyst was characterized by scanning electron microscopy (SEM), X-ray diffraction (XRD), Fouriertransform infrared (FTIR) spectroscopy, N2 sorptometry, and elemental analysis. Morphologically, it is composed of polyhedral particles with different, irregular shapes and sizes. The specific surface area (SSA) of the photocatalyst was 8.92 m2 g-1. It showed very good photocatalytic performance and reusability. Total decolorization of Reactive Blue 19 (RB19) was achieved in less than 10 minutes, which is much faster in comparison with TiO{sub 2} P25. Also, bismuth oxo citrate showed higher photocatalytic activity than other photocatalysts based on bismuth compounds reported by other authors. Optimal photocatalysis parameters were pH 2 and photocatalyst dose of 250 mg dm-3. The decolorization rate was found to decrease as initial dye concentration increased. The photocatalytic data best fitted to L-H kinetic model with pseudo-first order reaction rate. Chrastil diffusion model showed that diffusion has not influence on the process. Water Environ. Res., 90, 719 (2018). © 2018 Water Environment Federation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling O and N Alloying in Nb for SRF Applications

N and O-alloyed superconducting radio frequency cavities exhibit extraordinary quality factors. Developing diffusion models that describe interstitial N and O in Nb is important for optimizing alloyed cavity quality factors and accelerating gradients. N and O-alloyed Nb samples are examined with SEM AND SIMS and their diffusion profiles modeled.

Lechner, Eric↗

Systematic study of Reynolds stress closure models in the computations of plane channel flows

The roles of pressure-strain and turbulent diffusion models in the numerical calculation of turbulent plane channel flows with second-moment closure models are investigated. Three turbulent diffusion and five pressure-strain models are utilized in the computations. The main characteristics of the mean flow and the turbulent fields are compared against experimental data. All the features of the mean flow are correctly predicted by all but one of the Reynolds stress closure models. The Reynolds stress anisotropies in the log layer are predicted to varying degrees of accuracy (good to fair) by the models. None of the models could predict correctly the extent of relaxation towards isotropy in the wake region near the center of the channel. Results from the directional numerical simulation are used to further clarify this behavior of the models.

Demuren, A. O.↗

Perspective - Systematic study of Reynolds stress closure models in the computations of plane channel flows

The roles of pressure-strain and turbulent diffusion models in the numerical calculation of turbulent plane channel flows with second-moment closure models are investigated. Three turbulent diffusion and five pressure-strain models are utilized in the computations. The main characteristics of the mean flow and the turbulent fields are compared against experimental data. All the features of the mean flow are correctly predicted by all but one of the Reynolds stress closure models. The Reynolds stress anisotropies in the log layer are predicted to varying degrees of accuracy (good to fair) by the models. None of the models could predict correctly the extent of relaxation towards isotropy in the wake region near the center of the channel. Results from the directional numerical simulation are used to further clarify this behavior of the models.

Demuren, A. O.↗

Modeling Ni redistribution in the hydrogen electrode of solid oxide cells through Ni(OH) 2 diffusion and Ni-YSZ wettability change

Ni redistribution in the hydrogen electrodes of solid oxide cells is an important degradation mechanism. Its driving force is still under debate. This work focuses on the interplay between Ni(OH) 2 diffusion and Ni-YSZ wettability change. With a Ni(OH) 2 diffusion model developed in our previous work, we further employ three models for the Ni-YSZ contact angle to quantify the Ni-YSZ wettability change. The microstructure evolutions in a reconstructed Ni-YSZ electrode are simulated under selected experimental operating conditions. It is shown that the phenomenological model can capture the driving force of Ni spreading/detachment on YSZ surface and the Ni migration reported in experiments can be qualitatively reproduced through the competition between Ni(OH) 2 diffusion and Ni-YSZ wettability change. It is also found that the initial microstructure of the active layer and the microstructure characteristics in the support layer can strongly affect the distribution of steam partial pressure and overpotential, hence the Ni redistribution, in the active layer, which may explain the inconsistencies in experiments. Further, our results show that there are two missing pieces in current theory of Ni redistribution: a mechanism of fast Ni diffusion under humid atmosphere and the physical reason behind the Ni spreading/detachment.

25 ENERGY STORAGE↗

A stochastic extension of the explicit algebraic subgrid-scale models

The explicit algebraic subgrid-scale (SGS) stress model (EASM) of Marstorp et al. [“Explicit algebraic subgrid stress models with application to rotating channel flow,” J. Fluid Mech. 639, 403–432 (2009)] and explicit algebraic SGS scalar flux model (EASFM) of Rasam et al. [“An explicit algebraic model for the subgrid-scale passive scalar flux,” J. Fluid Mech. 721, 541–577 (2013)] are extended with stochastic terms based on the Langevin equation formalism for the subgrid-scales by Marstorp et al. [“A stochastic subgrid model with application to turbulent flow and scalar mixing,” Phys. Fluids 19, 035107 (2007)]. The EASM and EASFM are nonlinear mixed and tensor eddy-diffusivity models, which improve large eddy simulation (LES) predictions of the mean flow, Reynolds stresses, and scalar fluxes of wall-bounded flows compared to isotropic eddy-viscosity and eddy-diffusivity SGS models, especially at coarse resolutions. The purpose of the stochastic extension of the explicit algebraic SGS models is to further improve the characteristics of the kinetic energy and scalar variance SGS dissipation, which are key quantities that govern the small-scale mixing and dispersion dynamics. LES of turbulent channel flow with passive scalar transport shows that the stochastic terms enhance SGS dissipation statistics such as length scale, variance, and probability density functions and introduce a significant amount of backscatter of energy from the subgrid to the resolved scales without causing numerical stability problems. The improvements in the SGS dissipation predictions in turn enhances the predicted resolved statistics such as the mean scalar, scalar fluxes, Reynolds stresses, and correlation lengths. Moreover, the nonalignment between the SGS stress and resolved strain-rate tensors predicted by the EASM with stochastic extension is in much closer agreement with direct numerical simulation data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling the Contribution of Precipitation Loss to a Radiation Belt Electron Dropout Observed by Van Allen Probes

A drift-diffusion model is used to simulate the low-altitude electron distribution, accounting for azimuthal drift, pitch angle diffusion, and atmospheric backscattering effects during a rapid electron dropout event on 21st August 2013, at L = 4.5. Additional external loss effects are introduced during times when the low-altitude electron distribution cannot be reproduced by diffusion alone. Here, the model utilizes low-altitude electron count rate data from five POES/MetOp satellites to quantify pitch angle diffusion rates. Low-altitude data provides critical constraint on the model because it includes the drift loss cone region where the electron distribution in longitude is highly dependent on the balance between azimuthal drift and pitch angle diffusion. Furthermore, a newly derived angular response function for the detectors onboard POES/MetOp is employed to accurately incorporate the bounce loss cone measurements, which have been previously contaminated by electrons from outside the nominal field-of-view. While constrained by low-altitude data, the model also shows reasonable agreement with high-altitude data. Pitch angle diffusion rates during the event are quantified and are faster at lower energies. Precipitation is determined to account for all of the total loss observed for 450 keV electrons, 88% for 600 keV and 38% for 900 keV. Predictions made in the MeV range are deemed unreliable as the integral energy channels E3 and P6 fail to provide the necessary constraint at relativistic energies.

79 ASTRONOMY AND ASTROPHYSICS↗

Diffusion-controlled reactions modeling in Geant4-DNA

Context Under irradiation, a biological system undergoes a cascade of chemical reactions that can lead to an alteration of its normal operation. There are different types of radiation and many competing reactions. As a result the kinetics of chemical species is extremely complex. The simulation becomes then a powerful tool which, by describing the basic principles of chemical reactions, can reveal the dynamics of the macroscopic system. To understand the dynamics of biological systems under radiation, since the 80s there have been on-going efforts carried out by several research groups to establish a mechanistic model that consists in describing all the physical, chemical and biological phenomena following the irradiation of single cells. This approach is generally divided into a succession of stages that follow each other in time: (1) the physical stage, where the ionizing particles interact directly with the biological material; (2) the physico-chemical stage, where the targeted molecules release their energy by dissociating, creating new chemical species; (3) the chemical stage, where the new chemical species interact with each other or with the biomolecules; (4) the biological stage, where the repairing mechanisms of the cell come into play. This article focuses on the modeling of the chemical stage. Method This article presents a general method of speeding-up chemical reaction simulations in fluids based on the Smoluchowski equation and Monte-Carlo methods, where all molecules are explicitly simulated and the solvent is treated as a continuum. The model describes diffusion-controlled reactions. This method has been implemented in Geant4-DNA. The keys to the new algorithm include: (1) the combination of a method to compute time steps dynamically with a Brownian bridge process to account for chemical reactions, which avoids costly fixed time step simulations; (2) a k–d tree data structure for quickly locating, for a given molecule, its closest reactants. The performance advantage is presented in terms of complexity, and the accuracy of the new algorithm is demonstrated by simulating radiation chemistry in the context of the Geant4-DNA project. Application The time-dependent radiolytic yields of the main chemical species formed after irradiation are computed for incident protons at different energies (from 50 MeV to 500 keV). Both the time-evolution and energy dependency of the yields are discussed. The evolution, at one microsecond, of the yields of hydroxyls and solvated electrons with respect to the linear energy transfer is compared to theoretical and experimental data. According to our results, at high linear energy transfer, modeling radiation chemistry in the trading compartment representation might be adopted.

60 APPLIED LIFE SCIENCES↗

Fundamental Study on Sorption Behavior of Europium to Mica Flakes - 19116

It is well known that biotite exists with flake form in actual underground environment, dominating sorption of cationic radionuclides in a geological disposal system consisting of crystalline rock. The behavior of sorption of radionuclides will differ depending on the difference in shape, while most of the previous studies have used powder biotite. This study focused on the influence of flake structure of phyllosilicate minerals such as biotite on the diffusion and sorption processes of radionuclides, examining the sorption behaviors of Eu on biotite, phlogopite and muscovite in flake form. These flakes (5 mm x 6 mm) were prepared for the sorption experiments and the initial concentrations of the Eu ions were set to 1 mM. The sorption experiments were conducted by mixing 0.4 g of biotite flakes and 4 mL at pH 5 for 7 days. As a result, the sorption amount of Eu increased little by little with internal diffusion in all minerals, respectively. The diffusion process of Eu into the layer structure could be also observed by secondary ion-microprobe mass spectrometry (SIMS) analysis. The apparent diffusion coefficient of Eu into each flake was estimated by fitting a 2-D diffusion model to the sorption data. The estimated diffusion values were clearly smaller than that of matrix diffusion in granite. Hence it was suggested that the so-called Kd-model, which is a local equilibrium model, is not always applicable to the estimation of the retardation effect (i.e., migration rate) of radionuclides in rock matrix. In other words, the migration rate might depend also on the diffusion process of radionuclide in the inside of phyllosilicate minerals in addition of the macro diffusion process in rock matrix. Besides, sorption amount and diffusion coefficient of biotite, phlogopite and muscovite were largely different. Not only the tetrahedral structure adjacent to the potassium but also the octahedral structure might affect the restraint of the interlayer ion (K ions) and the sorption behavior. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A two-phase model that unifies and extends the classical models of membrane transport

Two models describe solvent transport through swollen, nonporous membranes. The pore-flow model, based on fluid mechanics, works for porous membranes, whereas the solution-diffusion model invokes molecular diffusion to treat nonporous membranes. Both approaches make valid arguments for swollen polymer membranes, but they disagree in their predictions of intramembrane pressure and concentration profiles. Using a fluid-solid model that treats the solvent and membrane matrix as separate phases, we show both classical models to be valid, to represent complementary approaches to the same phenomenon, and to make identical predictions. The fluid-solid model clarifies recent reverse osmosis measurements; provides a predictive and mechanistic basis for empirical high-pressure limiting flux phenomena, in quantitative agreement with classic measurements; and gives a framework to treat nonporous but mechanically heterogeneous membrane materials.

Science & Technology - Other Topics↗

Pickup protons and water ions at Comet Halley - Comparisons with Giotto observations

The cometary ion pickup process along the sun-comet line at Comet Halley is investigated using a quasi-linear diffusion model including both pitch angle and energy diffusion, adiabatic compression, and convective motion with the solar wind flow. The model results are compared with energetic ion distributions observed by instruments on board the Giotto spacecraft. The observed power spectrum index of magnetic turbulence (gamma) is 2-2.5. The present simulation shows that when gamma was 2, the calculated proton distributions were much more isotropic than the observed ones. The numerical solutions of the quasi-linear diffusion equations show that the isotropization of the pickup ion distribution, particularly at the pickup velocity, is not complete even close to the bow shock. Given the observed turbulence level, quasi-linear theory yields pickup ion energy distributions that agree with the observed ones quite well and easily produces energetic ions with energies up to hundreds of keV.

Ye, G.↗

Nonlinear optimization for a tumor invasion PDE model

In this work, we introduce a methodology to approximate unknown parameters that appear on a non-linear reaction–diffusion model of tumor invasion. These equations consider that tumor-induced alteration of micro-environmental pH furnishes a mechanism for cancer invasion. A coupled system reaction–diffusion explaining this model is given by three partial differential equations for the non-dimensional spatial distribution and temporal evolution of the density of normal tissue, the neoplastic tissue growth and the excess concentration of H{sup +} ions. The tumor model parameters have a corresponding biological meaning: the reabsorption rate, the destructive influence of H{sup +} ions in the healthy tissue, the growth rate of tumor tissue and the diffusion coefficient. We propose to solve the direct problem using the Finite Element Method (FEM) and minimize an appropriate functional including both the real data (obtained via in-vitro experiments and fluorescence ratio imaging microscopy) and the numerical solution. The gradient of the functional is computed by the adjoint method.

60 APPLIED LIFE SCIENCES↗