Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

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 757 records · Page 42

Approximate Turbulent Boundary-Layer Development in Plane Compressible Flow Along Thermally Insulated Surfaces with Application to Supersonic-Tunnel Contour Correction

Numerical solutions of the differential equation obtained from the momentum theorem for the development of a turbulent boundary layer along a thermally insulated surface in two-dimensional and in radial shock-free flow are presented in tabular form for a range of Mach numbers from 0.100 to 10. The solution can be used in a step-wise procedure with any given distribution of favorable pressure gradients and for zero pressure gradients. Solutions are also given for use with moderate adverse pressure gradients. The mean velocity in the boundary layer is approximated by a power-law profile. In view of the stepwise integration methods to be used, the exponent designated the profile shape can be varied along the surface between the integral fraction limits 1/5 and 1/11 through interpolation. Agreement obtained between theoretical and experimental boundary-layer development in a supersonic nozzle at a nominal Mach number of 2 indicates the general validity of the approximations used in the analysis - in particular, the method of extrapolating low-speed skin-friction relations to high Mach number flows. The extrapolation method used assumes that the skin-friction coefficient depend primarily on Reynolds number, provided that the density and the kinematic viscosity are evaluated at surface conditions.

BOUNDARY LAYER - WIND TUNNEL WALLS↗

The Gas Distribution in the Outer Regions of Galaxy Clusters

Aims. We present our analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We have exploited the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius, We stacked the density profiles to detect a signal beyond T200 and measured the typical density and scatter in cluster outskirts. We also computed the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compared our average density and scatter profiles with the results of numerical simulations. Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond approximately r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict density profiles that are too steep, whereas runs including additional physics and/ or treating gas clumping agree better with the observed gas distribution. We report high-confidence detection of a systematic difference between cool-core and non cool-core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only small differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the ENZO simulations. Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Although their general behavior agrees more closely with the observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, AGN feedback, and consideration of gas clumping is required to construct realistic models of the outer regions of clusters.

Eckert, D.↗

Blocking Losses With a Photon Counter

It was not known how to assess accurately losses in a communications link due to photodetector blocking, a phenomenon wherein a detector is rendered inactive for a short time after the detection of a photon. When used to detect a communications signal, blocking leads to losses relative to an ideal detector, which may be measured as a reduction in the communications rate for a given received signal power, or an increase in the signal power required to support the same communications rate. This work involved characterizing blocking losses for single detectors and arrays of detectors. Blocking may be mitigated by spreading the signal intensity over an array of detectors, reducing the count rate on any one detector. A simple approximation was made to the blocking loss as a function of the probability that a detector is unblocked at a given time, essentially treating the blocking probability as a scaling of the detection efficiency. An exact statistical characterization was derived for a single detector, and an approximation for multiple detectors. This allowed derivation of several accurate approximations to the loss. Methods were also derived to account for a rise time in recovery, and non-uniform illumination due to diffraction and atmospheric distortion of the phase front. It was assumed that the communications signal is intensity modulated and received by an array of photon-counting photodetectors. For the purpose of this analysis, it was assumed that the detectors are ideal, in that they produce a signal that allows one to reproduce the arrival times of electrons, produced either as photoelectrons or from dark noise, exactly. For single detectors, the performance of the maximum-likelihood (ML) receiver in blocking is illustrated, as well as a maximum-count (MC) receiver, that, when receiving a pulse-position-modulated (PPM) signal, selects the symbol corresponding to the slot with the largest electron count. Whereas the MC receiver saturates at high count rates, the ML receiver may not. The loss in capacity, symbol-error-rate (SER), and count-rate were numerically computed. It was shown that the capacity and symbol-error-rate losses track, whereas the count-rate loss does not generally reflect the SER or capacity loss, as the slot-statistics at the detector output are no longer Poisson. It is also shown that the MC receiver loss may be accurately predicted for dead times on the order of a slot.

Moision, Burce E.↗

Yawed-Rolling Tire Mechanical Properties Testing of the Navy T-45 Aircraft Tires

