SEARCH · Engineering Papers

Results for “diffuser”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Soot Formation in Hydrocarbon/Air Laminar Jet Diffusion Flames

Soot processes within hydrocarbon/air diffusion flames are important because they affect the durability and performance of propulsion systems, the hazards of unwanted fires, the pollutant and particulate emissions from combustion processes, and the potential for developing computational combustion. Motivated by these observations, this investigation involved an experimental study of the structure and soot properties of round laminar jet diffusion flames, seeking an improved understanding of soot formation (growth and nucleation) within diffusion flames. The present study extends earlier work in this laboratory concerning laminar smoke points (l) and soot formation in acetylene/air laminar jet diffusion flames (2), emphasizing soot formation in hydrocarbon/air laminar jet diffusion flames for fuels other than acetylene. In the flame system, acetylene is the dominant gas species in the soot formation region and both nucleation and growth were successfully attributed to first-order reactions of acetylene, with nucleation exhibiting an activation energy of 32 kcal/gmol while growth involved negligible activation energy and a collision efficiency of O.53%. In addition, soot growth in the acetylene diffusion flames was comparable to new soot in premixed flame (which also has been attributed to first-order acetylene reactions). In view of this status, a major issue is the nature of soot formation processes in diffusion flame involving hydrocarbon fuels other than acetylene. In particular, information is needed about th dominant gas species in the soot formation region and the impact of gas species other than acetylene on soot nucleation and growth.

Sunderland, P. B.↗

NASA NTRS1997

Spherical Diffusion Flames: Structure and Dynamics

The spherical geometry is the most suitable one to study fundamental issues concerning the structure and the dynamics of diffusion flames. From a theoretical point of view, it is the only geometry that permits the existence of a truly one dimensional stationary diffusion flame. A stationary planar diffusion flame with the fuel supplied upstream, at x approaches -infinity say, and the oxidant downstream, at x approaches +infinity is not possible. For a steady diffusion flame to exist, one must have nonzero fluxes of fuel and oxidant towards the flame. However, in the unlimited region behind the planar flame the only bounded solutions to the reaction-free convective-diffusive operator are constants. Hence the oxidant concentration behind the flame remains constant and there is no mechanism to generate the necessary flux towards the flame. A one-dimensional problem can be formulated if the reactants are supplied at finite locations; but the boundary conditions in this case introduce unnecessary complications and do not appropriately model the physical reality. Indeed, a planar diffusion flame can be established in the stagnation-point flow of two opposed jets but the flame in this case is stretched and the flow is essentially two-dimensional. The only stationary one-dimensional diffusion flame in an unlimited environment is therefore the spherical flame.

Matalon, Moshe↗

NASA NTRS1997

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↗

NASA NTRS1997

Effect of Crystal Defects on Minority Carrier Diffusion Length in 6H SiC Measured Using the Electron Beam Induced Current Method

We report values of minority carrier diffusion length in n-type 6H SiC measured using a planar Electron Beam Induced Current (EBIC) method. Values of hole diffusion length in defect free regions of n-type 6H SiC, with a doping concentration of 1.7El7 1/cu cm, ranged from 1.46 microns to 0.68 microns. We next introduce a novel variation of the planar method used above. This 'planar mapping' technique measured diffusion length along a linescan creating a map of diffusion length versus position. This map is then overlaid onto the EBIC image of the corresponding linescan, allowing direct visualization of the effect of defects on minority carrier diffusion length. Measurements of the above n-type 6H SiC resulted in values of hole diffusion length ranging from 1.2 micron in defect free regions to below 0.1 gm at the center of large defects. In addition, measurements on p-type 6H SiC resulted in electron diffusion lengths ranging from 1.42 micron to 0.8 micron.

Tabib-Azar, Massood↗

NASA NTRS1999

Ocean Turbulence: One-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 2

