Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric spectrum model”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Studies of optical and chemical properties of aged and fresh biomass burning absorbing aerosols for climate models

Particulate emissions from burning fuels from sub-Saharan Africa is currently poorly characterized. With increasing population and industrialization, emissions are expected to grow. Optical methods for measuring light absorption in the ultraviolet (UV) and visible (vis) regions of the spectrum are crucial for understanding the optical properties of BB aerosols (Bond and Bergstrom, 2006; Lack and Cappa, 2010; Moschos et al., 2021). One such instrument is the dual-spot AE33 model aethalometer is a common filter-based instrument that can be used to measure aerosol absorption properties (Hansen et al., 1984). This instrument has many figures of merit; however, the method's limitations include potential biases from scattering within the filter material, multiple scattering effects of the particles, filter loading (i.e. particle shadowing), and the need for wavelength-specific correction factors (Cλ), which depend on accurate intercomparison with other in situ techniques (Bond et al., 1999; Collaud Coen et al., 2010). This work aims to: 1. Optically characterize BB aerosol generated from African fuel 2. Correct the aethalometer at multiple wavelengths 3. Parametrize these correction factors. 4. Combusting the same fuels with two difference tube furnace systems assess the general applicability of this flexible means of producing BB aerosol. This was accomplished by conducting laboratory studies at Los Alamos National Laboratory (LANL) and North Carolina A&T State University (NCA&T), where aethalometer measurements were compared to in situ optical measurements. Mass cross sections of BB aerosol from these under-studied fuels will be determined along with their wavelength dependence.

09 BIOMASS FUELS↗
DOE OSTI2020

nPINNs: Nonlocal physics-informed neural networks for a parametrized nonlocal universal Laplacian operator. Algorithms and applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius $\delta$ goes to zero, and to the fractional Laplacian as $\delta$ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We also provide theoretical convergence rates with respect to $\delta$ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, $\delta$ and $\alpha$, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order $\alpha(y)$, which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as $\delta$. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

nPINNs: nonlocal Physics-Informed Neural Networks for a parametrized nonlocal universal Laplacian operator. Algorithms and Applications

Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integro-differential equations with sparse, noisy, unstructured, and multifidelity data. PINNs incorporate all available information, including governing equations (reflecting physical laws), initial-boundary conditions, and observations of quantities of interest, into a loss function to be minimized, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal Laplace operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius δ goes to zero, and to the fractional Laplacian as δ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We provide theoretical convergence rates with respect to δ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, δ and α, characterizing the kernel of the unified operator. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon, that is, different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order α(γ), which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as δ. More importantly, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

