Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coherent errors”

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.

431 records · Page 24

The Origin of the X-Ray and Ultraviolet Emission in NGC 7469

We present a spectral analysis of a approx 30 day, near-continuous observation of the Seyfert 1 galaxy NGC 7469 with RXTE. Daily integrations show strong spectral changes during the observation. Our main result is that we find the X-ray spectral index to be correlated with the UV flux. Furthermore, the broadband X-ray photon flux is also correlated with the UV continuum. These correlations point toward a model in which the X-rays originate via thermal Comptonization of UV seed photons. Furthermore, the UV is also correlated with the extrapolation of the X-ray power law into the soft X-ray/EUV region. Our data analysis therefore reopens the possibility that the UV photons and their variability arise from reprocessing, as long as the primary source of heating is photoelectric absorption in the re-processor, rather than Compton down-scattering. A coherent picture of the X-ray/UV variability can therefore be constructed whereby absorption and reprocessing of EUV/soft X-rays in a standard accretion disk produce a variable seed photon distribution, which is in turn up-scattered into the X-ray band. We also find a significant correlation between the 2-10 keV flux and the 6.4 keV iron K-alpha line, suggesting that at least some portion of the line originates within approx. 1 light day of the X-ray continuum source. Neither the power-law photon index nor the Compton reflection component are correlated with the 2-10 keV flux. The latter is not correlated with the iron K-alpha line flux either. We do find an apparent correlation between the X-ray spectral index and the strength of the Compton reflection component. In an Appendix we show, however, that this can be produced by a combination of statistical and systematic errors. We conclude that the apparent variations in the Compton reflection component may be an artifact of these effects.

Nandra, K.↗

Application of Tracking and Data Relay Satellite (TDRS) Differenced One-Way Doppler (DOWD) Tracking Data for Orbit Determination and Station Acquisition Support of User Spacecraft Without TDRS Compatible Transponders

Many spacecraft are launched today with only an omni-directional (omni) antenna and do not have an onboard Tracking and Data Relay Satellite (TDRS) transponder that is capable of coherently returning a carrier signal through TDRS. Therefore, other means of tracking need to be explored and used to adequately acquire the spacecraft. Differenced One-Way Doppler (DOWD) tracking data are very useful in eliminating the problems associated with the instability of the onboard oscillators when using strictly one-way Doppler data. This paper investigates the TDRS DOWD tracking data received by the Goddard Space Flight Center (GSFC) Flight Dynamics Facility (FDF) during the launch and early orbit phases for the the Interplanetary Physics Laboratory (WIND) and the National Oceanographic and Atmospheric Administration (NOAA)-J missions. In particular FDF personnel performed an investigation of the data residuals and made an assessment of the acquisition capabilities of DOWD-based solutions. Comparisons of DOWD solutions with existing data types were performed and analyzed in this study. The evaluation also includes atmospheric editing of the DOWD data and a study of the feasibility of solving for Doppler biases in an attempt to minimize error. Furthermore, by comparing the results from WIND and NOAA-J, an attempt is made to show the limitations involved in using DOWD data for the two different mission profiles. The techniques discussed in this paper benefit the launches of spacecraft that do not have TDRS transponders on board, particularly those launched into a low Earth orbit. The use of DOWD data is a valuable asset to missions which do not have a stable local oscillator to enable high-quality solutions from the one-way/return-link Doppler tracking data.

Olszewski, A. D., Jr.↗

Satellite Data Sets in the Polar Regions

We have generated about two decades of consistently derived geophysical parameters in the polar regions. The key parameters are sea ice concentration, surface temperature, albedo, and cloud cover statistics. Sea ice concentrations were derived from the Scanning Multichannel Microwave Radiometer (SMMR) data and the Special Scanning Cl Microwave Imager (SSM/I) data from several platforms using the enhanced Bootstrap Algorithm for the period 1978 through 1999. The new algorithm reduces the errors associated with spatial and temporal variations in the emissivity and surface temperatures of sea ice. Also, bad data at ocean/land interfaces are identified and deleted in an unsupervised manner. Surface ice temperature, albedo and cloud cover statistics are derived simultaneously from the Advanced Very High Resolution Radiometer (AVHRR) data from 1981 through 1999 and mapped at a higher resolution but the same format as the ice concentration data. The technique makes use these co-registered ice concentration maps to enable cloud masking to be done separately for open ocean, sea ice and land areas. The effect of inversion is minimized by taking into consideration the expected changes in the effect of inversion with altitude, especially in the Antarctic. A technique for ice type regional classification has also been developed using multichannel cluster analysis and a neural network. This provide a means to identify large areas of thin ice, first year ice, and older ice types. The data sets have been shown to be coherent with each other and provide a powerful tool for in depth studies of the currently changing Arctic and Antarctic environment.

