Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffused interface model”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗
NASA NTRS2002

C/sic Life Prediction for Propulsion Applications

Accurate life prediction is critical to successful use of ceramic matrix composites (CMC). The tools to accomplish this are immature and not oriented toward the behavior of carbon fiber reinforced silicon carbide (C/SiC), the primary system of interest for many reusable and single mission launch vehicle propulsion and airframe applications. This paper describes an approach and process made to satisfy the need to develop an integrated life prediction system that addresses mechanical durability and environmental degradation of C/SiC. Issues such as oxidation, steam and hydrogen effects on material behavior are discussed. Preliminary tests indicate that steam will aggressively remove SiC seal coat and matrix in line with past experience. The kinetics of water vapor reaction with carbon fibers is negligible at 600 C, but comparable to air attack at 1200 C. The mitigating effect of steam observed in fiber oxidation studies has also been observed in stress rupture tests. Detailed microscopy of oxidized specimens is being carried out to develop the oxidation model. Carbon oxidation kinetics are reaction controlled at intermediate temperatures and diffusion controlled at high temperatures (approximately 1000 C). Activation energies for T-300 and interface pyrolytic carbon were determined as key inputs to the oxidation model. Crack opening as a function of temperature and stress was calculated. Mechanical property tests to develop and verify the probabilistic life model are very encouraging except for residual strength prediction. Gage width is a key variable governing edge oxidation of seal coated specimens. Future efforts will include architectural effects, enhanced coatings, biaxial tests, and LCF. Modeling will need to account for combined effects.

Levine, Stanley R.↗
NASA NTRS2002

C/SIC Life Prediction for Propulsion Applications

Accurate life prediction is critical to successful use of ceramic matrix composites (CMC). The tools to accomplish this are immature and not oriented toward the behavior of carbon fiber reinforced silicon carbide (C/SiC), the primary system of interest for many reusable and single mission launch vehicle propulsion and airframe applications. This paper describes an approach and progress made to satisfy the need to develop an integrated life prediction system that addresses mechanical durability and environmental degradation of C/SiC. Issues such as oxidation, steam and hydrogen effects on material behavior are discussed. Preliminary tests indicate that stream will aggressively remove SiC seal coat and matrix in line with past experience. The kinetics of water vapor reaction with carbon fibers is negligible at 600 C, but comparable to air attack at 1200 C. The mitigating effect of steam observed in fiber oxidation studies has also been observed in stress rupture tests. Detailed microscopy of oxidized specimens is being carried out to develop the oxidation model. Carbon oxidation kinetics are reaction controlled at intermediate temperatures and diffusion controlled at high temperatures (approx. 1000 C). Activation energies for T-300 and interface pyrolytic carbon were determined as key inputs to the oxidation model. Crack opening as a function of temperature and stress was calculated. Mechanical property tests to develop and verify the probabilistic life model are very encouraging except for residual strength prediction. Gage width is a key variable governing edge oxidation of seal coated specimens. Future efforts will include architectural effects, enhanced coatings, biaxial tests, and LCF. Modeling will need to account for combined effects.

Levine, Stanley R.↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

AC impedance analysis of polypyrrole thin films

The AC impedance spectra of thin polypyrrole films were obtained at open circuit potentials from -0.4 to 0.4 V vs SCE. Two limiting cases are discussed for which simplified equivalent circuits are applicable. At very positive potentials, the predominantly nonfaradaic AC impedance of polypyrrole is very similar to that observed previously for finite porous metallic films. Modeling of the data with the appropriate equivalent circuit permits effective pore diameter and pore number densities of the oxidized film to be estimated. At potentials from -0.4 to -0.3 V, the polypyrrole film is essentially nonelectronically conductive and diffusion of polymer oxidized sites with their associated counterions can be assumed to be linear from the film/substrate electrode interface. The equivalent circuit for the polypyrrole film at these potentials is that previously described for metal oxide, lithium intercalation thin films. Using this model, counterion diffusion coefficients are determined for both semi-infinite and finite diffusion domains. In addition, the limiting low frequency resistance and capacitance of the polypyrrole thin fims was determined and compared to that obtained previously for thicker films of the polymer. The origin of the observed potential dependence of these low frequency circuit components is discussed.

Penner, Reginald M.↗
NASA NTRS2003

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗
DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

Some TEM observations of Al2O3 scales formed on NiCrAl alloys

