Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “isotherm model parameters”

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 163 records · Page 9

The interpretation of simultaneous soft X-ray spectroscopic and imaging observations of an active region

Simultaneous soft X-ray spectroscopic and broad-band imaging observations of an active region have been analyzed together to determine the parameters which describe the coronal plasma. From the spectroscopic data, models of temperature-emission measure-elemental abundance have been constructed which provide acceptable statistical fits. By folding these possible models through the imaging analysis, models which are not self-consistent can be rejected. In this way, only the oxygen, neon, and iron abundances of Pottasch (1967), combined with either an isothermal or exponential temperature-emission-measure model, are consistent with both sets of data. Contour maps of electron temperature and density for the active region have been constructed from the imaging data. The implications of the analysis for the determination of coronal abundances and for future satellite experiments are discussed.

Davis, J. M.↗

Computer Modeling of Non-Isothermal Crystallization

A realistic computer model for simulating isothermal and non-isothermal phase transformations proceeding by homogeneous and heterogeneous nucleation and interface-limited growth is presented. A new treatment for particle size effects on the crystallization kinetics is developed and is incorporated into the numerical model. Time-dependent nucleation rates, size-dependent growth rates, and surface crystallization are also included. Model predictions are compared with experimental measurements of DSC/DTA peak parameters for the crystallization of lithium disilicate glass as a function of particle size, Pt doping levels, and water content. The quantitative agreement that is demonstrated indicates that the numerical model can be used to extract key kinetic data from easily obtained calorimetric data. The model can also be used to probe nucleation and growth behavior in regimes that are otherwise inaccessible. Based on a fit to data, an earlier prediction that the time-dependent nucleation rate in a DSC/DTA scan can rise above the steady-state value at a temperature higher than the peak in the steady-state rate is demonstrated.

Kelton, K. F.↗

Application of BISON to UO 2 MiniFuel fission gas release analysis

There has been a recent push to accelerate fuel qualification by developing revolutionary capabilities to reduce irradiation periods, and thereby, reduce the time required to qualify a new fuel system. One such capability is the MiniFuel irradiation capsule designed to miniaturize fuel samples and irradiate “mini” fuel samples under isothermal temperature conditions. MiniFuel allows steady-state irradiations to decouple the traditionally coupled fission rate (i.e., power) and temperature parameters to understand and generate microstructures observed in fuel operated in a commercial reactor. Furthermore, this process offers the possibility to gather in situ data as well as postirradiation or transient data such as thermal conductivity, specific heat, fission gas diffusion and release, etc. However, accelerating fuel qualification is not solely reliant on generating large amounts of data but also on developing an informed test matrix designed to rapidly generate impactful data. Additionally, this process is reliant on fuel performance codes, such as BISON, to evaluate MiniFuel irradiations using existing material models. This process pinpoints model/data gaps, identifies desired irradiation conditions, and subsequently supports model validation and development. This work describes the use of BISON to perform a number of sensitivity studies designed to understand conditions that lead to fission gas release (FGR) under steady-state isothermal irradiation conditions and temperature transient conditions. The model is applied to a UO 2 MiniFuel example and shows an overall good qualitative agreement with experimental FGR annealing tests under different temperature conditions. It also accounts well for microstructural effects on FGR. When quantitatively compared with FGR data from previously irradiated 103 MWd/kgU UO 2 discs under thermal annealing, the model shows a less satisfactory agreement with the experimental data. Finally, a UO 2 MiniFuel test matrix is proposed to help to extend the model's operational range and validate the new FGR model capabilities to higher burnups and transient conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Calculation of atmospheric loss from microwave radiometric noise temperature measurements

Microwave propagation loss in the atmosphere can be inferred from microwave radiometric noise temperature measurements. The relevant equations are given and a derivation and calculation is made assuming various physical models. Comparison is made with the commonly used lumped element atmospheric model (isothermal and uniform loss) and the model with linear temperature and exponential loss distributions. The results are useful for estimating the integral inversion differences due to the model selection. This indicates that the commonly used lumped element atmospheric model is a very good approximation with judicious choice of the effective physical temperature. For the worst case comparison, the lumped element model agrees with the variable parameter model within 0.2 dB up to a propagation loss of 3 dB.

