Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pattern matching”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stochastic Optimization and Uncertainty Quantification of Natrium-based Nuclear-Renewable Energy Systems for Flexible Power Applications in Deregulated Markets

Rapid integration of variable renewable energy sources (VRES) has made modeling and stochastic optimization of hybrid energy systems crucial for studying their long-term performance and viability. However, most studies have focused on just historical data, which may be unreliable for capturing short-term fluctuations, rare events, and long-term patterns of energy demand, price, and the variability of renewable energy sources. For this study, optimal synthetic time series models were developed using Wasserstein distance. The models were validated by comparing the key statistical measures against those of the historical data. They were then used to optimize the integrated Natrium-style advanced energy systems and their long-term (30 years) economics. The stochastic model performs bi-level optimization to find the optimal sizes for the balance of plant and thermal energy storage, while also optimizing energy dispatch to achieve the maximum net present value. In studies of two deregulated markets (California ISO and the Electric Reliability Council of Texas), the integrated Natrium-style system performed better in CAISO than in ERCOT, given higher and more consistent electricity prices during peak-demand periods. The potentially enlarged cost associated with the variable operation and maintenance of the TES system also plays a significant role in driving the system sizing, thus its impacts on the system are investigated in detail through comparison against a baseline case. The study also finds that the bi-level optimization results based on stochastic gradient descent closely match the grid search results. The uncertainty quantification of the stochastic signals provides further NPV-related insights and probability distributions for the case studies. The normal standard error of the mean of NPV for the case with and without TES VOM for CAISO were found to be 7.73M (plus-minus sign) 1.09M USD and 104.99M (plus-minus sign) 1.25M USD, respectively based on a 95% confidence. Given the relatively small NPV variance based on 150 samples, the analysis affords the most robust possible prediction of the techno-economic performance of the integrated Natrium-style energy systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Optimization-based trip chain emulation for electrified ride-sourcing charging demand analyses

Range anxiety remains one of the key concerns for ride-sourcing drivers to adopt battery electric vehicles (BEVs). To investigate the feasibility of using BEVs for ride-sourcing services, we propose an optimization-based methodology to estimate the daily driving trip patterns of ride-sourcing vehicles based on widely available non-identifiable trip data. Furthermore, we investigate the charging needs of electrified ride-sourcing vehicles using agent-based simulation. The methodologies are illustrated through a case study in the city of Chicago. Through sensitivity analysis on driver working hours and initial charging status, we quantify the range of daily average vehicle miles traveled (VMT) per car and identify the hot spots of current public charging demand and potential unsatisfied charging demand. This study can be used to determine the priorities of future charging infrastructure investment to further mitigate range anxiety and promote the adoption of electrified ride-sourcing services

42 ENGINEERING↗
DOE OSTI2020

The SDSS/APOGEE catalogue of HgMn stars

ABSTRACT We report on H-band spectra of chemically peculiar Mercury–Manganese (HgMn) stars obtained via the SDSS/APOGEE survey. As opposed to other varieties of chemically peculiar stars such as classical Ap/Bp stars, HgMn stars lack strong magnetic fields and are defined by extreme overabundances of Mn, Hg, and other heavy elements. A satisfactory explanation for the abundance patterns remains to be determined, but low rotational velocity is a requirement and involvement in binary/multiple systems may be as well. The APOGEE HgMn sample currently consists of 269 stars that were identified among the telluric standard stars as those whose metallic absorption content is limited to or dominated by the H-band Mn ii lines. Due to the fainter magnitudes probed by the APOGEE survey as compared to past studies, only 9/269 stars in the sample were previously known as HgMn types. The 260 newly identified HgMn stars represents a more than doubling of the known sample. At least 32 per cent of the APOGEE sample are found to be binary or multiple systems, and from multi-epoch spectroscopy, we were able to determine orbital solutions for at least one component in 32 binaries. Many of the multilined systems include chemically peculiar companions, with noteworthy examples being the HgMn+Ap/Bp binary HD 5429, the HgMn+HgMn binary HD 298641, and the HgMn+Bp Mn + Am triple system HD 231263. As a further peculiarity, roughly half of the sample produces narrow emission in the C i 16895 Å line, with widths and radial velocities that match those of the Mn ii lines.

