Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

3-D Modeling of Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe Growth Under Microgravity Conditions

A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.

Bune, Andris V.↗
DOE OSTI2019

Rate constants in spatially inhomogeneous systems

Here, we introduce a theory and accompanying importance sampling method for computing rate constants in spatially inhomogeneous systems. Using the relationship between rate constants and path space partition functions, we illustrate that the relative change in the rate of a rare event through space is isomorphic to the calculation of a free energy difference, albeit in a trajectory ensemble. Like equilibrium free energies, relative rate constants can be estimated by importance sampling. An extension to transition path sampling is proposed that combines biased path ensembles and weighted histogram analysis to accomplish this estimate. We demonstrate that rate constants can also be decomposed into different contributions, including relative changes in stability, barrier height, and flux. This decomposition provides a means of interpretation and insight into rare processes in complex environments. We confirm these ideas with a simple model of diffusion with spatially varying diffusivity and illustrate their utility in a model of ion pair dissociation near an electrochemical interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Explaining persistent incomplete mixing in multicomponent reactive transport with Eulerian stochastic model

We present an Eulerian stochastic advection–diffusion–reaction (SADR) model and use it to explain incomplete mixing typically observed in transport experiments with bimolecular reactions. Unlike traditional advection– dispersion–reaction (ADR) models, the SADR model describes mechanical and diffusive mixing as two separate processes. In the SADR model, mechanical mixing is driven by random advective velocity whose variance is given by the coefficient of mechanical dispersion. The diffusive mixing is modeled as a Fickian diffusion process with the effective diffusion coefficient. We demonstrate that the sum of the two coefficients is equal to the dispersion coefficient, but only the effective diffusion coefficient contributes to the mixing-controlled reactions. We use experimental results of Gramling et al. * * to show that for transport and bimolecular reactions in porous media, the SADR model is significantly more accurate than the ADR model, which overestimates the concentration of the reaction product by as much as 60%. We also show that the SADR model predicts an experimentally observed bell- shaped spatial distribution of the reactive product concentration, while the ADR model results in a concentration distribution with an unphysical kink.

Bimolecular reactions, incomplete mixing, Stochast↗
DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗
NASA NTRS2018

Sedimentation Efficiency of Condensation Clouds in Substellar Atmospheres

Condensation clouds in substellar atmospheres have been widely inferred from spectra and photometric variability. Up until now, their horizontally averaged vertical distribution and mean particle size have been largely characterized using models, one of which is the eddy diffusion-sedimentation model from Ackerman and Marley that relies on a sedimentation efficiency parameter, f(sub sed), to determine the vertical extent of clouds in the atmosphere. However, the physical processes controlling the vertical structure of clouds in substellar atmospheres are not well understood. In this work, we derive trends in f(sub sed) across a large range of eddy diffusivities (K (sub zz) ), gravities, material properties, and cloud formation pathways by fitting cloud distributions calculated by a more detailed cloud microphysics model. We find that f(sub sed) is dependent on K (sub zz), but not gravity, when K (sub zz) is held constant. f(sub sed) is most sensitive to the nucleation rate of cloud particles, as determined by material properties like surface energy and molecular weight. High surface energy materials form fewer, larger cloud particles, leading to large f(sub sed) (>1), and vice versa for materials with low surface energy. For cloud formation via heterogeneous nucleation, f(sub sed) is sensitive to the condensation nuclei flux and radius, connecting cloud formation in substellar atmospheres to the objects' formation environments and other atmospheric aerosols. These insights could lead to improved cloud models that help us better understand substellar atmospheres. For example, we demonstrate that f(sub sed) could increase with increasing cloud base depth in an atmosphere, shedding light on the nature of the brown dwarf L/T transition.

brown dwarfs; planets and satellites: atmospheres↗
NASA NTRS2000

Postural control model interpretation of stabilogram diffusion analysis

