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,261 records · Page 70

Erosion, Transportation, and Deposition on Outer Solar System Satellites: Landform Evolution Modeling Studies

Mass movement and landform degradation reduces topographic relief by moving surface materials to a lower gravitational potential. In addition to the obvious role of gravity, abrasive mechanical erosion plays a role, often in combination with the lowering of cohesion, which allows disaggregation of the relief-forming material. The identification of specific landform types associated with mass movement and landform degradation provides information about local sediment particle size and abundance and transportation processes. Generally, mass movements can be classified in terms of the particle sizes of the transported material and the speed the material moved during transport. Most degradation on outer planet satellites appears consistent with sliding or slumping, impact erosion, and regolith evolution. Some satellites, such as Callisto and perhaps Hyperion and Iapetus, have an appearance that implies that some additional process is at work, most likely sublimation-driven landform modification and mass wasting. A variant on this process is thermally driven frost segregation as seen on all three icy Galilean satellites and perhaps elsewhere. Titan is unique among outer planet satellites in that Aeolian and fluvial processes also operate to erode, transport, and deposit material. We will evaluate the sequence and extent of various landform-modifying erosional and volatile redistribution processes that have shaped these icy satellites using a 3-D model that simulates the following surface and subsurface processes: 1) sublimation and re-condensation of volatiles; 2) development of refractory lag deposits; 3) disaggregation and downward sloughing of surficial material; 4) radiative heating/cooling of the surface (including reflection, emission, and shadowing by other surface elements); 5) thermal diffusion; and 6) vapor diffusion. The model will provide explicit simulations of landform development and thusly predicts the topographic and volatile evolution of the surface and final landscape form as constrained by DEMs. We have also simulated fluvial and lacustrine modification of icy satellites landscapes to evaluate the degree to which fluvial erosion of representative initial landscapes can replicate the present Titan landscape.

erosion↗

HYSPLIT/MACCS Atmospheric Dispersion Model Technical Documentation and Benchmark Analysis

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

97 MATHEMATICS AND COMPUTING↗

Illumination in diverse codimensions

This paper derives a model of diffuse and specular illumination in arbitrarily large dimensions, based on a few characteristics of material and light in three-space. It then describes how to adjust for the anomaly of excess brightness in large codimensions. If a surface is grooved or furry, it can be illuminated with a hybrid model that incorporates both the one dimensional geometry (the grooves or fur) and the two dimensional geometry (the surface).

Banks, David C.↗

Charge decay in the spatial afterglow of plasmas and its impact on diffusion regimes

Abstract The spatial afterglow is a region at the boundary of a non-equilibrium plasma where charged species relax into ambient equilibrium. In many applications, the spatial afterglow is the part of the plasma that interacts with surfaces, such as suspended particles or a material substrate. However, compared to the bulk plasma, there has been little effort devoted to studying the properties of the spatial afterglow, and a fundamental analysis has not yet been developed. Here, we apply double Langmuir probe measurements and develop an advection-diffusion-recombination model to provide a detailed description of charged species in the spatial afterglow over a wide range of pressures, temperatures, plasma dimensions, and flow rates. We find that the density of charged species in the spatial afterglow decays by orders of magnitude, which leads to a transition from ambipolar to free diffusion. These insights can be used to explain or predict experimental observations of phenomena, such as the charging of dust grains and the dose of charged species to a biomaterial.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cherenkov Telescope Array sensitivity to the putative millisecond pulsar population responsible for the Galactic Centre excess

ABSTRACT The leading explanation of the Fermi Galactic Centre γ-ray excess is the extended emission from an unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt γ-rays, also inject large quantities of electrons/positrons (e±) into the interstellar medium. These e± could potentially inverse-Compton (IC) scatter ambient photons into γ-rays that fall within the sensitivity range of the upcoming Cherenkov Telescope Array (CTA). In this article, we examine the detection potential of CTA to this signature by making a realistic estimation of the systematic uncertainties on the Galactic diffuse emission model at TeV-scale γ-ray energies. We forecast that, in the event that e± injection spectra are harder than E−2, CTA has the potential to robustly discover the IC signature of a putative Galactic bulge MSP population sufficient to explain the Galactic Centre excess for e± injection efficiencies in the range of ≈2.9–74.1 per cent, or higher, depending on the level of mismodelling of the Galactic diffuse emission components. On the other hand, for spectra softer than E−2.5, a reliable CTA detection would require an unphysically large e± injection efficiency of ${\gtrsim} 158{{\ \rm per\ cent}}$. However, even this pessimistic conclusion may be avoided in the plausible event that MSP observational and/or modelling uncertainties can be reduced. We further find that, in the event that an IC signal were detected, CTA can successfully discriminate between an MSP and a dark matter origin for the radiating e±.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamics of solar cosmic ray events - Processes at large heliocentric distances /much greater than 1 AU/