Stelzried, C.↗

High Temperature Degradation Mechanisms in Polymer Matrix Composites

Polymer matrix composites are increasingly used in demanding structural applications in which they may be exposed to harsh environments. The durability of such materials is a major concern, potentially limiting both the integrity of the structures and their useful lifetimes. The goal of the current investigation is to develop a mechanism-based model of the chemical degradation which occurs, such that given the external chemical environment and temperatures throughout the laminate, laminate geometry, and ply and/or constituent material properties, we can calculate the concentration of diffusing substances and extent of chemical degradation as functions of time and position throughout the laminate. This objective is met through the development and use of analytical models, coupled to an analysis-driven experimental program which offers both quantitative and qualitative information on the degradation mechanism. Preliminary analyses using a coupled diffusion/reaction model are used to gain insight into the physics of the degradation mechanisms and to identify crucial material parameters. An experimental program is defined based on the results of the preliminary analysis which allows the determination of the necessary material coefficients. Thermogravimetric analyses are carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. Powdered samples are used to eliminate diffusion effects. Tests in both inert and oxidative environments allow the separation of thermal and oxidative contributions to specimen mass loss. The concentration dependency of the oxidative reactions is determined from the tests in pure oxygen. Short term isothermal tests at different temperatures are carried out on neat resin and unidirectional macroscopic specimens to identify diffusion effects. Mass loss, specimen shrinkage, the formation of degraded surface layers and surface cracking are recorded as functions of exposure time. Geometry effects in the neat resin, and anisotropic diffusion effects in the composites, are identified through the use of specimens with different aspect ratios. The data is used with the model to determine reaction coefficients and effective diffusion coefficients. The empirical and analytical correlations confirm the preliminary model results which suggest that mass loss at lower temperatures is dominated by oxidative reactions and that these reaction are limited by diffusion of oxygen from the surface. The mechanism-based model is able to successfully capture the basic physics of the degradation phenomena under a wide range of test conditions. The analysis-based test design is successful in separating out oxidative, thermal, and diffusion effects to allow the determination of material coefficients. This success confirms the basic picture of the process; however, a more complete understanding of some aspects of the physics are required before truly predictive capability can be achieved.

Cunningham, Ronan A.↗

High Temperature Degradation Mechanisms in Polymer Matrix Composites

Polymer matrix composites are increasingly used in demanding structural applications in which they may be exposed to harsh environments. The durability of such materials is a major concern, potentially limiting both the integrity of the structures and their useful lifetimes. The goal of the current investigation is to develop a mechanism-based model of the chemical degradation which occurs, such that given the external chemical environment and temperatures throughout the laminate, laminate geometry, and ply and/or constituent material properties, we can calculate the concentration of diffusing substances and extent of chemical degradation as functions of time and position throughout the laminate. This objective is met through the development and use of analytical models, coupled to an analysis-driven experimental program which offers both quantitative and qualitative information on the degradation mechanism. Preliminary analyses using coupled diffusion/reaction model are used to gain insight into the physics of the degradation mechanisms and to identify crucial material parameters. An experimental program is defined based on the results of the preliminary analysis which allows the determination of the necessary material coefficients. Thermogravimetric analyses are carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. Powdered samples are used to eliminate diffusion effects. Tests in both inert and oxidative environments allow the separation of thermal and oxidative contributions to specimen mass loss. The concentration dependency of the oxidative reactions is determined from the tests in pure oxygen. Short term isothermal tests at different temperatures are carried out on neat resin and unidirectional macroscopic specimens to identify diffusion effects. Mass loss, specimen shrinkage, the formation of degraded surface layers and surface cracking are recorded as functions of exposure time. Geometry effects in the neat resin, and anisotropic diffusion effects in the composites, are identified through the use of specimens with different aspect ratios. The data is used with the model to determine reaction coefficients and effective diffusion coefficients. The empirical and analytical correlations confirm the preliminary model results which suggest that mass loss at lower temperatures is dominated by oxidative reactions and that these reaction are limited by diffusion of oxygen from the surface. The mechanism-based model is able to successfully capture the basic physics of the degradation phenomena under a wide range of test conditions. The analysis-based test design is successful in separating out oxidative, thermal, and diffusion effects to allow the determination of material coefficients. This success confirms the basic picture of the process; however, a more complete understanding of some aspects of the physics are required before truly predictive capability can be achieved.

