Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complete equation of state”

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 109 records · Page 6

Dynamics of heavy quarks in strongly coupled $\mathcal{N}$ = 4 SYM plasma

We calculate the probability distribution P(k) for a heavy quark with velocity v propagating through strongly coupled N = 4 SYM plasma in the ’t Hooft limit (N c → ∞, λ = g 2 N c → ∞) at a temperature T to acquire a momentum k due to interactions with the plasma. This distribution encodes the well-known drag coefficient η D and the transverse and longitudinal momentum diffusion coefficients κ T and κ L . The jet quenching parameter $\hat{q}$ can be extracted from P(k) for v = 1. Going beyond these known Gaussian characteristics of P(k), our calculation determines all of the higher order and mixed moments to leading order in 1/$\sqrt{λ}$ for the first time. These non-Gaussian features of P(k) include qualitatively novel correlations between longitudinal energy loss and transverse momentum broadening at nonzero v. We show that all higher moments scale characteristically with an effective temperature of the boosted plasma in the heavy quark rest frame, and we demonstrate that these non-Gaussian characteristics can be sizable in magnitude and even dominant in physically relevant situations. We use these results to derive a Kolmogorov equation for the evolution of the probability distribution for the total momentum of a heavy quark that propagates through strongly coupled plasma. This evolution equation accounts for all higher order correlations between transverse momentum broadening and longitudinal energy loss, which we have calculated from first principles. It reduces to a Fokker-Planck equation when truncated to only include the effects of η D , κ T and κ L . Remarkably, while heavy quarks do not reach kinetic equilibrium with the plasma if evolved with this Fokker-Planck equation, by showing that the Boltzmann distribution is a static solution of the all-order Kolmogorov equation that we have derived we demonstrate that heavy quarks do reach kinetic equilibrium if evolved with this equation. Our results thus provide a dynamically complete framework for understanding the thermalization of a heavy quark that may be initially far from equilibrium in the strongly coupled N = 4 SYM plasma — as well as new insight into heavy quark transport and equilibration in quark-gluon plasma.

Holography and Hydrodynamics↗

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Completion of Surrogate Testing for Low-Temperature TEX and a Look Towards the Future

To address the mounting need for below room temperature nuclear data validation, the Low-Temperature Thermal Epithermal eXperiments (LT-TEX) have been designed. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An optimization-based approach to tailor the mechanical response of soft metamaterials undergoing rate-dependent instabilities

An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metamaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. Here, these effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.

97 MATHEMATICS AND COMPUTING↗

Synthesis and Ultrahigh Pressure Compression of High-Entropy Boride (Hf0.2Mo0.2Nb0.2Ta0.2Zr0.2)B2 to 220 GPa

The high-entropy boride (Hf0.2Mo0.2Nb0.2Ta0.2Zr0.2)B2 material was synthesized under high-pressures and high-temperatures in a large-volume Paris-Edinburgh (PE) press from a ball-milled powder mix of HfO2, MoO3, Nb2O5, Ta2O5, ZrO2, carbon black, and boron carbide. The transformation process was monitored in situ by energy-dispersive x-ray diffraction with conversion starting at 1100 °C and completed by 2000 °C with the formation of a single hexagonal AlB2-type phase. The synthesized sample was recovered, powdered, and mixed with platinum pressure marker and studied under high pressure by angle-dispersive x-ray diffraction in a diamond anvil cell. The hexagonal AlB2-type phase of (Hf0.2Mo0.2Nb0.2Ta0.2Zr0.2)B2 was found to be stable up to the highest pressure of 220 GPa reached in this study (volume compression V/V0 = 0.70). The third order Birch-Murnaghan equation of state fit to the high-pressure data up to 220 GPa results in an ambient pressure unit cell volume V0=28.16±0.04 Å3, bulk modulusKo = 407 ± 6 GPa, pressure derivative of bulk-modulus K0′ = 2.73 ± 0.045 GPa. Our study indicates that this high-entropy boride (Hf0.2Mo0.2Nb0.2Ta0.2Zr0.2)B2 material is stable to ultrahigh pressures and temperatures and exhibit high bulk modulus similar to other incompressible transition metal borides like ReB2 and Os2B3.

36 MATERIALS SCIENCE↗

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Symmetry structure of a Riccati equation appearing in penetration mechanics

