Solar Wind Magnetic Field Correlation Length: Correlation Functions versus Cross-field Displacement Diffusivity Test
Explore the source record for details and available documents.
B. R. Ragot↗
SEARCH · Engineering Papers
Results for “diffusion length”
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.
A Proof of the Asymptotic Variance of Path Length Estimators for Single-Collision Monte Carlo Source Iteration in the Thick Diffusion Limit
Here, we prove a theorem relating the variance of path length estimators for single-collision Monte Carlo source iteration to a parameter that becomes infinitesimally small in an important physical regime arising in radiative transfer. In our usage, “single-collision Monte Carlo source iteration” refers to Monte Carlo Boltzmann transport methods in which each Monte Carlo particle history includes no more than a single collision, and the physics of multiple scattering is modeled by lagging the scattering source term and iterating until this term converges. Our theorem can be used to construct variance reduction techniques which improve the order of the estimator variance. This enables calculations that would otherwise require impractically large sample sizes to achieve practical estimator uncertainties. We believe this is the first postulation of a theorem relating estimator variance to a limiting case parameter for single-collision Monte Carlo source iteration, and the first proof of such a theorem. We illustrate the theorem’s value with an example in which the authors of a transport method used the theorem to design a variance reduction technique that improved the uncertainty of their solution by a factor of about 500 for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material.
Mathematics and Computing↗
The effects of the Asselin time filter on numerical solutions to the linearized shallow-water wave equations
In the present investigation, a one-dimensional linearized analysis is used to determine the effect of Asselin's (1972) time filter on both the computational stability and phase error of numerical solutions for the shallow water wave equations, in cases with diffusion but without rotation. An attempt has been made to establish the approximate optimal values of the filtering parameter nu for each of the 'lagged', Dufort-Frankel, and Crank-Nicholson diffusion schemes, suppressing the computational wave mode without materially altering the physical wave mode. It is determined that in the presence of diffusion, the optimum filter length depends on whether waves are undergoing significant propagation. When moderate propagation is present, with or without diffusion, the Asselin filter has little effect on the spatial phase lag of the physical mode for the leapfrog advection scheme of the three diffusion schemes considered.
Schlesinger, R. E.↗
A new method of optimal design for a two-dimensional diffuser by using dynamic programming
A new method for predicting the optimal velocity distribution on the wall of a two dimensional diffuser is presented. The method uses dynamic programming to solve the optimal control problem with inequality constraints of state variables. The physical model of optimization is designed to prevent the separation of the boundary layer while approaching the maximum pressure ratio in a diffuser of a specified length. The computational results are in fair agreement with the experimental ones. Optimal velocity distribution on a diffuser wall is said to occur when the flow decelerates quickly at first and then smoothly, while the flow is near separation, but always protected from it. The optimal velocity distribution can be used to design the contour of the diffuser.
Gu, Chuangang↗
Strain-Rate-Free Diffusion Flames: Initiation, Properties, and Quenching
For about a half century, the stabilization of a steady planar deflagration on a heat-sink-type flat-flame burner has been of extraordinary service for the theoretical modeling and diagnostic probing of combusting gaseous mixtures. However, most engineering devices and most unwanted fire involve the burning of initially unmixed reactants. The most vigorous burning of initially separated gaseous fuel and oxidizer is the diffusion flame. In this useful idealization (limiting case), the reactants are converted to product at a mathematically thin interface, so no interpenetration of fuel and oxidizer occurs. This limit is of practical importance because it often characterizes the condition of optimal performance (and sometimes environmentally objectionable operation) of a combustor. A steady planar diffusion flame is most closely approached in the laboratory in the counterflow apparatus. The utility of this simple-strain-rate flow for the modeling and probing of diffusion flames was noted by Pandya and Weinberg 35 years ago, though only in the last decade or so has its use become internationally common place. However, typically, as the strain rate a is reduced below about 20 cm(exp -1), and the diffusion-flame limit (reaction rate much faster than the flow rate) is approached, the burning is observed to become unstable in earth gravity. The advantageous steady planar flow is not available in the diffusion-flame limit in earth gravity. This is unfortunate because the typical spatial scale in a counterflow is (k/a)(sup 1/2), where k denotes a characteristic diffusion coefficient; thus, the length scale becomes large, and the reacting flow is particularly amenable to diagnostic probing, as the diffusion-flame limit is approached. The disruption of planar symmetry is owing the fact that, as the strain rate a decreases, the residence time (l/a) of the throughput in the counterflow burner increases. Observationally, when the residence time exceeds about 50 msec, the inevitably present convective (Rayleigh-Benard) instabilities, associated with hot-under-cold (flame-under-fresh-reactant) stratification of fluid in a gravitational field, have time to grow to finite amplitude during transit of the burner.
Fendell, Francis↗
Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)
In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.
03 NATURAL GAS↗
An experimental investigation of two large annular diffusers with swirling and distorted inflow
Two annular diffusers downstream of a nacelle-mounted fan were tested for aerodynamic performance, measured in terms of two static pressure recovery parameters (one near the diffuser exit plane and one about three diameters downstream in the settling duct) in the presence of several inflow conditions. The two diffusers each had an inlet diameter of 1.84 m, an area ratio of 2.3, and an equivalent cone angle of 11.5, but were distinguished by centerbodies of different lengths. The dependence of diffuser performance on various combinations of swirling, radially distorted, and/or azimuthally distorted inflow was examined. Swirling flow and distortions in the axial velocity profile in the annulus upstream of the diffuser inlet were caused by the intrinsic flow patterns downstream of a fan in a duct and by artificial intensification of the distortions. Azimuthal distortions or defects were generated by the addition of four artificial devices (screens and fences). Pressure recovery data indicated beneficial effects of both radial distortion (for a limited range of distortion levels) and inflow swirl. Small amounts of azimuthal distortion created by the artificial devices produced only small effects on diffuser performance. A large artificial distortion device was required to produce enough azimuthal flow distortion to significantly degrade the diffuser static pressure recovery.
Eckert, W. T.↗
Large Bandwidth and Low Noise in a Diffusion-Cooled Hot-Electron Bolometer Mixer
X Heterodyne measurements have been made at 533 GHz using a novel superconducting hot-electron bolometer in a waveguide mixer. The bolometer is a 0.3 ??ong niobium microbridge with T c =(approx. equal to)5K. The short length ensures that electron diffusion dominates over electron-phonon interactions as the electron cooling mechanism, which should allow heterodyne detection with intermediate frequencies (if{italics})) of several GHz. A Y-factor response of 1.15dB has been obtained at an if(italics) of 1.4 GHz with 77k and 295k loads, indicating a receiver noise temperature of 650 k DSB. The -3 dB rolloff in the if(italics) response occurs at 1.7 GHz.
Heterodyne↗
Large Bandwidth and Low Noise in a Diffusion-Cooled Hot-Electron Bolometer Mixer
Heterodyne measurements have been made at 533 GHz using a novel superconducting hot-electron bolometer in a waveguide mixer. The bolometer is a 0.3 micrometer long niobium microbridge with a superconducting transition temperature of 5 K. The short length ensures that electron diffusion dominates over electron-phonon interactions as the electron cooling mechanism, which should allow heterodyne detection with intermediate frequencies (if) of several GHz. A Y-factor response of 1.15 dB has been obtained at an if of 1.4 GHz with 77 and 295 K loads, indicating a receiver noise temperature of 650 K DSB. The -3 dB rolloff in the if response occurs at 1.7 GHz.
Skalare, A.↗
Role of a Multivalent Ion–Solvent Interaction on Restricted Mg 2+ Diffusion in Dimethoxyethane Electrolytes
The diffusion behavior of Mg 2+ in electrolytes is not as readily accessible as that from Li + or Na + utilizing PFG NMR, due to the low sensitivity, poor resolution, and rapid relaxation encountered when attempting 25 Mg NMR. In MgTFSI 2 /DME solutions, “bound” DME (coordinating to Mg 2+ ) and “free” DME (bulk) are distinguishable from 1 H NMR. With the exchange rates between them obtained from 2D 1 H EXSY NMR, we can extract the self-diffusivities of free DME and bound DME (which are equal to that of Mg 2+ ) before the exchange occurs using PFG diffusion NMR measurements coupled with analytical formulas describing diffusion under two-site exchange. Further, the high activation enthalpy for exhange (65–70 kJ/mol) can be explained by the structural change of bound DME as evidenced by its reduced C–H bond length. Comparison of the diffusion behaviors of Mg 2+ , TFSI – , DME, and Li + reveals a relative restriction to Mg 2+ diffusion that is caused by the long-range interaction between Mg 2+ and solvent molecules, especially those with suppressed motions at high concentrations and low temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Investigation of Perforated Convergent-divergent Diffusers with Initial Boundary Layer
An experimental investigation was made at Mach number 1.90 of the performance of a series of perforated convergent-divergent supersonic diffusers operating with initial boundary layer, which was induced and controlled by lengths of cylindrical inlets affixed to the diffusers. Supercritical mass-flow and peak total-pressure recoveries were decreased slightly by use of the longest inlets (4 inlet diameters in length). Combinations of cylindrical inlets, perforated diffusers, and subsonic diffuser were evaluated as simulated wind tunnels having second throats. Comparisons with noncontracted configurations of similar scale indicated conservatively computed power reductions of 25 percent.
DIFFUSERS, SUPERSONIC↗
Henry’s Solubility and Diffusion Coefficients for 29 Volatile Organic Compounds in Polydimethylsiloxane Sylgard 184 at 293 K
Two-dimensional (2D) inverse gas chromatography (IGC) enables simultaneous determination of Henry’s solubility and Fickian diffusion coefficients for volatile organic compounds (VOCs) in polymer films. This technique offers a significant advantage over traditional cylindrical column IGC by providing precise control and measurement of the film thickness (here, 0.064 ± 0.002 mm), which is the critical length scale for accurate diffusivity determination. We apply this methodology to characterize VOC transport in Sylgard 184, a widely used polydimethylsiloxane (PDMS)-based polymer containing substantial silica filler content. At room temperature (20 °C), we measured solubility and diffusion coefficients for 29 common VOCs spanning diverse chemical functionalities, including alkanes, aromatics, chlorinated solvents, ketones, esters, and alcohols. Comparison with literature data for pure PDMS reveals that VOC solubility in Sylgard 184 is generally higher; for most non-hydrogen-bonding compounds it remains within a factor of 2 of pure PDMS, whereas alcohols are enhanced by roughly 1.8 to 3.7 times, which we attribute to favorable interactions with residual silanol groups on the silanized silica filler. Diffusion coefficients range from 1.0 × 10 –6 cm 2 /s (n-undecane) to 8.9 × 10–5 cm 2 /s (acetonitrile) and align well with extrapolated literature values for PDMS. This study provides essential thermodynamic and transport data for predicting VOC permeation in Sylgard 184 while demonstrating the utility of 2D IGC as a robust technique for characterizing rubbery polymer membranes across diverse industrial applications.
organic↗
Crack diffusion coefficient - A candidate fracture toughness parameter for short fiber composites
In brittle matrix composites, crack propagation occurs along random trajectories reflecting the heterogeneous nature of the strength field. Considering the crack trajectory as a diffusive process, the 'crack diffusion coefficient' is introduced. From fatigue crack propagation experiments on a set of identical SEN polyester composite specimens, the variance of the crack tip position along the loading axis is found to be a linear function of the effective 'time'. The latter is taken as the effective crack length. The coefficient of proportionality between variance of the crack trajectory and the effective crack length defines the crack diffusion coefficient D which is found in the present study to be 0.165 mm. This parameter reflects the ability of the composite to deviate the crack from the energetically most efficient path and thus links fracture toughness to the microstructure.
Mull, M. A.↗
The propagation of solar flare particles in a turbulent coronal loop
Energetic solar flare particles, both electrons and protons, must survive the turbulent environment of a flaring loop and propogate to the lower corona or chromosphere in order to produce hard X-ray and gamma ray bursts. This plasma turbulence, often observed in soft X-ray line widths to be in excess of 100 km/s, is presumably capable of efficiently scattering the fast flare particles. This prevents to some degree the free streaming of accelerated particles and depending on the amplitude of the turbulence, restricts the particles to diffusive propagation along the length of the loop to the target chromosphere. In addition this turbulence is capable of performing additional acceleration of the fast particles by the second order Fermi mechanism. For compact flares with rise times 10s, the acceleration effect is small and the propagation of the particles is governed by spatial diffusion and energy loss in the ambient medium. The solution of the time dependent diffusion equation with velocity dependent diffusion and energy loss coefficients yields for the case of nonrelativistic protons particle precipitation rates which are necessary for calculating thick target gamma ray emission and also yields the total thin target emissivity.
Ryan, J. M.↗
Role of pressure diffusion in non-homogeneous shear flows
A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.
Demuren, A. O.↗
A Full Navier-Stokes Analysis of Subsonic Diffuser of a Bifurcated 70/30 Supersonic Inlet for High Speed Civil Transport Application
A full Navier-Stokes analysis was performed to evaluate the performance of the subsonic diffuser of a NASA Lewis Research Center 70/30 mixed-compression bifurcated supersonic inlet for high speed civil transport application. The PARC3D code was used in the present study. The computations were also performed when approximately 2.5 percent of the engine mass flow was allowed to bypass through the engine bypass doors. The computational results were compared with the available experimental data which consisted of detailed Mach number and total pressure distribution along the entire length of the subsonic diffuser. The total pressure recovery, flow distortion, and crossflow velocity at the engine face were also calculated. The computed surface ramp and cowl pressure distributions were compared with experiments. Overall, the computational results compared well with experimental data. The present CFD analysis demonstrated that the bypass flow improves the total pressure recovery and lessens flow distortions at the engine face.
Kapoor, Kamlesh↗
Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data
Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.
Reimer, C↗
Understanding the effect of irradiation and silicon carbide microstructure on silver diffusion in TRISO particles
The silicon carbide (SiC) layer in tristructural isotropic (TRISO) fuel particles serves as a barrier to prevent the escape of fission from the fuel kernel. The release of silver (Ag) is a concern due to the long half-life of the 110mAg isotope. In this study, a combination of atomistic and phase-field methods focuses on quantifying the effect of irradiation and SiC microstructure on Ag diffusion. Atomistic calculations provide the Ag diffusivity in SiC bulk and GBs, also providing an estimate of the effect of irradiation defects on Ag diffusivity in bulk SiC. These lower length scale efforts are leveraged to develop a mesoscale effective Ag diffusion coefficient (Deff) in SiC. Deff is defined as a function of temperature, microstructure variables, and fluence. Deff is implemented in the fuel performance code Bison to predict Ag release from AGR-1 and AGR-2 TRISO fuel particles. We hereby quantify the impact of SiC grain size and irradiation on Ag release and improve Bison's predictions.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