Comiso, Josefino C.↗

In situ midcircuit qubit measurement and reset in a single-species trapped-ion quantum computing system

We implement in situ midcircuit measurement and reset (MCMR) operations on a full-scale trapped-ion quantum computing system by using metastable qubit states in 171 Yb + ions. We compare two methods for isolating data qubits from measured qubits: one shelves the data qubits into the metastable state and the other drives the measured qubit to the metastable state without disturbing the other qubits. We experimentally demonstrate both methods on a crystal of two 171 Yb + ions using both the 𝑆 1/2 ground-state hyperfine clock qubit and the 𝑆 1/2 −𝐷 3/2 optical qubit. These MCMR methods result in errors on the data qubit of about 2% without degrading the measurement fidelity. With straightforward reductions in laser noise, these errors can be suppressed to less than 0.1%. The demonstrated methods allow MCMR to be performed in a single-species ion chain without shuttling or additional qubit-addressing optics, greatly simplifying the system architecture and allowing straightforward integration with existing trapped-ion quantum computers.

coherent control↗

Planar Doppler Velocimetry for Large-Scale Wind Tunnel Testing

Recently, Planar Doppler Velocimetry (PDV) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike other optical means of measuring flow velocities, PDV is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDV have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDV in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDV capabilities for applications in large-scale wind tunnels that are intended primarily for the production testing of subsonic aircraft. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes and separations. To accomplish the objectives indicated, the PDV measurement process is first modeled and its performance evaluated computationally. The noise sources considered include those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with coherent illumination from pulsed lasers. The signal noise estimates are incorporated into the PDV signal analysis process and combined with computed scattering signals using a Mie scattering theory for polydisperse smoke particles. The relevant parameters incorporate a range of practical aerodynamic test conditions and facility sizes. The results define the optimum instrument configurations, show that the expected signal levels from a practical PDV system are sufficiently large to allow its useful application in large facilities, and show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing. Experimental studies using several experimental bench-top setups are then described that validate the physics of the PDV model and to calibrate its computed results. The validated model allows estimates of the uncertainties of PDV measurements and a complete definition of the PDV capabilities to be made with sufficient confidence to decide the viability of PDV for large-scale wind tunnel applications.

McKenzie, Robert L.↗

Navigation Doppler Lidar Signal Processing Architecture and Algorithms

A coherent Doppler lidar has been developed to address the need for a high-performance, compact, and cost-effective velocity and altitude sensor onboard landing vehicles. One of the core subsystems within NDL is its signal processing architecture. The primary function of this subsystem is to digitize the analog output of the homodyned signal received from the photo-diode associated with each telescope, process this signal in the frequency domain, and then turn this information into velocity and range measurements. A secondary function of this subsystem is to setup and control the internal components within NDL and to enable the instrument to interface with external systems by receiving commands, sending telemetry, and synchronizing its measurements with external clocks and pulses. This paper will give an overview of the signal processing architecture within NDL. This will include a discussion on its fundamental processing components and the core signal processing algorithms. Command, telemetry and internal component controls will also be discussed as will how the signal processing architecture has also evolved over the years of development. We will discuss provisions for robust operation including Triple-Mode Redundancy (TMR) within NDL’s Field Programmable Gate Array (FPGA) and storage elements, a fault-tolerant boot process, and Error Detection and Correction (EDAC) on the Static Random-Access Memories (SRAMs). Finally, we conclude this paper with how this architecture could continue to evolve.

Doppler Lidar↗

Measurement Via Optical Near-Nulling and Subaperture Stitching