The T-45 Goshawk is a United States Navy Jet aircraft used primarily as a trainer. The aircraft design makes use of "off the shelf" hardware as much as possible and was found to have unusual directional control issues during around operations. The aircraft was involved in numerous pilot-induced-oscillation incidents as well as observed to have unusual directional control reactions to failed main gear tires, a condition that is normally handled relatively easily by conventional aircraft steering control techniques. The behavior of the aircraft's tires had previously been modeled in simulators as a result of approximations provided in 40-year-old reference publications. Since knowledge of the true tire cornering and braking behavior is essential to modeling, understanding, and fixing directional control problems, the United States Navy requested assistance from the NASA Langley Research Center's (LARC) Aircraft Landing Dynamics Facility (ALDF) to define the yawed-rolling mechanical properties of the T-45 aircraft tires. The purpose of this report is to document the results of testing the subject tires at the NASA LaRC ALDF in September 1998. Brief descriptions of the Instrumented Tire Test Vehicle (ITTV) are included to familiarize the reader with the ITTV capabilities, data acquisition system, test and measurement techniques, data accuracy, and analysis and presentation of the testing results.

Daugherty, Robert H.↗

The Gas Distribution in Galaxy Cluster Outer Regions

Aims. We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r200 and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. Results. As opposed to some recent Suzaku results, and confirming previous evidence from ROSAT and Chandra, we observe a steepening of the density profiles beyond approximately r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or treating gas clumping are in better agreement with the observed gas distribution. We report for the first time the high-confidence detection of a systematic difference between cool-core and non-cool core clusters beyond 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. Conclusions. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Although their general behavior is in better agreement with the observations, simulations including cooling and star formation convert a large amount of gas into stars, which results in a low gas fraction with respect to the observations. Consequently, a detailed treatment of gas cooling, star formation, AGN feedback, and taking into account gas clumping is required to construct realistic models of cluster outer regions.

Eckert, D.↗

Direct Numerical Simulation of an Airfoil with Sand Grain Roughness on the Leading Edge

As part of a computational study of acoustic radiation due to the passage of turbulent boundary layer eddies over the trailing edge of an airfoil, the Lattice-Boltzmann method is used to perform direct numerical simulations of compressible, low Mach number flow past an NACA 0012 airfoil at zero degrees angle of attack. The chord Reynolds number of approximately 0.657 million models one of the test conditions from a previous experiment by Brooks, Pope, and Marcolini at NASA Langley Research Center. A unique feature of these simulations involves direct modeling of the sand grain roughness on the leading edge, which was used in the abovementioned experiment to trip the boundary layer to fully turbulent flow. This report documents the findings of preliminary, proof-of-concept simulations based on a narrow spanwise domain and a limited time interval. The inclusion of fully-resolved leading edge roughness in this simulation leads to significantly earlier transition than that in the absence of any roughness. The simulation data is used in conjunction with both the Ffowcs Williams-Hawkings acoustic analogy and a semi-analytical model by Roger and Moreau to predict the farfield noise. The encouraging agreement between the computed noise spectrum and that measured in the experiment indicates the potential payoff from a full-fledged numerical investigation based on the current approach. Analysis of the computed data is used to identify the required improvements to the preliminary simulations described herein.

Ribeiro, Andre F. P.↗

Direct Numerical Simulation of Hypersonic Turbulent Boundary Layer inside an Axisymmetric Nozzle

As a first step toward a study of acoustic disturbance field within a conventional, hypersonic wind tunnel, direct numerical simulations (DNS) of a Mach 6 turbulent boundary layer on the inner wall of a straight axisymmetric nozzle are conducted and the results are compared with those for a flat plate. The DNS results for a nozzle radius to boundary-layer thickness ratio of 5:5 show that the turbulence statistics of the nozzle-wall boundary layer are nearly unaffected by the transverse curvature of the nozzle wall. Before the acoustic waves emanating from different parts of the nozzle surface can interfere with each other and undergo reflections from adjacent portions of the nozzle surface, the rms pressure fluctuation beyond the boundary layer edge increases toward the nozzle axis, apparently due to a focusing effect inside the axisymmetric configuration. Spectral analysis of pressure fluctuations at both the wall and the freestream indicates a similar distribution of energy content for both the nozzle and the flat plate, with the peak of the premultiplied frequency spectrum at a frequency of [(omega)(delta)]/U(sub infinity) approximately 6.0 inside the free stream and at [(omega)(delta)]/U(sub infinity) approximately 2.0 along the wall. The present results provide the basis for follow-on simulations involving reverberation effects inside the nozzle.

Huang, Junji↗

Simple approximations for the asymptotic description of the interaction between a normal shock wave and a turbulent boundary layer at transonic speeds

