Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,351 records · Page 75

Imaging properties and its improvements of scanning/imaging x-ray microscope

A scanning / imaging X-ray microscope (SIXM) system has been developed at SPring-8. The SIXM consists of a scanning X-ray microscope with a one-dimensional (1D) X-ray focusing device and an imaging (full-field) X-ray microscope with a 1D X-ray objective. The motivation of the SIXM system is to realize a quantitative and highly-sensitive multimodal 3D X-ray tomography by taking advantages of both the scanning X-ray microscope using multi-pixel detector and the imaging X-ray microscope. Data acquisition process of a 2D image is completely different between in the horizontal direction and in the vertical direction; a 1D signal is obtained with the linear-scanning while the other dimensional signal is obtained with the imaging optics. Such condition have caused a serious problem on the imaging properties that the imaging quality in the vertical direction has been much worse than that in the horizontal direction. In this paper, two approaches to solve this problem will be presented. One is introducing a Fourier transform method for phase retrieval from one phase derivative image, and the other to develop and employ a 1D diffuser to produce an asymmetrical coherent illumination.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On origin and destruction of relativistic dust and its implication for ultrahigh energy cosmic rays

Dust grains may be accelerated to relativistic speeds by radiation pressure, diffusive shocks, and other acceleration mechanisms. Such relativistic grains have been suggested as primary particles of ultrahigh energy cosmic rays (UHECRs). In this paper, we first revisit the problem of acceleration by radiation pressure and calculate maximum grain velocities achieved. We find that grains can be accelerated to relativistic speeds with Lorentz factor γ<2 by powerful radiation sources, which is lower than earlier estimates showing that γ could reach ∼10. We then investigate different destruction mechanisms for relativistic grains traversing a variety of environments. In solar radiation, we find that the destruction by thermal sublimation and Coulomb explosions is important. We also quantify grain destruction due to electronic sputtering by ions and grain–grain collisions. Electronic sputtering is found to be rather inefficient, whereas the evaporation following grain–grain collisions is shown to be an important mechanism for which the a ∼ 0.01–1 μm grains would be destroyed after sweeping a gas column N{sub coll}∼5×10{sup 19}–5×10{sup 20} cm{sup −2}. Relativistic dust in the interstellar medium and intergalactic medium (IGM) would be disrupted by Coulomb explosions due to collisional charging after traversing a gas column N{sub coll}∼10{sup 17} cm{sup −2} unless grain material is very strong. We show that photoelectric emission by optical and ultraviolet background radiation is also significant for the destruction of relativistic dust in the IGM. The obtained results indicate that relativistic dust from galaxies would be destroyed before reaching the Earth's atmosphere and unlikely to account for UHECRs.

79 ASTRONOMY AND ASTROPHYSICS↗

Resonance broadened quasi-linear (RBQ) model for fast ion distribution relaxation due to Alfvénic eigenmodes

The burning plasma performance is limited by the confinement of the super-alfvénic fusion products such as alpha particles and the auxiliary heating ions capable of exciting the Alfvénic eigenmodes (AEs) (Gorelenkov et al 2014 Nucl. Fusion 54 125001). Here in this work the effect of AEs on fast ions is formulated within the quasi-linear (QL) theory generalized for this problem recently (Duarte 2017 PhD Thesis University of São Paulo, Brazil). The generalization involves the resonance line broadened interaction of energetic particles (EP) with AEs supplemented by the diffusion coefficients depending on EP position in the velocity space. A new resonance broadened QL code (or RBQ1D) based on this formulation allowing for EP diffusion in radial direction is built and presented in details. In RBQ1D applications we reduce the wave particle interaction (WPI) dynamics to 1D case when the particle kinetic energy is nearly constant. The diffusion equation for EP distribution evolution is then solved simultaneously for all particles along the angular momentum direction. We make initial applications of the RBQ1D to a DIII-D plasma with elevated q-profile where the beam ions show stiff transport properties (Collins et al (The DIII-D Team) 2016 Phys. Rev. Lett. 116 095001). In conclusion, AE driven fast ion profile relaxation is studied for validations of the QL approach in realistic conditions of beam ion driven instabilities in DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multilevel Sequential Monte Carlo Samplers for Normalizing Constants