We develop and test a 1-point closure turbulence model with the following features: 1) we include the salinity field and derive the expression for the vertical turbulent diffusivities of momentum K(sub m) , heat K(sub h) and salt K(sub s) as a function of two stability parameters: the Richardson number R(sub i) (stratification vs. shear) and the Turner number R(sub rho) (salinity gradient vs. temperature gradient). 2) to describe turbulent mixing below the mixed layer (ML), all previous models have adopted three adjustable "background diffusivities" for momentum, heat and salt. We propose a model that avoids such adjustable diffusivities. We assume that below the ML, the three diffusivities have the same functional dependence on R( sub i) and R(sub rho) as derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking.measured by Gargett et al. The procedure frees the model from adjustable background diffusivities and indeed we employ the same model throughout the entire vertical extent of the ocean. 3) in the local model, the turbulent diffusivities K(sub m,h,s) are given as analytical functions of R(sub i) and R(sub rho). 5) the model is used in an O-GCM and several results are presented to exhibit the effect of double diffusion processes. 6) the code is available upon request.

Canuto, V. M.↗

NASA NTRS2004

MODIS Solar Diffuser Attenuation Screen Modeling Results

On-orbit calibration of the reflected solar bands on the EOS Moderate Resolution Imaging Spectroradiometer (MODIS) is accomplished by have the instrument view a high reflectance diffuse surface illuminated by the sun. For some of the spectral bands this proves to be much too bright a signal that results in the saturation of detectors designed for measuring low reflectance (ocean) surfaces signals. A mechanical attenuation device in the form of a pin hole screen is used to reduce the signals to calibrate these bands. The sensor response to solar illumination of the SD with and without the attenuation screen in place will be presented. The MODIS detector response to the solar diffuser is smooth when the attenuation screen is absent, but has structures up to a few percent when the attenuation screen is present. This structure corresponds to non-uniform illumination from the solar diffuser. Each pin hole produces a pin-hole image of the sun on the solar diffuser, and there are very many pin hole images of the sun on the solar diffuser for each MODIS detector. Even though there are very many pin-hole images of the sun on the solar diffuser, it is no longer perfectly uniformly illuminated. This non-uniformly illuminated solar diffuser produces intensity variation on the focal planes. The results of a very detailed simulation will be discussed which show how the illumination of the focal plane changes as a result of the attenuation, and the impacts on the calibration will be discussed.

Waluschka, Eugene↗

NASA NTRS2004

The Performance of a Subsonic Diffuser Designed for High Speed Turbojet-Propelled Flight

An initial-phase subsonic diffuser has been designed for the turbojet flowpath of the hypersonic x43B flight demonstrator vehicle. The diffuser fit into a proposed mixed-compression supersonic inlet system and featured a cross-sectional shape transitioning flowpath (high aspect ratio rectangular throat-to-circular engine face) and a centerline offset. This subsonic diffuser has been fabricated and tested at the W1B internal flow facility at NASA Glenn Research Center. At an operating throat Mach number of 0.79, baseline Pitot pressure recovery was found to be just under 0.9, and DH distortion intensity was about 0.4 percent. The diffuser internal flow stagnated, but did not separate on the offset surface of this initial-phase subsonic diffuser. Small improvements in recovery (+0.4 percent) and DH distortion (-32 percent) were obtained from using vane vortex generator flow control applied just downstream of the diffuser throat. The optimum vortex generator array patterns produced inflow boundary layer divergence (local downwash) on the offset surface centerline of the diffuser, and an inflow boundary layer convergence (local upwash) on the centerline of the opposite surface.

Biesiadny, Thomas J.↗

NASA NTRS1950

High-Subsonic Performance Characteristics and Boundary-Layer Investigations of a 12 10-Inch-Inlet-Diameter Conical Diffuser