The asymptotic description of the interaction between a normal shock wave and a turbulent boundary layer is reviewed. The layers necessary in a rational analysis of the interaction are discussed with emphasis on the differences from an interaction with a laminar boundary layer, the uncoupling of solutions for the distribution of pressure and skin friction at the wall, and the role of the Reynolds shear stress in these solutions. The accuracy of asymptotic solutions in flows at Reynolds numbers of technical interest is discussed. Solutions for the distribution of pressure and skin friction at the wall and the shape of the shock are considered for the case where the flow is near separation. For the pressure and skin friction, it is possible to write two simplified partial solutions, one valid at the beginning of the interaction and one valid somewhat downstream of the shock wave. A solution composed of these two parts and a linear interpolation between them appears to give good comparison with experiment; one unknown constant, independent of the parameters of the interaction, must be found from experiment. The simplified relations are presented. Comparison of numerical computations with experimental data indicates a possible value for the constant and shows quite satisfactory results.

Adamson, T. C., Jr.↗

Incorporation of a progressive failure analysis method in the CSM testbed software system

Analysis of the postbuckling behavior of composite shell structures pose many difficult and challenging problems in the field of structural mechanics. Current analysis methods perform well for most cases in predicting the postbuckling response of undamaged components. To predict component behavior accurately at higher load levels, the analysis must include the effects of local material failures. The CSM testbed software system is a highly modular structural analysis system currently under development at Langley Research Center. One of the primary goals of the CSM testbed is to provide a software environment for the development of advanced structural analysis methods and modern numerical methods which will exploit advanced computer architecture such as parallel-vector processors. Development of a progressive failure analysis method consists of the design and implementation of a processor which will perform the ply-level progressive failure analysis and the development of a geometrically nonlinear analysis procedure which incorporates the progressive failure processor. Regarding the development of the progressive failure processor, two components are required: failure criteria and a degradation model. For the initial implementation, the failure criteria of Hashin will be used. For a matrix failure which typically indicates the development of transverse matrix cracks, the ply properties will be degraded. Work to date includes the design of the progressive failure analysis processor and initial plans for the controlling geometrically nonlinear analysis procedure. The implementation of the progressive failure analysis has begun. Access to the model database and the Hashin failure criteria are completed. Work is in progress on the input/output operations for the processor related data and the finite element model updating procedures. In total the progressive failure processor is approximately one-third complete.

Arenburg, Robert T.↗

Effects of pressure gradients on turbulent premixed flames

The influence of a constant acceleration on a turbulent premixed flame is studied by direct numerical simulation. This acceleration induces a mean pressure gradient across the flame brush, leading to a modification of the turbulent flame structure due to differential buoyancy mechanisms between heavy cold fresh and light hot burnt gases. Such a pressure gradient may be encountered in practical applications in ducted flames. A favorable pressure gradient, i.e. the pressure decreases from unburnt to burnt gases, is found to decrease the flame wrinkling, the flame brush thickness, and the turbulent flame speed. A favorable pressure gradient also promotes counter-gradient turbulent transport. On the other hand, adverse pressure gradients tend to increase the flame brush thickness and turbulent flame speed, and promote classical gradient turbulent transport. The balance equation for the turbulent flux of the Favre averaged progress variable is also analyzed. The first results show that the fluctuating pressure term, cannot be neglected as generally assumed in models. Simple models assuming that a high mean pressure gradient may only be balanced by the cross-dissipation term seem too approximate. This analysis has to be continued to compare simulation data and closure schemes proposed for the transport equation. The analysis developed by Veynante et al.(1995) has been extended to imposed acceleration and mean pressure gradients. A simple model for the turbulent flux is proposed and validated from simulation data. Then, a modified criterion is derived to delineate between counter-gradient and gradient turbulent diffusion. In fact, counter-gradient diffusion may occur in most practical applications, especially for ducted flames.

Veynante, D.↗

Noise Sources in a Low-Reynolds-Number Turbulent Jet at Mach 0.9

The mechanisms of sound generation in a Mach 0.9, Reynolds number 3600 turbulent jet are investigated by direct numerical simulation. Details of the numerical method are briefly outlined and results are validated against an experiment at the same flow conditions. Lighthill's theory is used to define a nominal acoustic source in the jet, and a numerical solution of Lighthill's equation is compared to the simulation to verify the computational procedures. The acoustic source is Fourier transformed in the axial coordinate and time and then filtered in order to identify and separate components capable of radiating to the far field. This procedure indicates that the peak radiating component of the source is coincident with neither the peak of the full unfiltered source nor that of the turbulent kinetic energy. The phase velocities of significant components range from approximately 5% to 50% of the ambient sound speed which calls into question the commonly made assumption that the noise sources convect at a single velocity. Space-time correlations demonstrate that the sources are not acoustically compact in the streamwise direction and that the portion of the source that radiates at angles greater than 45 deg. is stationary. Filtering non-radiating wavenumber components of the source at single frequencies reveals that a simple modulated wave forms for the source, as might be predicted by linear stability analysis. At small angles from the jet axis the noise from these modes is highly directional, better described by an exponential than a standard Doppler factor.

