Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data intensive computing”

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 739 records · Page 41

Ames S-32 O-16 O-18 Line List for High-Resolution Experimental IR Analysis

By comparing to the most recent experimental data and spectra of the SO2 628 ν1/ν3 bands (see Ulenikov et al., JQSRT 168 (2016) 29-39), this study illustrates the reliability and accuracy of the Ames-296K SO2 line list, which is accurate enough to facilitate such high-resolution spectroscopic analysis. The SO2 628 IR line list is computed on a recently improved potential energy surface (PES) refinement, denoted Ames-Pre2, and the published purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface. Progress has been made in both energy level convergence and rovibrational quantum number assignments agreeing with laboratory analysis models. The accuracy of the computed 628 energy levels and line list is similar to what has been achieved and reported for SO2 626 and 646, i.e. 0.01-0.03 cm(exp −1) for bands up to 5500 cm(exp −1). During the comparison, we found some discrepancies in addition to overall good agreements. The three-IR-list based feature-by-feature analysis in a 0.25 cm(exp −1) spectral window clearly demonstrates the power of the current Ames line lists with new assignments, correction of some errors, and intensity contributions from varied sources including other isotopologues. We are inclined to attribute part of detected discrepancies to an incomplete experimental analysis and missing intensity in the model. With complete line position, intensity, and rovibrational quantum numbers determined at 296 K, spectroscopic analysis is significantly facilitated especially for a spectral range exhibiting such an unusually high density of lines. The computed 628 rovibrational levels and line list are accurate enough to provide alternatives for the missing bands or suspicious assignments, as well as helpful to identify these isotopologues in various celestial environments. The next step will be to revisit the SO2 828 and 646 spectral analyses.

Sulfur dioxide↗

Ab Initio Calculations of Water Line Strengths

We report on the determination of a high quality ab initiu potential energy surface (PES) and dipole moment function for water. This PES is empirically adjusted to improve the agreement between the computed line positions and those from the HITRAN 92 data base with J less than 6 for H2O. The changes in the PES are small, nonetheless including an estimate of core (oxygen 1s) electron correlation greatly improves the agreement with experiment. Using this adjusted PES, we can match 30,092 of the 30,117 transitions in the HITRAN 96 data base for H2O with theoretical lines. The 10,25,50,75, and 90 percentiles of the difference between the calculated and tabulated line positions are -0.11, -0.04, -0.01, 0.02, and 0.07 l/cm. Non-adiabatic effects are not explicitly included. About 3% of the tabulated line positions appear to be incorrect. Similar agreement using this adjusted PES is obtained for the oxygen 17 and oxygen 18 isotopes. For HDO, the agreement is not as good, with root-mean-square error of 0.25 l/cm for lines with J less than 6. This error is reduced to 0.02 l/cm by including a small asymmetric correction to the PES, which is parameterized by simultaneously fitting to HDO md D2O data. Scaling this correction by mass factors yields good results for T2O and HTO. The intensities summed over vibrational bands are usually in good agreement between the calculations and the tabulated results, but individual lines strengths can differ greatly. A high temperature list consisting of 307,721,352 lines is generated for H2O using our PES and dipole moment function.

Schwenke, David W.↗

Hybrid LES/RANS Simulation of Transverse Sonic Injection into a Mach 2 Flow

A computational study of transverse sonic injection of air and helium into a Mach 1.98 cross-flow is presented. A hybrid large-eddy simulation / Reynolds-averaged Navier-Stokes (LES/RANS) turbulence model is used, with the two-equation Menter baseline (Menter-BSL) closure for the RANS part of the flow and a Smagorinsky-type model for the LES part of the flow. A time-dependent blending function, dependent on modeled turbulence variables, is used to shift the closure from RANS to LES. Turbulent structures are initiated and sustained through the use of a recycling / rescaling technique. Two higher-order discretizations, the Piecewise Parabolic Method (PPM) of Colella and Woodward, and the SONIC-A ENO scheme of Suresh and Huyhn are used in the study. The results using the hybrid model show reasonably good agreement with time-averaged Mie scattering data and with experimental surface pressure distributions, even though the penetration of the jet into the cross-flow is slightly over-predicted. The LES/RANS results are used to examine the validity of commonly-used assumptions of constant Schmidt and Prandtl numbers in the intense mixing zone downstream of the injection location.