A subaperture stitching interferometer system provides near-nulling of a subaperture wavefront reflected from an object of interest over a portion of a surface of the object. A variable optical element located in the radiation path adjustably provides near-nulling to facilitate stitching of subaperture interferograms, creating an interferogram representative of the entire surface of interest. This enables testing of aspheric surfaces without null optics customized for each surface prescription. The surface shapes of objects such as lenses and other precision components are often measured with interferometry. However, interferometers have a limited capture range, and thus the test wavefront cannot be too different from the reference or the interference cannot be analyzed. Furthermore, the performance of the interferometer is usually best when the test and reference wavefronts are nearly identical (referred to as a null condition). Thus, it is necessary when performing such measurements to correct for known variations in shape to ensure that unintended variations are within the capture range of the interferometer and accurately measured. This invention is a system for nearnulling within a subaperture stitching interferometer, although in principle, the concept can be employed by wavefront measuring gauges other than interferometers. The system employs a light source for providing coherent radiation of a subaperture extent. An object of interest is placed to modify the radiation (e.g., to reflect or pass the radiation), and a variable optical element is located to interact with, and nearly null, the affected radiation. A detector or imaging device is situated to obtain interference patterns in the modified radiation. Multiple subaperture interferograms are taken and are stitched, or joined, to provide an interferogram representative of the entire surface of the object of interest. The primary aspect of the invention is the use of adjustable corrective optics in the context of subaperture stitching near-nulling interferometry, wherein a complex surface is analyzed via multiple, separate, overlapping interferograms. For complex surfaces, the problem of managing the identification and placement of corrective optics becomes even more pronounced, to the extent that in most cases the null corrector optics are specific to the particular asphere prescription and no others (i.e. another asphere requires completely different null correction optics). In principle, the near-nulling technique does not require subaperture stitching at all. Building a near-null system that is practically useful relies on two key features: simplicity and universality. If the system is too complex, it will be difficult to calibrate and model its manufacturing errors, rendering it useless as a precision metrology tool and/or prohibitively expensive. If the system is not applicable to a wide range of test parts, then it does not provide significant value over conventional null-correction technology. Subaperture stitching enables simpler and more universal near-null systems to be effective, because a fraction of a surface is necessarily less complex than the whole surface (excepting the extreme case of a fractal surface description). The technique of near-nulling can significantly enhance aspheric subaperture stitching capability by allowing the interferometer to capture a wider range of aspheres. More over, subaperture stitching is essential to a truly effective near-nulling system, since looking at a fraction of the surface keeps the wavefront complexity within the capability of a relatively simple nearnull apparatus. Furthermore, by reducing the subaperture size, the complexity of the measured wavefront can be reduced until it is within the capability of the near-null design.

Forbes, Greg↗

Measurement Capabilities of Planar Doppler Velocimetry in Large-Scale Wind Tunnels