Freund, Jonathan B.↗

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

Effect of enclosure shape on natural convection velocities

A numerical analysis was performed to compare natural convection velocities in two dimensional enclosures of various shape. The following shapes were investigated: circle, square, horizontal and upright 2 x 1 aspect ratio rectangles, horizontal and upright half circles, diamond. In all cases, the length scale in the various dimensionless parameters, such as Rayleigh number, is defined as the diameter of the equal area circle. Natural convection velocities were calculated for Rayleigh numbers of 1000 and 5000 with the temperature difference taken to be across (1) the maximum horizontal dimension, (2) the median horizontal line (line through centroid) and (3) the horizontal distance such that the temperature gradient is the same for shapes of equal area. For the class of shapes including the square, upright half circle and upright rectangle, the computed velocities were found to agree very closely with that of the equal area circle when the temperature difference is taken to be across the maximum horizontal dimension (condition (a)). The velocities for the horizontal rectangle and half circle were found to be approximately one half that of the equal area circle for the same condition. Better overall agreement among all shapes was obtained by setting the temperature difference across a distance such that the temperature gradients were equal for shapes of equal area.

Robertson, S. J.↗

Cloud atlas for the FIRE Cirrus Intensive Field Observation (IFO)

An Intensive Field Observation (IFO) of cirrus clouds was conducted over the mid-western U.S. during the period October 13 to November 2, 1986. This activity, part of the First ISCCP Regional Experiment (FIRE), included measurements made from specially deployed instruments on the ground, balloons, and aircraft as well as observations from existing operational and experimental satellites. One of the sets of satellite observations was the radiance measurements made with the 5-channel AVHRR radiometer on the NOAA 9 polar orbiting meteorological satellite. The ground resolution of the measurements at nadir is approx. 1 km. It is these measurements, made once each day at approximately 2:30 p.m. local time, that were used in determining the present cloud atlas. The area covered by the atlas is slightly larger than the area specified for the IFO, in order to be in alignment with the grid that will be used in a forthcoming atlas for the larger, ETO region. The atlas contains four pages of information for each satellite pass. The 1st page of each group shows the distribution of measured radiances in channel 1 (normalized to the incoming solar flux multiplied by the cosine of the solar zenith angle) and in channel 4 for the area as a whole and for each analysis box. The 2nd page shows the images in: channels 1 and 2, channel 3R; and channel 4. The 3rd page shows the retrieved parameters in graphical form for the region as a whole and for each analysis box, where cloud fraction appears as a contour plot with respect to optical thickness and cloudtop temperature. The 4th page provides a statistical summary of the retrieved parameters in numerical form for each analysis box.

Arking, Albert↗

Results of the Workshop on Impact Cratering: Bridging the Gap Between Modeling and Observations

On February 7-9,2003, approximately 60 scientists gathered at the Lunar and Planetary Institute in Houston, Texas, for a workshop devoted to improving knowledge of the impact cratering process. We (co-conveners Elisabetta Pierazzo and Robert Herrick) both focus research efforts on studying the impact cratering process, but the former specializes in numerical modeling while the latter draws inferences from observations of planetary craters. Significant work has been done in several key areas of impact studies over the past several years, but in many respects there seem to be a disconnect between the groups employing different approaches, in particular modeling versus observations. The goal in convening this workshop was to bring together these disparate groups to have an open dialogue for the purposes of answering outstanding questions about the impact process and setting future research directions. We were successful in getting participation from most of the major research groups studying the impact process. Participants gathered from five continents with research specialties ranging from numerical modeling to field geology, and from small-scale experimentation and geochemical sample analysis to seismology and remote sensing.With the assistance of the scientific advisory committee (Bevan French, Kevin Housen, Bill McKinnon, Jay Melosh, and Mike Zolensky), the workshop was divided into a series of sessions devoted to different aspects of the cratering process. Each session was opened by two invited t a b , one given by a specialist in numerical or experimental modeling approaches, and the other by a specialist in geological, geophysical, or geochemical observations. Shorter invited and contributed talks filled out the sessions, which were then concluded with an open discussion time. All modelers were requested to address the question of what observations would better constrain their models, and all observationists were requested to discuss how their observations can constrain modeling efforts.

Herrick, Robert↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s Global Earth Observing System

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy, Jr.↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and increases in number from approximately 6 thousand RO profiles per day during 2019 to approximately 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2022. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