Chojnowski, S. Drew↗
DOE OSTI2020

Unstable Invasion of Sedimenting Granular Suspensions

Here we investigate the development of mobility inversion and fingering when a granular suspension is injected radially between horizontal parallel plates of a cell filled with a miscible fluid. While the suspension spreads uniformly when the suspension and the displaced fluid densities are exactly matched, even a small density difference is found to result in a dense granular front which develops fingers with angular spacing that increase with granular volume fraction and decrease with injection rate. We show that the timescale over which the instability develops is given by the volume fraction dependent settling timescale of the grains in the cell. We then show that the mobility inversion and the nonequilibrium Korteweg surface tension due to granular volume fraction gradients determine the number of fingers at the onset of the instability in these miscible suspensions.

02 PETROLEUM↗
DOE OSTI2022

A deep learning-based direct forecasting of CO 2 plume migration

Accurate and timely forecasts of CO 2 plume evolution in geological reservoirs are crucial for CO 2 migration detection, leakage risk assessment, and operation decision support. Conventional forecasting usually adopts a two-step strategy, first calibrating reservoir model parameters against observations using iterative inverse modeling (or history matching) and then applying the calibrated model for predictions. This method impedes real-time forecasts due to the heavy computational demand in inverse modeling and may suffer from poor prediction accuracy because of the limited observation data. In this work, we propose a deep learning-based latent space mapping framework to forecast CO 2 plume migration directly by avoiding the inverse modeling. We first use the convolutional autoencoder to map the high-dimensional complex plume extents onto low-dimensional latent space. Next, we use neural networks to learn the relationship between the observation variables and the prediction latent variables. And then for given observation data, we infer the prediction values directly. This one-step direct forecasting is computationally efficient which requires a few number of parallelizable reservoir simulations and it can provide accurate predictions with limited observations by learning the observation-prediction relationship in the reduced dimension. Therefore, our proposed method enables an in-time forecast of dynamic CO 2 plume distributions. In this work, we demonstrate the effectiveness and accuracy of our method in predicting the CO 2 plume migration using four metrics such as plume area, centroid movement distance, and plume spreading in the primary and secondary directions. And the spatio-temporal evolution patterns of plume migration under diverse geological complexities are also accurately quantified.

15 GEOTHERMAL ENERGY↗
NASA NTRS2018

Chain Formation as a Mechanism for Mass-Independent Fractionation of Sulfur Isotopes in the Archean Atmosphere

The anomalous abundances of sulfur isotopes in ancient sediments provide the strongest evidence for an anoxic atmosphere prior to ∼2.45 Ga, but the mechanism for producing this 'mass-independent' fractionation pattern remains in question. The prevailing hypothesis has been that it is created by differences in the UV photolysis rates of different SO2 isotopologues. We investigate here a recently proposed additional source of fractionation during gas-phase formation of elemental sulfur (S4 and S8). Because two minor S isotopes rarely occur in the same chain, the longer S4 and S8 chains should be strongly, and roughly equally, depleted in all minor isotopes. This gives rise to negative ∆(sup 33) S values and positive ∆(sub 36) S values in elemental sulfur-just the opposite of (and much larger than) what is predicted from SO2 photolysis itself. Back-reactions during chain formation, specifically photolysis of S2 and S3, pass sulfur having the opposite fractionation back to atomic S, and thence to other sulfur species, causing H2S, SO2, sulfate, and short-chain elemental sulfur to have positive ∆(sup 33) S and negative ∆(sup 36) S. Positive ∆(sup 33) S values in elemental sulfur produced in laboratory SO2 photolysis experiments could be caused by the initial fractionation during photolysis, combined with rapid condensation of short-chain sulfur species on the walls of the reaction chamber, along with a scarcity of back-reactions. The simulated fractionations produced by the chain formation mechanism do not directly match fractionations from the rock record. The mismatch might be explained if the isotopic signals leaving the atmosphere were significantly modulated by life, by uncertainties in the rates of reactions of both major and minor isotopic sulfur species, or by the relatively large potential range of atmospheric parameters. Further work is needed to better constrain these uncertainties, but this novel mechanism suggests new avenues to explore in our search for a explanation for the S-MIF record.

Great Oxidation Event↗
DOE OSTI2022