Boles, John A.↗

Further studies with data collected by NASA's airborne Doppler lidar in Oklahoma in 1981

Continued study of the lidar data collected in 1981 has resulted in significant new improvements in the analysis techniques reported by Bluestein et al. (1985) and McCaul (1985). Through comparison of fore- and aft-derived scalar fields of intensity and spectral width, the self-consistency of the lidar moment estimates was assessed. Reflectivity estimates were found to be quite stable and reliable, while spectral widths were prone to become noisy if signal to noise ratio (SNR) fell below 12 dB. In addition, spectral widths contained a significant component due to radial velocity gradients in areas along gust fronts, and these components were different along the fore and aft lines of sight. Significant improvement in agreement between the fore and aft fields of spectral width was obtained by estimating the radial velocity gradient component and then removing it from the raw measured widths to yield only the turbulent portion of the contribution to width. Additional analyses showed that lidar-derived vorticity estimates were consistent with several approximate models of vorticity growth along gust front zones, and with the hypothesis that Helmholtz instability could have been responsible for vortices seen along part of the gust front of 30 June 1981. Computations of divergence transverse to axes through an isolated cumulus congestus indicated that the strongest convergence tended to lie along an axis parallel to the congestus. This and the results of other additional analyses seem to suggest that the lidar winds do indeed accurately reflected the basic features of the real wind field.

Bluestein, H. B.↗

Convection Parametrization and Multi-Nesting Dependence of a Heavy Rainfall Event over Namibia with Weather Research and Forecasting (WRF) Model

Namibia is considered as one of the countries that are most vulnerable to climate change due to its generally dry climate and the percentage of its population that rely on subsistence agriculture for their livelihoods. Early warning systems are an important aspect of adapting to climate change. Weather forecasting relies on the use of numerical weather prediction models and these need to be configured properly. In this study we investigate the effects of using multi-nests and a convection scheme on the simulation of a heavy rainfall event over the Northwestern region of Kunene, Namibia. The event which was associated with a cut-off low system, was short-lived and resulted in over 45 mm of rainfall in one hour. For the multi-nest, a 9 km grid length parent domain is nested within the Global Forecast System (GFS) simulations, which in turn forces a 3 km grid spacing child domain. A different set of simulations are produced using a single nest of 3 km grid spacing, nested directly inside the GFS data. The simulations are produced with the convection scheme switched on and off. The impact of a single versus multi nest is found to be small in general, with slight differences in the location of high rainfall intensity. Switching off the convection schemes results in high rainfall intensity and increased detail in the simulations, including when a grid spacing of 9 km mis used. Using a grid spacing of 3 km with the convection scheme on, results in a loss of detail in the simulations as well as lower rainfall amounts. The study shows a need for different configurations to be tested before an optimum configuration can be selected for operational forecasting. We recommend further tests with different synoptic forcing to be conducted to identify a suitable configuration for Namibia.

flooding↗

Multiphase Simulations of the SLS Launch Environment

