Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHASE VELOCITY”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Shock compression of silver to 300 GPa: Wave profile measurements and melting transition

Recent in-situ x-ray diffraction (XRD) measurements on laser-shock compressed Ag foils demonstrated a face-centered-cubic to body-centered-cubic transformation at ~150 GPa and melting between 172 and 197 GPa [Phys. Rev. Lett. 124, 235701 (2020)]. As a complement to the XRD work, we conducted plate impact experiments to obtain shock velocity and wave profile measurements on thicker Ag samples shock-compressed to peak stresses between 30 and 300 GPa. The shock velocity–particle velocity results were fitted very well by a linear relation over the entire stress range, providing an accurate determination of the Ag Hugoniot (locus of the stress-volume states achieved under shock compression). For peak stresses below 187 GPa and above 210 GPa—corresponding to the solid and liquid phases, respectively—the wave profiles show clean single waves. No wave profile features related to the fcc-bcc transformation at ~150 GPa were observed, implying minimal volume change for the transformation. For stresses between 187 and 210 GPa, an initial jump was followed by a time-dependent increase in the particle velocity (20–80 ns risetime) to the peak state—corresponding to the solid-liquid mixed phase response. Unlike the solid and liquid response, the mixed-phase response cannot be readily analyzed analytically. Instead, numerical simulations incorporating an accurate multiphase equation of state for Ag—not currently available—are required to analyze the wave profiles measured at 187–210 GPa stresses. Furthermore, the present work shows the potential for using wave profile measurements to examine the melting transition under shock compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Anisotropic thermal-conductivity degradation in the phase-field method accounting for crack directionality

Dynamic loads applied to metals may lead to brittle or ductile fracture depending on the imposed strain rates, material properties, and specimen geometry. With an increase in the applied velocity, brittle to ductile failure mode transition may also be observed. Furthermore, at high strain rates, ductile fracture may be preceded by shear bands, which are narrow bands of intense plastic deformation, typically accompanied by a significant rise in temperature. Reliable models are needed to predict the response of metals subject to dynamic loads. In particular, capturing the interplay between heat conduction and crack propagation using the phase-field fracture method is still an open research field. To accurately capture the heat transfer physics across crack surfaces, damage models degrading thermal-conductivity are necessary. While isotropic thermal-conductivity degradation models were proposed, they may lead to errors as they are indifferent to crack directionality. In our current work, a new anisotropic approach is proposed in which thermal-conductivity, which depends on the phase-field gradient, is degraded solely across the crack. Here, it is shown that this approach improves the near-field approximation of temperature and heat flux compared with isotropic degradation, when taking the discontinuous crack solutions as reference. Additionally, the anisotropic degradation approach is implemented in a unified dynamic fracture model and a couple of examples demonstrate the viability of this approach.

42 ENGINEERING↗
DOE OSTI2018

Axial dispersion and mixing of coolant gas within a separate-effect prismatic modular reactor

Multiphase Reactors Engineering and Applications Laboratory performed gas phase dispersion experiments in a separate-effect cold-flow experimental setup for coolant flow within heated channels of the prismatic modular reactor under accident scenario using gaseous tracer technique. The separate-effect experimental setup was designed on light of local velocity measurements obtained by using hot wire anemometry. The measurements consist of pulse-response of gas tracer that is flowing through the mimicked riser channel using air as a carrier. The dispersion of the gas phase within the separate-effect riser channel was described using one-dimensional axial dispersion model. The axial dispersion coefficient and Peclet number of the coolant gas phase and their residence time distribution within were measured. Effect of heating intensities in terms of heat fluxes on the coolant gas dispersion along riser channels were mimicked in the current study by a certain range of volumetric air flow rate ranging from 0.0015 to 0.0034 m 3 /s which corresponding to heating intensity range from 200 to 1400 W/m 2 . Results confirm a reduction in the response curve spreads is achieved by increasing the volumetric air velocity (representing heating intensity). Also, the results reveal a reduction in values of axial dispersion coefficient with increasing the air volumetric flow rate.