Cunningham, Ronan A.↗

The X-ray spectrm of the Coma Cluster of galaxies

The HEAO 1 A-2 and OSO 8 observations of the Coma Cluster of galaxies in the 2-60 keV energy band are presented. The spectrum is not well fitted by an isothermal bremsstrahlung model with an Fe line. To model the spectrum more accurately, we have performed fits of polytropic models and find that the gas is well fitted by an equation of state with polytropic index ranging between 1.38 and 1.6. This analysis reduces the total cluster mass by a factor of 3-5 below the optically derived virial mass. Using the best-fitting parameters, the gas is found to be a significant fraction (25 percent) of the total cluster mass. The Fe K-alpha line blend is strongly detected, from which we infer 0.18-0.36 solar abundance in the X-ray emitting gas. The Fe K-beta line blend is also weakly detected. The spectrum is not consistent with the existence of a hard X-ray tail at the flux level reported by Bazzano el al. (1984). In addition, the HEAO I scanning mode data place strong constraints on possible contaminating sources in the region surrounding Coma.

Henriksen, M. J.↗

Influence of porous aluminosilicate grain size materials in experimental and modelling Cs + adsorption kinetics and wastewater column process

This paper focuses on the influence of the grain size of a geopolymer based adsorbent on its Cs + adsorption performances both in batch and fixed-bed process. The geopolymer phase was used as a binder to support NaY zeolite particle in a 20 wt% charged porous composite with 160 m 2 .g –1 of porous surface area. These samples were shaped with three grain sizes (50 /100/500 µm) to remove 80–90 mg/g of Cs + in batch and column operations. After their microstructural and porous characterizations, their efficiency and adsorption characteristics were investigated through adsorption isotherms and kinetic in the two processes. While the grain size has no influence on the maximal extraction capacity of the adsorbent, it strongly affects the sorption kinetic. By coupling experimental data and a modelling approach, the complex sorption mechanism was highlighted, suggesting a new insight of the contaminant sorption kinetic. Then, comparison of batch and column adsorption experiments illustrates the detailed explanation of various process parameters for column study. The results show challenges for fixed-bed column utilization by the choice of the appropriate grain size as a compromise between the material sorption kinetic and hydrodynamic considerations. Furthermore, this is of high importance to more accurately optimize the design of column adsorption to assess the transport of Cs+ in multi-porous tailored grain size materials.

36 MATERIALS SCIENCE↗

Phase Field Modeling of Microstructure Development in Microgravity

This newly funded project seeks to extend our NASA-sponsored project on modeling of dendritic microstructures to facilitate collaboration between our research group and those of other NASA investigators. In our ongoing program, we have applied advanced computational techniques to study microstructural evolution in dendritic solidification, for both pure isolated dendrites and directionally solidified alloys. This work has enabled us to compute dendritic microstructures using both realistic material parameters and experimentally relevant processing conditions, thus allowing for the first time direct comparison of phase field computations with laboratory observations. This work has been well received by the materials science and physics communities, and has led to several opportunities for collaboration with scientists working on experimental investigations of pattern selection and segregation in solidification. While we have been able to pursue these collaborations to a limited extent, with some important findings, this project focuses specifically on those collaborations. We have two target collaborations: with Prof. Glicksman's group working on the Isothermal Dendritic Growth Experiment (IDGE), and with Prof. Poirier's group studying directional solidification in Pb-Sb alloys. These two space experiments match well with our two thrusts in modeling, one for pure materials, as in the IDGE, and the other directional solidification. Such collaboration will benefit all of the research groups involved, and will provide for rapid dissemination of the results of our work where it will have significant impact.

Dantzig, Jonathan A.↗

A nonlinear model for gas chromatograph systems

Fundamental engineering design techniques and concepts were studied for the optimization of a gas chromatograph-mass spectrometer chemical analysis system suitable for use on an unmanned, Martian roving vehicle. Previously developed mathematical models of the gas chromatograph are found to be inadequate for predicting peak heights and spreading for some experimental conditions and chemical systems. A modification to the existing equilibrium adsorption model is required; the Langmuir isotherm replaces the linear isotherm. The numerical technique of Crank-Nicolson was studied for use with the linear isotherm to determine the utility of the method. Modifications are made to the method eliminate unnecessary calculations which result in an overall reduction of the computation time of about 42 percent. The Langmuir isotherm is considered which takes into account the composition-dependent effects on the thermodynamic parameter, mRo.