In the design of projectile penetration experiments a matter of considerable interest is scaling: that is, the potential relevance of small-scale experiments to their full-scale counterparts, in a manner analogous to that most often encountered in the context of fluid mechanics. From the theoretical standpoint, phenomena associated with scaling and scalability can be assessed using the well-established tools of dimensional analysis and the Buckingham-Pi Theorem. However, the familiar precepts of dimensional analysis are themselves a specific manifestation of the broader group invariance properties or symmetries of a mathematical model. Here, this work explores these notions – that is, dimensional analysis, the conditions for realizing complete similarity, and any additional symmetry structures – in the context of a Riccati differential equation appearing in the context of penetration mechanics. The aim of the investigation is twofold: 1) to complement existing empirical considerations with a concrete theoretical basis, and 2) to provide a deeper theoretical understanding of the projectile penetration model and its many implications.

97 MATHEMATICS AND COMPUTING↗

Non-Abelian electric field correlator at NLO for dark matter relic abundance and quarkonium transport

We perform a complete next-to-leading order calculation of the non-Abelian electric field correlator in a SU(N c ) plasma, which encodes properties of the plasma relevant for heavy particle bound state formation and dissociation, and is different from the correlator for the heavy quark diffusion coefficient. The calculation is carried out in the real-time formalism of thermal field theory and includes both vacuum and finite temperature contributions. By working in the R ξ gauge, we explicitly show the results are gauge independent, infrared and collinear safe. The renormalization group equation of this electric field correlator is determined by that of the strong coupling constant. Our next-to-leading order calculation can be directly applied to any dipole singlet-adjoint transition of heavy particle pairs. For example, it can be used to describe dissociation and (re)generation of heavy quarkonia inside the quark-gluon plasma well below the melting temperature, as well as heavy dark matter pairs (or charged co-annihilating partners) in the early universe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling and ZVS Operation of the Isolated Modular Multilevel DC–DC Converter With a Unified Trapezoidal Wave Modulation

Here, this article introduces a unified trapezoidal wave (UTW) modulation scheme for the isolated modular multilevel dc–dc (IM2dc) converter, which consolidates multiple existing modulation strategies for the IM2dc converter, including the quasi-square wave (QSW) modulation, trapezoidal wave modulation, and sinusoidal wave modulation, into a unified framework. Furthermore, this article introduces the harmonic state-space (HSS) equations to model the IM2dc converter based on the UTW modulation method. The HSS model operates in the frequency domain, enabling it to circumvent the complexities associated with time-domain analysis and seamlessly integrate with the UTW modulation. This article proceeds to analyze the real and reactive power transfer characteristics of the IM2dc converter, as well as the power factor, considering the influence of multiple modulation parameters. Subsequently, it delves into the examination of zero voltage switching (ZVS) conditions for the IM2dc converter based on the UTW modulation and the HSS model. The complete ZVS boundaries of the IM2dc converter, taking various voltage ratios into account, are derived. This article also illustrates the effects of including harmonic orders in the modeling process, modulation parameters, and internal harmonic ripples on the ZVS boundaries. Finally, experimental validation of the analyses is conducted on a down-scaled prototype.

42 ENGINEERING↗

What if GW190425 did not produce a black hole promptly?

ABSTRACT It is widely believed that the binary neutron star merger GW190425 produced a black hole promptly upon merger. Motivated by the potential association with the fast radio burst FRB 20190425A, which took place 2.5 h after the merger, we revisit the question of the outcome of GW190425 by means of numerical relativity simulations. We show that current laboratory and astrophysical constraints on the equation of state of dense matter do not rule out the formation of a long-lived remnant. However, the formation of a stable remnant would have produced a bright kilonova, in tension with upper limits by ZTF at the location and time of FRB 20190425A. Moreover, the ejecta would have been optically thick to radio emission for days to months, preventing a putative FRB from propagating out. The predicted dispersion measure is also several orders of magnitude larger than that observed for FRB 20190425A. Our results indicate that FRB 20190425A and GW190425 are not associated. However, we cannot completely rule out the formation of a long-lived remnant, due to the incomplete coverage of the relevant sky regions. More observations of GW190425-like events, including potential upper limit, have the potential to constrain nuclear physics. To this aim, it is important that follow-up observational campaigns of gravitational wave events are informed by the properties of the source, such as their chirp mass, and we urge the LIGO-Virgo-KAGRA collaboration to promptly release them publicly.