Said, Ibrahim↗
DOE OSTI2017

Field evaluation of remote wind sensing technologies: Shore-based and buoy mounted LIDAR systems

In developing a national energy strategy, the United States has a number of objectives, including increasing economic growth, improving environmental quality, and enhancing national energy security. Wind power contributes to these objectives through the deployment of clean, affordable and reliable domestic energy. To achieve U.S. wind generation objectives, the Wind and Water Power Program within the Department of Energy’s (DOE) Office of Energy Efficiency and Renewable Energy (EERE) instituted the U.S. Offshore Wind: Removing Market Barriers Program in FY 2011. Accurate and comprehensive information on offshore wind resource characteristics across a range of spatial and temporal scales is one market barrier that needs to be addressed through advanced research in remote sensing technologies. There is a pressing need for reliable offshore wind-speed measurements to assess the availability of the potential wind energy resource in terms of power production and to identify any frequently occurring spatial variability in the offshore wind resource that may impact the operational reliability and lifetime of wind turbines and their components and to provide a verification program to validate the “bankability” of the output of these alternative technologies for use by finance institutions for the financing of offshore wind farm construction. The application of emerging remote sensing technologies is viewed as a means to cost-effectively meet the data needs of the offshore wind industry. In particular, scanning and buoy mounted LIDAR have been proposed as a means to obtain accurate offshore wind data at multiple locations without the high cost and regulatory hurdles associated with the construction of offshore meteorological towers. However; before these remote sensing technologies can be accepted the validity of the measured data must be evaluated to ensure their accuracy. The proposed research will establish a unique coastal ocean test-bed in the Mid-Atlantic for the evaluation of LIDAR-based wind measurement systems to validate the accuracy of remotely measured wind data in marine applications. Specifically, the test-bed will be utilized to systematically evaluate the capability of emerging scanning LIDAR and buoy mounted vertically profiling LIDAR by: (1) Evaluating a fixed scanning LIDAR against land-based 50 and 60 meter high meteorological masts fitted with research quality cup-vane and/or sonic anemometers; (2) Evaluating a buoy mounted vertically profiling LIDAR fixed on land and floating in a sheltered bay against a co-located 60 meter high meteorological mast fitted with a research quality cup-vane and/or sonic anemometers and the fixed scanning LIDAR; and (3) Offshore field evaluation of both LIDAR platforms through a comparison of the fixed scanning LIDAR data and data obtained by the buoy mounted LIDAR located 10 miles offshore. The proposed research will systematically validate Light Detection and Ranging (LIDAR) based wind measurement systems and assess the temporal and spatial variability of the offshore wind resource in the Mid-Atlantic east of New Jersey. The goal of the proposed project is to address the technical and commercial challenges of the offshore wind energy industry by validating and assessing cost-effective, over ocean wind resource characterization technologies. The objective is to systematically evaluate the capability of both scanning and vertically profiling LIDARs to accurately measure 3D wind fields through comparison with fixed met masts and intercomparison among LIDAR platforms. Once validated, data collected by both buoy mounted vertically profiling LIDARs and shore-based, pulsed horizontally scanning LIDARs can be used to accurately assess offshore wind resources and to quantify the spatial and temporal variability in the offshore wind fields. One of the fundamental research questions to be addressed in phase 1 is the assessment of various measurement and data processing schemes to retrieve accurate wind vectors in the marine environment over large sampling ranges (10 to 12 km) and varying atmospheric aerosol levels. Atmospheric conditions and aerosol content within the coastal ocean region of the Mid-Atlantic seaboard of the US can vary significantly over short time periods in response to frontal passages and extratropical and tropical low pressure system passage offshore of the coast. Since aerosols provide the scattering medium for the determination of LIDAR Doppler shifts in the atmosphere the accuracy and range of LIDAR derived velocity measurements as a function of variation in aerosol content in the marine environment is a key research question to be addressed. In phase 1, it is desired to capture as much variation in atmospheric conditions and aerosol content as possible. To this end, collocated measurements of LIDAR and standard anemometer wind fields will be captured by the project PIs over all four seasons and during specific events (e.g., coastal low pressure system passage) in year 1. Additionally, since the meteorological masts are permanent structures, additional events can be captured over the three year duration of the field research project. All research instruments are owned by Fishermen’s Energy and made available to the PIs though a lease agreement as part of the DOE grant. Energy Fishermen’s Energy will be responsible for the operation and maintenance of the scanning LIDAR and met mast anemometers. On a daily basis, environmental data and systems performance indicators will be transmitted from each measurement station to the Fishermen’s project team consisting of both in-house personnel and equipment manufacturer engineers. Data sets include compiled LIDAR files as well as data sets from ancillary sensors. Diagnostic parameters to be monitored include standard deviations of measured values, battery levels and charging systems output, and the operational status. Once data have been confirmed as complete and reliable, files will be transferred to the Garrad Hassan (a subcontractor to Fishermen’s Energy) for incorporation in to the validation database, which is accessible to other scientific team members. Data collection times and durations will be determined by the PI and Co-PIs in consultation with instrument engineers to ensure the capture of data representative of the expected range of mid-Atlantic atmospheric conditions (e.g., temperature, moisture, coastal low pressure systems, tropical systems, rain, snow, fog). The collection and processing of the data is a function of site specific measurement requirements (Kelley et.al. 2007; Hannon et.al. 2008). To determine the optimal profiles of wind speed and direction from the LIDAR radial velocities as a function of azimuth angle, rigorous estimates of the bias and random error of each radial velocity estimate are required. Lockheed Martin Coherent Technologies, Inc., under contract with Fishermen’s Energy, will provide analyses of raw and processed data using various scan patterns to determine optimal performance settings for the pulsed scanning LIDAR. Once optimized, appropriate processing and analyses techniques will be evaluated by Garrad Hassan for use in validating the accuracy of the LIDAR wind field measurements against the standard anemometer measurements from the meteorological masts. The most attractive capability of the scanning LIDAR is the ability to provide high spatial resolution observations in a three-dimensional volume which provides superior statistical accuracy due to the large number of samples obtained. Each radial scan provides measurements in 100 range gates over a distance of 10 to 12 km at an update rate of 5 to 10 Hz and rotation of 2.5° per second. Each rotation at a fixed azimuth requires 2.4 minutes. Depending on the number of azimuths desired a complete scan can take up to 10 minutes or longer to complete. Once collected the radial velocities are processed to produce vector wind velocity estimates based on a set of data distributed in angle and range around points of interest, typically a standard grid within the radial wind map. To calculate wind vectors over a limited spatial area of interest for the comparison of data with other measurement platforms a localized least-squares approach has been applied by Hannon et.al. (2008) and a Variation Assimilation (VAR) processing technique has been applied by Chan and al Assimilation (VAR) processing technique has been applied by Chan and Shao (2006). Additionally, Kelley et.al. (2007) applied a “stare” technique that fixed a scanning LIDAR in both azimuth and elevation angles to measure over collocated sampling volumes of the LIDAR and a 3D sonic anemometer mounted to a fixed mast. Although the stare technique is limited to wind directions aligned with the sampling radial of the LIDAR, it does provide a direct comparison of sampling volumes. Each of the processing techniques described above (and possibly others) will be evaluated to determine the validity of the LIDAR derived wind fields in the marine environment. Numerical methods such as linear regression and comparison of probability density functions of wind fields measured by each instrument platform will be used to assess the processing techniques. Linear regression has the advantage of directly evaluating corresponding pairs of wind data measured by each instrument and can lend insight into deviations and bias between instruments as a function of wind speed. Assessment of the coherency between the probability density function of the wind measured by each instrument provides insight into processes that may not be accurately resolved by each instrument at specific frequencies. Once processed and assessed the most appropriate technique will be utilized to provide valid wind measurements from the pulsed scanning LIDAR. A detailed analysis of the measurement data from the LIDAR and the three meteorological towers and comparison of the coastal wind characteristics from the different systems will be performed by the CO_PI at NREL. This analysis will evaluate how the WT LIDAR performance and measurement of the wind characteristics vary with distance from the LIDAR and by atmospheric conditions; using the tower measurements at different distances (4.8, 9.6, and 19 km) from the LIDAR as a reference. The comparative analysis will include, to the extent possible, evaluation of parameters such as wind speed and direction distributions, wind shear, turbulence intensity and their variations by atmospheric conditions, month or season, and time of day.