Collins and De Luca [Collins JJ. De Luca CJ (1993) Exp Brain Res 95: 308-318] introduced a new method known as stabilogram diffusion analysis that provides a quantitative statistical measure of the apparently random variations of center-of-pressure (COP) trajectories recorded during quiet upright stance in humans. This analysis generates a stabilogram diffusion function (SDF) that summarizes the mean square COP displacement as a function of the time interval between COP comparisons. SDFs have a characteristic two-part form that suggests the presence of two different control regimes: a short-term open-loop control behavior and a longer-term closed-loop behavior. This paper demonstrates that a very simple closed-loop control model of upright stance can generate realistic SDFs. The model consists of an inverted pendulum body with torque applied at the ankle joint. This torque includes a random disturbance torque and a control torque. The control torque is a function of the deviation (error signal) between the desired upright body position and the actual body position, and is generated in proportion to the error signal, the derivative of the error signal, and the integral of the error signal [i.e. a proportional, integral and derivative (PID) neural controller]. The control torque is applied with a time delay representing conduction, processing, and muscle activation delays. Variations in the PID parameters and the time delay generate variations in SDFs that mimic real experimental SDFs. This model analysis allows one to interpret experimentally observed changes in SDFs in terms of variations in neural controller and time delay parameters rather than in terms of open-loop versus closed-loop behavior.

NASA Program Biomedical Research and Countermeasur↗
NASA NTRS2019

Modeling transmittance of MODIS solar diffuser stability monitor sun view screen

Earlier work by Sun et al. [2005] modeled the Sun view screen transmittance of the solar diffuser stability monitor (SDSM) on Moderate Resolution Imaging Spectroradiometer (MODIS) which explains the oscillatory pattern in the Terra/Aqua MODIS SDSM Sun view response collected during the post-launch yaw maneuvers. The SDSM screen consists of small size pinholes which limit the screen transmittance. The SDSM Sun view transmittance pattern is affected by the size and distribution of pinholes, screen thickness, and orientation between the aperture stop and the Sun view screen. Planning of post-launch yaw maneuver test often requires comprehensive angular coverage to resolve the angular variation of the screen transmittance. This paper developed a Gauss-Kronrod quadrature method-based model to simulate the SDSM Sun view transmittance. This method performs iterative integration of the interception area between pinhole and aperture stop which is critical in the screen transmittance calculation. Parametric dependences of the MODIS SDSM Sun view transmittance on various geometric factors such as the pinhole size or spacing and misalignment between the SDSM Sun view screen and aperture stop are studied. The implications for the post-launch yaw maneuver planning to resolve the angular variation of the SDSM Sun view transmittance due to the SDSM Sun view screen misalignment are discussed.

MODIS↗
NASA NTRS1996

Diffusion, Viscosity and Crystal Growth in Microgravity

The diffusivity of TriGlycine Sulfate (TGS), Potassium Dihydrogen Phosphate (KDP), Ammonium Dihydrogen Phosphate (ADF) and other compounds of interest to microgravity crystal growth, in supersaturated solutions as a function of solution concentration, 'age' and 'history was studied experimentally. The factors that affect the growth of crystals from water solutions in microgravity have been examined. Three non-linear optical materials have been studied, potassium dihydrogen phosphate (KDP), ammonium dihydrogen phosphate (ADP) and triglycine sulfate (TGC). The diffusion coefficient and viscosity of supersaturated water solutions were measured. Also theoretical model of diffusivity and viscosity in a metastable state, model of crystal growth from solution including non-linear time dependent diffusivity and viscosity effect and computer simulation of the crystal growth process which allows simulation of the microgravity crystal growth were developed.

Myerson, Allan S.↗
DOE OSTI2020

Modeling gas–shell mixing in ICF with separated reactants

Mixing between fuel and shell materials in ICF implosions can affect implosion dynamics and even prevent ignition. We use data from a series of separated reactant experiments on the National Ignition Facility to calibrate and test the predictive power of gas–shell mix models. Two models are used to estimate fuel–shell mix: a Reynolds-averaged turbulence model and molecular diffusion. Minor uncertainties in capsule manufacture, experimental conditions, and values for mix model parameters produce significant variation in simulation results. Using input/output pairs from 1D simulations, we train Gaussian process surrogate models to predict experimental quantities of interest. The surrogates are used to construct posteriors for mix model parameters by marginalizing over uncertainties in capsule manufacture and experimental conditions. Mix models are calibrated with a subset of experimental data (neutron yields, ion temperature, and bang time) and tested using the remaining data. In general, both the diffusion and turbulence model correctly predict experimental DT and TT neutron yields. Despite having more free parameters, the turbulence model underpredicts ion temperature at high convergence ratio. Furthermore, the simpler diffusion model correctly predicts these temperatures, suggesting nonhydrodynamic gas–shell mix. The computational model consistently overpredicts DD neutron yield, indicating possible shortcomings outside of the mix model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Supernova PTF 12glz: A Possible Shock Breakout Driven through an Aspherical Wind