Feinberg, M. P.↗

Direct electric field heating and acceleration of electrons in solar flares

We show that the observed properties of solar flare X-ray and microwave emission can be explained through the Joule heating and electric field acceleration of runaway electrons in current channels. The global properties of the flaring region required for this are presented. We have fit a hybrid thermal/nonthermal electron distribution, consisting of hot, isothermal electrons with a nonthermal tail of runaway electrons, to high-resolution hard X-ray and microwave spectra and have obtained excellent fits to both. The hybrid model relaxes the electron number and energy flux requirements for the hard X-ray emission over those of a purely nonthermal model. The model also provides explanations for several previously unexplained aspects of the high-resolution microwave spectra. The fit parameters can be related to physical properties (such as the electric field strength in the current channels) of the acceleration region.

Holman, Gordon D.↗

Globular cluster winds with central accretion by a massive compact object or subcluster

Steady-state isothermal wind flows are computed with accretion by a central mass concentration in spherically symmetric systems of gas-losing stars. Solutions are uniquely determined by two dimensionless parameters lambda and m, where lambda is approximately equal to the square of the ratio of the central escape velocity to the sound speed and m is approximately equal to the ratio of the central point mass to the cluster core mass. The models are applied to globular clusters for two types of central mass concentration - a single massive black hole and a subcluster of massive stars. A hard (greater than a keV) X-ray source due to a central black hole fed by clusterwide inflow does not provide a consistent model for globular cluster X-ray sources. However, in the case of hot (T being greater than or approximately equal to 2 x 10 to the 5th K) winds, gas trapping by a central subcluster of neutron stars or of binaries containing white dwarfs could explain recent UV and H-alpha observations. Similar applications to elliptical galaxies are discussed.

Durisen, R. H.↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Anisotropic magneto-resistivity and magnetocaloric effect in DyAl 2

Here, in this study, we experimentally and theoretically investigate the correlation between the anisotropic magnetic resistivity and magnetocaloric effect in DyAl 2 single crystals. DyAl 2 crystallizes in the C15 Laves phase structure with cubic symmetry. While some earlier studies revealed that the isothermal entropy change in dialuminides follows a similar trend as the electrical resistivity change as a function of temperature and magnetic field, the correlation between anisotropic behavior of these two physical parameters is yet to be explored to the best of our knowledge. To address this gap, we measured electrical resistivity of DyAl 2 single crystals along two different directions in zero and applied magnetic fields and compared the experimental results with our theoretical model. For the theoretical analysis, we employed a model Hamiltonian in the mean-field approximation, considering the exchange interaction, Zeeman effect, and crystalline electric field associated with the cubic symmetry. Our findings reveal a significant dependence of DyAl 2 resistivity on the direction of the applied magnetic field, in agreement with our theoretical results. These outcomes emphasize the interplay between magnetocaloric effect and magneto-resistivity, underscoring the potential of the magnetocaloric effect as a valuable tool for gaining deeper insights into basic physical properties including electron transport behavior.

36 MATERIALS SCIENCE↗

Raman and X-ray diffraction study of pressure-induced phase transition in synthetic Mg 2 TiO 4

Synthetic Mg 2 TiO 4 qandilite was investigated to 50 and 40.4 GPa at room temperature using Raman spectroscopy and X-ray diffraction, respectively. The Raman measurements showed that cubic Mg 2 TiO 4 spinel transforms to a high pressure tetragonal (I4 1 /amd, No.141) phase at 14.7 GPa. Owing to sluggish kinetics at room temperature, the spinel phase coexists with the tetragonal phase between 14.7 and 24.3 GPa. In the X-ray diffraction experiment, transformation of the cubic Mg 2 TiO 4 to the tetragonal structure was complete by 29.2 GPa, ~5 GPa higher than the transition pressure obtained by Raman measurements, owing to slow kinetics. The obtained isothermal bulk modulus of Mg 2 TiO 4 spinel is K T0 = 148(3) GPa when K T0 ’ = 6.6, or K T0 = 166(1) GPa when K T0 ’ is fixed at 4. The isothermal bulk modulus of the high-pressure tetragonal phase is calculated to be 209(2) GPa and V 0 = 270(2) Å 3 when K T0 ’ is fixed at 4, and the volume reduction on change from cubic to tetragonal phase is about 9%. The calculated thermal Grüneisen parameters (γ th ) of cubic and tetragonal Mg 2 TiO 4 phases are 1.01 and 0.63. Based on the radii ratio of spinel cations, a simple model is proposed to predict post-spinel structures.