17 WIND ENERGY↗
DOE OSTI2015

Dynamic density functional theory with hydrodynamic interactions: Theoretical development and application in the study of phase separation in gas-liquid systems

Building on recent developments in dynamic density functional theory, we have developed a version of the theory that includes hydrodynamic interactions. This is achieved by combining the continuity and momentum equations eliminating velocity fields, so the resulting model equation contains only terms related to the fluid density and its time and spatial derivatives. The new model satisfies simultaneously continuity and momentum equations under the assumptions of constant dynamic or kinematic viscosity and small velocities and/or density gradients. We present applications of the theory to spinodal decomposition of subcritical temperatures for one-dimensional and three-dimensional density perturbations for both a van der Waals fluid and for a lattice gas model in mean field theory. In the latter case, the theory provides a hydrodynamic extension to the recently studied dynamic mean field theory. We find that the theory correctly describes the transition from diffusive phase separation at short times to hydrodynamic behaviour at long times.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Phase shift method to estimate solids circulation rate in circulating fluidized beds

While solids circulation rate is a critical design and control parameter in circulating fluidized bed (CFB) reactor systems, there are no available techniques to measure it directly at conditions of industrial interest. Cold flow tests have been conducted at NETL in an industrial scale CFB unit where the solids flow has been the topic of research in order to develop an independent method which could be applied to CFBs operating under the erosive and corrosive high temperatures and pressures of a coal fired boiler or gasifier. The dynamic responses of the CFB loop to modest modulated aeration flows in the return leg or standpipe were imposed to establish a periodic response in the unit without causing upset in the process performance. The resulting periodic behavior could then be analyzed with a dynamic model and the average solids circulation rate could be established. This method was applied to the CFB unit operated under a wide range of operating conditions including fast fluidization, core annular flow, dilute and dense transport, and dense suspension upflow. In addition, the system was operated in both low and high total solids inventories to explore the influence of inventory limiting cases on the estimated results. The technique was able to estimate the solids circulation rate for all transport circulating fluidized beds when operating above upper transport velocity, U{sub tr2}. For CFB operating in the fast fluidized bed regime (i.e., U{sub g}< U{sub tr2}), the phase shift technique was not successful. The riser pressure drop becomes independent of the solids circulation rate and the mass flow rate out of the riser does not show modulated behavior even when the riser pressure drop does.