Observations of solar cosmic ray events at heliocentric distances up to 6 AU and beyond obtained by the Pioneer 10 and 11 spacecraft are discussed in terms of the propagation of energetic charged particles in interplanetary space. Following a review of the diffusion propagation model used to explain cosmic ray transport, the statistical studies of McCarthy and O'Gallagher (1976) and Zwickl and Webber (1977) of the relations of event parameters with radial distance and the simultaneous observation studies of Hamilton (1977) are reviewed, and it is noted that the results imply a slowly increasing radial diffusion coefficient out to about 6 AU. More recent analyses of data obtained at heliocentric distances greater than 10 AU are then presented which indicate that the coefficient of radial diffusion may actually be decreasing with radial distance beyond 5 AU. Finally, theoretical predictions of the radial variation of the diffusion coefficient are presented which take into account the background interplanetary medium and are shown to be in agreement with observations.

Hamilton, D. C.↗

How much hydrogen is there in a white dwarf?

Stratified hydrogen/helium envelope models in diffusive equilibrium are calculated for a 0.6-solar-mass white dwarf for effective temperatures between 10,000 and 80,000 K in order to investigate the observational constraints placed on the total hydrogen mass. Convective mixing is included ab initio in the calculations, and synthetic spectra are used for comparing these models with observational materials. It is shown that evolutionary changes in the surface composition of white dwarfs cannot be explained by a model in which a small amount of hydrogen floats to the surface from initially being mixed in the outer parts of a helium envelope. It is pointed out that the shape of the hydrogen lines can be used for constraining theories of convective overshoot.

Macdonald, James↗

Diffusion-driven growth of nanowires by low-temperature molecular beam epitaxy

With ZnTe as an example, we use two different methods to unravel the characteristics of the growth of nanowires (NWs) by gold-catalyzed molecular beam epitaxy at low temperature. In the first approach, CdTe insertions have been used as markers, and the nanowires have been characterized by scanning transmission electron microscopy, including geometrical phase analysis and energy dispersive electron spectrometry; the second approach uses scanning electron microscopy and the statistics of the relationship between the length of the tapered nanowires and their base diameter. Axial and radial growth are quantified using a diffusion-limited model adapted to the growth conditions; analytical expressions describe well the relationship between the NW length and the total molecular flux (taking into account the orientation of the effusion cells), and the catalyst-nanowire contact area. A long incubation time is observed. This analysis allows us to assess the evolution of the diffusion lengths on the substrate and along the nanowire sidewalls, as a function of temperature and deviation from stoichiometric flux.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffusive mass transport in agglomerated glassy fallout from a near-surface nuclear test

Aerodynamically-shaped glassy fallout is formed when vapor phase constituents from the nuclear device are incorporated into molten carriers (i.e. fallout precursor materials derived from soil or other near-field environmental debris). The effects of speciation and diffusive transport of condensing constituents are not well defined in models of fallout formation. Previously we reported observations of diffuse micrometer scale layers enriched in Na, Fe, Ca, and 235 U, and depleted in Al and Ti, at the interfaces of agglomerated fallout objects. Here in this paper, we derive the timescales of uranium mass transport in such fallout as it cools from 2500 K to 1500 K by applying a 1-dimensional planar diffusion model to the observed 235 U/ 30 Si variation at the interfaces. By modeling the thermal transport between the fireball and the carrier materials, the time of mass transport is calculated to be <0.6 s, <1 s, <2 s, and <3.5 s for fireball yields of 0.1 kt, 1 kt, 10 kt, and 100 kt respectively. Based on the calculated times of mass transport, a maximum temperature of deposition of uranium onto the carrier material of ~2200 K is inferred (1σ uncertainty of ~200 K). We also determine that the occurrence of micrometer scale layers of material enriched in relatively volatile Na-species as well as more refractory Ca-species provides evidence for an oxygen-rich fireball based on the vapor pressure of the two species under oxidizing conditions. These results represent the first application of diffusion-based modeling to derive material transport, thermal environments, and oxidation-speciation in near-surface nuclear detonation environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A search with Copernicus for interstellar N2 in diffuse clouds

