SEARCH · Engineering Papers
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.
The 'flooded agglomerate' model of the Teflon-bonded gas diffusion electrode is discussed. A mathematical treatment of the 'flooded agglomerate' model is given; it can be used to predict the performance of the electrode as a function of measurable physical parameters.
High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.
The ramp-reversal memory (RRM) effect in metal–insulator transition metal oxides (TMOs), a non-volatile resistance change induced by repeated temperature cycling, has attracted considerable interest in neuromorphic computing and non-volatile memory devices. Our previous defect motion model successfully explained RRM in vanadium dioxide (VO 2 ), capturing observed critical temperature shifts and memory accumulation throughout the sample. However, this approach lacked interactions between metallic and insulating domains. Here, we extend our model by combining a correlated Random Field Ising Model with defect diffusion-segregation, enabling accurate hysteresis modeling while predicting the relationship between RRM and domain interactions. Our simulations demonstrate that the maximum RRM occurs when the turnaround temperature approaches the inflection point. This peak in RRM vs. turnaround temperature is consistent with prior transport measurements, as well as our own optical measurements reported here. Significantly, we find that increasing nearest-neighbor interactions enhances the maximum memory effect, thus providing a clear mechanism for optimizing RRM performance. Since our model employs minimal assumptions, we predict that RRM should be a widespread phenomenon in materials exhibiting patterned phase coexistence of electronic domains. This work not only advances fundamental understanding of memory behavior in TMOs but also establishes a much-needed theoretical framework for optimizing device applications.
The diffusion and dissociation mechanisms governing the evolution of oxygen and vacancy-oxygen defects in Czochralski-grown Si samples have been studied. The samples were irradiated at (i) room temperature or (ii) elevated temperature (350 °C) by MeV electrons and then isothermally annealed at 8 different temperatures in the range of 300 °C to 500 °C. The evolution of the concentrations of oxygen complexes (O{sub n}, n ≤ 3) and mono-vacancy-oxygen defects (VO{sub n}, n ≤ 4) have been followed by infrared absorption measurements of local vibrational modes originating from the individual defects. The experimental kinetics data have been compared with simulation results based on the theory for diffusion limited reactions, assuming a model where sequential build-up of the VO{sub n} defects is a key ingredient. A close quantitative agreement is obtained for both sets of samples despite quite different initial conditions prior to the annealing, which adds evidence to the validity of the model. Values for the diffusivity and dissociation rates of VO{sub n} (n ≤ 4) and O{sub n} (n ≤ 3) have been deduced and in general, the mobility and stability of VO{sub n} decrease and increase with n, respectively. For all the defects, partial dissociation appears as a prevailing process during diffusion, while full dissociation of VO{sub n} is limited by an energy barrier identical to that of interstitial oxygen (O{sub i}) diffusion (∼2.55 eV). The oxygen dimer and trimer are fast diffusers but slower than substitutional oxygen, i.e., VO; VO is found to be the most mobile species, whilst O{sub i} is the slowest one with a difference in diffusivity of up to 7 orders of magnitude in the studied temperature range.
The properties of water under confinement are of practical and fundamental interest. Here in this work we study the properties of water in the self-assembled lyotropic phases of gemini surfactants with a focus on testing the standard analysis of quasi-elastic neutron scattering (QENS) experiments. In QENS experiments the dynamic structure factor is measured and fit to models to extract the translational diffusion constant, D T , and rotational relaxation time, τ R . We test this procedure by using simulation results for the dynamic structure factor, extracting the dynamic parameters from the fit as is typically done in experiments, and comparing the values to those directly measured in the simulations. We find that the decoupling approximation, where the intermediate scattering function is assumed to be a product of translational and rotational contributions, is quite accurate. The jump-diffusion and isotropic rotation models, however, are not accurate when the degree of confinement is high. In particular, the exponential approximations for the intermediate scattering function fail for highly confined water and the values of D T and τ R can differ from the measured value by as much as a factor of two. Other models have more fit parameters, however, and with the range of energies and wave-vectors accessible to QENS, the typical analysis appears to be the best choice. In the most confined lamellar phase, the dynamics are sufficiently slow that QENS does not access a large enough time scale and neutron spin echo measurements would be a valuable technique in addition to QENS.
A mathematical model is presented for the description of axisymmetric laminar-jet diffusion flames. The analysis includes the effects of inertia, viscosity, diffusion, gravity and combustion. These mechanisms are coupled in a boundary layer type formulation and solutions are obtained by an explicit finite difference technique. A dimensional analysis shows that the maximum flame width radius, velocity and thermodynamic state characterize the flame structure. Comparisons with experimental data showed excellent agreement for normal gravity flames and fair agreement for steady state low Reynolds number zero gravity flames. Kinetics effects and radiation are shown to be the primary mechanisms responsible for this discrepancy. Additional factors are discussed including elipticity and transient effects.
Complex heat treatment operations and advanced manufacturing processes such as laser or electron-beam welding will see the metallic workpiece experience a considerable range of temperatures and heating/cooling rates. These intrinsic conditions will have a significant bearing upon the microstructure of the material, and in turn upon the thermo-mechanical properties. In this work, a diffusion-based approach to model the growth and shrinkage of precipitates in the alpha + beta field of the common titanium alloy Ti-6Al-4V is established. Further, the numerical model is extended using a JMA-type approach to explore the dependency of the beta-transus temperature on extremely high heating rates, whereby dissolution alone is too slow to accurately describe the alpha to beta-phase transformation. Experimental heat treatments at rates of 5, 50, and 500 °C/s were performed, and metallographic analysis of the samples was used to validate the numerical modeling framework predictions for lamellar shrinkage, while data from the literature has been used to evaluate the numerical modeling framework for the growth of equiaxed microstructures. The agreement between measurements and numerical predictions was found to be good.
Two-dimensional Euler and Navier-Stokes solutions of the flow through three inlet/diffuser configurations with terminal shock systems are reported. Calculations without bleed indicate that the terminal shock location is very sensitive to the outflow back pressure. For cases where there are no available experimental results, it becomes difficult to estimate the back pressure that will result in a terminal shock. Estimates based on quasi-one-dimensional analysis are not found adequate for complex two-dimensional flows. It is found that since the flow downstream of the terminal shock is subsonic, and what happens at the outflow boundary affects the flow inside the inlet, enough of the subsonic diffuser must be modeled to accurately predict the terminal shock region. The diffuser portion should be fairly long with the outflow boundary occurring in a region of more or less uniform flow to be able to prescribe a uniform back pressure. The third configuration studied was investigated with and without incorporating bleed in the code. It is found that the use of bleed stabilizes the shock location and allows solutions which without bleed result in unstarting of the inlet. Comparisons are made with available experimental data.
The influence of breaking gravity waves on the dynamics and chemical composition of the 60 to 110 km region is investigated with a two dimensional model that includes a parameterization of gravity wave momentum deposition and diffusion. The dynamical model is described by Garcia and Solomon (1983) and Solomon and Garcia (1983) and includes a complete chemical scheme for the mesosphere and lower thermosphere. The parameterization of Lindzen (1981) is used to calculate the momentum deposited and the turbulent diffusion produced by the gravity waves. It is found that wave momentum deposition drives a very vigorous mean meridional circulation, produces a very cold summer mesopause and reverse the zonal wind jets above about 85 km. The seasonal variation of the turbulent diffusion coefficient is consistent with the behavior of mesospheric turbulences inferred from MST radar echoes. The large degree of consistency between model results and various types of dynamical and chemical data supports very strongly the hypothesis that breaking gravity waves play a major role in determining the zonally-averaged dynamical and chemical structure of the 60 to 110 km region of the atmosphere.
Fire and turbulent flow flame propagation using gas jets - model study
A simple theoretical model which can be used as an aid in the design of the baffle aperture region of a hollow cathode equipped ion thruster was developed. An analysis of the ion and electron currents in both the main and cathode discharge chambers is presented. From this analysis a model of current flow through the aperture, which is required as an input to the design model, was developed. This model was verified experimentally. The dominant force driving electrons through the aperture was the force due to the electrical potential gradient. The diffusion process was modeled according to the Bolm diffusion theory. A number of simplifications were made to limit the amount of detailed plasma information required as input to the model to facilitate the use of the model in thruster design. This simplified model gave remarkably consistant results with experimental results obtained with a given thruster geometry over substantial changes in operating conditions. The model was uncertain to about a factor of two for different thruster cathode region geometries. The design usefulness was limited by this factor of two uncertainty and by the accuracy to which the plasma parameters required as inputs to the model were specified.
A theoretical model for the trap-assisted tunneling process in diffused n-on-p and implanted n(+)-on-p HgCdTe photodiodes is presented. The model describes the traps and the trap characteristics: concentration, energy level, and capture cross sections. We have observed that the above two types of diodes differ in the voltage dependence of the trap-assisted tunneling current and dynamic resistance. Our model takes this difference into account and offers an explanation of the phenomenon. The good fit between measured and calculated DC characteristics of the photodiodes (for medium and high reverse bias and for temperatures from 65 to 140 K) supports the validity of the model.
Si self-diffusion in the presence of end-of-range (EOR) defects is investigated using {sup nat}Si/{sup 28}Si isotope multilayers. The isotope multilayers were amorphized by Ge ion implantation, and then annealed at 800–950 °C. The behavior of Si self-interstitials is investigated through the {sup 30}Si self-diffusion. The experimental {sup 30}Si profiles show further enhancement of Si self-diffusion at the EOR defect region, in addition to the transient enhanced diffusion via excess Si self-interstitials by EOR defects. To explain this additional enhanced diffusion, we propose a model which takes into account enhanced diffusion by tensile strain originated from EOR defects. The calculation results based on this model have well reproduced the experimental {sup 30}Si profiles.
Pore pressure diffusion along faults influences induced seismicity and rupture mechanics. In situ hydraulic diffusivity measurements along faults are rare and generally lower than inferred from seismicity migration. Here we use the tidal response of deep geothermal boreholes to measure fault diffusivity and permeability. Initial interpretations of the observation with a homogeneous confined aquifer model result in diffusivities of 10 -3 –10 -1 m 2 /s. However, this model mixes signals from both the conduit and the host rock. We develop a model for tidal response with a fault passing through the aquifer based on the fault-guided fracture network and solve for hydraulic properties in both the fault and the host rock. The resulting fault permeability is 2 × 10 -14 –7 × 10 -14 m 2 (90% confidence interval[CI]) and fault diffusivity is 0.08–0.33 m 2 /s (90% CI), which is 2 orders of magnitude higher than the host rock diffusivity in some wells, thus highlighting the role of faults as fluid conduits.
Not Available
The need to understand potential climate impacts and feedbacks in Arctic regions has prompted recent interest in modeling of permafrost dynamics in a warming climate. A new fine-scale integrated surface/subsurface thermal hydrology modeling capability is described and demonstrated in proof-of-concept simulations. The new modeling capability combines a surface energy balance model with recently developed three-dimensional subsurface thermal hydrology models and new models for nonisothermal surface water flows and snow distribution in the microtopography. Surface water flows are modeled using the diffusion wave equation extended to include energy transport and phase change of ponded water. Variation of snow depth in the microtopography, physically the result of wind scour, is also modeled heuristically with a diffusion wave equation. The multiple surface and subsurface processes are implemented by leveraging highly parallel community software. Fully integrated thermal hydrology simulations on the tilted open book catchment, an important test case for integrated surface/subsurface flow modeling, are presented. Fine-scale 100-year projections of the integrated permafrost thermal hydrological system on an ice wedge polygon at Barrow Alaska in a warming climate are also presented. Finally, these simulations demonstrate the feasibility of microtopography-resolving, process-rich simulations as a tool to help understand possible future evolution of the carbon-rich Arctic tundra in a warming climate.
The MODIS instrument relies on solar calibration to achieve the required radiometric accuracy. This solar calibration occurs as the TERRA spacecraft comes up over the North Pole. The earth underneath the spacecraft is still dark for approximately one minute and the sun is just rising over the earth's north polar regions. During this time the sun moves through about 4 degrees, the scan mirror rotates about 19 times and about 50 frames (exposures) are made of the white solar diffuser. For some of MODIS'S bands the brightness of the screen is reduced, to prevent detector saturation, by means of a pinhole screen, which produces approximately 600 pinhole images of the sun within the field of view of any one detector. Previous attempts at creating a detailed radiometric model of this calibration scenario produced intensity variations on the focal planes with insufficient detail to be useful. The current computational approach produces results, which take into account the motion of the sun and the scan mirror and produces variations, which strongly resemble the observed focal plane intensity variations. The computational approach and results and a comparison with actual observational data are presented.
Not provided.