SEARCH · Search NASA

Results for “Numerical diffusion”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Wide bandwidth frequency-domain thermoreflectance: Volumetric heat capacity, anisotropic thermal conductivity, and thickness measurements

A frequency-domain thermoreflectance (FDTR) system with a frequency range of 1 Hz to 75 MHz is presented. The wide bandwidth of pump modulation frequency enables enhanced simultaneous thermophysical property characterization of bulk and thin-film materials. The wide bandwidth FDTR system is demonstrated with simultaneous measurements of in-plane and cross-plane thermal conductivities of sapphire and muscovite mica, thickness and thermal conductivity of gold/titanium thin films, and isotropic thermal conductivity and volumetric heat capacity of lithium niobate and silicon. Thermophysical properties measured with FDTR are compared to literature values or independent measurements for verification. Finally, at low frequencies, a negligible error was achieved in the numerical integration of the heat diffusion equation with a Gauss–Legendre quadrature method solved with 400 integration points and an upper integration limit of 16weff, where weff is the effective pump–probe 1/e2 radius.

Ziade, Elbara (ORCID:0000000249020497)↗

DOE OSTI2021

Evidence for a thermally driven charge-density-wave transition in 1T-TaS 2 thin-film devices: Prospects for GHz switching speed

We report on the room-temperature switching of 1T-TaS 2 thin-film charge-density-wave devices, using nanosecond-duration electrical pulsing to construct their time-resolved current–voltage characteristics. The switching action is based upon the nearly commensurate to incommensurate charge-density-wave phase transition in this material, which has a characteristic temperature of 350 K at thermal equilibrium. For sufficiently short pulses, with rise times in the nanosecond range, self-heating of the devices is suppressed, and their current–voltage characteristics are weakly nonlinear and free of hysteresis. This changes as the pulse duration is increased to ~200 ns, where the current develops pronounced hysteresis that evolves nonmonotonically with the pulse duration. By combining the results of our experiments with a numerical analysis of transient heat diffusion in these devices, we clearly reveal the thermal origins of their switching. In spite of this thermal character, our modeling suggests that suitable reduction of the size of these devices should allow their operation at GHz frequencies.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Third-order photon correlations extract single-nanocrystal multiexciton properties in solution

Colloidal semiconductor nanocrystals are considered promising materials for high-flux optical applications, including lasing, light-emitting diodes, biological imaging, and quantum optics. In high-flux applications, multiexcitons can significantly contribute to emission, influencing its brightness, spectral purity, and kinetics. As a result, understanding and controlling multiexciton emission in colloidal nanocrystal materials is of the utmost importance. In the past, single-nanocrystal photon correlation methods have been applied to understand biexciton and triexciton efficiencies, lifetimes, and spectra. While powerful, such methods suffer from user selection bias and require stable emission from single nanocrystals. To compensate for this shortcoming, second-order correlation methods were developed to extract sample-averaged biexciton properties from a solution of nanocrystals. Until now, however, the analogous third-order solution photon correlation methods remained unexplored. In this work, we present a pair of third-order photon correlation techniques to obtain the sample-averaged single-nanocrystal triexciton quantum yield and lifetime in a solution-phase experiment. These techniques derive from the relationship between the Poisson probability of nanocrystal photon absorption and the intrinsic probability of nanocrystal photon emission. We validate the theoretical background of these techniques by creating a numerical model to simulate the diffusion and emission of many nanocrystals in solution. Our simulations confirm that the average triexciton quantum yield and triexciton lifetime can be extracted from a solution of nanocrystals. These techniques will enable researchers to gain a better understanding of the fundamental multiexciton properties of colloidal nanocrystals.