Improving deep learning performance for predicting large-scale geological ${{CO}_{2}}$ sequestration modeling through feature coarsening

Physics-based reservoir simulation for fluid flow in porous media is a numerical simulation method to predict the temporal-spatial patterns of state variables (e.g. pressure p) in porous media, and usually requires prohibitively high computational expense due to its non-linearity and the large number of degrees of freedom (DoF). This work describes a deep learning (DL) workflow to predict the pressure evolution as fluid flows in large-scale 3-dimensional(3D) heterogeneous porous media. In particular, we develop an efficient feature coarsening technique to extract the most representative information and perform the training and prediction of DL at the coarse scale, and further recover the resolution at the fine scale by spatial interpolation. We validate the DL approach to predict pressure field against physics-based simulation data for a field-scale 3D geologic CO 2 sequestration reservoir model. We evaluate the impact of feature coarsening on DL performance, and observe that the feature coarsening not only decreases the training time by >74% and reduces the memory consumption by >75%, but also maintains temporal error 0.63% on average. Besides, the DL workflow provides predictive efficiency with 1406 times speedup compared to physics-based numerical simulation. The key findings from this research significantly improve the training and prediction efficiency of deep learning model to deal with large-scale heterogeneous reservoir models, and thus it can also be further applied to accelerate workflows of history matching and reservoir optimization for close-loop reservoir management.

58 GEOSCIENCES↗
NASA NTRS2018

Chain Formation as a Mechanism for Mass-Independent Fractionation of Sulfur Isotopes in the Archean Atmosphere

The anomalous abundances of sulfur isotopes in ancient sediments provide the strongest evidence for an anoxic atmosphere prior to ∼2.45 Ga, but the mechanism for producing this ‘mass-independent’ fractionation pattern remains in question. The prevailing hypothesis has been that it is created by differences in the UV photolysis rates of different SO2 isotopologues. We investigate here a recently proposed additional source of fractionation during gas-phase formation of elemental sulfur (S4 and S8). Because two minor S isotopes rarely occur in the same chain, the longer S4 and S8 chains should be strongly, and roughly equally, depleted in all minor isotopes. This gives rise to negative Δ(exp 33)S values and positive Δ(exp 33)S values in elemental sulfur—just the opposite of (and much larger than) what is predicted from SO2 photolysis itself. Back-reactions during chain formation, specifically photolysis of S2 and S3, pass sulfur having the opposite fractionation back to atomic S, and thence to other sulfur species, causing H2S, SO2, sulfate, and short-chain elemental sulfur to have positive Δ(exp 33)S and negative Δ(exp 33)S. Positive Δ(exp 33)S values in elemental sulfur produced in laboratory SO2 photolysis experiments could be caused by the initial fractionation during photolysis, combined with rapid condensation of short-chain sulfur species on the walls of the reaction chamber, along with a scarcity of back-reactions. The simulated fractionations produced by the chain formation mechanism do not directly match fractionations from the rock record. The mismatch might be explained if the isotopic signals leaving the atmosphere were significantly modulated by life, by uncertainties in the rates of reactions of both major and minor isotopic sulfur species, or by the relatively large potential range of atmospheric parameters. Further work is needed to better constrain these uncertainties, but this novel mechanism suggests new avenues to explore in our search for a explanation for the S-MIF record.

C. E. Harman↗
NASA NTRS2019

Assessment of Precipitation Anomalies in California Using TRMM and MERRA Data

After more than a decade of moderate seasonal deviations from the expected climate, it is easy to forget that California is actually prone to instabilities in precipitation patterns that occur on various scales. Using modern satellite and reanalysis data we reassess certain aspects of the precipitation climate in California from the past three decades. California has a well-pronounced rain season that peaks in December-February. However, the 95% confidence interval around the climatological precipitation during these months imply that deviations on the order of 60% of the expected amounts are very likely during the most important period of the rain season. While these positive and negative anomalies alternate almost every year and tend to cancel each other, severe multi-year declines of precipitation in California seem to appear on decadal scales. The 1986-1994 decline of precipitation was similar to the current one that started in 2011, and is apparent in the reanalysis data. In terms of accumulated deficits of precipitation, that episode was no less severe than the current one. While El Niño (the warm phase of the El Niño Southern Oscillation, ENSO) is frequently cited as the natural forcing expected to bring a relief, our assessment is that ENSO has been driving at best only 6% of precipitation variability in California in the past three decades. It means El Niño needs to be stronger and longer, in order to have a higher likelihood of a positive impact, and the current one does not match these criteria. Using fractional risk analysis of precipitation populations during normal and dry periods, we show that the likelihood of losing the most intensive precipitation events drastically increases during the multi-year drying events. Since storms delivering up to 50% of precipitation in California are driven by atmospheric rivers making landfall, thus the importance of their suppression and blockage by persistent ridges of atmospheric pressure in the northeast Pacific.