Over the past few years, Planar Doppler Velocimetry (PDT) has been shown by several laboratories to offer an attractive means for measuring three-dimensional velocity vectors everywhere in a light sheet placed in a flow. Unlike some other optical means of measuring flow velocities, PDT is particularly attractive for use in large wind tunnels where distances to the sample region may be several meters, because it does not require the spatial resolution and tracking of individual scattering particles or the alignment of crossed beams at large distances. To date, demonstrations of PDT (also called Doppler Global Velocimetry by some authors) have been made either in low speed flows without quantitative comparison to other measurements, or in supersonic flows where the Doppler shift is large and its measurement is relatively insensitive to instrumental errors. Moreover, most reported applications have relied on the use of continuous-wave lasers, which limit the measurement to time-averaged velocity fields. This work summarizes the results of two previous studies of PDT in which the use of pulsed lasers to obtain instantaneous velocity vector fields is evaluated. The objective has been to quantitatively define and demonstrate PDT capabilities for applications in large-scale wind tunnels that are intended primarily for the testing of rotorcraft and subsonic aircraft at speeds typically less than 100 m/s. For such applications, the adequate resolution of low-speed flow fields requires accurate measurements of small Doppler shifts that are obtained at distances of several meters from the sample region and with a field of view that is sufficient to encompass the entire region of interest. The use of pulsed lasers provides the unique capability to obtain not only time-averaged fields, but also their statistical fluctuation amplitudes and the spatial excursions of unsteady flow regions such as wakes, separations, and rotor-tip vortices. To accomplish the objectives of these studies, the PDT measurement process was first modeled and its performance evaluated computationally. The noise sources considered included those related to the optical and electronic properties of Charge-Coupled Device (CCD) arrays and to speckle effects associated with the coherent illumination of aerosols from pulsed lasers. The signal noise estimates were incorporated into the PDT signal analysis process and combined with computed scattering signals using a Mae scattering theory for a distributed range of aerosol particle sizes. The results are used to define the necessary instrument configuration, to show that the expected signal levels from a practical PDT system are sufficiently large to allow its useful application in large facilities, and to show that the expected velocity measurement uncertainties are small compared to the mean velocities of interest for most subsonic, large-scale wind tunnel testing, Experimental studies using several bench-top setups are then described that validate the physics of the PDT model and demonstrate the ability to obtain accurate PDT measurements using procedures that are compatible with large-scale wind tunnel operations. The PDT measurement capabilities found in this study lead to the conclusion that PDT offers significant advantages compared to other means of measuring velocity fields in large-scale wind tunnels, including those operating at low speeds.

McKenzie, Robert L.↗

Sensitive detection of structural dynamics using a statistical framework for comparative crystallography

Chemical and conformational changes are crucial to protein function and its pharmacological control. X-ray crystallography can reveal these changes in atomic detail, but standard analysis methods, which refine separate datasets, often overlook differences that are subtle or arise in only a subset of molecules. Direct comparison of crystallographic datasets is, in principle, more powerful, but systematic errors (“scales”) often mask changes in the crystallographic observables (“structure factors”). Machine learning algorithms that jointly estimate scales and structure factors can address this limitation. Here, we augment this approach with multivariate, structured priors derived from crystallographic theory, implemented in the variational deep learning framework Careless. Doing so strongly improves the detection of protein dynamics, element-specific anomalous signals, and the binding of drug candidates, offering a robust approach to comparative crystallography and, potentially, to detection of protein dynamics by other structure determination methods.