NASA’s Space Launch System (SLS), which will send astronauts back to the Moon in the next few years, is powered by four RS-25 engines and two RSRMV solid rocket boosters (SRBs). During launch the SLS propulsion system generates intense acoustics and other powerful waves, such as ignition overpressure (IOP) which, if unmitigated, have the potential to damage the vehicle and possibly cause loss of mission or crew. To protect the vehicle from these powerful waves, the SLS launch pad design includes an ignition overpressure/sound suppression (IOP/SS) system which sprays 270,000 gallons per minute of water very close to the SRB and RS-25 nozzles. The SRB and RS-25 engine plumes, and the proximity of the IOP/SS water, create a complex multiphase (gas and liquid) environment during the SLS ignition sequence. The interplay among these systems creates challenges related to water spray into/onto engine nozzles, potential debris transport, and additional transient loads due to strong plume-water interactions - all of which the SLS vehicle must be able to withstand. Prior to the Artemis I launch, the SLS multiphase liftoff environment was largely unknown due to differences from the Space Shuttle and other programs. Some data was available from tests of individual systems, but no integrated testing or analysis was available. Even post-launch analysis of Artemis I cannot provide a full understanding of the complex physics involved due to limited (or obstructed) camera views and instrumentation. Computational fluid dynamics (CFD) is being used to investigate the details of the multiphase environment which could not be measured, help comprehend the data gathered from the launch, and ultimately identify phenomena that are a concern for future flights. Project Details Engineers at NASA’s Marshall Space Flight Center (MSFC) have executed simulations using the Loci/STREAM-Volume of Fluid (VoF) multiphase CFD solver to understand this environment. Initial efforts successfully validated the CFD solver on various tests, giving confidence to simulate the SLS multiphase liftoff environment prior to the Artemis I launch. The CFD simulation of the SLS ignition sequence was conducted in three phases. First the IOP/SS water system was simulated for approximately 6 seconds to reach a quasi-steady state. Next, the RS-25 engine plumes were activated and held at full power for 1 second. Lastly, the SRB booster was activated and the simulation was carried out until just prior to vehicle motion. This simulation process mimics the conditions that exist at launch. Results and Impact The SLS ignition sequence simulation results provide deep understanding of the underlying physics occuring during launch. Observations from the simulation include reduction of water splashing into/onto the engine nozzles, change in angling of the dense water sheets, and the origin of the powerful ignition overpressure (IOP) wave. These observations directly inform the SLS program on subjects including plume-water induced side loads, debris transport, and the acoustic launch environment. Additionally, with post launch comparison of CFD observations to flight data, these tools can be applied to launch vehicles and environments other than SLS with confidence. Why HPC Matters The SLS ignition sequence CFD simulations are conducted on meshes up to hundreds of millions of cells on thousands of processors for weeks at a time. These simulations generate terabytes of data that must also be stored and archived for future use on HPC systems. Simply put, the CFD simulations would not be possible without NASA HPC resources. What’s Next Comparisons between the Artemis I flight data and the CFD simulations will be continued to both improve confidence in the CFD results and provide deeper understanding into the SLS multiphase launch environment. This will be used to provide insight for decision making for the first manned SLS flight, Artemis II. Future simulations will target new configurations of the SLS IOP/SS water required to support the more powerful variants of the SLS vehicle, such as Block 1B. Additionally, this capability provides NASA the ability to investigate launch environments for vehicles other than SLS to support other missions.

Travis Rivord↗

Fortran for the nineties

Fortran has largely enjoyed prominence for the past few decades as the computer programming language of choice for numerically intensive scientific, engineering, and process control applications. Fortran's well understood static language syntax has allowed resulting parsers and compiler optimizing technologies to often generate among the most efficient and fastest run-time executables, particularly on high-end scalar and vector supercomputers. Computing architectures and paradigms have changed considerably since the last ANSI/ISO Fortran release in 1978, and while FORTRAN 77 has more than survived, it's aged features provide only partial functionality for today's demanding computing environments. The simple block procedural languages have been necessarily evolving, or giving way, to specialized supercomputing, network resource, and object-oriented paradigms. To address these new computing demands, ANSI has worked for the last 12-years with three international public reviews to deliver Fortran 90. Fortran 90 has superseded and replaced ISO FORTRAN 77 internationally as the sole Fortran standard; while in the US, Fortran 90 is expected to be adopted as the ANSI standard this summer, coexisting with ANSI FORTRAN 77 until at least 1996. The development path and current state of Fortran will be briefly described highlighting the many new Fortran 90 syntactic and semantic additions which support (among others): free form source; array syntax; new control structures; modules and interfaces; pointers; derived data types; dynamic memory; enhanced I/O; operator overloading; data abstraction; user optional arguments; new intrinsics for array, bit manipulation, and system inquiry; and enhanced portability through better generic control of underlying system arithmetic models. Examples from dynamical astronomy, signal and image processing will attempt to illustrate Fortran 90's applicability to today's general scalar, vector, and parallel scientific and engineering requirements and object oriented programming paradigms. Time permitting, current work proceeding on the future development of Fortran 2000 and collateral standards will be introduced.

Himer, J. T.↗

Optimization of Car Body under Constraints of Noise, Vibration, and Harshness (NVH), and Crash