The microstructural development of Al2O3 scales on NiCrAl alloys has been examined by transmission electron microscopy. Voids were observed within grains in scales formed on a pure NiCrAl alloy. Both voids and oxide grains grew measurably with oxidation time at 1100 C. The size and amount of porosity decreased towards the oxide-metal growth interface. The voids resulted from an excess number of oxygen vacancies near the oxidemetal interface. Short-circuit diffusion paths were discussed in reference to current growth stress models for oxide scales. Transient oxidation of pure, Y-doped, and Zr-doped NiCrAl was also examined. Oriented alpha-(Al, Cr)2O3 and Ni(Al, Cr)2O4 scales often coexisted in layered structures on all three alloys. Close-packed oxygen planes and directions in the corundum and spinel layers were parallel. The close relationship between oxide layers provided a gradual transition from initial transient scales to steady state Al2O3 growth.

Smialek, J.↗
NASA NTRS1979

Some TEM observations of Al2O3 scales formed on NiCrAl alloys

The microstructural development of Al2O3 scales on NiCrAl alloys has been examined by transmission electron microscopy. Voids have been observed within grains in scales formed on a pure NiCrAl alloy. Both voids and oxide grains grew measurably with oxidation time at 1100 C. The size and amount of porosity decreased towards the oxide-metal growth interface. It was postulated that the voids resulted from an excess number of oxygen vacancies near the oxide-metal interface. Short-circuit diffusion paths were discussed in reference to current growth stress models for oxide scales. Transient oxidations of pure, Y-doped, and Zr-doped NiCrAl was also examined. Oriented alpha-(Al,Cr)2O3 and Ni(Al,Cr)2O4 scales often coexisted in layered structures on all three alloys. Close-packed oxygen planes and directions in the corundum and spinel layers were parallel. The close relationships between oxide layers provided a gradual transition from initial transient scales to steady state Al2O3 growth.

Smialek, J.↗
NASA NTRS1996

Ignition Delay Associated with a Strained Strip

Ignition processes associated with two adjacent fuel-oxidizer interferences bounding a strained fuel strip are explored here using single-step activation energy asymptotics. Calculations are made for constant as well as temporally decaying strain fields. There possible models of ignition are determined: one in which the two interfaces ignite independently as diffusion flames; one in which the two interfaces ignite dependently and in which ignition occurs to form a single , premixed flame at very high strain rates before ignition is completely prevented. In contrast to a single, isolated interface in which ignition can be prevented by overmatching heat production with heat convection due to strain, ignition of a strained fuel strip can also be prevented if the finite extend of fuel is diluted by oxidizer more quickly than heat production can cause a positive feedback thermal runaway. These behaviors are dependent on the relative sizes of timescales associated with species and heat diffusion, with convection due to strain, and with the chemical reaction. The result here indicate that adjacent, strained species interfaces may ignite quite differently in nature from ignition of a single, strained intrface and that their interdependence should be considered as the interfaces are brought closer together in complex strain fields. Critical strain rates leading to complete ignition delay are found to be considerably smaller for the fuel strip than those for single interfaces as the fuel strip is made thin in comparison to diffusion and chemical length scales.

Gerk, T. J.↗
NASA NTRS2016

Materials Outgassing Rate Decay in Vacuum at Isothermal Conditions

As a laboratory for scientific research, the International Space Station has been in Low Earth Orbit for nearly 20 years and is expected to be on-orbit for another 10 years. The ISS has been maintaining a relatively pristine contamination environment for science payloads. Materials outgassing induced contamination is currently the dominant source for sensitive surfaces on ISS and modeling the outgassing rate decay over a 20 to 30 year period is challenging. Materials outgassing is described herein as a diffusion-reaction process using ASTM E 1559 rate data. The observation of -1/2 (diffusion) or non-integers (reaction limited) as rate decay exponents for common ISS materials indicate classical reaction kinetics is unsatisfactory in modeling materials outgassing. Non-randomness of reactant concentrations at the interface is the source of this deviation from classical reaction kinetics. A diffusion limited decay was adopted as the result of the correlation of the contaminant layer thicknesses on returned ISS hardware, the existence of high outgassing silicone exhibiting near diffusion limited decay, and the confirmation of non-depleted material after ten years in the Low Earth Orbit.Keywords: Materials Outgassing, ASTM E 1559, Reaction Kinetics, Diffusion, Space Environments Effects, Contamination

Huang, Alvin Y.↗
NASA NTRS1995

Stream interfaces and energetic ions II: Ulysses test of Pioneer results

Ulysses measurements of energetic and solar wind particles taken near 5 AU between 20 and 30 degrees south latitude during a well-developed recurring corotating interaction region (CIR) show that the CIR's corotating energetic ion population (CEIP) associated with the trailing reverse shock starts within the CIR at the stream interface. This is consistent with an earlier result obtained by Pioneers 10 and 11 in the ecliptic plane between 4 and 6 AU. The Ulysses/Pioneer finding is noteworthy since the stream interface is not magnetically connected to the reverse shock, but lies 12-17 corotation hours from it. Thus, the finding seems to be inconsistent with the basic model that generates CEIP particles at the reverse shock and propagates them along field lines. Eliminating the inconsistency probably entails an extension of the standard model such as cross-field diffusion or a non-shock energization process operating near the stream interface closer to the sun.