We study systematics associated with estimating simple stellar population (SSP) parameters – age, metallicity [M/H], α-enhancement [α/Fe], and initial mass function (IMF) shape – and associated M * /L gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs), and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H β , 〈Fe〉, TiO 2SDSS , and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius R e : (a) ages are younger and [α/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about R e /2; (c) this makes M * /L about 2 × larger in the central regions than beyond R e /2. While the models of Conroy et al. return similar [M/H] and [α/Fe] profiles, the age and (hence) M * /L profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs, and S0s: younger ages and less enhanced [α/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred M * /L gradients significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Polyakov model in ’t Hooft flux background: a quantum mechanical reduction with memory

We construct a compactification of Polyakov model on T 2 × $\mathbb{R}$ down to quantum mechanics which remembers non-perturbative aspects of field theory even at an arbitrarily small area. Standard compactification on small T 2 × $\mathbb{R}$ possesses a unique perturbative vacuum (zero magnetic flux state), separated parametrically from higher flux states, and the instanton effects do not survive in the Born-Oppenheimer approximation. By turning on a background magnetic GNO flux in co-weight lattice corresponding to a non-zero ’t Hooft flux, we show that N-degenerate vacua appear at small torus, and there are N - 1 types of flux changing instantons between them. We construct QM instantons starting with QFT instantons using the method of replicas. For example, SU(2) gauge theory with flux reduces to the double-well potential where each well is a fractional flux state. Despite the absence of a mixed anomaly, the vacuum structure of QFT and the one of QM are continuously connected. We also compare the quantum mechanical reduction of the Polyakov model with the deformed Yang-Mills, by coupling both theories to TQFTs. In particular, we compare the mass spectrum for dual photons and energy spectrum in the QM limit. We give a detailed description of critical points at infinity in the semi-classical expansion, and their role in resurgence structure.

nonperturbative effects↗
DOE OSTI2025

SASHIMI-SIDM: semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales

We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses. While the dependence of the collapse time on mass and concentration implies such behaviour, our semi-analytical approach allows us to quantify and illustrate this trend clearly across the full mass spectrum of subhalos, including for subhalo masses below the resolution limit of any current cosmological SIDM simulation. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si.

cosmological simulations↗
DOE OSTI2021

Evaluation of gamma-ray transmission through rectangular collimator slits for application in nuclear fuel spectrometry

Gamma-ray spectrometry is widely applied in several science fields, and in particular in non-destructive gamma scanning and gamma emission tomography of irradiated nuclear fuel. Usually, a collimator is used in the experimental setup, to selectively interrogate a region of interest in the fuel. For the optimization of instrument design, as well as for planning measurement campaigns, predictive models for the transmitted gamma-ray intensity through the collimator are needed. Commonly, Monte Carlo Radiation Transport tools are used for accurate prediction of gamma-ray transport, however, the long computation time requirements when used in low-efficiency experimental setups present challenges. Here, the full-energy peak intensity transmitted through a rectangular collimator slit was examined. A uniform planar surface source emitting isotropically was considered, and the rate of photons reaching an ideal counter plane on the opposite side of the collimator was evaluated by analytical integration. To find a closed-form primitive function, some idealizations were required, and thereby parametric models were obtained for the optical field of view, dependent on slit dimensions (length, height and width) and source-to-collimator distance. For contributions from outside the optical field of view, where a closed-form expression cannot be found, fast numerical integral methods were instead used. The results were validated using the Monte Carlo code MCNP6 and show an agreement within three percent for the numerical method. For the analytical method, deviations up to tens of percent were obtained, which is deemed to still be sufficient for instrument design and measurement planning, where often the order of magnitude of the count rate is not a priori known. The method is planned for use in iterative optimization routines in the design of Gamma Emission Tomography devices, as well as for the prediction of gamma spectra obtained in the planning of fuel inspections. An application of the proposed method was demonstrated in spectrum prediction for a short cooling-time fuel rod test from the Halden reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2000

SAO and Kelvin Waves in the EuroGRIPS GCMS and the UK Meteorological Offices Analyses

This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.

Amodei, M.↗
NASA NTRS2009

Probabilistic Estimation of Critical Flaw Sizes in the Primary Structure Welds of the Ares I-X Launch Vehicle

The primary structure of the Ares I-X Upper Stage Simulator (USS) launch vehicle is constructed of welded mild steel plates. There is some concern over the possibility of structural failure due to welding flaws. It was considered critical to quantify the impact of uncertainties in residual stress, material porosity, applied loads, and material and crack growth properties on the reliability of the welds during its pre-flight and flight. A criterion--an existing maximum size crack at the weld toe must be smaller than the maximum allowable flaw size--was established to estimate the reliability of the welds. A spectrum of maximum allowable flaw sizes was developed for different possible combinations of all of the above listed variables by performing probabilistic crack growth analyses using the ANSYS finite element analysis code in conjunction with the NASGRO crack growth code. Two alternative methods were used to account for residual stresses: (1) The mean residual stress was assumed to be 41 ksi and a limit was set on the net section flow stress during crack propagation. The critical flaw size was determined by parametrically increasing the initial flaw size and detecting if this limit was exceeded during four complete flight cycles, and (2) The mean residual stress was assumed to be 49.6 ksi (the parent material s yield strength) and the net section flow stress limit was ignored. The critical flaw size was determined by parametrically increasing the initial flaw size and detecting if catastrophic crack growth occurred during four complete flight cycles. Both surface-crack models and through-crack models were utilized to characterize cracks in the weld toe.

Pai, Shantaram S.↗
DOE OSTI2025

Plasma phosphorylated tau217 strongly associates with memory deficits in the Alzheimer’s disease spectrum

Abstract Plasma phosphorylated tau (p-tau) biomarkers open unprecedented opportunities for identifying carriers of Alzheimer’s disease pathophysiology in early disease stages using minimally invasive techniques. Plasma p-tau biomarkers are believed to reflect tau phosphorylation and secretion. However, it remains unclear to what extent the magnitude of plasma p-tau abnormalities reflects neuronal network disturbance in the form of cognitive impairment. To address this question, we included 103 cognitively unimpaired elderly and 40 cognitively impaired, amyloid-β-positive individuals from the TRIAD cohort, in addition to 336 cognitively unimpaired and 216 cognitively impaired, amyloid-β-positive older adults from the BioFINDER-2 cohort. Participants had tau PET scans, amyloid PET scans or amyloid CSF, p-tau217, p-tau181 and p-tau231 blood measures, structural T1-MRI and cognitive assessments. In this cross-sectional study, we used regression models and correlation analyses to assess the relationship between plasma biomarkers and cognitive scores. Furthermore, we applied receiver operating characteristic curves to assess cognitive impairment across plasma biomarkers. Finally, we categorized participants into amyloid (A), p-tau (T1) and tau PET (T2) positive (+) or negative (−) profiles and ran non-parametric comparisons to assess differences across cognitive domains. We found that plasma p-tau217 was more associated with cognitive performance than p-tau181 and p-tau231 and that this relationship was particularly strong for memory scores (TRIAD: βp-tau217 = −0.53, βp-tau181 = −0.35 and βp-tau231 = −0.24; BioFINDER-2: βp-tau217 = −0.52, βp-tau181 = −0.24 and βp-tau231 = −0.29). Associations in amyloid-β-positive participants resembled these results, but other cognitive scores also showed strong associations in cognitively impaired individuals. Moreover, plasma p-tau217 outperformed plasma p-tau181 and plasma p-tau231 in identifying memory impairment (area under the curve values for TRIAD: p-tau217 = 0.86, p-tau181 = 0.77 and p-tau231 = 0.75; and for BioFINDER-2: p-tau217 = 0.86, p-tau181 = 0.76 and p-tau231 = 0.81) and in identifying executive function impairment only in the BioFINDER-2 cohort (p-tau217 = 0.82, p-tau181 = 0.76 and p-tau231 = 0.76). Lastly, we showed that subtle memory deficits were present in A+T1+T2− participants for plasma p-tau217 (P = 0.007) and plasma p-tau181 (P = 0.01) in the TRIAD cohort and for all biomarkers across cognitive domains in A+T1+T2− and A+T1+T2− individuals (P < 0.001 in all) in the BioFINDER-2 cohort. The A+T1+T2− individuals showed cognitive deficits in both cohorts (P < 0.001 in all). Together, our results suggest that plasma p-tau217 stands out as a biomarker capable of identifying memory deficits attributable to Alzheimer’s disease and that memory impairment certainly occurs in amyloid-β- and plasma p-tau-positive individuals who have no significant amounts of tau in the neocortex.

Neurosciences & Neurology↗
NASA NTRS2011

Saturn Ring Data Analysis and Thermal Modeling

CIRS, VIMS, UVIS, and ISS (Cassini's Composite Infrared Specrtometer, Visual and Infrared Mapping Spectrometer, Ultra Violet Imaging Spectrometer and Imaging Science Subsystem, respectively), have each operated in a multidimensional observation space and have acquired scans of the lit and unlit rings at multiple phase angles. To better understand physical and dynamical ring particle parametric dependence, we co-registered profiles from these three instruments, taken at a wide range of wavelengths, from ultraviolet through the thermal infrared, to associate changes in ring particle temperature with changes in observed brightness, specifically with albedos inferred by ISS, UVIS and VIMS. We work in a parameter space where the solar elevation range is constrained to 12 deg - 14 deg and the chosen radial region is the B3 region of the B ring; this region is the most optically thick region in Saturn's rings. From this compilation of multiple wavelength data, we construct and fit phase curves and color ratios using independent dynamical thermal models for ring structure and overplot Saturn, Saturn ring, and Solar spectra. Analysis of phase curve construction and color ratios reveals thermal emission to fall within the extrema of the ISS bandwidth and a geometrical dependence of reddening on phase angle, respectively. Analysis of spectra reveals Cassini CIRS Saturn spectra dominate Cassini CIRS B3 Ring Spectra from 19 to 1000 microns, while Earth-based B Ring Spectrum dominates Earth-based Saturn Spectrum from 0.4 to 4 microns. From our fits we test out dynamical thermal models; from the phase curves we derive ring albedos and non-lambertian properties of the ring particle surfaces; and from the color ratios we examine multiple scattering within the regolith of ring particles.

Saturn's rings↗
NASA NTRS2002

A Parametric Study of Fine-scale Turbulence Mixing Noise

The present paper is a study of aerodynamic noise spectra from model functions that describe the source. The study is motivated by the need to improve the spectral shape of the MGBK jet noise prediction methodology at high frequency. The predicted spectral shape usually appears less broadband than measurements and faster decaying at high frequency. Theoretical representation of the source is based on Lilley's equation. Numerical simulations of high-speed subsonic jets as well as some recent turbulence measurements reveal a number of interesting statistical properties of turbulence correlation functions that may have a bearing on radiated noise. These studies indicate that an exponential spatial function may be a more appropriate representation of a two-point correlation compared to its Gaussian counterpart. The effect of source non-compactness on spectral shape is discussed. It is shown that source non-compactness could well be the differentiating factor between the Gaussian and exponential model functions. In particular, the fall-off of the noise spectra at high frequency is studied and it is shown that a non-compact source with an exponential model function results in a broader spectrum and better agreement with data. An alternate source model that represents the source as a covariance of the convective derivative of fine-scale turbulence kinetic energy is also examined.

Khavaran, Abbas↗
DOE OSTI2025

Atacama Cosmology Telescope: Multiprobe cosmology with unWISE galaxies and ACT DR6 CMB lensing

We present a joint analysis of the cosmic microwave background (CMB) lensing power spectra measured from the Data Release 6 of the Atacama Cosmology Telescope (ACT) and Planck PR4, cross-correlations between the ACT and Planck lensing reconstruction and galaxy clustering from unWISE, and the unWISE clustering auto-spectrum. We obtain 1.5% constraints on the matter density fluctuations at late times parametrized by the best constrained parameter combination 𝑆$^{3⁢x⁢2⁢pt}_{8}$ ≡ 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.4 = 0.815 ± 0.012. The commonly used 𝑆 8 ≡ 𝜎 8⁢ (Ω 𝑚 /0.3) 0.5 parameter is constrained to 𝑆 8 = 0.816 ± 0.015. In combination with baryon acoustic oscillation (BAO) measurements we find 𝜎 8 = 0.815 ± 0.012. We also present sound-horizon-independent estimates of the present day Hubble rate of 𝐻 0 = 66.4$^{+3.2}_{−3.7}$ km s −1 Mpc −1 from our large scale structure data alone and 𝐻 0 = 64.3$^{+2.1}_{−2.4}$ km s −1 Mpc −1 in combination with uncalibrated supernovae from Pantheon+. Using parametric estimates of the evolution of matter density fluctuations, we place constraints on cosmic structure in a range of high redshifts typically inaccessible with cross-correlation analyses. Combining lensing cross- and autocorrelations, we derive a 3.3% constraint on the integrated matter density fluctuations above 𝑧 = 2.4, one of the tightest constraints in this redshift range and fully consistent with a Λ cold dark matter (Λ⁢CDM) model fit to the primary CMB from Planck. Finally, combining with primary CMB observations and using the extended low redshift coverage of these combined datasets we derive constraints on a variety of extensions to the Λ⁢CDM model including massive neutrinos, spatial curvature, and dark energy. We find in flat Λ⁢CDM⁢ ∑𝑚 𝜈 < 0.12 eV at 95% confidence using the large scale structure data, BAO measurements from Sloan Digital Sky Survey, and primary CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

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

Multichannel spectral estimation in acoustics: A state-space approach

Spectral estimation is a necessary methodology to analyze the frequency content of noisy data sets especially in acoustic applications. Many spectral techniques have evolved starting with the classical Fourier transform methods based on the well-known Wiener-Khintchine relationship relating the covariance-to-spectral density as a transform pair culminating with more elegant model-based parametric techniques that apply prior knowledge of the data to produce a high-resolution spectral estimate. Multichannel spectral representations are a class of both nonparametric, as well as parametric, estimators that provide improved spectral estimates. In any case, classical nonparametric multichannel techniques can provide reasonable estimates when coupled with peak-peaking methods as long as the signal levels are reasonably high. Furthermore, parametric multichannel methods can perform quite well in low signal level environments even when applying simple peak-picking techniques. In this paper, the performance of both nonparametric (periodogram) and parametric (state-space) multichannel spectral estimation methods are investigated when applied to both synthesized noisy structural vibration data as well as data obtained from a sounding rocket flight. It is demonstrated that for the multichannel problem, state-space techniques provide improved performance, offering a parametric alternative compared to classical methods.

42 ENGINEERING↗
NASA NTRS2010

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗
NASA NTRS1992

The kinematics of buoyancy induced mixing

A model problem to study mixing driven by buoyancy flow fields due to a time dependent body force which consists of a steady and an oscillatory component is considered. Flow fields generated by a time dependent body force, which consists of a steady and an oscillatory component from vibrations or g-jitter, mix two fluids inside a cavity by stretching and folding the interface into various morphological patterns. These patterns show two basic structures depending on the dominant components of the body force. These structures are whirls and tendrils, corresponding to the dominance of either the steady or oscillatory component, respectively. A combination of these structures also occurs in the interaction region. In this region the effect of the steady component opposes the oscillatory component and the result is to smooth the interface, and prevent breakup of the interface. For the parametric range considered, two basic flow regimes occur: convective and chaotic. In the convective regime, the morphological patterns are topological. In the chaotic regime, mixing occurs by a repetitive sequence of bubble formation, necking, and breakup which spreads the concentration field throughout the cavity. The length stretch of the interface increases exponentially with time. The dynamical system characteristics show that in the convective regime, at a point in space, the flow field is oscillatory with elliptical phase space trajectories, and its power spectrum indicates that the flow field responds to the corresponding input frequency. Whereas, in the chaotic regime, the flow field is aperiodic, the phase space trajectories show irregular patterns, and the power spectrum shows a broadband distribution.

Duval, Walter M. B.↗
NASA NTRS2000

Marangoni Effects on Near-Bubble Microscale Transport During Boiling of Binary Fluid Mixtures

In earlier investigations, Marangoni effects were observed to be the dominant mechanism of boiling transport in 2-propanol/water mixtures under reduced gravity conditions. In this investigation we have examined the mechanisms of binary mixture boiling by exploring the transport near a single bubble generated in a binary mixture between a heated surface and cold surface. The temperature field created in the liquid around the bubble produces vaporization over the portion of its interface near the heated surface and condensation over portions of its interface near the cold surface. Experiments were conducted using different mixtures of water and 2-propanol under 1g conditions and under reduced gravity conditions aboard the KC135 aircraft. Since 2-propanol is more volatile than water, there is a lower concentration of 2-propanol near the hot surface and a higher concentration of 2-propanol near the cold plate relative to the bulk quantity. This difference in interface concentration gives rise to strong Marangoni effects that move liquid toward the hot plate in the near bubble region for 2-propanol and water mixtures. In the experiments in this study, the pressure of the test system was maintained at about 5 kPa to achieve the full spectrum of boiling behavior (nucleate boiling, critical heat flux and film boiling) at low temperature and heat flux levels. Heat transfer data and visual documentation of the bubble shape were extracted from the experimental results. In the 1-g experiments at moderate to high heat flux levels, the bubble was observed to grow into a mushroom shape with a larger top portion near the cold plate due to the buoyancy effect. The shape of the bubble was somewhat affected by the cold plate subcooling and the superheat of the heated surface. At low superheat levels for the heated surface, several active nucleation sites were observed, and the vapor stems from them merged to form a larger bubble. The generation rate of vapor is moderate in this regime and the bubble shape is cylindrical in appearance. In some instances, the bubble interface appeared to oscillate. At higher applied heat flux levels, the top of the bubble became larger, apparently to provide more condensing interface area adjacent to the cold plate. Increasing the applied heat flux ultimately led to dry-out of the heated surface, with conditions just prior to dryout corresponding to the maximum heat flux (CHF). A more stable bubble was observed when the system attained the minimum heat flux (for film boiling). In this regime, most of the surface under the bottom of the bubble was dry with nucleate boiling sometimes occuring around the contact perimeter of the bubble at heated surface. Different variations (e.g. gap between two plates, molar concentration of the liquid mixture) of the experiments were examined to determine parametric effects on the boiling process and to determine the best conditions for the KC135 reduced gravity tests. Variation of the gap was found to have a minor impact on the CHF. However, reducing the gap between the hot and cold surface was observed to significantly reduce the minimum heat flux for fixed molar concentration of 2-propanol. In the reduced gravity experiments aboard the KC135 aircraft, the bubble formed in the 6.4 mm gap was generally cylindrical or barrel shaped and it increased its extent laterally as the surface superheat increased. In reduced gravity experiments, dryout of the heated surface under the bubble was observed to occur at a lower superheated temperature than for 1g conditions. Observed features of the boiling process and heat transfer data under reduced gravity will be discussed in detail. The results of the reduced gravity experiments will also be compared to those obtained in comparable 1g experiments. In tandem with the experiments we are also developing a computational model of the transport in the liquid surrounding the bubble during the boiling process. The computational model uses a level set method to model motion of the interface. It will incorporate a macroscale treatment of the transport in the liquid gap between the surfaces and a microscale treatment of transport in the regions between the bubble interface and the solid surfaces. The features of the model will be described in detail. Future research directions suggested by the results to date will also be discussed.

V. Carey↗