Savtchenko, Andrey K.↗
DOE OSTI2022

Blanco DECam Bulge Survey (BDBS)

The Blanco DECam Bulge Survey (BDBS) has imaged more than 200 square degrees of the southern Galactic bulge, providing photometry in the ugrizy filters for ~250 million unique stars. The presence of a strong foreground disk population, along with complex reddening and extreme image crowding, has made it difficult to constrain the presence of young and intermediate age stars in the bulge population. We employed an accurate cross-match of BDBS with the latest data release (EDR3) from the Gaia mission, matching more than 140 million sources with BDBS photometry and Gaia EDR3 photometry and astrometry. We relied on Gaia EDR3 astrometry, without any photometric selection, to produce clean BDBS bulge colour-magnitude diagrams (CMDs). Gaia parallaxes were used to filter out bright foreground sources, and a Gaussian mixture model fit to Galactic proper motions could identify stars kinematically consistent with bulge membership. We applied this method to 127 different bulge fields of 1 deg 2 each, with |$\ell$| ≤ 9.5 ° and –9.5 ° ≤ b ≤ –2.5 °. The astrometric cleaning procedure removes the majority of blue stars in each field, especially near the Galactic plane, where the ratio of blue to red stars is .10%, increasing to values ~20% at higher Galactic latitudes. We rule out the presence of a widespread population of stars younger than 2 Gyr. The vast majority of blue stars brighter than the turnoff belong to the foreground population, according to their measured astrometry. We introduce the distance between the observed red giant branch bump and the red clump as a simple age proxy for the dominant population in the field, and we confirm the picture of a predominantly old bulge. Further work is needed to apply the method to estimate ages to fields at higher latitudes, and to model the complex morphology of the Galactic bulge. We also produce transverse kinematic maps, recovering expected patterns related to the presence of the bar and of the X-shaped nature of the bulge.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2009

Spatial Cognitive Performance During Adaptation to Conflicting Tilt-Translation Stimuli as a Sensorimotor Spaceflight Analog

The need to resolve new patterns of sensory feedback in altered gravitoinertial environments requires cognitive processes to develop appropriate reference frames for spatial orientation awareness. The purpose of this study was to examine deficits in spatial cognitive performance during adaptation to conflicting tilt-translation stimuli. Fourteen subjects were tilted within a lighted enclosure that simultaneously translated at one of 3 frequencies. Tilt and translation motion was synchronized to maintain the resultant gravitoinertial force aligned with the longitudinal body axis, resulting in a mismatch analogous to spaceflight in which the canals and vision signal tilt while the otoliths do not. Changes in performance on different spatial cognitive tasks were compared 1) without motion, 2) with tilt motion alone (pitch at 0.15, 0.3 and 0.6 Hz or roll at 0.3 Hz), and 3) with conflicting tilt-translation motion. The adaptation paradigm was continued for up to 30 min or until the onset of nausea. The order of the adaptation conditions were counter-balanced across 4 different test sessions. There was a significant effect of stimulus frequency on both motion sickness and spatial cognitive performance. Only 3 of 14 were able to complete the full 30 min protocol at 0.15 Hz, while 7 of 14 completed 0.3 Hz and 13 of 14 completed 0.6 Hz. There were no changes in simple visual-spatial cognitive tests, e.g., mental rotation or match-to-sample. There were significant deficits during 0.15 Hz adaptation in both accuracy and reaction time during a spatial reference task in which subjects are asked to identify a match of a 3D reoriented cube assemblage. Our results are consistent with antidotal reports of cognitive impairment that are common during sensorimotor adaptation with G-transitions. We conclude that these cognitive deficits stem from the ambiguity of spatial reference frames for central processing of inertial motion cues.

