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 37 records · Page 2

The physical implications of an isothermal model for the hot intracluster medium

X-ray fluxes from HEAO-1 A2 and Einstein Imaging Proportional Counter (IPC) observations of clusters of galaxies were used to constrain the parameter beta in the isothermal surface brightness profile. Beta is found primarily to have values between .50 and .75 for 15 clusters. Eight of these objects have values of beta previously measured using imaging observations. For these clusters good agreement is found with the values reported here implying that this profile is a good description of the surface brightness out to 8 to 10 core radii. The total gas mass and radial distribution (assuming spherical symmetry) within the cluster resulting from the isothermal model imply an extended halo of hot gas which has 30 to 60% of the virial mass for some clusters.

Hendriksen, M. J.↗

A robust three-parameter reference curve for condensed phase materials

Many isothermal equations of state lead to numerical problems under tension conditions due to negative bulk modulus. We summarize the properties of four existing models and devise a continuous non-piecewise reference curve to address this issue. The final functional form contains three fitting parameters that are directly related to the bulk modulus and its pressure derivative at the reference state and at infinite pressures. We also include bounds that guarantee the monotonicity of its pressure derivative in all volume space. To demonstrate the qualitative features of each isotherm, we conduct calibrations against experimental data for cyclotetramethylene tetranitramine, magnesium oxide, gold, and platinum. Furthermore, pressure residuals and root mean square deviations show that the proposed functional form is capable of adequately fitting the data compared to other three-parameter models.

36 MATERIALS SCIENCE↗

Standing shocks in a two-fluid solar wind

We present a numerical study of the formation of standing shocks in the solar wind using a two-fluid time-dependent model in the presence of Alfven waves. Included in this model is the adiabatic cooling and thermal conduction of both electrons and protons. In this study, standing shocks develop in the flow when additional critical points form as a result of either localized momentum addition or rapid expansion of the flow tube below the existing sonic point. While the flow speed and density exhibit the same characteristics as found in earlier studies of the formation of standing shocks, the inclusion of electron and proton heat conduction produces different signatures in the electron and proton temperature profiles across the shock layer. Owing to the strong heat conduction, the electron temperature is nearly continuous across the shock, but its gradient has a negative jump across it, thus producing a net heat flux out of the shock layer. The proton temperature exhibits the same characteristics for shocks produced by momentum addition but behaves differently when the shock is formed by the rapid divergence of the flow tube. The adiabatic cooling in a rapidly diverging flow tube reduces the proton temperature so substantially that the proton heat conduction becomes negligible in the vicinity of the shock. As a result, protons experience a positive jump in temperature across the shock. While Alfven waves do not affect the formation of standing shocks, they contribute to the change of the mmomentum and energy balance across them. We also find that for this solar wind model the inclusion of thermal conduction and adiabatic cooling for the elctrons and protons increases significantly the range of parameters characterizing the formation of standing shocks over those previously found for isothermal and polytropic models.

Habbal, Shadia R.↗

Development of a Rational Modeling Approach for the Design, and Optimization of the Multifiltration Unit

This thesis includes the development and verification of an adsorption model for analysis and optimization of the adsorption processes within the International Space Station multifiltration beds. The fixed bed adsorption model includes multicomponent equilibrium and both external and intraparticle mass transfer resistances. Single solute isotherm parameters were used in the multicomponent equilibrium description to predict the competitive adsorption interactions occurring during the adsorption process. The multicomponent equilibrium description used the Fictive Component Analysis to describe adsorption in unknown background matrices. Multicomponent isotherms were used to validate the multicomponent equilibrium description. Column studies were used to develop and validate external and intraparticle mass transfer parameter correlations for compounds of interest. The fixed bed model was verified using a shower and handwash ersatz water which served as a surrogate to the actual shower and handwash wastewater.

Hand, David W.↗

Development of a near-isothermal transcritical CO 2 compression system with a liquid piston compressor

Compressors are critical components in vapor compression cycle systems, significantly contributing to energy consumption. As global demand for HVAC&R systems rises, enhancing compressor efficiency becomes increasingly vital. This paper introduces a novel liquid piston compressor integrated with a gas cooler for the transcritical CO 2 refrigeration cycle. The liquid piston enables CO 2 refrigerant compression within various types of heat exchangers, facilitating the transfer of compression heat to the heat transfer fluid. By releasing significant heat, this compressor allows for removing or downsizing the traditional gas cooler in the refrigeration system. This paper presents the experimental performance of the first near-isothermal compressor utilizing a liquid piston in a vapor compression cycle. The critical parameters affecting heat transfer are analyzed by using a 1-D simulation model to achieve a near-isothermal compression process. The results show that the developed prototype successfully reduced the compression temperature increase from 95 K to 10 K, achieving 90 % isothermal efficiency. Furthermore, the 1-D simulation results suggest the smaller internal diameter tubes benefit the isothermal efficiency the most.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Thermally conductive magnetohydrodynamic flows in helmet-streamer coronal structures