42 ENGINEERING↗
DOE OSTI2016

Energy spread minimization in a cascaded laser wakefield accelerator via velocity bunching

We propose a scheme to minimize the energy spread of an electron beam (e-beam) in a cascaded laser wakefield accelerator to the one-thousandth-level by inserting a stage to compress its longitudinal spatial distribution. In this scheme, three-segment plasma stages are designed for electron injection, e-beam length compression, and e-beam acceleration, respectively. The trapped e-beam in the injection stage is transferred to the zero-phase region at the center of one wakefield period in the compression stage where the length of the e-beam can be greatly shortened owing to the velocity bunching. After being seeded into the third stage for acceleration, the e-beam can be accelerated to a much higher energy before its energy chirp is compensated owing to the shortened e-beam length. A one-dimensional theory and two-dimensional particle-in-cell simulations have demonstrated this scheme and an e-beam with 0.2% rms energy spread and low transverse emittance could be generated without loss of charge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electrodynamics of Thin Sheets of Twisted Material

We construct a minimal theory describing the optical activity of a thin sheet of a twisted material, the simplest example of which is twisted bilayer graphene. We introduce the notion of "twisted electrical conductivity", which parametrizes the parity-odd response of a thin film to a perpendicularly falling electromagnetic waves with wavelength larger than the thickness of the sheet. Here, we show that the low-frequency Faraday rotation angle has different behaviors in different phases. For an insulator, the Faraday angle behaves as ω 2 at low frequencies, with the coefficient being determined by the linear relationship between a component of the electric quadrupole moment and the external electric field. For superconductors, the Faraday rotation angle is constant when the frequency of the incoming EM waves is below the superconducting gap and is determined by the coefficient of the Lifshitz invariant in the Ginzburg-Landau functional describing the superconducting state. In the metallic state, we show that the twisted conductivity is proportional to the "magnetic helicity" (scalar product of the velocity and the magnetic moment) of the quasiparticle, averaged around the Fermi surface. The theory is general and is applicable to strongly correlated phases.