We present visible-light and ultraviolet (UV) observations of the supernova PTF 12glz. The SN was discovered and monitored in the near-UV and R bands as part of a joint GALEX and Palomar Transient Factory campaign. It is among the most energetic SNe IIn observed to date (≈1051 erg). If the radiated energy mainly came from the thermalization of the shock kinetic energy, we show that PTF 12glz was surrounded by ~1 M ⊙ of circumstellar material (CSM) prior to its explosive death. PTF 12glz shows a puzzling peculiarity: at early times, while the freely expanding ejecta are presumably masked by the optically thick CSM, the radius of the blackbody that best fits the observations grows at ≈7000 km s –1 . Such a velocity is characteristic of fast moving ejecta rather than optically thick CSM. This phase of radial expansion takes place before any spectroscopic signature of expanding ejecta appears in the spectrum and while both the spectroscopic data and the bolometric luminosity seem to indicate that the CSM is optically thick. We propose a geometrical solution to this puzzle, involving an aspherical structure of the CSM around PTF 12glz. By modeling radiative diffusion through a slab of CSM, we show that an aspherical geometry of the CSM can result in a growing effective radius. Furthermore, this simple model also allows us to recover the decreasing blackbody temperature of PTF 12glz. SLAB-Diffusion, the code we wrote to model the radiative diffusion of photons through a slab of CSM and evaluate the observed radius and temperature, is made available online.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Structure of turbulent non-premixed flames modeled with two-step chemistry

Direct numerical simulations of turbulent diffusion flames modeled with finite-rate, two-step chemistry, A + B yields I, A + I yields P, were carried out. A detailed analysis of the turbulent flame structure reveals the complex nature of the penetration of various reactive species across two reaction zones in mixture fraction space. Due to this two zone structure, these flames were found to be robust, resisting extinction over the parameter ranges investigated. As in single-step computations, mixture fraction dissipation rate and the mixture fraction were found to be statistically correlated. Simulations involving unequal molecular diffusivities suggest that the small scale mixing process and, hence, the turbulent flame structure is sensitive to the Schmidt number.

Chen, J. H.↗
DOE OSTI2023

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗
DOE OSTI2017

The Fermi Galactic Center GeV excess and implications for dark matter

Here, the region around the Galactic Center (GC) is now well established to be brighter at energies of a few GeV than what is expected from conventional models of diffuse gamma-ray emission and catalogs of known gamma-ray sources. We study the GeV excess using 6.5 yr of data from the Fermi Large Area Telescope. We characterize the uncertainty of the GC excess spectrum and morphology due to uncertainties in cosmic-ray source distributions and propagation, uncertainties in the distribution of interstellar gas in the Milky Way, and uncertainties due to a potential contribution from the Fermi bubbles. We also evaluate uncertainties in the excess properties due to resolved point sources of gamma rays. The GC is of particular interest, as it would be expected to have the brightest signal from annihilation of weakly interacting massive dark matter (DM) particles. However, control regions along the Galactic plane, where a DM signal is not expected, show excesses of similar amplitude relative to the local background. Furthermore, based on the magnitude of the systematic uncertainties, we conservatively report upper limits for the annihilation cross-section as a function of particle mass and annihilation channel.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Predictive simulation of non-steady-state transport of gases through rubbery polymer membranes

In this paper, a multiscale, physically-based, reaction-diffusion kinetics model is developed for non-steady-state transport of simple gases through a rubbery polymer. Experimental data from the literature, new measurements of non-steady-state permeation and a molecular dynamics simulation of a gas-polymer sticking probability for a typical system are used to construct and validate the model framework. Using no adjustable parameters, the model successfully reproduces time-dependent experimental data for two distinct systems: (1) O 2 quenching of a phosphorescent dye embedded in poly(n-butyl(amino) thionylphosphazene), and (2) O 2 , N 2 , CH 4 and CO 2 transport through poly(dimethyl siloxane). The calculations show that in the pre-steady-state regime, permeation is only correctly described if the sorbed gas concentration in the polymer is dynamically determined by the rise in pressure. Finally, the framework is used to predict selectivity targets for two applications involving rubbery membranes: CO 2 capture from air and blocking of methane cross-over in an aged solar fuels device.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Calculation of the Influence of the Ion Current Density and Temperature on the Accumulation Kinetics of Point Defects under the Irradiation of Si with Light Ions