79 ASTRONOMY AND ASTROPHYSICS↗

Electromagnetic Transient Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (Distribution Grid Models)

The distribution and transmission grids are observing an increased penetration of power electronics in loads and generations. For example, there is increasing interest in integrating in extreme fast charging (XFC) systems for fast charging of electrical vehicles. As these systems are integrated, developing high-fidelity electromagnetic transient model of XFC systems in distribution grids and evaluating their interactions with the power grid would be of significant interest. This model will be utilized for design of XFC systems, to identify upgrades in distribution and/or transmission grids, for planning purposes by transmission planners or operators or owners, among others. It can also be utilized in operations for improved reliable performance of the grid and/or XFC station. The challenge with simulating these models is the high computational complexity introduced by the large number of states present in the system and the time-step needed to simulate the system. In this paper, advanced simulations algorithms are applied to reduce the computational complexity of simulating large-scale XFC systems. The algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. While the first three algorithms split the matrix that needs to be inverted from a large matrix to much smaller matrices, the final algorithm reduces the computational burden of applying higher-order integration approaches in the complete system. The comparison made in the previous sentence is with respect to use of homogeneous integration approaches used in conventional electromagnetic transient simulators like power systems computer aided design (PSCAD). The approaches mentioned here have resulted in speed-up of 36x in the simulation of a single distribution system with 15 XFCs.

Debnath, Suman↗

A limit on the total lepton number in the Universe from BBN and the CMB