This article considers the sequential Monte Carlo (SMC) approximation of ratios of normalizing constants associated to posterior distributions which in principle rely on continuum models. Therefore, the Monte Carlo estimation error and the discrete approximation error must be balanced. A multilevel strategy is utilized to substantially reduce the cost to obtain a given error level in the approximation as compared to standard estimators. Two estimators are considered and relative variance bounds are given. The theoretical results are numerically illustrated for two Bayesian inverse problems arising from elliptic partial differential equations (PDEs). The examples involve the inversion of observations of the solution of (i) a 1-dimensional Poisson equation to infer the diffusion coefficient, and (ii) a 2-dimensional Poisson equation to infer the external forcing.

97 MATHEMATICS AND COMPUTING↗

Diffusion bonding of tungsten-vanadium-zirconium using vacuum hot pressing for the development of a low decay heat cladding solution for tungsten spallation targets

Tantalum has been used as cladding material for water-cooled solid tungsten targets at many leading spallation neutron production facilities thanks to its high neutron yield, manageable radiation damage behavior, and excellent corrosion/erosion resistance in radiation environments. However, from a safety hazard perspective, thermal neutron capture of tantalum in spallation environments causes a high specific decay heat in the target volume, which often becomes a limiting factor in increasing the beam power on the target. In this paper, we studied vacuum hot pressing (VHP) parameters to diffusion bond zirconium to tungsten to explore the feasibility of using zirconium alloys as an alternative cladding material to tantalum. Zirconium alloys have long been used as cladding material for early generation solid spallation targets, and nuclear fuel rods. In spallation environments zirconium has significantly lower decay heat with shorter decay time compared to tantalum. The hot isostatic pressing (HIP) of zirconium and tungsten is known to produce limited bonding quality due to the formation of the brittle ZrW2 intermetallic layer. To overcome this problem, placing a vanadium interlayer between tungsten and zirconium has been proposed by exploring parameter space in binary alloy phase diagrams. Under the VHP conditions, 860 ◦C at 70 MPa for 4 hours, Zr-V and V-W showed good diffusion bonding, which demonstrates the feasibility of a single step HIP process to make the zirconium alloy clad tungsten spallation volumes.

Lee, Yong Joong↗

Formulation analysis and computation of an optimization-based local-to-nonlocal coupling method

We present an optimization-based coupling method for local and nonlocal continuum models. Our approach couches the coupling of the models into a control problem where the states are the solutions of the nonlocal and local equations, the objective is to minimize their mismatch on the overlap of the local and nonlocal problem domains, and the virtual controls are the nonlocal volume constraint and the local boundary condition. We present the method in the context of Local-to-Nonlocal diffusion coupling. Numerical examples illustrate the theoretical properties of the approach.

97 MATHEMATICS AND COMPUTING↗

Predictive Large-Eddy Simulation of Supercritical-Pressure Reactive Flows in the Cold Ignition Regime

This report describes a one-year study performed under DOE sponsorship, continuing the investigation of high-pressure turbulent reactive flows. The interest is in the effect of the chemical species distribution in high-pressure turbulent flows in the presence of strong temperature gradients as would occur during reactions in realistic flows. The prime example of such flows are boundary layers in which the wall is at a lower temperature than that of the fluid, as would be the case in Diesel engines. Previous DOE BES work in the program further highlighted the importance of the boundary layer as a configuration for fundamental studies: soot formation in boundary layers is still a problem poorly understood and depends on the availability of particular chemical species at that location, and when boundary layer Large Eddy Simulations results were compared with measurements, the agreement was unfavorable, showing that this important ‘unit’ problem is not well understood. To understand this unit problem, modeling and Direct Numerical Simulations of this unit problem were conducted for the simplest possible multispecies system, that is a binary-species system. The results discovered a new phenomenon, that is, Soret effect induced uphill diffusion. The far-reaching implication is that through an imposed wall/fluid temperature difference it is possible to control the distribution of the species in the boundary layer. Because the results have been documented in one paper published in the refereed literature, and also in conference papers, this final report is in the form of an Executive Summary succinctly describing the results and putting them in perspective with respect to existing information. The refereed and conference papers published are individually listed as Appendices and attached to this report. One manuscript is still in progress and is thus not listed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Space-Time Finite Element Tensor Network Approach for the Time-Dependent Convection–Diffusion–Reaction Equation with Variable Coefficients