36 MATERIALS SCIENCE↗

A hybrid thermal/nonthermal model for the energetic emissions from solar flares

We present an alternative interpretation of the high-resolution solar flare spectra obtained in 1981 by Lin et al. In our interpretation electron heating and acceleration are simultaneous and physically linked, rather than heating being a secondary effect of particle acceleration. We show that the observed properties of solar flare X-ray emission can in general be explained through the Joule heating and electric field acceleration of runaway electrons in current channels. We have fitted a hybrid thermal/nonthermal electron distribution, consisting of hot isothermal electrons with a nonthermal tail of runaway electrons, to representative spectra obtained by Lin et al. The hybrid model relaxes the electron number and energy flux requirements for the hard X-ray emission over those of a purely nonthermal model. The low-energy 'cutoff' to the nonthermal part of the X-ray spectrum is determined by the critical velocity in the electron distribution above which runaway acceleration occurs. We relate the fit parameters to the physical properties (such as the electric field strength in the current channels) of the acceleration region.

Holman, Gordon D.↗

Lateral Shock of the R Aquarii Jet

The R Aqr jet was observed with the VLA B-configuration at two epochs separated by approximately 13.2 yr. Comparison of the resulting 6 cm continuum images show that the radio jet has undergone a lateral counterclockwise rotation of approximately 6 deg-12 deg on the plane of the sky. The model of jet parcels on independent trajectories is difficult to reconcile with these observations and leads us to consider a path-oriented jet (i.e., younger parcels follow the same path as older parcels). Comparison of the most recent radio image with a nearly contemporaneous HST/FOC ultraviolet image at approximately 2330 Angstroms suggests that the ultraviolet emission lies along the leading side of the rotating radio jet. In conjunction with a proper motion analysis of the jet material that yields empirical space-velocity and resulting acceleration-magnitude relationships as a function of distance from the central source, we evaluate the observational results in terms of a schematic model in which the jet emission consists of plane-parallel isothermal shocks along the leading edge of rotation. In such a radiating shock, the ultraviolet-emitting region is consistent with the adiabatic region in the form of a high-temperature, low-density sheath that surrounds the cooled postshock radio-emitting region. Within the context of the schematic model, we obtain the temperatures, densities, and pressures within the preshock, adiabatic, and postshock regions as a function of distance from the central source; the physical parameters so derived compare favorably to previously published estimates. We obtain a total jet mass of 3.1 x 10(exp -5) solar mass and an age of approximately 115 yr. We evaluate the model in the context of its density-boundary condition, its applicability to an episodic or quasi-continuous jet, and angular momentum considerations.

Hollis, J. M.↗

Bound oscillations on thin magnetic flux tubes - Convective instability and umbral oscillations

The possibility that 'tube waves' can be trapped on slender solar magnetic flux tubes is investigated. For rigid isothermal flux tubes, it is found that the flux tube geometry can by itself lead to waves which are trapped on the part of the tube that expands with height. Some geometries lead to trapped modes with eigenperiods near 180 s, if parameters appropriate to sunspot umbrae are chosen. It is possible that the umbral oscillations are a manifestation of such trapped waves, if sunspot umbrae consist of an assembly of slender flux tubes, as in the spaghetti model of Parker (1979). For flux tubes which have a constant ratio of Alfven speed to sound speed, it is found that it is primarily the variation of temperature with height which determines whether trapped waves can exist. Certain temperature profiles lead to disturbances for which omega squared is less than zero, corresponding to convective instability or Rayleigh-Taylor instability.

Hollweg, J. V.↗