Kayanickupuram, A. J.↗
NASA NTRS1989

Infrared spectra of crystalline and glassy silicates and application to interstellar dust

The infrared spectra of crystalline minerals predicted in theoretical condensation sequences do not match the astronomical observations. Since the astronomical spectra are a closer match to glassy silicates, the authors undertook a study to measure the infrared spectra of glassy silicates that have compositions similar to silicate minerals predicted in theoretical condensation sequences. The data should support observations aimed at elucidating condensation chemistry in dust forming regions. The authors measured the mass absorption coefficients, from 2.5 to 25 microns, of ground samples of olivine, diopside, and serpentine and also smoke samples that were prepared from these minerals. The smoke samples prepared in this way are predominantly glassy with nearly the same composition as the parent minerals. The crystalline samples consisted of pure olivine ((Fe(0.1)Mg(0.9))(2)SiO(4)), serpentine, diopside. Sample purity was confirmed by x ray diffraction. Each mineral was ground for 10 hours and a measured mass of the powder was mixed with KBr powder for absorption measurements using the method of Borghesi et a. (1985). The smoke samples were prepared from the same samples used for grinding by vaporizing the minerals using pulsed laser radiation in air. The smoke samples formed by condensation of the resulting vapor. The smoke settled onto infrared transparent KRS-5 substrates and onto a quartz crystal microbalance used to obtain mass measurements. A description of the preparation method is given in Stephens (1980). The glassy diopside showed only diffuse electron diffraction peaks and hence was nearly amorphous, while the serpentine smoke showed a weak diffraction pattern corresponding to MgO. The smoke from olivine showed a weak diffraction pattern corresponding to Fe2O3 and/or Fe3O4. The mass absorption coefficients, from 2.5 to 25 microns, of crystalline diopside, olivine, and serpentine and their corresponding smoke samples are shown in figures.

Stephens, John R.↗
NASA NTRS1992

Slot antenna SIS mixers for submillimeter wavelengths

We are developing improved versions of a slot antenna SIS mixer which we have previously described. The initial work demonstrated a double sideband noise temperature of 420 K for a 500 GHz quasi-optical SIS mixer employing a twin-slot antenna on a quartz dielectric substrate. A quartz hyperhemispherical lens is used to focus the incoming radiation onto the twin-slot antenna. The advantages of a twin-slot antenna include a low impedance (35 omega) and a clean, symmetric beam pattern into the dielectric with a 70 percent efficiency. In our original mixer, the radiation was coupled from the two slot antennas onto superconducting microstrip lines which fed the SIS junction. By performing an impedance transformation using tapered lines and by feeding the radiation from the two slots to the junction in parallel, the effective (real) impedance seen by the junction was reduced to just 4 omega. This very low impedance allowed a junction area of 2.3 sq micron to be used at 500 GHz, which was manufactured using optical lithography. However, no attempt was made to tune out the junction capacitance. We estimate that this capacitance reduces the impedance coupling efficiency to eta(sub Z) approx. equals 0.23, for our junction with omega R(sub N) C = 5.3 at 500 GHz. The recent development of techniques using electron-beam lithography to manufacture junctions with very small areas (approx. equals 0.1 sq microns) now allows considerably more flexibility in the design of SIS mixer circuits. We have redesigned the slot-antenna mixer to take advantage of this possibility. In particular, we have included a novel circuit which allows the junction capacitance to be tune out over a broad bandwidth. For instance, mixers designed for 800 GHz using NbN/MgO/NbN junctions with realistic parameters achieve a 3 dB impedance bandwidth of nearly 400 GHz. Furthermore, our circuit uses only short lengths of microstrip and should be less sensitive to RF losses than other designs. The improved impedance match should give a large reduction in noise temperature as compared to our previous mixer. The new devices are currently under fabrication. Further details of the design and any available experimental results are discussed.

Zmuidzinas, Jonas↗
NASA NTRS2011

Dielectric Covered Planar Antennas at Submillimeter Wavelengths for Terahertz Imaging