To be competitive on the today's market, cars have to be as light as possible while meeting the Noise, Vibration, and Harshness (NVH) requirements and conforming to Government-man dated crash survival regulations. The latter are difficult to meet because they involve very compute-intensive, nonlinear analysis, e.g., the code RADIOSS capable of simulation of the dynamics, and the geometrical and material nonlinearities of a thin-walled car structure in crash, would require over 12 days of elapsed time for a single design of a 390K elastic degrees of freedom model, if executed on a single processor of the state-of-the-art SGI Origin2000 computer. Of course, in optimization that crash analysis would have to be invoked many times. Needless to say, that has rendered such optimization intractable until now. The car finite element model is shown. The advent of computers that comprise large numbers of concurrently operating processors has created a new environment wherein the above optimization, and other engineering problems heretofore regarded as intractable may be solved. The procedure, shown, is a piecewise approximation based method and involves using a sensitivity based Taylor series approximation model for NVH and a polynomial response surface model for Crash. In that method the NVH constraints are evaluated using a finite element code (MSC/NASTRAN) that yields the constraint values and their derivatives with respect to design variables. The crash constraints are evaluated using the explicit code RADIOSS on the Origin 2000 operating on 256 processors simultaneously to generate data for a polynomial response surface in the design variable domain. The NVH constraints and their derivatives combined with the response surface for the crash constraints form an approximation to the system analysis (surrogate analysis) that enables a cycle of multidisciplinary optimization within move limits. In the inner loop, the NVH sensitivities are recomputed to update the NVH approximation model while keeping the Crash response surface constant. In every outer loop, the Crash response surface approximation is updated, including a gradual increase in the order of the response surface and the response surface extension in the direction of the search. In this optimization task, the NVH discipline has 30 design variables while the crash discipline has 20 design variables. A subset of these design variables (10) are common to both the NVH and crash disciplines. In order to construct a linear response surface for the Crash discipline constraints, a minimum of 21 design points would have to be analyzed using the RADIOSS code. On a single processor in Origin 2000 that amount of computing would require over 9 months! In this work, these runs were carried out concurrently on the Origin 2000 using multiple processors, ranging from 8 to 16, for each crash (RADIOSS) analysis. Another figure shows the wall time required for a single RADIOSS analysis using varying number of processors, as well as provides a comparison of 2 different common data placement procedures within the allotted memories for each analysis. The initial design is an infeasible design with NVH discipline Static Torsion constraint violations of over 10%. The final optimized design is a feasible design with a weight reduction of 15 kg compared to the initial design. This work demonstrates how advanced methodology for optimization combined with the technology of concurrent processing enables applications that until now were out of reach because of very long time-to-solution.

Kodiyalam, Srinivas↗

Rapid Chemometric Filtering of Spectral Data

A method of rapid, programmable filtering of spectral transmittance, reflectance, or fluorescence data to measure the concentrations of chemical species has been proposed. By programmable is meant that a variety of spectral analyses can readily be performed and modified in software, firmware, and/or electronic hardware, without need to change optical filters or other optical hardware of the associated spectrometers. The method is intended to enable real-time identification of single or multiple target chemical species in applications that involve high-throughput screening of multiple samples. Examples of such applications include (but are not limited to) combinatorial chemistry, flow cytometry, bead assays, testing drugs, remote sensing, and identification of targets. The basic concept of the proposed method is to perform real-time crosscorrelations of a measured spectrum with one or more analytical function(s) of wavelength that could be, for example, the known spectra of target species. Assuming that measured spectral intensities are proportional to concentrations of target species plus background spectral intensities, then after subtraction of background levels, it should be possible to determine target species concentrations from cross-correlation values. Of course, the problem of determining the concentrations is more complex when spectra of different species overlap, but the problem can be solved by use of multiple analytical functions in combination with computational techniques that have been developed previously for analyses of this type. The method is applicable to the design and operation of a spectrometer in which spectrally dispersed light is measured by means of an active-pixel sensor (APS) array. The row or column dimension of such an array is generally chosen to be aligned along the spectral-dispersion dimension, so that each pixel intercepts light in a narrow spectral band centered on a wavelength that is a known function of the pixel position. The proposed method admits of two hardware implementations for computing cross-correlations in real time.