Multiple Copernicus scans of two N2 band regions (near 958.5 and 960.2A) of Delta Sco and Epsilon Per are reported. The observations indicate upper limits for the number of N2 molecules equal to 1.0-3.8 times 10 to the -12th/sq cm and 1.2-4.4 times 10 to the -12th/sq cm, respectively; the limits depend on the cloud temperature. It is suggested that the limits are consistent with the column densities predicted by chemical models for diffuse interstellar clouds, and the predicted relative abundances are presented in terms of the ratio of N(N2)/(2N(H2) + N(Hl)).

Lutz, B. L.↗

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional computational fluid dynamic (CFD) flow solver with the Monte Carlo simulation of the pdf evolution equation was developed. The algorithm was validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames was carried out using this algorithm. Numerical results compared favorably with experimental data. The computations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed in this work. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional CFD flow solver with the Monte Carlo simulation of the pdf evolution equation has been developed. The algorithm has been validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames has been carried out using this algorithm. Numerical results compared favorably with experimental data. The computuations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗

An Investigation of Thermal Properties of 2D Materials [Dissertation]

Studying the thermal conductivity of 2D materials is important due to the applications of 2D materials in fields such as thermal management, thermoelectricity, renewable energy, and sensors. As such, measurements of the thermal conductivity of these 2D materials become important to measure. Thermal conductivity is often difficult to measure for 2D materials due to their atomically thin nature and many experimental methods for doing so requiring contact with the sample, which can alter the thermal properties. A non-contact method for calculating the thermal conductivity of 2D materials supported on substrates in order to model the thermal conductivity of 2D materials for devices, is proposed and experimentally performed in this dissertation. The optothermal Raman technique is a useful non-contact diagnostic technique useful in determining the thermal conductivity of 2D materials. The optothermal Raman typically does not account for heat losses due to convection or radiation or substrate resistance, which are shown to be important factors to consider when developing an optothermal Raman model. Additionally, the calculation of the interfacial thermal conductance between the bottom surface of the sample and the top surface of the substrate, plays an important role in determining the final value of the thermal conductivity of a supported sample, and will yield differing results based on whether or not the conductance is calculated using an approach such as the Diffuse Mismatch Model (DMM) or calculated directly by varying the laser heating profile (usually done by changing the laser objective). This is shown to be the case for both graphene on Ni, graphene on Cu, and SnSe 2 on Cu. In addition to experimentally calculating the thermal conductivity of a 2D material with the optothermal Raman technique, the thermal conductivity of 2D materials can also be calculated using computational methods. The three-phonon method is a method which can be used to simulate phonon scattering processes and determine the thermal conductivity of semiconductors, wherein phonon scattering is the dominant mechanism which determines the thermal conductivity. The three-phonon method uses relaxation times for phonon scattering with other phonons, electrons, and other material system elements, such as isotopes or material defects, in order to create a single-mode relaxation time approximation (SMRTA), which is used to calculate the final value of the thermal conductivity. An important consideration when determining the thermal conductivity of a 2D material using this method is the device geometry, which is reflected in this work as the phonon-boundary scattering relaxation time. This inclusion is important along with the inclusion of phonon-electron scattering in accurately determining the thermal conductivity of a 2D material. In both the optothermal Raman experiments and the three-phonon method computations, strain is shown to have a demonstrable effect on the thermal conductivity of 2D materials. When a 1.1% strain was applied to the mechanical properties of SnSe, the three-phonon processes yielded a lower thermal conductivity than the no-strain case. For the optothermal Raman experiments, the strain induced in the Cu substrate and transferred to a single-layer graphene (SLG) sample yields a trend where the thermal conductivity of the SLG decreases with respect to strain applied. In the case where the interfacial thermal conductance was calculated directly, the conductance increased with respect to strain applied. This presents strain as a reliable and viable method for tuning the thermal properties of 2D materials for device applications.