At temperatures below the QCD phase transition, any substantial lepton number in the Universe can only be present within the neutrino sector. In this work, we systematically explore the impact of a non-vanishing lepton number on Big Bang Nucleosynthesis (BBN) and the Cosmic Microwave Background (CMB). Relying on our recently developed framework based on momentum averaged quantum kinetic equations for the neutrino density matrix, we solve the full BBN reaction network to obtain the abundances of primordial elements. We find that the maximal primordial total lepton number L allowed by BBN and the CMB is -0.12 (-0.10) ≤ L ≤ 0.13 (0.12) for NH (IH), while specific flavor directions can be even more constrained. This bound is complementary to the limits obtained from avoiding baryon overproduction through sphaleron processes at the electroweak phase transition since, although numerically weaker, it applies at lower temperatures and is obtained completely independently. We publicly release the C++ code COFLASY-C on GitHub (https://github.com/mariofnavarro/COFLASY/tree/COFLASY-C) which solves for the evolution of the neutrino quantum kinetic equations numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of Strain on the Resistivity and Thermal Conductivity of High Purity Niobium

High purity niobium (Nb) is a technologically important material for large-scale accelerator and nano-scale quantum computing applications in microwave frequency range. The high thermal conductivity and low resistivity of Nb are critical to the high performance at temperature range of 0.01–4.0 K. The presence of interstitials such as O, N, H, and C act as scattering centers and alter the mean free path, reducing resistivity and thermal conductivity, and contributing significantly to Nb's thermal performance for temperatures of 2.0 K and above. The residual resistivity ratio (RRR), defined as the ratio of the normal state resistivity at 300 K to that at 4.2 K (Nb, T c = 9.2 K), is an accepted direct estimate of the impurity content of fully recrystallized Nb. Complete re-crystallization of Nb is challenging unless very high temperatures are employed, which is often impractical, hence, in practice, dislocation and dislocation structures impact the thermal performance of Nb due to strong phonon scattering contributions. Furthermore, this paper reports on the degradation of thermal conductivity and RRR of high purity Nb large grain, single crystal with fixed impurity, varying strain, and dislocation content levels. Experimental thermal conductivity data fits the Boltzmann transport equation incorporating dislocation density.

36 MATERIALS SCIENCE↗

Interpretation of wide-angle x-ray diffraction patterns from mesophase pitch-based carbon fibers – a simulation and experimental study

Carbon fiber is a critical material in a wide range of industries, where it is highly valued for its high specific strength/stiffness, excellent wear resistance, efficient electrical and thermal transport properties, chemical resistance, and low coefficient of thermal expansion. The properties of a specific carbon fiber are closely tied to its structural characteristics at all length scales. Here, in this work, we applied wide-angle x-ray diffraction to a set of heat-treated mesophase pitch-based carbon fibers, with the goal of elucidating the crystalline structures as a function of fiber orientation. To assist with analysis and interpretation of the experimental data, we employed diffraction pattern simulations using the scalar and vector forms of the Debye scattering equation to determine the influence of basal plane orientation, crystalline ordering (turbostratic-graphitic), and basal plane asymmetry on the diffraction patterns. The results presented here suggest that growth of the transverse crystallites in mesophase pitch-based carbon fiber is fixed until graphitization temperatures are reached. The work completed here provides a framework for the analysis of carbon fiber and other oriented carbon-based materials via diffraction.

36 MATERIALS SCIENCE↗

Transport model comparison studies of intermediate-energy heavy-ion collisions

Transport models are the main method to obtain physics information on the nuclear equation of state and in-medium properties of particles from low to relativistic-energy heavy-ion collisions. The Transport Model Evaluation Project (TMEP) has been pursued to test the robustness of transport model predictions in reaching consistent conclusions from the same type of physical model. To this end, calculations under controlled conditions of physical input and set-up were performed with various participating codes. These included both calculations of nuclear matter in a box with periodic boundary conditions, which test separately selected ingredients of a transport code, and more realistic calculations of heavy-ion collisions. Over the years, six studies have been performed within this project. In this intermediate review, we summarize and discuss the present status of the project. We also provide condensed descriptions of the 26 participating codes, which contributed to some part of the project. These include the major codes in use today. After a compact description of the underlying transport approaches, we review the main results of the studies completed so far. They show, that in box calculations the differences between the codes can be well understood and a convergence of the results can be reached. These studies also highlight the systematic differences between the two families of transport codes, known under the names of Boltzmann–Uehling–Uhlenbeck (BUU) and Quantum Molecular Dynamics (QMD) type codes. However, when the codes were compared in full heavy-ion collisions using different physical models, as recently for pion production, they still yielded substantially different results. This calls for further comparisons of heavy-ion collisions with controlled models and of box comparisons of important ingredients, like momentum-dependent fields, which are currently underway. Furthermore, our evaluation studies often indicate improved strategies in performing transport simulations and thus can provide guidance to code developers. Results of transport simulations of heavy-ion collisions from a given code will have more significance if the code can be validated against benchmark calculations such as the ones summarized in this review.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fluctuations in Hill’s equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton’s energy into Standard Model particles. This conversion often proceeds through nonperturbative mechanisms such as parametric resonance, which is described by Hill’s equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill’s equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. Furthermore, these results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.

Cosmology↗

Fostering a Guiding Multiscale Model for the Development of Advanced MgB 2 Hydrogen Storage Materials (Final Technical Report)

Project Goal and Objective. The demand for energy and for an upgraded energy infrastructure has steadily grown, as have the needs for energy independence and alternatives to our reliance on petroleum. Hydrogen is considered the most viable fuels for wide-scale implementation in the near future as it is less-polluting, non-toxic, and has more stored energy than petroleum. It is envisioned that hydrogen can eventually become the prime energy carrier, integrating the transportation, grid, and chemical sectors in a way that improves resiliency, diversifies feedstocks, and affords new economic opportunities. A key remaining challenge is the development materials with enhanced gravimetric and volumetric hydrogen storage capacities that offer a higher performance than compressed gas. These materials would eliminate the need for large-scale compression, thereby dramatically reducing the footprint and cost of gas storage. The high gravimetric and volumetric hydrogen capacities of complex hydrides has prompted an intensive investigation of the potential of this class of materials as hydrogen storage media over the past 25 years. Among the many complex hydrides that have been explored, magnesium borohydride, Mg(BH 4 ) 2 , has been found to possess the best combination of practical thermodynamic properties. These include a gravimetric H 2 density of 14.9 wt% H 2 and thermodynamics for the dehydrogenation of Mg(BH 4 ) 2 to MgB 2 (equation 1) (ΔH° = 39 kJ/mol H 2 , ΔS = 112 J/K mol H 2 ) which lie in the narrow window required Mg(BH 4 ) 2 $\Leftrightarrow$ MgB 2 + 4 H 2 (1) for reversibility under moderate pressure and temperature. However, overcoming the extremely slow kinetics of the reversible release of hydrogen by this material in the solid state is a daunting challenge. At temperatures greater than 400 °C, the borohydride releases up to 14 wt% hydrogen giving MgB 2 . We discovered that the direct re-hydrogenation of MgB 2 to Mg(BH 4 ) 2 can be accomplished under 950 bar H 2 at 400 °C. While this demonstrated that complete reversibility can be achieved, the conditions employed are far too extreme for commercial hydrogen storage applications. More recently, we found through US DOE funded research projects (EERE HyMARC and HySCOR), that hydrogen cycling, can be accomplish at much milder conditions upon modification of the borohydride or boride. Guided by these discoveries these discoveries, the objective of this research project was to obtain key information that will enable the development of a model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 . The ultimate goal of our efforts is to attain a model of this transformation that can be utilized to accelerate development further advanced materials. This project directly follows on discoveries that were made over the course of a US DOE, EERE HyMARC project that was focused on improvement of the hydrogen cycling kinetics of modified MgB 2 . We found that that mechanical milling with graphene results the desired, pronounced kinetic enhancement. The dramatic lowering of the conditions required for the hydrogenation of MgB 2 is a significant step towards overcoming its chemical inertness allowing its development as a practical onboard hydrogen storage material. However, the exact nature of the modification(s) of MgB 2 that is responsible for its activation towards hydrogenation is completely unknown. This situation is not unique, as efforts to develop hydrogen storage materials typically have a narrow focus rather than a comprehensive approach that takes atomic level bonding and structure; molecular dynamics; long range, nano- and mesoscale-structure and their interconnection all into account. The goal of this project was the development of a comprehensive, multi-scale computational model of reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 that can be utilized for development of higher performance versions of the modified material. Development of the model requires determination of: 1) the bulk, nano-scale, and meso-scale structural changes occurring at elevated pressure following mechano-chemical modification of MgB 2 ; 2) the reaction pathway of the reversible hydrogenation of MgB 2 to Mg(BH 4 ) 2 ; 3) the effect of elevated pressure and mechano-chemical modification on the chemical reaction pathways; 4) the interactions at solid-gas interfaces; and particle surfaces; and 5) the kinetics and thermodynamic parameters associated with each step of the hydrogenation reaction pathway. This investigation required advanced techniques as preliminary, standard XRD, 11 B NMR, and FTIR analysis showed no signs of material modification. In order to gain this level of understanding of modified MgB 2 , required the teaming of a diverse group of experts and state-of-the art experimental capabilities at the University of Hawaii at Manoa (UHM) and collaborating National Laboratories: Craig Jensen , Department of Chemistry (PI and Project Director), solid state, solution, and high pressure NMR spectroscopy; solid-state synthesis; and high pressure hydrogenation (collaboration with SNL); Godwin Severa , Hawaii Natural Energy Institute (co-PI) calorimetry; infrared and Raman spectroscopy (collaboration with NREL); Dera , high pressure X-ray diffraction including in situ experiments (collaboration with ANL); Hope Ishii , Hawaii Institute of Geophysics electron microscopy investigations (collaboration with LBNL); and Joe Brown , Mechanical Engineering , material electronic structure and electric field effects.

08 HYDROGEN↗

A Semi-Empirical Correlation for the Onset of Density Wave Oscillations in a Helical Coil Steam Generator

NuScale Power has completed full length helical coil steam generator (HCSG) tests at SIET Laboratories in Piacenza, Italy. The test program included an investigation of the onset and evolution of density wave oscillations (DWO). This paper describes the mechanisms leading to the onset of DWO as inferred from tests performed in the full-height TF-2 counter-flow test facility. This research supports the concept that the onset of DWO is caused by liquid bridging inside the boiling length of the tubes. Using a separate cylinders approximation, a DWO onset correlation is obtained by performing a steady-state momentum balance on the density front created by liquid bridging. This yields a momentum balance equation in terms of a set of dimensionless superficial velocities, a corresponding dimensionless pressure drop term and an average void fraction. Envelope theory is then used to obtain a critical void fraction which is substituted into the dimensionless momentum balance to obtain a DWO onset correlation similar to the Wallis counter-current flooding correlation. A comparison of the DWO onset correlation to the TF-2 DWO experimental data is presented and shows excellent agreement when the average equilibrium vapor quality is greater than 0.2.

Density Wave Oscillations, Helical Coil Steam Gene↗