Performance and boundary-layer data were taken in a 12 degree 10-inch inlet-diameter conical diffuser of 2:1 exit- to inlet-area ratio. These data were taken for two inlet-boundary-layer conditions. The first condition was that of a thinner inlet boundary later (boundary-layer displacement thickness, delta* approximately equal to 0.034) produced by an inlet section approximately 1 inlet diameter in length between the entrance bell and the diffuser. The second condition was a thicker inlet boundary layer (delta* approximately equal to 0.120) produced by an additional inlet section length of approximately 6 diameters. Longitudinal static-pressure distributions were measured fro wall static orifices. Transverse total- and static-pressure surveys were made at the inlet and exit stations. Boundary-layer velocity distributions were measured at seven stations between the inlet and exit. These data were obtained for a Reynolds number (based on inlet diameter) range of 1 x 10(exp 6) to 3.9 x 10(exp 6). The corresponding Mach number range was from M = 0.2 to choking. At the maximum-power-available condition supersonic flow was obtained as far as 4.5 inches downstream from the diffuser inlet with a maximum Mach number of M approximately equal to 1.5. The total-pressure loss through the diffuser in percentage of inlet dynamic pressure was approximately 2.5 percent for the thinner inlet boundary later and 5.5 percent for the thicker inlet boundary later over the lower subsonic range. These valued increased with increasing flow rate- the values for the thicker inlet boundary later more than those for the thinner inlet boundary layer. The diffuser effectiveness, expressed as the ratio of the actual static-pressure rise to the ideal static-pressure rise, was about 85 percent for the thinner inlet boundary layer and about 67 percent for the thicker inlet boundary later in the lower subsonic range. These values decrease with increasing flow rate. Separated flow was observed for both inlet-boundary-layer conditions in the region of adverse pressure gradient just downstream of the transition curvature from inlet section to diffuser. The flow for the thinner-inlet-boundary-layer condition did not fully re-establish itself along the diffuser walls. The thicker inlet-boundary-layer flow, while not completely re-establishing the normal flow pattern downstream of the separated region, did re-establish more successfully than the thinner inlet boundary layer.

Little, B. H., Jr.↗

NASA NTRS1952

Amplitude of Supersonic Diffuser Flow Pulsations

A theoretical method for evaluating the stability characteristics and the amplitude and the frequency of pulsation of ram-jet engines without heat addition is presented herein. Experimental verification of the theoretical results are included where data were available. Theory and experiment show that the pulsation amplitude of a high mass-flow-ratio diffuser having no cone surface flow separation increases with decreasing mass flow. The theoretical trends for changes in amplitude, frequency, and mean-pressure recovery with changes in plenum-chamber volume were experimentally confirmed. For perforated convergent-divergent-type diffusers, a stability hysteresis loop was predicted on the pressure-recovery mass-flow-ratio curve. At a given mean mass-flow ratio, the higher.value of mean pressure recovery corresponded to oscillatory flow in the diffuser while the lower branch was stable. This hysteresis has been observed experimentally. The theory indicates that for a ram-jet engine of given diameter, the amplitude of pulsation of a supersonic diffuser is increased by decreasing the relative size of the plenum chamber with respect to the diffuser volume down to a critical value at which oscillations cease. In the region of these critical values, the stable mass-flow range of the diffuser may be increased either by decreasing the combustion chamber volume or by increasing the length of the diffuser.

Sterbentz, William H.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

NASA NTRS2013

Kinetic Monte Carlo Simulation of Cation Diffusion in Low-K Ceramics

Low thermal conductivity (low-K) ceramic materials are of interest to the aerospace community for use as the thermal barrier component of coating systems for turbine engine components. In particular, zirconia-based materials exhibit both low thermal conductivity and structural stability at high temperature, making them suitable for such applications. Because creep is one of the potential failure modes, and because diffusion is a mechanism by which creep takes place, we have performed computer simulations of cation diffusion in a variety of zirconia-based low-K materials. The kinetic Monte Carlo simulation method is an alternative to the more widely known molecular dynamics (MD) method. It is designed to study "infrequent-event" processes, such as diffusion, for which MD simulation can be highly inefficient. We describe the results of kinetic Monte Carlo computer simulations of cation diffusion in several zirconia-based materials, specifically, zirconia doped with Y, Gd, Nb and Yb. Diffusion paths are identified, and migration energy barriers are obtained from density functional calculations and from the literature. We present results on the temperature dependence of the diffusivity, and on the effects of the presence of oxygen vacancies in cation diffusion barrier complexes as well.

Good, Brian↗

NASA NTRS2016

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗

NASA NTRS2016

Generating a Simulated Fluid Flow over a Surface Using Anisotropic Diffusion