36 MATERIALS SCIENCE↗

A microstructure sensitive grain boundary sliding and slip based constitutive model for machining of $\mathrm{Ti}$-$\mathrm{6Al}$-$\mathrm{4V}$

We report a composite dual phase internal state variable constitutive model was developed for Ti-6Al-4V. The proposed model includes diffusion assisted grain boundary sliding based physics in addition to a traditional slip-based plasticity. Influence of microstructure on the flow stress is introduced via dislocation density and mean grain size internal state variables. The dislocation density evolves according to a physics based law that considers dislocation nucleation and annihilation processes. Grain refinement is driven by dynamic recrystallization, which is modeled phenomenologically. The model is calibrated with uniaxial stress–strain data that ranges between quasi-static and dynamic rates across a wide range of temperatures. Validation against machining data shows that the model predicts chip segmentation frequency, machining forces, and tool temperatures reasonably well. The newly introduced grain boundary sliding physics was found to dominate deformation following sufficient grain refinement. This deformation mode provides softening at the constitutive level without the need for invoking damage based softening mechanisms. This physical interpretation is something that has not previously been explored in the machining literature.

36 MATERIALS SCIENCE↗

ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES) - Merging Observations and Laboratory Data

The Diffuse Interstellar Bands (DIBs) are a set of ~500 absorption bands that are detected in the spectra of stars with interstellar clouds in the line of sight. DIBs are found from the NUV to the NIR in the spectra of reddened stars spanning different interstellar environments in our local, and in other galaxies. DIB carriers are a significant part of the interstellar chemical inventory. They are stable and ubiquitous in a broad variety of environments and play a unique role in interstellar physics/chemistry. It has long been realized that the solving of the DIB problem requires a strong synergy between astronomical observations, laboratory astrophysics, and astrophysical modeling of line-of-sights. PAHs are among the molecular species that have been proposed as DIB carriers. We will present an assessment of the PAH-DIB model in view of the progress and the advances that have been achieved over the past years through a series of studies involving astronomical observations of DIBs, laboratory simulation of interstellar analogs for neutrals and ionized PAHs, theoretical calculations of PAH spectra and the modelization of diffuse and translucent interstellar clouds. We will present a summary of what has been learned from these complementary studies, the constraints that can now be derived for the PAHs as DIB carriers in the context of the PAH-DIB model and how these constraints can be applied to the EDIBLES project. The spectra of several neutral and ionized PAHs isolated in the gas phase at low temperature have been measured in the laboratory under experimental conditions that mimic interstellar conditions and are compared with an extensive set of astronomical spectra of reddened, early type stars. The comparisons of astronomical and laboratory data provide upper limits for the abundances of specific neutral PAH molecules and ions along specific lines-of-sight. Something that is not attainable from infrared observations alone. We present the characteristics of the laboratory facilities, MIS and COSmIC, that have been developed for this study and discuss the findings resulting from the comparison of the laboratory data with high resolution, high S/N ratio astronomical observations. MIS stands for Matrix Isolation Spectroscopy, a well-proven technique for isolating cold molecular species in inert solid environments. COSmIC stands for Cosmic Simulation Chamber. It combines a supersonic free jet expansion with discharge plasma and high-sensitivity cavity ringdown spectroscopy and time-of-flight mass spectrometry detection tools for the generation and the detection of cold, isolated gas-phase molecules and ions under experimental conditions that closely mimic interstellar conditions. The column densities of the individual neutral PAH molecules and ions probed in these surveys are derived from the comparison of these unique laboratory data with high resolution, high S/N ratio astronomical observations. The comparisons of astronomical and laboratory data lead to clear and unambiguous conclusions regarding the expected abundances for PAHs of various sizes and charge states in the interstellar environments probed in the surveys. Band profile comparisons between laboratory and astronomical spectra lead to information regarding the molecular structures and characteristics associated with the DIB carriers in the corresponding lines-of-sight. These quantitative surveys of neutral and ionized PAHs in the optical range open the way for unambiguous quantitative searches of PAHs and complex organics in a variety of interstellar and circumstellar environments.