The behavior of thermally conductive plasma flows in helmet-streamer coronal structures is investigated within the framework of the axisymmetric nonrotating one-fluid MHD model. Continuous subsonic-supersonic solutions satisfying observed boundary conditions at the sun as well as the vanishing of the temperature at infinity are obtained and presented. Special attention is paid to the combined effects of conductive flow (and corresponding thermal force) and rapidly diverging magnetic field on the critical points. In this, the heliocentric distance of the neutral point determining the separation between closed and open field lines (cusp) is treated as a free parameter. These thermally conductive solutions are contrasted with those provided by corresponding isothermal models.

Cuperman, S.↗

Shape Parameter for a Non-Axisymmetric Isothermal Dendrite

In previous work, we found approximate solutions for paraboloids having perturbations with four-fold axial symmetry in order to model dendritic growth in cubic materials. These solutions provide self-consistent corrections through second order in a shape parameter e to the Peclet number-supercooling relation of the Ivantsov solution. The parameter e is proportional to the amplitude of the four-fold correction to the dendrite shape, as measured from the Ivantsov paraboloid of revolution. We calculate e by comparing the dendrite tip shape to the portion of the equilibrium shape near the growth direction, (001), for anisotropic surface free energy, where the ni are components of the unit normal of the crystal surface. This comparison results in epsilon = -2(epsilon 4), independent of the Peclet number. From the experimental value of epsilon 4, we find epsilon approximately 0.011, in good agreement with the measured value epsilon approximately 0.008 of LaCombe et al.

McFadden, G. B.↗

Numerical Model for Isothermal and Non-Isothermal Crystallization of Liquids and Glasses

A new numerical model of isothermal and non-isothermal first order phase transformations, such as the crystallization of liquids and glasses, is presented. This model computes directly the volume fraction transformed, taking into account time-dependent nucleation rates and cluster-size-dependent growth velocities. The model is applied to the crystallization of lithium disilicate glass, using the appropriate kinetic and thermodynamic parameters. The model is used (1) to determine the validity of common methods for computing the volume fraction transformed as a function of time in isothermal experiments when a time-dependent nucleation rate is expected, (2) to simulate non-isothermal differential scanning calorimetry (DSC) studies of glass devitrification as a function of scan rate, and (3) to compute the effect of preannealing on the DSC peak parameters. A novel behavior of the nucleation rate with scan rate is predicted, arising because the relaxation of the cluster distribution cannot be described by a single relaxation time. Comparisons of the calculations with experimental data on this glass demonstrate the validity of the model.

Kelton, K. F.↗

Optimizing Process Parameters to Produce Single Phase YBa2Cu3O7 Powder

Process parameters such as temperature and time are varied to obtain single phase (pure) YBa2Cu307 powders. X-ray diffraction patterns of superconducting and nonsuperconducting powders are observed in situ in normal and elevated temperature environments to confirm known strong peaks that are unique to YBa2Cu307. Peaks vary in magnitude as a result of superconducting - nonsuperconducting phase transitions. Known strong peaks serve as our basis in monitoring transitional changes. Peak isotherms obtained experimentally are correlated with theoretical models to identify the weighting factor that characterizes the optimum isotherm. By isolating the parameters under which 100% peak growth occurs most rapidly, we obtain the optimum isotherm. The choices made for firing temperature and atmosphere were determined to be the primary factors necessary to ensure sample purity.

Hurley, J. S.↗

Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity (Quarterly Research Performance Progress Report: Q3Y3)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. To start, we have been working towards adjusting numerical parameters to reproduce DFIT tests with different fracture geometries. Then, we have focused modeling the DFIT conducted at well 58-32 (cycle 4 in zone 2). Specifically, we have carefully analyzed the influence of the aperture-contact stress relationship and investigated how two different laws affect model results. We have also introduced thermal effects in the model and rerun the simulations to estimate the effect of cooling on the DFIT response. Due to the small amount of fluid injected during a DFIT, thermal effects seem to be moderate.

15 GEOTHERMAL ENERGY↗

Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity: Quarterly Research Performance Progress Report

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. To start, we have been working towards adjusting numerical parameters to reproduce DFIT tests with different fracture geometries. Then, we have focused modeling the DFIT conducted at well 58-32 (cycle 4 in zone 2). Specifically, we have carefully analyzed the influence of the aperture-contact stress relationship and investigated how two different laws affect model results.

58 GEOSCIENCES↗

A Numerical Method for Calculating Stellar Occultation Light Curves from an Arbitrary Atmospheric Model