Hekstra, Doeke R. [Harvard Univ., Cambridge, MA (U↗

Cosmic axions revealed via amplified modulation of the ellipticity of a laser

We propose a new axion dark matter detection strategy that employs optical readout of laser beam ellipticity modulations caused by axion-induced electric fields in a microwave cavity, using electro-optic (EO) crystals, enhanced by externally injected radio-frequency (rf) power. Building upon the variance-based probing method [Phys. Rev. D 107, 103005 (2023).], we extend this concept to the optical domain: A weak probe laser interacts with an EO crystal coupled to the resonant microwave cavity field at cryogenic temperatures, and the axion-induced electric field is revealed through induced ellipticity. The injected rf signal coherently interferes with that of the axion field, amplifying the optical response and significantly improving sensitivity. While our EO-based method employs a Fabry-Pérot resonator, we do not require Michelson interferometers. Our method, hence, enables compact, high-frequency axion searches, across the 0.5–50 GHz range. Operating at cryogenic temperatures not only suppresses thermal backgrounds, but, critically, allows the probing method to mitigate the quantum noise. This approach offers a scalable path forward for axion detection over the ∼(few−200) μ⁢eV mass range—covering the preferred parameter space for the postinflationary Peccei-Quinn axion dark matter—using compact, tunable systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Absolute Position of Targets Measured Through a Chamber Window Using Lidar Metrology Systems

Lidar is a useful tool for taking metrology measurements without the need for physical contact with the parts under test. Lidar instruments are aimed at a target using azimuth and elevation stages, then focus a beam of coherent, frequency modulated laser energy onto the target, such as the surface of a mechanical structure. Energy from the reflected beam is mixed with an optical reference signal that travels in a fiber path internal to the instrument, and the range to the target is calculated based on the difference in the frequency of the returned and reference signals. In cases when the parts are in extreme environments, additional steps need to be taken to separate the operator and lidar from that environment. A model has been developed that accurately reduces the lidar data to an absolute position and accounts for the three media in the testbed air, fused silica, and vacuum but the approach can be adapted for any environment or material. The accuracy of laser metrology measurements depends upon knowing the parameters of the media through which the measurement beam travels. Under normal conditions, this means knowledge of the temperature, pressure, and humidity of the air in the measurement volume. In the past, chamber windows have been used to separate the measuring device from the extreme environment within the chamber and still permit optical measurement, but, so far, only relative changes have been diagnosed. The ability to make accurate measurements through a window presents a challenge as there are a number of factors to consider. In the case of the lidar, the window will increase the time-of-flight of the laser beam causing a ranging error, and refract the direction of the beam causing angular positioning errors. In addition, differences in pressure, temperature, and humidity on each side of the window will cause slight atmospheric index changes and induce deformation and a refractive index gradient within the window. Also, since the window is a dispersive media, the effect of both phase and group indices have to be considered. Taking all these factors into account, a method was developed to measure targets through multiple regions of different materials and produce results that are absolute measurements of target position in three-dimensional space, rather than simply relative position. The environment in which the lidar measurements are taken must be broken down into separate regions of interest and each region solved for separately. In this case, there were three regions of interest: air, fused silica, and vacuum. The angular position of the target inside the chamber is solved using only phase index and phase velocity, while the ranging effects due to travel from air to glass to vacuum/air are solved with group index and group velocity. When all parameters are solved simultaneously, an absolute knowledge of the position of each target within an environmental chamber can be derived. Novel features of this innovation include measuring absolute position of targets through multiple dispersive and non-dispersive media, deconstruction of lidar raw data from a commercial off-the-shelf unit into reworkable parameters, and use of group velocities to reduce range data. Measurement of structures within a vacuum chamber or other harsh environment, such as a furnace, may now be measured as easily as if they were in an ambient laboratory. This analysis permits transformation of the raw data into absolute spatial units (e.g., mm). This technique has also been extended to laser tracker, theodolite, and cathetometer measurements through refractive media.

Kubalak, David↗

The New Sun-Sky-Lunar Cimel CE318-T Multiband Photometer - A Comprehensive Performance Evaluation

This paper presents the new photometer CE318-T, able to perform daytime and night-time photometric measurements using the sun and the moon as light source. Therefore,this new device permits a complete cycle of diurnal aerosol and water vapour measurements valuable to enhance atmospheric monitoring to be extracted. In this study wehave found significantly higher precision of triplets when comparing the CE318-T master instrument and the Cimel AErosol RObotic NET work (AERONET) master (CE318-AERONET) triplets as a result of the new CE318-T tracking system. Regarding the instrument calibration, two new methodologies to transfer the calibration from a reference instrument using only daytime measurements (Sun Ratio and Sun-Moon gain factor techniques) are presented and discussed. These methods allow the reduction of the previous complexities inherent to nocturnal calibration. A quantitative estimation of CE318-T AOD uncertainty by means of error propagation theory during daytime revealed AOD uncertainties (u(sup D)(sub AOD)) for Langley-calibrated instruments similar to the expected values for other reference instruments (0.002-0.009). We have also found u(sup D)(sub AOD) values similar to the values reported in sun photometry for field instruments (approximately 0.015). In the case of the night-time period, the CE318-T-estimated standard combined uncertainty (u(sup N)(sub AOD)) is dependent not only on the calibration technique but also on illumination conditions and the instrumental noise. These values range from 0.011-0.018 for Lunar Langley-calibrated instruments to 0.012-0.021 for instruments calibrated using the Sun Ratio technique. In the case of moon-calibrated instruments using the Sun-Moon gain factor method and sun calibrated using the Langley technique, we found u(sup N)(sub AOD) ranging from 0.016 to 0.017 (up to 0.019 in 440 nm channel), not dependent on any lunar irradiance model. A subsequent performance evaluation including CE318-T and collocated measurements from independent reference instruments has served to assess the CE318-T performance as well as to confirm its estimated uncertainty. Daytime AOD evaluation, performed at Izana station from March to June 2014, encompassed measurements from a reference CE318-T, a CE318-AERONET master instrument, a Precision Filter Radiometer (PFR) and a Precision Spectroradiometer (PSR) prototype, reporting low AOD discrepancies between the four instruments (up to 0.006). The nocturnal AOD evaluation was performed using CE318-T- and starphotometer-collocated measurements and also by means of a day/night coherence transition test using the CE318-T master instrument and the CE318 daytime data from the CE318-AERONET master instrument. Results showed low discrepancies with the star photometer at 870 and 500 nm channels(less tna or equal to 0.013) and differences with AERONET daytime data (1 h after and before sunset and sunrise) in agreement with the estimated u(sup N)(sub AOD) values at all illumination conditions in the case of channels within the visible spectral range, and only for high moon's illumination conditions in the case of near infrared channels. Perceptible water vapour (PWV) validation showed a good agreement between CE318-T and Global Navigation Satellite System (GNSS) PWV values for all illumination conditions, within the expected precision for sun photometry. Finally, two case studies have been included to highlight the ability of the new CE318-T to capture the diurnal cycle of aerosols and water vapour as well as short-term atmospheric variations, critical for climate studies.

Barreto, Africa↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Ultra-Low Noise Laser and Optical Frequency Comb-Based Timing System for the Event Horizon Explorer (EHE) Mission

Very Long Baseline Interferometry (VLBI) requires a highly stable time reference and synchronization system to maintain coherence between observations recorded independently at multiple stations – for example, the array of ground-based radio telescopes that comprise the Event Horizon Telescope (EHT). Instability in the time reference or inaccuracy in the synchronization impacts phase conference and reduces the performance of the VLBI system. For space-based telescopes like the Event Horizon Explorer (EHE) mission, the precision timing system must meet the size, weight, and power (SWaP) constraints of the spacecraft platform and be capable of operating in the space environment. The EHE mission is a proposed extension to the EHT that can improve the angular and time resolution of observations by enabling longer interferometric baselines than are possible on the Earth and sampling a wide range of Fourier spatial frequencies throughout the orbital motion. In this effort, we demonstrate the performance of one option for the EHE mission: the use of a space-qualified, ultra-low noise laser developed as part of the Laser Interferometer Space Antenna (LISA) mission as the timing reference, and an optical frequency comb to transfer the stability of this laser to the microwave regime for instrumentation use. We describe the implementation of the microwave down-conversion, in which the LISA cavity-stabilized laser is locked to a self-referenced optical frequency comb to divide the optical frequency down to 100 MHz. The phase noise of the 100 MHz signal is measured and validated using a phase noise analyzer that is referenced to a separate laboratory ultra-stable laser system. We present the results of this experiment, which demonstrates that the performance of this system meets the EHE requirement: a relative frequency error of 1e-14 from 1 to 30 seconds.

Hannah Tomio↗

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗

Search for Pulsations from a Nearby Millisecond Pulsar and Wasilewski 49: Mirror for a Hidden Seyfert 1 Nucleus

Five studies are reported in this final report. The recently discovered 5.3 ms pulsar J1012+5307 at a distance of 520 pc is in an area of the sky which is particularly deficient in absorbing gas. The column density along the line of sight is less than 7.5 x 10(exp 19) CM(exp -2) which facilitates soft X-ray observations. Halpern reported a possible ROSAT Position Sensitive Proportional Counter (PSPC) detection of the pulsar in a serendipitous, off-axis observation. We have now confirmed the X-ray emission of PSR J1012+,5307 in a 23 ksec observation with the ROSAT High Resolution Imager (HRI). A point source is detected within 3" of the radio position. Its count rate of 1.6 +/- 0.3 x 10(exp -3) s(exp -1) corresponds to an unabsorbed 0. 1-2.4 keV flux of 6.4 x 10(exp -14) ergs cm(exp -2) s(exp -1), similar to that reported previously. This counts-to-flux conversion is valid for N(sub H) = 5 x 10(exp 19) cm(exp -2), and either a power-law spectrum of photon index 2.5 or a blackbody of kT = 0.1 keV. The implied X-ray luminosity of 2.0 x 10(exp 30) ergs s(exp -1) is 5 X 10(exp -4) of the pulsar's spin-down power dot-E, and similar to that of the nearest millisecond pulsar J0437-4715, which is nearly a twin of J1012+5307 in P dot-E. We subjected the 37 photons (and 13 background counts) within the source region to a pulsar search, but no evidence for pulsation was found. The pulsar apparently emits over a large fraction of its rotation cycle, and the absence of sharp modulation can be taken as evidence for surface thermal emission, as is favored for PSR J0437-4715, rather than magnetospheric X-ray emission which is apparent in the sharp pulses of the much more energetic millisecond pulsar B1821-24. A further test of this interpretation will be made with a longer ROSAT observation, which will increase the number of photons collected by a factor of 5, and permit a more sensitive examination of the light curve for modulation due to emission from heated polar caps. If found, such modulation will be further evidence that surface reheating by the impact of particles accelerated along open fiel;d lines operates in these approximately 10(exp -9) year old pulsars. In a second study, a new AM Her star serendipitously in a 25 day observation was detected with the EUVE satellite. A coherent period of 85.82 min is present in the EUVE Deep Survey imager light curve of this source. A spectroscopic optical identification is made with a 19th magnitude blue star that has H and He emission lines, and broad cyclotron humps typical of a magnetic cataclysmic variable. A lower limit to the polar magnetic field of 50 MG is estimated from the spacing of the cyclotron harmonics. EUVE J0425.6-5714 is also detected in archival ROSAT HRI observations spanning two months, and its stable and highly structured light curve permits us to fit a coherent ephemeris linking the ROSAT and EUVE data over a 1.3 yr gap. The derived period is 85.82107 +/- 0.00020 min, and the ephemeris should be accurate to 0.1 cycles until the year 2005. A narrow but partial X-ray eclipse suggests that this object belongs to the group of AM Her stars whose viewing geometry is such that the accretion stream periodically occults the soft X-ray emitting accretion spot on the surface of the white dwarf. A non-detection of hard X-rays from ASCA observations that are contemporaneous with the ROSAT HRI shows that the soft X-rays must dominate by at least an order of magnitude, which is consistent with a known trend among AM Her stars with large magnetic field. This object should not be confused with the Seyfert galaxy IH 0419-577 (= LB 1727), another X-ray/EUV source which lies only 3'95 away, and was the principal target of these monitoring observations. In a third report; the identity of the persistent high-energy (> 100 MeV) gamma-ray sources in the Galaxy, still largely a mystery is investigated. The second installment of the EGRET (2EG) lists a total of 128 sources, of which 51 are likely or possibly identified with AGNs, five with rotation-powered pulsars, and one is the LMC. There are 71 unidentified sources, of which 33, or almost half, lie in the narrow band of absolute value of b < 10 degrees along the Galactic plane. This excess of low-latitude sources must, therefore, constitute a Galactic population that is either similar to the already identified gamma-ray pulsars, or an entirely new class of gamma-ray emitters associated with the disk population. We are continuing our program, begun in A06, that is aimed at intermediate-latitude sources, arguing that X-ray detection of them is the most plausible method of identifying the Galactic population. The sources at high latitude must statistically be mostly AGNS, and are more straightforwardly identified through radio and optical means. For the fourth report, identification of the EGRET Source the identity of the persistent, high-energy gamma-ray sources in the Galactic plane is a mystery. The most likely scenario is a population of middle-aged pulsars, many of which could be radio quiet like Geminga. We have an ongoing program of ROSAT, VLA, and optical observations of selected EGRET error circles at intermediate Galactic latitude. For one of these fields, at (l, b) = (106 degrees, +3), our complete census of X-ray and radio sources, reveals a remarkable association between a radio shell with unique properties, and a complete census of X-ray and radio sources reveals a remarkable association between a radio shell with unique properties, and a compact X-ray source. Further observations are needed to determine whether or not this source has a hard X-ray spectrum like that of other gamma ray pulsars and, ideally, to find its pulsations. A fifth report reviews, Wasilewski 49, an interacting pair of Seyfert galaxies at z = 0.063, one of which contains a hidden Seyfert I nucleus as evidenced by broad wings on its Balmer lines.

Halpern, Jules P.↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