Most optical systems require antennas with directive patterns. This means that the physical area of the antenna will be large in terms of the wavelength. When non-cooled systems are used, the losses of microstrip or coplanar waveguide lines impede the use of standard patch or slot antennas for a large number of elements in a phased array format. Traditionally, this problem has been solved by using silicon lenses. However, if an array of such highly directive antennas is to be used for imaging applications, the fabrication of many closely spaced lenses becomes a problem. Moreover, planar antennas are usually fed by microstrip or coplanar waveguides while the mixer or the detector elements (usually Schottky diodes) are coupled in a waveguide environment. The coupling between the antenna and the detector/ mixer can be a fabrication challenge in an imaging array at submillimeter wavelengths. Antennas excited by a waveguide (TE10) mode makes use of dielectric superlayers to increase the directivity. These antennas create a kind of Fabry- Perot cavity between the ground plane and the first layer of dielectric. In reality, the antenna operates as a leaky wave mode where a leaky wave pole propagates along the cavity while it radiates. Thanks to this pole, the directivity of a small antenna is considerably enhanced. The antenna consists of a waveguide feed, which can be coupled to a mixer or detector such as a Schottky diode via a standard probe design. The waveguide is loaded with a double-slot iris to perform an impedance match and to suppress undesired modes that can propagate on the cavity. On top of the slot there is an air cavity and on top, a small portion of a hemispherical lens. The fractional bandwidth of such antennas is around 10 percent, which is good enough for heterodyne imaging applications.The new geometry makes use of a silicon lens instead of dielectric quarter wavelength substrates. This design presents several advantages when used in the submillimeter-wave and terahertz bands: a) Antenna fabrication compatible with lithographic techniques. b) Much simpler fabrication of the lens. c) A simple quarter-wavelength matching layer of the lens will be more efficient if a smaller portion of the lens is used. d) The directivity is given by the lens diameter instead of the leaky pole (the bandwidth will not depend anymore on the directivity but just on the initial cavity). The feed is a standard waveguide, which is compatible with proven Schottky diode mixer/detector technologies. The development of such technology will benefit applications where submillimeter- wave heterodyne array designs are required. The main fields are national security, planetary exploration, and biomedicine. For national security, wideband submillimeter radars could be an effective tool for the standoff detection of hidden weapons or bombs concealed by clothing or packaging. In the field of planetary exploration, wideband submillimeter radars can be used as a spectrometer to detect trace concentrations of chemicals in atmospheres that are too cold to rely on thermal imaging techniques. In biomedicine, an imaging heterodyne system could be helpful in detecting skin diseases.

Chattopadhyay, Goutam↗
NASA NTRS1993

Compressibility effects on rotor forces in the leakage path between a shrouded pump impeller and its housing

A modified approach to Childs' previous work on fluid-structure interaction forces in the leakage path between an impeller shroud and its housing is presented in this paper. Three governing equations consisting of continuity, path-momentum, and circumferential-momentum equations were developed to describe the leakage path inside a pump impeller. Radial displacement perturbations were used to solve for radial and circumferential force coefficients. In addition, impeller-discharge pressure disturbances were used to obtain pressure oscillation responses due to precessing impeller pressure wave pattern. Childs' model was modified from an incompressible model to a compressible barotropic-fluid model (the density of the working fluid is a function of the pressure and a constant temperature only). Results obtained from this model yielded interaction forces for radial and circumferential force coefficients. Radial and circumferential forces define reaction forces within the impeller leakage path. An acoustic model for the same leakage path was also developed. The convective, Coriolis, and centrifugal acceleration terms are removed from the compressible model to obtain the acoustics model. A solution due to impeller discharge pressure disturbances model was also developed for the compressible and acoustics models. The results from these modifications are used to determine what effects additional perturbation terms in the compressible model have on the acoustic model. The results show that the additional fluid mechanics terms in the compressible model cause resonances (peaks) in the force coefficient response curves. However, these peaks only occurred at high values of inlet circumferential velocity ratios greater than 0.7. The peak pressure oscillation was shown to occur at the wearing ring seal. Introduction of impeller discharge disturbances with n = 11 diametral nodes showed that maximum peak pressure oscillations occurred at nondimensional precession frequencies of f = 6.4 and f = 7.8 for this particular pump. Bolleter's results suggest that for peak pressure oscillations to occur at the wearing ring seal, the nondimensional excitation frequency should be on the order of f = 2.182 for n = 11. The resonances found in this research do not match the excitation frequencies predicted by Bolleter. At the predicted peak excitation frequencies given by Bolleter, the compressible model shows an attenuation of the pressure oscillations at the seal exit. The compressibility of the fluid does not have a significant influence on the model at low values of nondimensional excitation frequency. At high values of nondimensional frequency, the effects of compressibility become more significant. For the acoustic analysis, the convective, Coriolis, and centrifugal acceleration terms do affect the results to a limited extent for precession excitation and to a large extent for a pressure excitation when the fluid operates at relatively high Mach numbers.