Horowitz, Jonah R. [Massachusetts Institute of Tec↗

NASA NTRS1974

Experimental determination of turbulence in a GH2-GOX rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a gaseous rocket combustion chamber have been determined from the experimental measurements of the tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and Spalding's numerical method for solving the conservation equations of fluid mechanics was used to calculate these quantities. Taylor's theory was extended to consider the inhomogeneity of the turbulence field in the axial direction of the combustion chamber. An exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the maximum value of the intensity of turbulence is about 15% and the Lagrangian correlation coefficient drops to about 0.12 in one inch of the chamber length.

Tou, P.↗

NASA NTRS1978

Measurement of the turbulence in a gaseous rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient in a gaseous hydrogen (GH2)-gaseous oxygen (GOX) rocket combustion chamber are studied. A helium tracer gas was injected at an upstream point on the combustion chamber centerline. Samples collected downstream were analyzed for the spreading of the tracer gas concentration. Using Taylor's turbulent diffusion theory and Gosman's numerical method, turbulence parameters were found. In order to evaluate turbulence field inhomogeneity in the axial direction, Taylor's theory was extended. An exponential function was used to approximate the Lagrangian correlation coefficient. It was found that turbulence intensity varies from 15% near the injector to 4% at the nozzle entrance, and that the Lagrangian coefficient rapidly decreases when the dispersion time difference increases.

Tou, P. P.↗

NASA NTRS1979

Interfacial geometry and D-variation effects in two-phase systems

Numerical solutions of the governing diffusion equation for two-phase concentration dependent diffusion coefficients are examined. Solutions were also calculated for planar, cylindrical, and spherical geometries to compare the effect of interface geometries with those caused by concentration-dependent diffusion coefficients, and two methods of averaging D were considered to determine the best averaging method for different types of D-variations. The effects of interface-location criteria on mass conservation and convergence of interface location, diffusion coefficient variation in the alpha and beta-phases of a two-phase binary alloy system, effect of D(alpha) variation in a cylindrical couple on beta-phase thickness, and geometry and D-variation effects on the degree of homogenization were determined. It is concluded that typical D(alpha)-variations can have a greater influence on the kinetics of interdiffusion than the geometry.

Tenney, D. R.↗

NASA NTRS1981

Visible light scatter measurements of the Advanced X-ray Astronomical Facility /AXAF/ mirror samples

NASA is studying the properties of mirror surfaces for X-ray telescopes, the data of which will be used to develop the telescope system for the Advanced X-ray Astronomical Facility. Visible light scatter measurements, using a computer controlled scanner, are made of various mirror samples to determine surface roughness. Total diffuse scatter is calculated using numerical integration techniques and used to estimate the rms surface roughness. The data measurements are then compared with X-ray scatter measurements of the same samples. A summary of the data generated is presented, along with graphs showing changes in scatter on samples before and after cleaning. Results show that very smooth surfaces can be polished on the common substrate materials (from 2 to 10 Angstroms), and nickel appears to give the lowest visible light scatter.

Griner, D. B.↗

NASA NTRS1985

The expected cosmic ray density and stream distributions at the heliolatitudinal asymmetry of solar wind

The results of the spatial distribution of cosmic ray density, gradients, and anisotropy obtained on the basis of the numerical solution of the anisotropic diffusion equation with an account of solar wind velocity change depending on the latitudinal angle theta of the form U=u sub oe sup alpha theta and the diffusion coefficient depending on the spatial coordinates and the particle rigidity are presented. It is shown that the increase of the solar wind velocity and the diffusion coefficient with heliolatitude leads to gradient distributions that are in accord with experimental data observed in space. The results of the energetic spectrum of 11 and 22-year cosmic ray variations obtained with an account of direction of the general magnetic field of the Sun are presented are given.

Alania, M. V.↗

NASA NTRS1985

The effect of the neutral sheet structure of the interplanetary magnetic field on cosmic ray distribution in space

Results of the numerical solution of the anistoropic diffusion equation are presented. The modulation depth of galactic cosmic rays is defined by the degree of curvature of the neutral current sheet in the heliosphere. The effect of the regular interplanetary magnetic field (IMF) on cosmic ray anisotropy in the period of solar activity minimum (in 1976) is analyzed by the data of the neutron super-monitors of the world network, and the heliolatitudinal gradient and cosmic ray diffusion coefficient are defined.

Alania, M. V.↗

NASA NTRS1984

Effect of growth parameters on compositional variations in directionally solidified HgCdTe alloys

A series of Hg(1-x)Cd(x)Te alloy crystals was grown by directional solidification with compositions ranging from x = 0.2 to x = 0.4. The measured axial compositional profiles were interpreted in terms of an exact numerical solution of the appropriate diffusion equation that takes into account both the variations of the segregation coefficient and solidification rate with composition. The solutions for all growth rates agree generally well with the experimental data, and confirm earlier observations that employed approximate analytical solutions. The effective mass diffusion coefficients showed no significant composition or melt temperature dependence. Slightly higher diffusion coefficients were obtained, however, for the highest growth rates.

Szofran, F. R.↗

NASA NTRS1990

Dynamical effects of annihilation in pair-dominated winds

The steady, spherically symmetric flow of an ideal fluid dominated by photons and ultrarelativistic electron-positron pairs is analyzed. A new wind equation and a set of critical point conditions are obtained which describe the relativistic flow of an annihilation gas in which the flow velocity exceeds the diffusion velocity of the photons. Numerical results are reported which suggest the possible existence of trapped, pure-pair winds driven by the combined pressure of the pairs and the photons. Most of the annihilation occurs below the critical radius in trapped flows, and a substantial fraction of the total energy of the injected pairs is converted into kinetic energy and radiation. Accurate numerical solutions for the flow velocity and the positron loss rate in optically thin, Newtonian winds are obtained, and a useful approximate analytic relation between the positron loss rate and the flow velocity is derived which suggests that a large number of pairs may survive the annihilation region, ultimately escaping the potential well.

Becker, Peter A.↗

NASA NTRS1990

Effects of radiation of InP cells epitaxially grown on Si and GaAs substrates

The properties of heteroepitaxial InP cells were determined both before and after 10-MeV proton irradiations. Numerical values, obtained for the diffusion and recombination components of the reverse saturation currents, were found to be consistent with the distribution of dislocations. The radiation resistance of the heteroepitaxial cells was significantly greater than that observed for n/p homoepitaxial InP cells. The carrier removal rate, obtained by C-V measurements, was 1800/cm for 10-MeV protons compared with 2.2/cm for 1-MeV electrons. The high carrier removal rate was found to have no significant effect on the cell's series resistance. It was concluded that the heteroepitaxial cell performance is dominated by the high dislocation density attributable to lattice constant mismatch. Although the efficiencies of the present cells are low, the recent achievement of 13.7 percent AM0 efficiencies using a GaAs substrate demonstrates the marked improvement that can be attained using more appropriate transition layers.

Weinberg, I.↗

NASA NTRS1992

Upper bounds for convergence rates of vector extrapolation methods on linear systems with initial iterations

The application of the minimal polynomial extrapolation (MPE) and the reduced rank extrapolation (RRE) to a vector sequence obtained by the linear iterative technique x(sub j) + 1 = Ax(sub j) = b,j = 1,2,..., is considered. Both methods produce a two dimensional array of approximations s(sub n,k) to the solution of the system (I - A)x = b. Here, s(sub n,k) is obtained from the vectors x(sub j), n is less than or equal to j is less than or equal to n + k + 1. It was observed in an earlier publication by the first author that the sequence s(sub n,k), k = 1,2,..., for n greater than 0, but fixed, possesses better convergence properties than the sequence s(sub 0,k), k = 1,2,.... A detailed theoretical explanation for this phenomenon is provided in the present work. This explanation is heavily based on approximations by incomplete polynomials. It is demonstrated by numerical examples when the matrix A is sparse that cycling with s(sub n,k) for n greater than 0, but fixed, produces better convergence rates and costs less computationally than cycling with s(sub 0,k). It is also illustrated numerically with a convection-diffusion problem that the former may produce excellent results where the latter may fail completely. As has been shown in an earlier publication, the results produced by s(sub 0,k) are identical to the corresponding results obtained by applying the Arnoldi method or generalized minimal residual scheme (GMRES) to the system (I - A)x = b.

Sidi, Avram↗

NASA NTRS1991

Characterization of directionally solidified Hg(1-x)Zn(x)Se semiconducting alloys

Hg(1-x)Zn(x)Se alloys with composition between x = 0.08 and 0.115 were synthesized from elemental constituents and were resolidified using a Bridgman-Stockbarger growth technique. By performing precision mass-density measurements on selected wafers cut perpendicular to the growth axis, it was shown that the axial compositional profiles fit a numerical solution to the 1D diffusion equation which takes into account the variation of interface velocity with time. Infrared transmission-edge measurements performed on selected transverse slices from each ingot showed that the relative radial variations in composition decreased with decreasing growth rate. Van der Pauw measurements on selected wafers showed that the 10 exp 18/cu cm electron concentration, typical of as-grown crystals, could be reduced by approximately an order of magnitude by annealing in Se vapor.

Cobb, S. D.↗

NASA NTRS1994

On the dynamics of some grid adaption schemes

The dynamics of a one-parameter family of mesh equidistribution schemes coupled with finite difference discretisations of linear and nonlinear convection-diffusion model equations is studied numerically. It is shown that, when time marched to steady state, the grid adaption not only influences the stability and convergence rate of the overall scheme, but can also introduce spurious dynamics to the numerical solution procedure.

Sweby, Peter K.↗

NASA NTRS2002

Numerical Simulations of Vortex Generator Vanes and Jets on a Flat Plate

Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Allan, Brian G.↗