In this paper, we present a new space-time Galerkin-like method, where we treat the discretization of spatial and temporal domains simultaneously. This method utilizes a mixed formulation of the tensor-train (TT) and quantized tensor-train (QTT) (please see Section Tensor-Train Decomposition), designed for the finite element discretization (Q1-FEM) of the time-dependent convection–diffusion–reaction (CDR) equation. We reformulate the assembly process of the finite element discretized CDR to enhance its compatibility with tensor operations and introduce a low-rank tensor structure for the finite element operators. Recognizing the banded structure inherent in the finite element framework’s discrete operators, we further exploit the QTT format of the CDR to achieve greater speed and compression. Additionally, we present a comprehensive approach for integrating variable coefficients of CDR into the global discrete operators within the TT/QTT framework. The effectiveness of the proposed method, in terms of memory efficiency and computational complexity, is demonstrated through a series of numerical experiments, including a semi-linear example.

convection–diffusion–reaction equation↗

ELECTRON HEATING IN MAGNETOROTATIONAL INSTABILITY: IMPLICATIONS FOR TURBULENCE STRENGTH IN THE OUTER REGIONS OF PROTOPLANETARY DISKS

The magnetorotational instability (MRI) drives vigorous turbulence in a region of protoplanetary disks where the ionization fraction is sufficiently high. It has recently been shown that the electric field induced by the MRI can heat up electrons and thereby affect the ionization balance in the gas. In particular, in a disk where abundant dust grains are present, the electron heating causes a reduction of the electron abundance, thereby preventing further growth of the MRI. By using the nonlinear Ohm's law that takes into account electron heating, we investigate where in protoplanetary disks this negative feedback between the MRI and ionization chemistry becomes important. We find that the “e-heating zone,” the region where the electron heating limits the saturation of the MRI, extends out up to 80 AU in the minimum-mass solar nebula with abundant submicron-sized grains. This region is considerably larger than the conventional dead zone whose radial extent is ∼20 AU in the same disk model. Scaling arguments show that the MRI turbulence in the e-heating zone should have a significantly lower saturation level. Submicron-sized grains in the e-heating zone are so negatively charged that their collisional growth is unlikely to occur. Our present model neglects ambipolar and Hall diffusion, but our estimate shows that ambipolar diffusion would also affect the MRI in the e-heating zone.

79 ASTRONOMY AND ASTROPHYSICS↗

A new method for in situ electron temperature determinations from plasma wave phenomena

A discrepancy has been reported between the values of the electron temperature t sub e deduced from satellite electrostatic probe measurements and ground based radar backscatter measurements. This discrepancy (radar backscatter temperature less than the probe temperature) is not present when the probe experiment is flown on a rocket to lower altitudes but reappears when the rocket probe attains higher altitudes. There is a need for an independent radio wave method for making in-situ t sub e measurements in order to help resolve this problem. A new method of determining t sub e from the satellite resonant phenomena is presented. It is based on the splitting (observed at high latitudes) of the diffuse resonance which occurs at the frequency f sub d1 between f sub h and 2 f sub h. The advantage of this method over the other methods involving ionospheric resonances is the simplicity of the required calculations. There is, however, the limitation of plasma conditions where the f sub d1 resonance can be observed and the limitation to mid-to-high latitudes where the splitting is observed.

Oya, H.↗

Fiber distributed feedback laser

Utilizing round optical fibers as communication channels in optical communication networks presents the problem of obtaining a high efficiency coupling between the optical fiber and the laser. A laser is made an integral part of the optical fiber channel by either diffusing active material into the optical fiber or surrounding the optical fiber with the active material. Oscillation within the active medium to produce lasing action is established by grating the optical fiber so that distributed feedback occurs.

Elachi, C.↗

Inexpensive cross-linked polymeric separators made from water soluble polymers