Cao, Nhai The↗
DOE OSTI2021

New insights into the drainage of inundated ice-wedge polygons using fundamental hydrologic principles

Abstract. The pathways and timing of drainage from the inundated centers of ice-wedge polygons in a warming climate have important implications for carbon flushing, advective heat transport, and transitions from methane to carbon dioxide dominated emissions. Here, we expand on previous research using a recently developed analytical model of drainage from a low-centered polygon. Specifically, we perform (1) a calibration to field data identifying necessary model refinements and (2) a rigorous model sensitivity analysis that expands on previously published indications of polygon drainage characteristics. This research provides intuition on inundated polygon drainage by presenting the first in-depth analysis of drainage within a polygon based on hydrogeological first principles. We verify a recently developed analytical solution of polygon drainage through a calibration to a season of field measurements. Due to the parsimony of the model, providing the potential that it could fail, we identify the minimum necessary refinements that allow the model to match water levels measured in a low-centered polygon. We find that (1) the measured precipitation must be increased by a factor of around 2.2, and (2) the vertical soil hydraulic conductivity must decrease with increasing thaw depth. Model refinement (1) accounts for runoff from rims into the ice-wedge polygon pond during precipitation events and possible rain gauge undercatch, while refinement (2) accounts for the decreasing permeability of deeper soil layers. The calibration to field measurements supports the validity of the model, indicating that it is able to represent ice-wedge polygon drainage dynamics. We then use the analytical solution in non-dimensional form to provide a baseline for the effects of polygon aspect ratios (radius to thaw depth) and coefficient of hydraulic conductivity anisotropy (horizontal to vertical hydraulic conductivity) on drainage pathways and temporal depletion of ponded water from inundated ice-wedge polygon centers. By varying the polygon aspect ratio, we evaluate the relative effect of polygon size (width), inter-annual increases in active-layer thickness, and seasonal increases in thaw depth on drainage. The results of our sensitivity analysis rigorously confirm a previous analysis indicating that most drainage through the active layer occurs along an annular region of the polygon center near the rims. This has important implications for transport of nutrients (such as dissolved organic carbon) and advection of heat towards ice-wedge tops. We also provide a comprehensive investigation of the effect of polygon aspect ratio and anisotropy on drainage timing and patterns, expanding on previously published research. Our results indicate that polygons with large aspect ratios and high anisotropy will have the most distributed drainage, while polygons with large aspect ratios and low anisotropy will have their drainage most focused near their periphery and will drain most slowly. Polygons with small aspect ratios and high anisotropy will drain most quickly. These results, based on parametric investigation of idealized scenarios, provide a baseline for further research considering the geometric and hydraulic complexities of ice-wedge polygons.

Harp, Dylan R. (ORCID:0000000197778000)↗
DOE OSTI2023

Familial Associations of Prevalence and Cause-Specific Mortality for Thoracic Aortic Disease and Bicuspid Aortic Valve in a Large-Population Database