Son, Dam Thanh↗
DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗
DOE OSTI2016

Pore-scale multiphase flow modeling and imaging of CO 2 exsolution in Sandstone

This study utilizes synchrotron X-ray micro-tomography and pore scale modeling to investigate the process of gas exsolution and how it affects non-wetting phase relative permeability. Exsolved gas distributions are measured on Domengine and Boise sandstone samples using synchrotron X-ray micro-tomography. Observed gas phase distributions are compared to a new model that simulates the growth and distribution of exsolved gas phase at the pore-scale. Water relative permeability curves are calculated using a Stokes flow simulator with modeled and observed gas distributions, under various conditions, such as rock geometry, and pressure depletion rates. By comparing the actual bubble distributions with modeled distributions, we conclude that exsolved gas is more likely to form and accumulate at locations with higher water velocities. This suggests that convective delivery of CO2 to the gas bubble is a primary mechanism for bubble growth, as compared to diffusive transport through the aqueous phase. For carbonated brine flowing up a fault at half a meter per day, with 5% exsolved gas, the water relative permeability is estimated to be 0.6~0.8 for various sandstones. The reduction of water mobility reduces upward brine migration when even a small amount of exsolution occurs.

58 GEOSCIENCES↗
DOE OSTI2024

Interactions of laser-driven tin ejecta microjets over phase transition boundaries

Ejecta microjets offer an experimental methodology to study high-speed particle laden-flow interactions, as microjets consist of millions of particulates traveling at velocities of several kilometers per second and are easily generated by most common shock drives. Previous experiments on the OMEGA Extended Performance laser found that collisions between two counter-propagating laser-driven tin ejecta microjets varied as a function of drive pressure; jets generated near shock pressures of 10 GPa passed through each other without interacting, whereas jets generated at shock pressures of over 100 GPa interacted strongly, forming a cloud around the center interaction point. In this paper, we present a more systematic scan of tin ejecta microjet collisions over intermediate pressure regimes to identify how and at what shock pressure interaction behavior onsets. Radiographs of interacting microjets at five different laser drive energies qualitatively demonstrate that interaction behavior onsets slowly as a function of laser drive energy. Quantitative mass and density metrics from each radiograph provide trends on jet characteristics and collisional mass dispersion. It is observed that jetting mass, jet densities, and mass dispersion increase with increasing drive pressures and that the increased jet density at the higher drive energies may account for the increased mass dispersion. This work provides an important step in the understanding of tin ejecta microjet collisions and paves the way for future studies on the physics dominating high-speed particle-laden flow interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A 17 year climatology of the macrophysical properties of convection in Darwin