Polyvinyl alcohol (PVA) crosslinked chemically with aldehyde reagents produces membranes which demonstrate oxidation resistance, dimensional stability, low ionic resistivity, low zincate diffusivity, and low zinc dendrite penetration rate which make them suitable for use as alkaline battery separators. They are intrinsically low in cost and environmental health and safety problems associated with commercial production appear minimal. Preparation, property measurements, and cell test results in Ni/Zn and Ag/Zn cells are described and discussed.

Hsu, L. C.↗

Thermal response of large area high temperature superconducting YBaCuO infrared bolometer

Thermal analysis of large area high temperature superconducting infrared detector operating in the equilibrium mode (bolometer) was performed. An expression for the temperature coefficient beta=1/R(dR/dT) in terms of the thermal conductance and the thermal time constant of the detector were derived. A superconducting transition edge bolometer is a thermistor consisting of a thin film superconducting YBaCuO evaporated into a suitable thermally isolated substrate. The operating temperature of the bolometer is maintained close to the midpoint of the superconducting transition region where the resistance R has a maximum dynamic range. A detector with a strip configuration was analyzed and an expression for the temperature rise (delta T) above the ambient due to a uniform illumination with a source of power density P(sub i) was calculated. An expression for the thermal responsivity of the detector was derived using the thermal diffusion analysis with appropriate boundary conditions. It was found that the thermal responsibility depends upon the spatial modulation frequency and the angular frequency of the incoming radiation. The problem of the thermal cross talk between different detector elements was addressed. In the case of monolithic HTS detector array with a row of square elements of dimensions 2a and CCD or CID readout electronics the thermal spread function was derived for different spacing between elements. This analysis can be critical for future design and applications of large area focal plane arrays as broad band optical detectors made of granular thin films HTS YBaCuO.

Khalil, Ali E.↗

Catalytic ignition model in a monolithic reactor with in-depth reaction

Two transient models have been developed to study the catalytic ignition in a monolithic catalytic reactor. The special feature in these models is the inclusion of thermal and species structures in the porous catalytic layer. There are many time scales involved in the catalytic ignition problem, and these two models are developed with different time scales. In the full transient model, the equations are non-dimensionalized by the shortest time scale (mass diffusion across the catalytic layer). It is therefore accurate but is computationally costly. In the energy-integral model, only the slowest process (solid heat-up) is taken as nonsteady. It is approximate but computationally efficient. In the computations performed, the catalyst is platinum and the reactants are rich mixtures of hydrogen and oxygen. One-step global chemical reaction rates are used for both gas-phase homogeneous reaction and catalytic heterogeneous reaction. The computed results reveal the transient ignition processes in detail, including the structure variation with time in the reactive catalytic layer. An ignition map using reactor length and catalyst loading is constructed. The comparison of computed results between the two transient models verifies the applicability of the energy-integral model when the time is greater than the second largest time scale of the system. It also suggests that a proper combined use of the two models can catch all the transient phenomena while minimizing the computational cost.

Tien, Ta-Ching↗

Hot HB Stars in Globular Clusters: Physical Parameters and Consequences for Theory: Radiative Levitation Versus Helium Mixing - 5

Atmospheric parameters (T(sub eff), log g), masses and helium abundances are derived for 42 hot horizontal branch (HB) stars in the globular cluster NGC6752. For 19 stars we derive magnesium and iron abundances as well and find that iron is enriched by a factor of 50 on average with respect to the cluster abundance whereas the magnesium abundances are consistent with the cluster abundance. Radiation pressure may levitate heavy elements like iron to the surface of the star in a diffusive process. Taking into account the enrichment of heavy elements in our spectroscopic analyses we find that high iron abundances can explain part, but not all, of the problem of anomalously low gravities along the blue HB. The blue HB stars cooler than about 15,100 K and the sdB stars (T(sub eff) greater than or = 20,000 K) agree well with canonical theory when analysed with metal-rich ([M/H] = +0.5) model atmospheres, but the stars in between these two groups remain offset towards lower gravities and masses. Deep Mixing in the red giant progenitor phase is discussed as another mechanism that may influence the position of the blue HB stars in the (T(sub eff), log g)-plane but not their masses.

Moehler, S.↗

Structural Concepts Study of Non-circular Fuselage Configurations