Thoracic aortic disease and bicuspid aortic valve (BAV) likely have a heritable component, but large population-based studies are lacking. This study characterizes familial associations of thoracic aortic disease and BAV, as well as cardiovascular and aortic-specific mortality, among relatives of these individuals in a large-population database. In this observational case-control study of the Utah Population Database, we identified probands with a diagnosis of BAV, thoracic aortic aneurysm, or thoracic aortic dissection. Age- and sex-matched controls (10:1 ratio) were identified for each proband. First-degree relatives, second-degree relatives, and first cousins of probands and controls were identified through linked genealogical information. Cox proportional hazard models were used to quantify the familial associations for each diagnosis. We used a competing-risk model to determine the risk of cardiovascular-specific and aortic-specific mortality for relatives of probands. The study population included 3,812,588 unique individuals. Familial hazard risk of a concordant diagnosis was elevated in the following populations compared with controls: first-degree relatives of patients with BAV (hazard ratio [HR], 6.88 [95% CI, 5.62–8.43]); first-degree relatives of patients with thoracic aortic aneurysm (HR, 5.09 [95% CI, 3.80–6.82]); and first-degree relatives of patients with thoracic aortic dissection (HR, 4.15 [95% CI, 3.25–5.31]). In addition, the risk of aortic dissection was higher in first-degree relatives of patients with BAV (HR, 3.63 [95% CI, 2.68–4.91]) and in first-degree relatives of patients with thoracic aneurysm (HR, 3.89 [95% CI, 2.93–5.18]) compared with controls. Dissection risk was highest in first-degree relatives of patients who carried a diagnosis of both BAV and aneurysm (HR, 6.13 [95% CI, 2.82–13.33]). First-degree relatives of patients with BAV, thoracic aneurysm, or aortic dissection had a higher risk of aortic-specific mortality (HR, 2.83 [95% CI, 2.44–3.29]) compared with controls. Our results indicate that BAV and thoracic aortic disease carry a significant familial association for concordant disease and aortic dissection. The pattern of familiality is consistent with a genetic cause of disease. Furthermore, we observed higher risk of aortic-specific mortality in relatives of individuals with these diagnoses. In conclusion, this study provides supportive evidence for screening in relatives of patients with BAV, thoracic aneurysm, or dissection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

EI_MS_ML

The unambiguous identification of compounds from their electron ionization mass (EI-MS) spectra remains a significant unsolved problem in the field of metabolomics and analytical chemistry as a whole. Typically EI-MS spectra are compared using various mathematical operations that convert the spectral similarity or differences into a distance-like metric that roughly approximates the similarity of any two spectra. A commonly used metric for this is the cosine similarity metric which has values close to one for very similar spectra and a value of zero for very dissimilar spectra; however, no metric is perfect. Due to the prevalence of structurally-similar compounds such as isomers and the prevalence of certain fragmentation patterns across structurally-dissimilar compounds, the unambiguous assignment of EI-MS spectra compounds remains difficult. Frequently, querying an observed EI-MS spectrum against a large database such as the NIST17 library yields multiple possible assignments requiring the end user to distinguish between multiple high scoring hits, or multiple low scoring hits while keeping in mind that the correct hit may not be in the database at all. Although techniques such as orthogonal information from techniques such as chromatography can greatly aid in unambiguous assignment, this also requires more complicated experimental designs and access to more complicated analytical instrumentation. Substructures can be trivially detected and represented as strings using a previously published technique called node coloring from a known chemical structure. However, for experimentally-derived EI-MS spectra this information must be derived from the spectra itself (i.e., because we do not know what compound it represents). To achieve this, the software uses techniques from the field of machine learning and a large training dataset of EI-MS spectra corresponding to known structures annotated with substructure strings, to build models that can predict the presence of a given chemical substructure from an EI-MS spectrum directly.If these predictions are of high-quality (i.e., are unlikely to be false positives), the presence of one or more predicted substructures can be used to constrain the number of possible hits for a query spectrum. Mathematically, this restriction could be expressed in many forms, but the most straight-forward implementation is to weight the cosine similarity of a query spectrum and a plausible database match with a Tanimoto-like coefficient based on the ratio of the number of substructures predicted to the number of substructures present in the potential database hit. Determining which combination of models best reduces assignment ambiguity will be achieved using a combination of manual curation and optimization techniques such as genetic algorithms. This software will perform all the steps necessary to construct said models from a training dataset and evaluate them using a holdout dataset. Various statistical analyses can be performed to determine if this approach does decrease assignment ambiguity. For example, if this approach works, on average, the rank-order of the correct assignment for the holdout set of EI-MS spectra should decrease and the weighted cosine similarities for most of the possible matches in the database should be better than the unweighted cosine similarities. Furthermore, this same pipeline can be used on real experimental data to generate less ambiguous assignments.

Mitchell, Joshua↗