We present a method for speeding up numerical calculations of a light curve for a stellar occultation by a planetary atmosphere with an arbitrary atmospheric model that has spherical symmetry. This improved speed makes least-squares fitting for model parameters practical. Our method takes as input several sets of values for the first two radial derivatives of the refractivity at different values of model parameters, and interpolates to obtain the light curve at intermediate values of one or more model parameters. It was developed for small occulting bodies such as Pluto and Triton, but is applicable to planets of all sizes. We also present the results of a series of tests showing that our method calculates light curves that are correct to an accuracy of 10(exp -4) of the unocculted stellar flux. The test benchmarks are (i) an atmosphere with a l/r dependence of temperature, which yields an analytic solution for the light curve, (ii) an atmosphere that produces an exponential refraction angle, and (iii) a small-planet isothermal model. With our method, least-squares fits to noiseless data also converge to values of parameters with fractional errors of no more than 10(exp -4), with the largest errors occurring in small planets. These errors are well below the precision of the best stellar occultation data available. Fits to noisy data had formal errors consistent with the level of synthetic noise added to the light curve. We conclude: (i) one should interpolate refractivity derivatives and then form light curves from the interpolated values, rather than interpolating the light curves themselves; (ii) for the most accuracy, one must specify the atmospheric model for radii many scale heights above half light; and (iii) for atmospheres with smoothly varying refractivity with altitude, light curves can be sampled as coarsely as two points per scale height.

Chamberlain, D. M.↗

Simultaneous concentration and purification through gradient deformation chromatography

Mobile-phase additives, commonly used to modulate absorbate retention in gradient elution chromatography, are usually assumed to be either linearly retained or unretained. Previous theoretical work from our laboratory has shown that these modulators, such as salts in ion-exchange and hydrophobic interaction chromatography and organic modifiers in reversed-phase chromatography, can absorb nonlinearly, giving rise to gradient deformation. Consequently, adsorbate peaks that elute in the vicinity of the head of the deformed gradient may exhibit unusual shapes, form shoulders, and/or be concentrated. These effects for a reversed-phase sorbent with aqueous acetonitrile (ACN) as the modulator are verified experimentally. Gradient deformation is demonstrated experimentally and agrees with simulations based on ACN isotherm parameters that are independently determined from batch equilibrium studies using the layer model. Unusual absorbate peak shapes were found experimentally for single-component injections of phenylalanine, similar to those calculated by the simulations. A binary mixture of tryptophan and phenylalanine is used to demonstrate simultaneous concentration and separation, again in agreement with simulations. The possibility of gradient deformation in ion-exchange and hydrophobic interaction chromatography is discussed.

NASA Discipline Number 93-10↗

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗

Models for Galaxy halos in an open universe

The properties of simple, self-consistent infall models which embed an 'isothermal' dark halo in an open universe are presented. These models depend on one scaling parameter and two shape parameters. The first may be taken to be the circular velocity in the inner regions, and the others to be the density of the universe, Omega, and a measure, e, of the shape of orbits. These models provide a useful description of ensembles of satellite galaxies, or of binary galaxies, because they allow explicitly for the fact that such systems have orbital periods approaching the age of the universe and so cannot have a random orbital phase.

White, Simon D. M.↗

Tidal disruption of inviscid protoplanets

Roche showed that equilibrium is impossible for a small fluid body synchronously orbiting a primary within a critical radius now termed the Roche limit. Tidal disruption of orbitally unbound bodies is a potentially important process for planetary formation through collisional accumulation, because the area of the Roche limit is considerably larger then the physical cross section of a protoplanet. Several previous studies were made of dynamical tidal disruption and different models of disruption were proposed. Because of the limitation of these analytical models, we have used a smoothed particle hydrodynamics (SPH) code to model the tidal disruption process. The code is basically the same as the one used to model giant impacts; we simply choose impact parameters large enough to avoid collisions. The primary and secondary both have iron cores and silicate mantles, and are initially isothermal at a molten temperature. The conclusions based on the analytical and numerical models are summarized.

Boss, Alan P.↗

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗

Chemoviscosity modeling for thermosetting resin systems, part 3

A new analytical model for simulating chemoviscosity resin has been formulated. The model is developed by modifying the well established Williams-Landel-Ferry (WLF) theory in polymer rheology for thermoplastic materials. By introducing a relationship between the glass transition temperature (T sub g (t)) and the degree of cure alpha(t) of the resin system under cure, the WLF theory can be modified to account for the factor of reaction time. Temperature-dependent functions of the modified WLF theory parameters C sub 1 (T) and C sub 2 (T) were determined from the isothermal cure data. Theoretical predictions of the model for the resin under dynamic heating cure cycles were shown to compare favorably with the experimental data. This work represents a progress toward establishing a chemoviscosity model which is capable of not only describing viscosity profiles accurately under various cure cycles, but also correlating viscosity data to the changes of physical properties associated with the structural transformations of the thermosetting resin systems during cure.

Hou, T. H.↗