astrophysics↗

Physics-Informed Deep Neural Networks for Learning Parameters and Constitutive Relationships in Subsurface Flow Problems

In this paper, we present a physics informed deep neural network (DNN) method for estimating parameters and unknown physics (constitutive relationships) in partial differential equation (PDE) models. We use PDEs in addition to measurements to train DNNs to approximate unknown parameters and constitutive relationships as well as states. The proposed approach increases the accuracy of DNN approximations of partially known functions when a limited number of measurements is available and allows for training DNNs when no direct measurements of the functions of interest are available. We employ physics informed DNNs to estimate the unknown space-dependent diffusion coefficient in a linear diffusion equation and an unknown constitutive relationship in a non-linear diffusion equation. For the parameter estimation problem, we assume that partial measurements of the coefficient and states are available and demonstrate that under these conditions, the proposed method is more accurate than state-of-the-art methods. For the non-linear diffusion PDE model with a fully unknown constitutive relationship (i.e., no measurements of constitutive relationship are available), the physics informed DNN method can accurately estimate the non-linear constitutive relationship based on state measurements only. Finally, we demonstrate that the proposed method remains accurate in the presence of measurement noise.

42 ENGINEERING↗

Ensemble Grouping Strategies for Embedded Stochastic Collocation Methods Applied to Anisotropic Diffusion Problems

Previous work has demonstrated that propagating groups of samples, called ensembles, together through forward simulations can dramatically reduce the aggregate cost of sampling-based uncertainty propagation methods [E. Phipps, M. D'Elia, H. C. Edwards, M. Hoemmen, J. Hu, and S. Rajamanickam, SIAM J. Sci. Comput., 39 (2017), pp. C162--C193]. However, critical to the success of this approach when applied to challenging problems of scientific interest is the grouping of samples into ensembles to minimize the total computational work. For example, the total number of linear solver iterations for ensemble systems may be strongly influenced by which samples form the ensemble when applying iterative linear solvers to parameterized and stochastic linear systems. In this paper we explore sample grouping strategies for local adaptive stochastic collocation methods applied to PDEs with uncertain input data, in particular canonical anisotropic diffusion problems where the diffusion coefficient is modeled by truncated Karhunen--Loève expansions. Finally, we demonstrate that a measure of the total anisotropy of the diffusion coefficient is a good surrogate for the number of linear solver iterations for each sample and therefore provides a simple and effective metric for grouping samples.

97 MATHEMATICS AND COMPUTING↗

FermI large area telescope third source catalog

We present the third Fermi Large Area Telescope (LAT) source catalog (3FGL) of sources in the 100 MeV–300 GeV range. Based on the first 4 yr of science data from the Fermi Gamma-ray Space Telescope mission, it is the deepest yet in this energy range. Relative to the Second Fermi LAT catalog, the 3FGL catalog incorporates twice as much data, as well as a number of analysis improvements, including improved calibrations at the event reconstruction level, an updated model for Galactic diffuse γ-ray emission, a refined procedure for source detection, and improved methods for associating LAT sources with potential counterparts at other wavelengths. The 3FGL catalog includes 3033 sources above 4σ significance, with source location regions, spectral properties, and monthly light curves for each. Of these, 78 are flagged as potentially being due to imperfections in the model for Galactic diffuse emission. Twenty-five sources are modeled explicitly as spatially extended, and overall 238 sources are considered as identified based on angular extent or correlated variability (periodic or otherwise) observed at other wavelengths. For 1010 sources we have not found plausible counterparts at other wavelengths. More than 1100 of the identified or associated sources are active galaxies of the blazar class; several other classes of non-blazar active galaxies are also represented in the 3FGL. Pulsars represent the largest Galactic source class. From source counts of Galactic sources we estimate that the contribution of unresolved sources to the Galactic diffuse emission is ∼3% at 1 GeV.

79 ASTRONOMY AND ASTROPHYSICS↗