A preliminary study of structural concepts for noncircular fuselage configurations is presented. For an unconventional flying-wing type aircraft, in which the fuselage is inside the wing, multiple fuselage bays with non-circular sections need to be considered. In a conventional circular fuselage section, internal pressure is carried efficiently by a thin skin via hoop tension. If the section is non-circular, internal pressure loads also induce large bending stresses. The structure must also withstand additional bending and compression loads from aerodynamic and gravitational forces. Flat and vaulted shell structural configurations for such an unconventional, non-circular pressurized fuselage of a large flying-wing were studied. A deep honeycomb sandwich-shell and a ribbed double-wall shell construction were considered. Combinations of these structural concepts were analyzed using both analytical and simple finite element models of isolated sections for a comparative conceptual study. Weight, stress, and deflection results were compared to identify a suitable configuration for detailed analyses. The flat sandwich-shell concept was found preferable to the vaulted shell concept due to its superior buckling stiffness. Vaulted double-skin ribbed shell configurations were found to be superior due to their weight savings, load diffusion, and fail-safe features. The vaulted double-skin ribbed shell structure concept was also analyzed for an integrated wing-fuselage finite element model. Additional problem areas such as wing-fuselage junction and pressure-bearing spar were identified.

Mukhopadhyay, Vivel↗

Position-Dependency of Fuel Assembly Homogenization in a Pressurized Water Reactor

2-step nuclear reactor core analysis has long been established as the conventional calculational procedure for the pressurized water reactor (PWR) due to its efficiency. This is because the procedure significantly avoids the potentially huge computational burden arising from an explicit geometrical modeling and consideration of over ∼ 500 energy groups. It instead approaches the problem with spatial homogenization and group condensation: a 2-D lattice transport calculation to homogenize the fuel assembly cross-sections and a 3-D whole core nodal diffusion. One notes that the few-group constant generation using deterministic method could be complicated due to its limited flexibilities in geometrical modeling and inherent drawbacks of the multi-group approximation. Furthermore, accuracy of the deterministic few-group constant generation can especially be compromised during the resonance treatment calculation. As such, an alternative lattice calculation using the Monte Carlo (MC) approach, in place of the deterministic fuel assembly homogenization (few-group constant generation), has previously been proposed. During the assembly homogenization, directions and magnitudes of the current at boundaries are unknowns. Therefore, zero net current boundary conditions (infinite lattice) are instead used as the assembly node. The assumption is quite reasonable considering heterogeneities inside the node are more dominant than heterogeneities due to nodal location in the core. In fact, this approximation has successfully been applied to many LWR core analyses to reasonable accuracy. It should nonetheless be noted that the net current zero boundary assumption effectively ignores position-dependency of the fuel assembly homogenization in the general 2-step method. As such, there are a number of well-known researches on the extended assembly homogenization techniques, which introduce correction factors to compensate effects of the assembly environment, such as Generalized Equivalence Theory (GET), Super-homogenization method (SPH), and color set calculation method. Even though the assembly environment effect has long been acknowledged, its effect has not properly been investigated due to the limited computing power in the past and due to the reasonable accuracy of the extended assembly homogenization techniques. Nonetheless, quantification of position-dependency of the assembly homogenization in the core could still potentially offer additional value in terms of calculation accuracy to the 2-step reactor analysis. It is upon this motivation that the paper is prepared. In this research, position-dependent few-group constants were calculated and compared against those of infinite lattices. All calculations were completed using the MC Serpent2 code with ENDF/B7.0 nuclear data library. In this study, position-dependency of the fuel assembly homogenized cross sections in a PWR core has been investigated through comparison of the conventional 2-step reactor analysis with the reference high fidelity Monte Carlo simulation. It was found that the homogenized two-group constants are clearly position-dependent as the discrepancies between infinite lattice and whole-core calculations are quite significant especially in a rodded assembly. As such, the corresponding correction required due to the assembly environment effect is not negligible compared to the correction required due to assembly heterogeneity itself. It is therefore quite clear that position-dependency of the group cross sections should be considered in the 2-step procedure for a more accurate PWR core analysis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

54 ENVIRONMENTAL SCIENCES↗