A fluid-flow simulation over a computer-generated surface is generated using a diffusion technique. The surface is comprised of a surface mesh of polygons. A boundary-layer fluid property is obtained for a subset of the polygons of the surface mesh. A gradient vector is determined for a selected polygon, the selected polygon belonging to the surface mesh but not one of the subset of polygons. A maximum and minimum diffusion rate is determined along directions determined using the gradient vector corresponding to the selected polygon. A diffusion-path vector is defined between a point in the selected polygon and a neighboring point in a neighboring polygon. An updated fluid property is determined for the selected polygon using a variable diffusion rate, the variable diffusion rate based on the minimum diffusion rate, maximum diffusion rate, and the gradient vector.

Rodriguez, David L.↗

NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Accurate and precise age determination of impact cratering events remains challenging and often contentious; less than half of all known craters are regarded as accurately and precisely dated. Zircon (U-Th)/He (ZHe) dating of impactites can be employed to date medium to large impact structures as ZHe ages can be fully reset in minutes at T >1000°C, a plausible scenario in the central melt pool. In contrast, complete resetting of ZHe at 200-300°C, encountered near the crater margins or due to post-impact hydrothermal overprinting, may take >103-6 years. To test the reliability of ZHe impact dating, we have quantified the effects of shock-induced microstructures on helium diffusion kinetics in well-characterized variably shocked zircon. We investigated samples from two impact structures, the Chicxulub multi-ring crater and Ries complex crater, to compare diffusion kinetics from structures with different size, age, and hydrothermal system longevity. Shock microstructures were characterized by backscattered-electron imaging prior to determining the He diffusion kinetics by prograde and retrograde fractional-release experiments via light-bulb furnace with incremental step-heating (10°C) from 300°C to 600°C. Next, we examine the internal interconnectivity and sizes of the diffusion domains by EBSD. While we found that zircon with few shock microstructures exhibited no marked deviation from helium diffusion kinetics of undamaged zircon, zircon grains with planar microstructures and granular textures are characterized by a dramatic decrease in helium retentivity due to the reduction in the effective domain size and the introduction of interconnected fast diffusion pathways. A subset of grains were ZHe dated and showed that less deformed grains yielded a weighted mean age within error of the accepted impact ages, while the grains with planar microstructures or granular textures gave systematically younger ages. These new diffusion data and ZHe ages demonstrate that highly shocked grains are unsuitable for ZHe impact crater dating. Therefore, detailed characterization of impact-induced microstructures is critical for determining accurate ZHe impact ages and offers the possibility of investigating post-impact hydrothermal circulation.

Zircon↗

NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Absolute age determination of impact cratering events remains difficult and often controversial; a challenge that has resulted in < 50% of known impact craters regarded as accurately and precisely dated. Besides conventional 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has been applied to large- to medium-sized impact structures. ZHe dates can be fully reset in minutes at 1000°C, which is commonly reached in central sections of the melt sheet, whereas resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-6 years. There is a critical need to quantify the effects of shock-induced microstructures and impact metamorphism on helium diffusion kinetics in well-characterized, variably shocked zircon to further establish the reliability of (U-Th)/He for dating impacts. For this purpose, we investigated suevite and impact melt samples from two impact structures, the Chicxulub multi-ring basin and the Ries complex crater, which enables us to compare zircon helium diffusion kinetics from impact structures with differing sizes, ages, and hydrothermal system longevities. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incremental 10°C steps from 250°C to 600°C. Afterward, we characterize the diffusion domain sizes and their interconnectivity within the shocked zircon grains using electron backscatter diffraction (EBSD). We find that zircon with few shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features (PDFs) and granular textures classified by BSE and EBSD are characterized by a dramatic decrease in helium retentivity, similar to radiation damaged grains, due to a reduction in the effective domain size and the introduction of interconnected fast diffusion pathways created by shock microstructures. A subset of grains were dated by ZHe after the external morphology of the grains was determined by BSE imaging. The euhedral grains yielded a weighted mean age within the uncertainty of the accepted impact ages, whereas the grains with PDFs or granular textures exhibited younger ages. Thus, these diffusion experiments and ZHe dates suggest that the dramatic decrease in domain size likely renders shocked grains more susceptible to impact-induced hydrothermal resetting and subsequent overprinting. Hence, characterization of shock microstructures is critical for determining accurate impact ages using ZHe methods especially when applied to previously unconstrained craters. The thermochronometer also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

zircon↗

DOE OSTI2016