Numerical calculations of the accumulation kinetics of point defects—vacancies and divacancies—under the irradiation of Si with ions having masses of M{sub 1} ≤ 31 amu and energies of E ≤ 100 keV under various irradiation conditions are performed. The previously proposed diffusion-coagulation model is used without application of the “weak diffusion” approximation, which was performed during its analytical implementation. The main peculiarities of the dependences of the concentrations of vacancies and divacancies on the dose, ion-current density, and temperature under irradiation are analyzed. A physical interpretation of these results is given. The developed computing complex is rather flexible and makes it possible to analyze the influence of model input parameters by means of their variation and include additional processes into consideration if necessary.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Unifying Advective and Diffusive Descriptions of Bedform Pumping in the Benthic Biolayer of Streams

Abstract Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment‐water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow‐dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments.

Grant, Stanley B.↗
DOE OSTI2015

THE DETECTION OF A SN IIn IN OPTICAL FOLLOW-UP OBSERVATIONS OF ICECUBE NEUTRINO EVENTS

The IceCube neutrino observatory pursues a follow-up program selecting interesting neutrino events in real-time and issuing alerts for electromagnetic follow-up observations. In 2012 March, the most significant neutrino alert during the first three years of operation was issued by IceCube. In the follow-up observations performed by the Palomar Transient Factory (PTF), a Type IIn supernova (SN IIn) PTF12csy was found 0.°2 away from the neutrino alert direction, with an error radius of 0.°54. It has a redshift of z = 0.0684, corresponding to a luminosity distance of about 300 Mpc and the Pan-STARRS1 survey shows that its explosion time was at least 158 days (in host galaxy rest frame) before the neutrino alert, so that a causal connection is unlikely. The a posteriori significance of the chance detection of both the neutrinos and the SN at any epoch is 2.2σ within IceCube's 2011/12 data acquisition season. Also, a complementary neutrino analysis reveals no long-term signal over the course of one year. Therefore, we consider the SN detection coincidental and the neutrinos uncorrelated to the SN. However, the SN is unusual and interesting by itself: it is luminous and energetic, bearing strong resemblance to the SN IIn 2010jl, and shows signs of interaction of the SN ejecta with a dense circumstellar medium. High-energy neutrino emission is expected in models of diffusive shock acceleration, but at a low, non-detectable level for this specific SN. In this paper, we describe the SN PTF12csy and present both the neutrino and electromagnetic data, as well as their analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Adsorption Performance and Mechanisms of Methylene Blue Removal by Non-magnetic and Magnetic Particles Derived from the Vallisneria natans Waste

Vallisneria natans (V. natans) is a submerged aquatic plant, and its dried biomass waste was used to prepare magnetic particle biomass adsorbent by chemical co-precipitation via FeSO{sub 4}·7H{sub 2}O and FeCl{sub 3}·6H{sub 2}O as iron precursors. Brunauer–Emmett–Teller specific surface area, total pore volume, and average pore width were 7.434 m{sup 2} g{sup −1}, 0.013 cm{sup 3} g{sup −1}, and 6.992 nm for magnetic V. natans, correspondingly, and 3.575 m{sup 2} g{sup −1}, 0.004 cm{sup 3} g{sup −1}, and 4.476 nm for non-magnetic V. natans, respectively. The adsorption kinetics of methylene blue (MB) is a more appropriate fit for pseudo-second-order kinetic model than pseudo-first-order and intraparticle diffusion kinetic models. Isotherm results showed that the adsorption of MB onto the non-magnetic and magnetic adsorbents agrees well with Dubinin–Radushkevich isotherm model. Moreover, the thermodynamics study implied that the adsorption of MB onto both adsorbents is a spontaneous and exothermic process. The maximum adsorption capacities of MB on non-magnetic and magnetic V. natans were 657.90 and 473.93 mg g{sup −1} at 303 K, correspondingly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