Intriligator, Devrie S.↗
DOE OSTI2021

Carrier Diffusion Lengths in Continuously Grown and Etched-and-Regrown GaN Pin Diodes

Advanced GaN power devices are promising for many applications in high power electronics but performance limitations due to material quality in etched-and-regrown junctions prevent their widespread use. Carrier diffusion length is a critical parameter that not only determines device performance but is also a diagnostic of material quality. Here we present the use of electron-beam induced current to measure carrier diffusion lengths in continuously grown and etched-and-regrown GaN pin diodes as models for interfaces in more complex devices. Variations in the quality of the etched-and-regrown junctions are observed and shown to be due to the degradation of the n-type material. Finally, we observe an etched-and-regrown junction with properties comparable to a continuously grown junction.

42 ENGINEERING↗
NASA NTRS2015

Thermodynamic Derivation of the Activation Energy for Ice Nucleation

Cirrus clouds play a key role in the radiative and hydrological balance of the upper troposphere. Their correct representation in atmospheric models requires an understanding of the microscopic processes leading to ice nucleation. A key parameter in the theoretical description of ice nucleation is the activation energy, which controls the flux of water molecules from the bulk of the liquid to the solid during the early stages of ice formation. In most studies it is estimated by direct association with the bulk properties of water, typically viscosity and self-diffusivity. As the environment in the ice-liquid interface may differ from that of the bulk, this approach may introduce bias in calculated nucleation rates. In this work a theoretical model is proposed to describe the transfer of water molecules across the ice-liquid interface. Within this framework the activation energy naturally emerges from the combination of the energy required to break hydrogen bonds in the liquid, i.e., the bulk diffusion process, and the work dissipated from the molecular rearrangement of water molecules within the ice-liquid interface. The new expression is introduced into a generalized form of classical nucleation theory. Even though no nucleation rate measurements are used to fit any of the parameters of the theory the predicted nucleation rate is in good agreement with experimental results, even at temperature as low as 190 K, where it tends to be underestimated by most models. It is shown that the activation energy has a strong dependency on temperature and a weak dependency on water activity. Such dependencies are masked by thermodynamic effects at temperatures typical of homogeneous freezing of cloud droplets; however, they may affect the formation of ice in haze aerosol particles. The new model provides an independent estimation of the activation energy and the homogeneous ice nucleation rate, and it may help to improve the interpretation of experimental results and the development of parameterizations for cloud formation.

Cirrus clouds↗
NASA NTRS2009

Thermal Transport Model for Heat Sink Design

A document discusses the development of a finite element model for describing thermal transport through microcalorimeter arrays in order to assist in heat-sinking design. A fabricated multi-absorber transition edge sensor (PoST) was designed in order to reduce device wiring density by a factor of four. The finite element model consists of breaking the microcalorimeter array into separate elements, including the transition edge sensor (TES) and the silicon substrate on which the sensor is deposited. Each element is then broken up into subelements, whose surface area subtends 10 10 microns. The heat capacity per unit temperature, thermal conductance, and thermal diffusivity of each subelement are the model inputs, as are the temperatures of each subelement. Numerical integration using the Finite in Time Centered in Space algorithm of the thermal diffusion equation is then performed in order to obtain a temporal evolution of the subelement temperature. Thermal transport across interfaces is modeled using a thermal boundary resistance obtained using the acoustic mismatch model. The document concludes with a discussion of the PoST fabrication. PoSTs are novel because they enable incident x-ray position sensitivity with good energy resolution and low wiring density.

Chervenak, James A.↗
NASA NTRS2015

Influences of Models on the Unsteady Pressure Characteristics of the NASA National Transonic Facility

Pressure fluctuations have been measured over the course of several tests in the National Transonic Facility to study unsteady phenomenon both with and without the influence of a model. Broadband spectral analysis will be used to characterize the length scales of the tunnel. Special attention will be given to the large-scale, low frequency data that influences the Mach number and force and moment variability. This paper will also discuss the significance of the vorticity and sound fields that can be related to the Common Research Model and will also highlight the comparisons to an empty tunnel configuration. The effectiveness of vortex generators placed at the interface of the test section and wind tunnel diffuser showed promise in reducing the empty tunnel unsteadiness, however, the vortex generators were ineffective in the presence of a model.

Jones, Gregory↗