Abstract. The validation of convective processes in global climate models (GCMs) could benefit from the use of large datasets that provide long-term climatologies of the spatial statistics of convection. To that regard, echo top heights (ETHs), convective areas, and frequencies of mesoscale convective systems (MCSs) from 17 years of data from a C-band polarization (CPOL) radar are analyzed in varying phases of the Madden–Julian Oscillation (MJO) and northern Australian monsoon in order to provide ample validation statistics for GCM validation. The ETHs calculated using velocity texture and reflectivity provide similar results, showing that the ETHs are insensitive to various techniques that can be used. Retrieved ETHs are correlated with those from cloud top heights retrieved by Multifunctional Transport Satellites (MTSATs), showing that the ETHs capture the relative variability in cloud top heights over seasonal scales. Bimodal distributions of ETH, likely attributable to the cumulus congestus clouds and mature stages of convection, are more commonly observed when the active phase of the MJO is over Australia due to greater mid-level moisture during the active phase of the MJO. The presence of a convectively stable layer at around 5 km altitude over Darwin inhibiting convection past this level can explain the position of the modes at around 2–4 km and 7–9 km. Larger cells were observed during break conditions compared to monsoon conditions, but only during the inactive phase of the MJO. The spatial distributions show that Hector, a deep convective system that occurs almost daily during the wet season over the Tiwi Islands, and sea-breeze convergence lines are likely more common in break conditions. Oceanic MCSs are more common during the night over Darwin. Convective areas were generally smaller and MCSs more frequent during active monsoon conditions. In general, the MJO is a greater control on the ETHs in the deep convective mode observed over Darwin, with higher distributions of ETH when the MJO is active over Darwin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

RADIAL VELOCITY OBSERVATIONS AND LIGHT CURVE NOISE MODELING CONFIRM THAT KEPLER-91b IS A GIANT PLANET ORBITING A GIANT STAR

Kepler-91b is a rare example of a transiting hot Jupiter around a red giant star, providing the possibility to study the formation and composition of hot Jupiters under different conditions compared to main-sequence stars. However, the planetary nature of Kepler-91b, which was confirmed using phase-curve variations by Lillo-Box et al., was recently called into question based on a re-analysis of Kepler data. We have obtained ground-based radial velocity observations from the Hobby-Eberly Telescope and unambiguously confirm the planetary nature of Kepler-91b by simultaneously modeling the Kepler and radial velocity data. The star exhibits temporally correlated noise due to stellar granulation which we model as a Gaussian Process. We hypothesize that it is this noise component that led previous studies to suspect Kepler-91b to be a false positive. Our work confirms the conclusions presented by Lillo-Box et al. that Kepler-91b is a 0.73 ± 0.13 M {sub Jup} planet orbiting a red giant star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Physics based model for online fault detection in autonomous cryogenic loading system

We report the progress in the development of the chilldown model for a rapid cryogenic loading system developed at NASA-Kennedy Space Center. The nontrivial characteristic feature of the analyzed chilldown regime is its active control by dump valves. The two-phase flow model of the chilldown is approximated as one-dimensional homogeneous fluid flow with no slip condition for the interphase velocity. The model is built using commercial SINDA/FLUINT software. The results of numerical predictions are in good agreement with the experimental time traces. The obtained results pave the way to the application of the SINDA/FLUINT model as a verification tool for the design and algorithm development required for autonomous loading operation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Photoelectron spectroscopy of anthracene and fluoranthene radical anions

Here, we report the slow electron velocity map imaging spectroscopy of cryogenically cooled anthracene and fluoranthene radical anions, two similarly sized polycyclic aromatic hydrocarbon molecules. The results allow us to examine the lowest energy singlet and triplet states in the neutral molecules on equal footing from the anionic ground state. The analysis of the experimental spectra is aided by harmonic calculations and Franck-Condon simulations, which generally show good agreement with experimental values and spectra. The electron affinity of fluoranthene is measured to be 0.757(2) eV, which is larger than that of anthracene at 0.532(3) eV. The lowest energy triplet state in anthracene is observed at 1.872(3) eV above the singlet ground state, while that of fluoranthene is observed at 2.321(2) eV above its singlet ground state. Comparisons of experimental and calculated spectra show that in addition to the Franck-Condon active modes, there is a clear presence of vibrational modes that gain intensity via vibronic coupling in both the singlet and triplet states in both molecules. In addition, the triplet state generally exhibits increased vibronic coupling compared to the singlet state, with the fluoranthene triplet state exhibiting evidence of distortion from C 2v symmetry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise

SUMMARY We present a new high-resolution 3-D shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broad-band seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain empirical Green functions (EGFs) for 2016 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh-wave group velocities in the period range 7–30 s and associated uncertainties. The tomographic inversion is performed in two steps: First, we determine group-velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3-D shear wave velocity model. The final 3-D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of ∼50 km. In the upper and middle crusts, the Vs anomaly pattern matches the tectonic units of the Variscan Massif and Alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate- and high-velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely, the Manteigas–Vilariça–Bragança fault and the Porto–Tomar–Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3-D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.

Silveira, Graça (ORCID:0000000221102554)↗
DOE OSTI2012

Numerical calculation of neoclassical distribution functions and current profiles in low collisionality, axisymmetric plasmas

A new code, the Neoclassical Ion-Electron Solver (NIES), has been written to solve for stationary, axisymmetric distribution functions (f) in the conventional banana regime for both ions and electrons using a set of drift-kinetic equations (DKEs) with linearized Fokker-Planck-Landau collision operators. Solvability conditions on the DKEs determine the relevant non-adiabatic pieces of f (called h). We work in a 4D phase space in which {psi} defines a flux surface, {theta} is the poloidal angle, v is the magnitude of the velocity referenced to the mean flow velocity, and {lambda} is the dimensionless magnetic moment parameter. We expand h in finite elements in both v and {lambda}. The Rosenbluth potentials, {Phi} and {Psi}, which define the integral part of the collision operator, are expanded in Legendre series in cos{chi}, where {chi} is the pitch angle, Fourier series in cos{theta}, and finite elements in v. At each {psi}, we solve a block tridiagonal system for h{sub i} (independent of f{sub e}), then solve another block tridiagonal system for h{sub e} (dependent on f{sub i}). We demonstrate that such a formulation can be accurately and efficiently solved. NIES is coupled to the MHD equilibrium code JSOLVER [J. DeLucia et al., J. Comput. Phys. 37, 183-204 (1980)] allowing us to work with realistic magnetic geometries. The bootstrap current is calculated as a simple moment of the distribution function. Results are benchmarked against the Sauter analytic formulas and can be used as a kinetic closure for an MHD code (e.g., M3D-C{sup 1}[S. C. Jardin et al., Comput. Sci. Discovery 5, 014002 (2012)]).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Improvements in the reactor core isolation cooling (RCIC) pump model

As part of the efforts to understand the unexpected “self-regulating” mode of the RCIC (Reactor Core Isolation Cooling) systems that was observed during the Fukushima accidents and to extend BWR RCIC and PWR AFW (Auxiliary Feed Water) operational range and flexibility, mechanistic models for the Terry turbine are being developed and tested with the RELAP-7 code to simulate the RCIC system. In the previous work, we developed a set of analytical models for the normal working conditions, as well as the two-phase off-design conditions, of the RCIC system, based on Sandia’s original work. An under-expanded jet model was developed to replace the CFD-based reduced-order model in the Sandia work to obtain the velocity and thermodynamic conditions of the turbine stator inlet. The models include both an adiabatic expansion process inside the nozzle and a free expansion process outside of the nozzle to ambient pressure. For the two-phase off-design conditions, two well-established choking models were used - the Isentropic Homogenous Equilibrium Model (IHEM) and Moody’s model. The new twophase Terry turbine model used the choking models to calculate the mass flow rate, the critical pressure at the nozzle throat, and steam quality. The two-phase nozzle expansion model was also developed by considering the vapor phase with a similar model for the single-phase case with the assumption that the liquid phase would slip along the wall with a much slower speed and will not contribute to the impulse on the rotor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