Linking pore-scale and basin-scale effects on diffusive methane transport in hydrate bearing environments through multi-scale reservoir simulations

We explore the gas hydrate-generating capacity of diffusive methane transport induced by solubility gradients due to pore size contrasts in lithologically heterogeneous marine sediments. Through the use of 1D, 2D, and 3D reactive transport simulations, we investigate scale-dependent processes in diffusion-dominated gas hydrate systems. These simulations all track a sand body, or series of sands, surrounded by clays as they are buried through the gas hydrate stability zone. Methane is sourced by microbial methanogenesis in the clays surrounding the sand layers. In 1D, simulations performed in a Lagrangian reference frame demonstrate that gas hydrate in thin sands (3.6 m thick) can occur in high saturations (upward of 70%) at the edges of sand bodies within the upper 400 meters below the seafloor. Diffusion of methane toward the center of the sand layer depends on the concentration gradient within the sand: broader sand pore size distributions with smaller median pore sizes enhance diffusive action toward the sand’s center. Incorporating downhole log- and laboratory-derived sand pore size distributions, gas hydrate saturations in the center of the sand can reach 20% of the hydrate saturations at the sand’s edges. Furthermore, we show that hydrate-free zones exist immediately above and below the sand and are approximately 5 m thick, depending on the sand-clay solubility contrast. A moving reference frame is also adopted in 2D, and the angle of gravity is rotated relative to the grid system to simulate a dipping sand layer. This is important to minimize diffusive edge effects or numerical diffusion that might be associated with a dipping sand in an Eulerian grid system oriented orthogonal to gravity. Two-dimensional simulations demonstrate the tendency for gas hydrate to accumulate downdip in a sand body because of greater methane transport at depth due to larger sand-clay solubility contrasts. In 3D, basin-scale simulations illuminate how convergent sand layers in a multilayered system can compete for diffusion from clays between them, resulting in relatively low hydrate saturations. All simulations suggest that when hydrate present in clays dissociates with burial, the additional dissolved methane is soaked up by nearby sands preserving high hydrate saturations.

Nole, Michael↗

DOE OSTI2015

A High-Order Nonlinear Diffusion Acceleration for the SN Equations Discretized with the Discontinuous FEM I: Theory and Numerical Results

Recently, nonlinear diffusion acceleration (NDA) methods have received widespread attention in the radiation transport community because of their inherent efficiency (Willert 2014, Wang 2014, Park 2012). However, NDA is in general only conditionally stable \cite{Larsen2014}, while diffusion synthetic acceleration methods (DSA) are unconditionally stable if suitable discretization methods for the transport and diffusion systems are chosen \cite{Larsen1982}. The presented work introduces a new generalized NDA method for solving the transport equation discretized with $S_N$ (discrete ordinates) and discontinuous Finite Elements. In contrast to common coarse-mesh finite difference methods (CMFD), where the diffusion system is discretized using constant shape functions, the presented NDA method discretizes the diffusion system using the same (higher order) finite element basis functions as those for the transport system. This paper focuses on the solution of the steady-state transport equation with the low-order diffusion equation employing the same mesh and the same spatial shape functions, but the NDA methods currently under development for the Rattle$S_N$ake radiation transport solver\cite{Wang2013} based on the MOOSE\cite{Gaston2009} framework is more general: it allows for the diffusion system to be discretized using an arbitrary group structure and arbitrary mesh nested in their $S_N$ system counterparts, and arbitrary combinations of FEM basis functions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Hydraulically opened cone vertical tube diffuser with slanted anti-siphon hole

A diffuser for a water heater fill tubes having a tube wall with an outside diameter. The diffuser includes an elongated flexible diffuser body for positioning at the outlet end of the fill tube. The diffuser body has a sealing end for sealing with an outside surface of the fill tube wall closest to the inlet end of the fill tube, and an open end for positioning at an end closest to the outlet end of the fill tube. The open end has a diameter larger than an outside diameter of the tube wall, thereby creating an diffuser water outlet opening between the diffuser body and the tube wall for redirecting radial water flow emanating from the fill tube water outlet opening(s) toward the diffuser outlet opening. A fill tube assembly for a water heater, a water heater, and a method for heating water are also disclosed.

Rendall, Joseph D.↗