Beaman, Gregory↗

Stable isotope laser spectroscopy

Recent advances in semiconductor laser technology have produced a reliable lightweight device ideally suited for a spacecraft high resolution molecular spectrometer. Lead-salt tunable diode lasers (TDL) emit in several spectral modes, each with a very narrow linewidth of -0.0003/cm. This spectral resolution is much narrower than typical Doppler broadened molecular linewidths in the mid-IR range. Thus it is possible to detect individual rotational lines within the vibrational band and measure their intensity, which can be used to determine gas concentration. The narrow spectral lines of any impurity gas tend to lie between the narrow lines of the gas of interest. This represents a major advantage over the accepted gas chromatograph mass spectrometer (GCMS) technique for measuring gas concentrations and isotope ratios. The careful and extensive gas purification procedures required to remove impurities for reliable GCMS measurements will not be required for an IR laser gas analysis. The infrared laser gas analysis technique is being developed to measure stable isotopic ratios of gases such as CO2, CH4, N2O, and NH3. This will eventually lead to development of instruments capable of in situ istopic measurements on planets such as Mars. The carbon (C-12, C-13) isotope ratio is indicative of the type of carbon fixation mechanisms (e.g., photosynthesis, respiration) in operation on a planet, while the nitrogen (N-14, N-15) isotope ratio can probably be used to date nitrogen-bearing Martian samples. The absorbance ratio of two adjacent lines of CO2 in the 2300/cm (4.3 micron) region of the spectrum was measured. The precision of the measurement is presently better than 1 percent and significant improvement is anticipated as rapid sweep-integration techniques and computer controlled data acquistion capabilities are incorporated.

Becker, J. F.↗

X-ray Mapping of Terrestrial and Extraterrestrial Materials Using the Electron Microprobe

Lunar samples returned from the Apollo program motivated development of the Bence-Albee algorithm for the rapid and accurate analysis of lunar materials, and established interlaboratory comparability through its common use. In the analysis of mineral and rock fragments it became necessary to combine micro- and macroscopic analysis by coupling electron-probe microanalysis (EPMA) with automated stage point counting. A coarse grid that included several thousand points was used, and initially wavelength-dispersive (WDS) and later energydispersive (EDS) data were acquired at discrete stage points using approx. 5 sec count times. A approx 50 micrometer beam diameter was used for WDS and up to 500 micrometer beam diameter for EDS analysis. Average analyses of discretely sampled phases were coupled with the point count data to calculate the bulk composition using matrix algebra. Use of a defocused beam resulted in a contribution from multiple phases to each analytical point, and the analytical data were deconvolved relative to end-member phase chemistry on the fly. Impressive agreement was obtained between WDS and EDS measurements as well as comparison with bulk chemistry obtained by other methods. In the 30 years since these methods were developed, significant improvements in EPMA automation and computer processing have taken place. Digital beam control allows routine collection of x-ray maps by EDS, and stage mapping for WDS is conducted continuously at slew speed and incrementally by sampling at discrete points. Digital pulse processing in EDS systems has significantly increased the throughput for EDS mapping, and the ongoing development of Si-drift detector systems promises mapping capabilities rivaling WDS systems. Spectrum imaging allows a data cube of EDS spectra to be acquired and sophisticated processing of the original data is possible using matrix algebra techniques. The study of lunar and meteoritic materials includes the need to conveniently: (1) Characterize the sample at microscopic and macroscopic scales with relatively high sensitivity, (2) Determine the modal abundance of minerals, and (3) Identify and relocate discrete features of interest in terms of size and chemistry. The coupled substitution of cations in minerals can result in significant variation in mineral chemistry, but at similar average Z, leading to poor backscattered-electron (BSE) contrast discrimination of mineralogy. It is necessary to discriminate phase chemistry at both the trace element level and the major element level. To date, the WDS of microprobe systems is preferred for mapping due to high throughput and the ability to obtain the necessary intensity to discriminate phases at both trace and major element concentrations. It is desirable to produce fully quantitative compositional maps of geological materials, which requires the acquisition of k-ratio maps that are background and dead-time corrected, and which have been corrected by phi(delta z> or an equivalent algorithm at each pixel. To date, turnkey systems do not allow the acquisition of k-ratio maps and the rigorous correction in this manner. X-ray maps of a chondrule from the Ourique meteorite, and a comb-layered xenolith from the San Francisco volcanic field, have been analyzed and processed to extract phase information. The Ourique meteorite presents a challenge due to relatively low BSE contrast, and has been studied using spectrum imaging. X-ray maps for Si, Mg, and FeK(alpha) were used to produce RGB images. The xenolith sample contains sector-zoned augite, olivine, plagioclase, and basaltic glass. X-ray maps were processed using Lispix and ImageJ software to produce mineral phase maps. The x-ray maps for Mg, Ca, and Ti were used with traceback to generate binary images that were converted to RGB images. These approaches are successful in discriminating phases, but it is desirable to achieve the methods that were used on lunar samples 30 years ago on current microprobe systems. Curnt research includes x-ray mapping analysis of the Dalgety Downs chondrite by micro x-ray fluorescence and spectrum imaging, in collaboration with Kenny Witherspoon of IXRF Systems and Dale Newbury of NIST.

Carpenter, P.↗

Robotic Vision-Based Localization in an Urban Environment

A system of electronic hardware and software, now undergoing development, automatically estimates the location of a robotic land vehicle in an urban environment using a somewhat imprecise map, which has been generated in advance from aerial imagery. This system does not utilize the Global Positioning System and does not include any odometry, inertial measurement units, or any other sensors except a stereoscopic pair of black-and-white digital video cameras mounted on the vehicle. Of course, the system also includes a computer running software that processes the video image data. The software consists mostly of three components corresponding to the three major image-data-processing functions: Visual Odometry This component automatically tracks point features in the imagery and computes the relative motion of the cameras between sequential image frames. This component incorporates a modified version of a visual-odometry algorithm originally published in 1989. The algorithm selects point features, performs multiresolution area-correlation computations to match the features in stereoscopic images, tracks the features through the sequence of images, and uses the tracking results to estimate the six-degree-of-freedom motion of the camera between consecutive stereoscopic pairs of images (see figure). Urban Feature Detection and Ranging Using the same data as those processed by the visual-odometry component, this component strives to determine the three-dimensional (3D) coordinates of vertical and horizontal lines that are likely to be parts of, or close to, the exterior surfaces of buildings. The basic sequence of processes performed by this component is the following: 1. An edge-detection algorithm is applied, yielding a set of linked lists of edge pixels, a horizontal-gradient image, and a vertical-gradient image. 2. Straight-line segments of edges are extracted from the linked lists generated in step 1. Any straight-line segments longer than an arbitrary threshold (e.g., 30 pixels) are assumed to belong to buildings or other artificial objects. 3. A gradient-filter algorithm is used to test straight-line segments longer than the threshold to determine whether they represent edges of natural or artificial objects. In somewhat oversimplified terms, the test is based on the assumption that the gradient of image intensity varies little along a segment that represents the edge of an artificial object.

Mchenry, Michael↗

Accurate Infrared Line Lists for 20 Isotopologues of CS 2 at Room Temperature

To facilitate atmospheric and spectroscopic studies of carbon disulfide, or CS 2 , in both planetary and exoplanetary atmospheres, we adopt the “Best Theory + Reliable High-resolution Experiment” algorithm to generate semiempirical IR line lists for the 20 most abundant CS 2 isotopologues, denoted as Ames-296K. The IR lists are computed using the Ames-1 potential energy surface, refined using the experimental transition set and an ab initio dipole moment surface fitted from CCSD(T)/aug-cc-pV(T/Q/5+d)Z dipoles extrapolated to a one-particle basis set limit. The IR lists cover the range of 0–10,000 cm -1 , with an S 296K cutoff at 10 -31 cm -1 /molecule·cm -2 (abundance included). A “natural” IR line list at 296 K includes about 10 million lines of the 20 isotopologues, with their intensities scaled by the corresponding abundances. The zero-point energy, partition functions, and abundances are reported for each isotopologue. The energy levels in the global effective Hamiltonian model for 12 C 32 S 2 are adopted to improve the line position accuracy. This new IR list for the main isotopologue is denoted as A+I.296K. Reliable HITRAN2020 line positions are also utilized to improve the accuracy of the 32 S 12 C 34 S, 32 S 12 C 33 S, and 32 S 13 C 32 S isotopologue line lists. The final composite line list is validated against Pacific Northwest National Laboratory experimental cross sections, showing excellent agreement. The agreement supports the quality of the composite line list and the power of synergy between experiment and theory. The new data are proposed for use in updating and expanding the CS 2 data in HITRAN and other high-resolution IR databases. Supplementary files are available in Zenodo and AHED.

Exoplanet atmospheric composition↗

Reconstructed Historical Land Cover and Biophysical Parameters for Studies of Land-Atmosphere Interactions within the Eastern United States

The local environment where we live within the Earth's biosphere is often taken for granted. This environment can vary depending on whether the land cover is a forest, grassland, wetland, water body, bare soil, pastureland, agricultural field, village, residential suburb, or an urban complex with concrete, asphalt, and large buildings. In general, the type and characteristics of land cover influence surface temperatures, sunlight exposure and duration, relative humidity, wind speed and direction, soil moisture amount, plant life, birds, and other wildlife in our backyards. The physical and biological properties (biophysical characteristics) of land cover help to determine our surface environment because they directly affect surface radiation, heat, and soil moisture processes, and also feedback to regional weather and climate. Depending on the spatial scale and land use intensity, land cover changes can have profound impacts on our local and regional environment. Over the past 350 years, the eastern half of the United States, an area extending from the grassland prairies of the Great Plains to the Gulf and Atlantic coasts, has experienced extensive land cover and land use changes that began with land clearing in the 1600s, led to extensive deforestation and intensive land use practices by 1920, and then evolved to the present-day landscape. Determining the consequences of such land cover changes on regional and global climate is a major research issue. Such research requires detailed historical land cover data and modeling experiments simulating historical climates. Given the need to understand the effects of historical land cover changes in the eastern United States, some questions include: - What were the most important land cover transformations and how did they alter biophysical characteristics of the land cover at key points in time since the mid-1600s? - How have land cover and land use changes over the past 350 years affected the land surface environment including surface weather, hydrologic, and climatic variability? - How do the potential effects of regional human-induced land cover change on the environment compare to similar changes that are caused by the natural variations of the Earth's climate system? To help answer these questions, we reconstructed a fractional land cover and biophysical parameter dataset for the eastern United States at 1650, 1850, 1920, and 1992 time-slices. Each land cover fraction is associated with a biophysical parameter class, a suite of parameters defining the biophysical characteristics of that kind of land cover. This new dataset is designed for use in computer models of land-atmosphere interactions, to understand and quantify the effects of historical land cover changes on the water, energy, and carbon cycles

Steyaert, Louis T.↗

Predicting Real-Time Safety of the National Airspace System

To gain the situational awareness necessary for informed decision making regarding avoidance of airspace hazards, each operator must consolidate operations-relevant information from disparate sources and apply extensive domain knowledge to correctly interpret not just the current state of the NAS but forecast its (combined) evolution over the duration of the operation. This time- and workload-intensive process is periodically repeated throughout the operation so that changes can be managed in a timely manner.The imprecision, inaccuracies, inconsistency, and incompleteness of the incoming data further challenges the process. To facilitate informed decision making, this paper presents a model-based framework for the textitautomated real-time monitoring and prediction of possible effects of airspace hazards on the safety of the National Airspace System (NAS). First, hazards to flight are identified and transformed into sms, that is, quantities of interest that could be evaluated based on available data and are predictive of an unsafe event. The sms and associated thresholds that specify when an event transitions from emphsafe to emphunsafe are combined with models of airspace operations and aircraft dynamics. The framework can include any hazard to flight that can be modeled quantitatively. Models can be detailed and complex, or they can be considerably simplifed, as appropriate to the application. Real-time NAS safety monitoring and prediction begins with an estimate of the state of the NAS using the dynamic models. Given the state estimate and a probability distribution of future inputs to the NAS, we can then predict the evolution of the NAS - the future state - and the occurrence of hazards and unsafe events. The entire probability distribution of airspace sms is computed, not just point estimates, without significant assumptions regarding the distribution type andor parameters. We demonstrate our overall approach through a simulated scenario in which we predict the occurrence of some unsafe events and show how these predictions evolve in time as flight operations progress. Predictions accounting for common sources of uncertainty are included and it is shown how the predictions improve in time, become more confident, and change dynamically as new information is made available to the prediction algorithm.

safety↗

A Stochastic Model of Space-Time Variability of Tropical Rainfall: I. Statistics of Spatial Averages

Global maps of rainfall are of great importance in connection with modeling of the earth s climate. Comparison between the maps of rainfall predicted by computer-generated climate models with observation provides a sensitive test for these models. To make such a comparison, one typically needs the total precipitation amount over a large area, which could be hundreds of kilometers in size over extended periods of time of order days or months. This presents a difficult problem since rain varies greatly from place to place as well as in time. Remote sensing methods using ground radar or satellites detect rain over a large area by essentially taking a series of snapshots at infrequent intervals and indirectly deriving the average rain intensity within a collection of pixels , usually several kilometers in size. They measure area average of rain at a particular instant. Rain gauges, on the other hand, record rain accumulation continuously in time but only over a very small area tens of centimeters across, say, the size of a dinner plate. They measure only a time average at a single location. In making use of either method one needs to fill in the gaps in the observation - either the gaps in the area covered or the gaps in time of observation. This involves using statistical models to obtain information about the rain that is missed from what is actually detected. This paper investigates such a statistical model and validates it with rain data collected over the tropical Western Pacific from ship borne radars during TOGA COARE (Tropical Oceans Global Atmosphere Coupled Ocean-Atmosphere Response Experiment). The model incorporates a number of commonly observed features of rain. While rain varies rapidly with location and time, the variability diminishes when averaged over larger areas or longer periods of time. Moreover, rain is patchy in nature - at any instant on the average only a certain fraction of the observed pixels contain rain. The fraction of area covered by rain decreases, as the size of a pixel becomes smaller. This means that within what looks like a patch of rainy area in a coarse resolution view with larger pixel size, one finds clusters of rainy and dry patches when viewed on a finer scale. The model makes definite predictions about how these and other related statistics depend on the pixel size. These predictions were found to agree well with data. In a subsequent second part of the work we plan to test the model with rain gauge data collected during the TRMM (Tropical Rainfall Measuring Mission) ground validation campaign.

Kundu, Prasun K.↗

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗

Enhancements to Linear Stability-Based, CFD-integrated Transition Prediction for High-Speed Flows

Combining linear stability calculations with computational fluid dynamics (CFD) simulations has great potential for the automated modeling of high-speed flows, especially when adequate information about the configuration and the disturbance environment is available. However, a significant impediment to the applicability of this technique is the lack of an efficient method to calculate the crucial amplification ratio corresponding to the onset of transition in hypersonic flows. This ratio, also known as the "transition N-factor," is dependent upon the freestream disturbance environment as well as the surface properties of the test article. In response to the need for an engineering solution to predict the transition N-factor within conventional hypersonic wind tunnels, this paper presents a data-driven correlation that expands the existing correlations from straight circular cones with a narrow range of half angles to a broader array of axisymmetric configurations. Furthermore, when tested against a chosen dataset that was not used in its calibration, the suggested correlation shows good predictive accuracy with an RMS error of only 6.9%. Although similar accuracy may also be achieved via existing correlations based on similar datasets, predictions based on the proposed correlation have the advantage of not requiring an extensive amount of configuration-specific data. Practical applications often have access to the input parameters needed for this correlation, such as the freestream disturbance intensity, Mach number, and body-based slenderness Reynolds number. Additionally, this correlation outperforms the traditional assumption of a constant N-factor, particularly for configurations with blunted nose geometries. The development of this correlation is grounded in an extensive dataset encompassing conical models with body half-angles varying between 5 degrees and 16 degrees, Mach numbers ranging from 5 to 14, and nosetip-based Reynolds numbers approaching the transition reversal limit for blunt-nosed cones.

CFD↗