Search NASA⌕ Search

SEARCH · Search NASA

Results for “fully linear 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 577 records · Page 32

Testbed Results for Scalar and Vector Radiative Transfer Computations of Light in Atmosphere-Ocean Systems

We generate and tabulate reflectance values of the Stokes parameters I, Q, and U of upwelling radiance just above a rough ocean surface and at the top of the atmosphere (TOA) for 100 scattering geometries, four atmosphere-ocean systems, and four wavelengths. The atmosphere-ocean systems increase in complexity from (a) a molecular atmosphere above a rough ocean surface (AOS-I model); to (b) a pure water body below a rough ocean surface (AOS-II model); to (c) a fully-coupled simple atmosphere-ocean system (AOS-III model) containing a molecular atmosphere, rough ocean surface, and pure water; to (d) a fully-coupled complex atmosphere-ocean system (AOS-IV model) that includes scattering by molecules, rough ocean surface, pure water, and hydrosols. Our wavelengths (350, 450, 550, and 650 nm) capture the ultraviolet-visible range. Our tables provide radiative transfer (RT) testbed results for atmosphere-ocean systems with an accuracy that surpasses the measurement accuracy of state-of-the-art polarimeters. To validate the accuracy of these tables we performed computations using three independent RT codes that provide deterministic numerical solutions for the RT equation. The agreement is 10(exp –5) for AOS-IV model, and 10(exp –6) for the other models. The degree of linear polarization computed by these RT codes differs by ≤0.2% for 15 isolated cases of tabulated reflectance values, and by ≤0.1% for all remaining cases. We also provide comparisons with results obtained by a stochastic RT code for AOS-I model. The agreement between the deterministic and stochastic results for this model is 10(exp –5) at TOA, and 10(exp –6) above the ocean surface.

Polarimeters↗

Seasonality of the Migrating Semidiurnal Tide in the Tropical Upper Mesosphere and Lower Thermosphere and its Thermodynamic and Momentum Budget

This work uses the Specified Dynamics-Whole Atmosphere Community Climate Model with Ionosphere/Thermosphere eXtension (SD-WACCM-X) to determine and explain the seasonality of the migrating semidiurnal tide (SW2) components of tropical upper mesosphere and lower thermosphere (UMLT) temperature, zonal wind, and meridional wind. This work also quantifies aliasing due to SW2 in satellite-based tidal estimates. Results show that during equinox seasons, the vertical profiles of tropical UMLT temperature SW2 and zonal-wind SW2’s amplitudes have a double-peak structure while they, along with meridional-wind SW2, have a single-peak structure in their amplitudes in June solstice. Hough mode reconstruction reveals that a linear combination of five SW2 Hough modes cannot fully reproduce these features. Tendency analysis reveals that for temperature, the adiabatic term, nonlinear advection term, and linear advection term are important. For the winds, the classical terms, nonlinear advection term, linear advection term, and gravity wave drag are important. Results of our alias analysis then indicate that SW2 can induce an ∼60% alias in zonal-mean and DW1 components calculated from sampling like that of the Thermosphere–Ionosphere–Mesosphere Energetics and Dynamics satellite and the Aura satellite. This work concludes that in situ generation by wave–wave interaction and/or by gravity waves plays significant roles in the seasonality of tropical UMLT temperature SW2, zonal-wind SW2, and meridional-wind SW2. The alias analysis further adds that one cannot simply assume that SW2 in the tropical UMLT is negligible.

Tides↗

Poster — Thur Eve — 21: Off-axis dose perturbation effects in water in a 5 × 5 cm{sup 2} 18 MV photon beam for the PTW microLion and Exradin A1SL ionization chambers

A PTW microLion liquid ionization chamber and an Exradin A1SL air-filled ionization chamber have been modeled using the egs-chamber user code of the EGSnrc system to determine their perturbation effects in water in a 5 × 5 cm{sup 2} 18 MV photon beam. A model of the Varian CL21EX linear accelerator was constructed using the BEAMnrc Monte Carlo code, and was validated by comparing measured PDDs and profiles from the microLion and A1SL chambers to calculated results that included chamber models. Measured PDDs for a 5 × 5 cm{sup 2} field for the microLion chamber agreed with calculations to within 1.5% beyond a depth of 0.5 cm, and the A1SL PDDs agreed within 1.0% beyond 1.0 cm. Measured and calculated profiles at 10 cm depth agreed within 1.0% for both chambers inside the field, and within 4.0% near the field edge. Local percent differences increased up to 15% at 4 cm outside the field. The ratio of dose to water in the absence of the chamber relative to dose in the chamber's active volume as a function of off-axis distance was calculated using the egs-chamber correlated sampling technique. The dose ratio was nearly constant inside the field and consistent with the stopping power ratios of water to detector material, but varied up to 3.3% near the field edge and 5.2% at 4 cm outside the field. Once these perturbation effects are fully characterized for more field sizes and detectors, they could be applied to clinical water tank measurements for improved dosimetric accuracy.

07 ISOTOPE AND RADIATION SOURCES↗

Thermal traction contact performance evaluation under fully flooded and starved conditions

Ultra high speed traction tests were performed on two traction fluids commonly employed. Traction data on these fluids is required for purposes of traction drive design optimization techniques. To obtain the traction data, an existing twin disc traction test machine was employed. This machine was modified to accommodate the range of test variables. All the data reported was obtained under conditions of side slip, a technique whereby only low power levels are required to simulate real traction drive contacts. Theoretical traction predictions were performed for a representative number of curves that showed the influence of rolling velocity, of contact pressure and of aspect ratio. To establish the accuracy of the thermal model the predictions were performed ith increasing levels of independence of experimentally determined parameters. In the final resulting prediction only two non linear thermal parameters were used for the prediction of 15 different traction curves covering the entire range of variables as used in the investigation, with the exception of the influence of asperity traction. Comparison of these theoretical curves and corresponding experimental traces show very good agreement.

Tevaarwerk, J. L.↗

Simulating the role of surface forcing on observed multidecadal upper-ocean salinity changes

The ocean’s surface salinity field has changed over the observed record, driven by an intensification of the water cycle in response to global warming. However, the origin and causes of the coincident subsurface salinity changes are not fully understood. The relationship between imposed surface salinity and temperature changes and their corresponding subsurface changes is investigated using idealized ocean model experiments. The ocean’s surface has warmed by about 0.5°C (50 yr) –1 while the surface salinity pattern has amplified by about 8% per 50 years. The idealized experiments are constructed for a 50-yr period, allowing a qualitative comparison to the observed salinity and temperature changes previously reported. The comparison suggests that changes in both modeled surface salinity and temperature are required to replicate the three-dimensional pattern of observed salinity change. The results also show that the effects of surface changes in temperature and salinity act linearly on the changes in subsurface salinity. In addition, surface salinity pattern amplification appears to be the leading driver of subsurface salinity change on depth surfaces; however, surface warming is also required to replicate the observed patterns of change on density surfaces. This is the result of isopycnal migration modified by the ocean surface warming, which produces significant salinity changes on density surfaces.

54 ENVIRONMENTAL SCIENCES↗

Targeted brain proteomics uncover multiple pathways to Alzheimer's dementia

Objective Previous gene expression analysis identified a network of coexpressed genes that is associated with β‐amyloid neuropathology and cognitive decline in older adults. The current work targeted influential genes in this network with quantitative proteomics to identify potential novel therapeutic targets. Methods Data came from 834 community‐based older persons who were followed annually, died, and underwent brain autopsy. Uniform structured postmortem evaluations assessed the burden of β‐amyloid and other common age‐related neuropathologies. Selected reaction monitoring quantified cortical protein abundance of 12 genes prioritized from a molecular network of aging human brain that is implicated in Alzheimer's dementia. Regression and linear mixed models examined the protein associations with β‐amyloid load and other neuropathological indices as well as cognitive decline over multiple years preceding death. Results Average age at death was 88.6 years. Overall, 349 participants (41.9%) had Alzheimer's dementia at death. A higher level of PLXNB1 abundance was associated with more β‐amyloid load ( p = 1.0 × 10 −7 ) and higher PHFtau tangle density ( p = 2.3 × 10 −7 ), and the association of PLXNB1 with cognitive decline is mediated by these known Alzheimer's disease pathologies. On the other hand, higher IGFBP5, HSPB2, and AK4 and lower ITPK1 levels were associated with faster cognitive decline, and, unlike PLXNB1, these associations were not fully explained by common neuropathological indices, suggesting novel mechanisms leading to cognitive decline. Interpretation Using targeted proteomics, this work identified cortical proteins involved in Alzheimer's dementia and begins to dissect two different molecular pathways: one affecting β‐amyloid deposition and another affecting resilience without a known pathological footprint. Ann Neurol 2018;83:78–88

Yu, Lei↗

The impact of wind, solar, and other factors on the decline in wholesale power prices in the United States

Across multiple organized wholesale power markets in the United States, annual average prices declined by 19–64 dollars per MWh between 2008 and 2017 while retirements of thermal power plants accelerated. Several prominent changes over the last decade are often discussed as contributors to this decline in prices. These include growth in wind and solar, a reduction in the price of natural gas, and weakened load growth. Here we construct a fundamental supply curve model for each of seven organized wholesale market regions and use counterfactual simulations to assess the degree to which wind and solar—among other factors—have influenced wholesale electricity prices. We find that growth in wind and solar since 2008 reduced average annual wholesale electricity prices by less than 3 dollars per MWh. In contrast the decline in natural gas prices reduced wholesale prices by 7–53 dollars per MWh, depending on the region. This suggests that recent thermal-plant retirements in the U.S. are primarily due to low natural gas prices, not growth in wind and solar. Fully isolating the impact of individual factors, however, is limited by non-linear interactions between factors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Conservative tightly-coupled simulations of stochastic multiscale systems

Multiphysics problems often involve components whose macroscopic dynamics is driven by microscopic random fluctuations. The fidelity of simulations of such systems depends on their ability to propagate these random fluctuations throughout a computational domain, including subdomains represented by deterministic solvers. When the constituent processes take place in nonoverlapping subdomains, system behavior can be modeled via a domain-decomposition approach that couples separate components at the interfaces between these subdomains. Its coupling algorithm has to maintain a stable and efficient numerical time integration even at high noise strength. We propose a conservative domain-decomposition algorithm in which tight coupling is achieved by employing either Picard's or Newton's iterative method. Coupled diffusion equations, one of which has a Gaussian white-noise source term, provide a computational testbed for analysis of these two coupling strategies. Fully-converged (“implicit”) coupling with Newton's method typically outperforms its Picard counterpart, especially at high noise levels. This is because the number of Newton iterations scales linearly with the amplitude of the Gaussian noise, while the number of Picard iterations can scale superlinearly. At large time intervals between two subsequent inter-solver communications, the solution error for single-iteration (“explicit”) Picard's coupling can be several orders of magnitude higher than that for implicit coupling. Increasing the explicit coupling's communication frequency reduces this difference, but the resulting increase in computational cost can make it less efficient than implicit coupling at similar levels of solution error, depending on the communication frequency of the latter and the noise strength. This trend carries over into higher dimensions, although at high noise strength explicit coupling may be the only computationally viable option.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Polarimetric Radar Variables in Snowfall at Ka- and W-Band Frequency Bands: A Comparative Analysis

Dual-frequency millimeter-wavelength radar observations in snowfall are analyzed in order to evaluate differences in conventional polarimetric radar variables such as differential reflectivity (ZDR) specific differential phase shift (K DP ) and linear depolarization ratio (LDR) at traditional cloud radar frequencies at Ka and W bands (~35 and ~94 GHz, correspondingly). Low radar beam elevation (~5°) measurements were performed at Oliktok Point, Alaska, with a scanning fully polarimetric radar operating in the horizontal–vertical polarization basis. This radar has the same gate spacing and very close beam widths at both frequencies, which largely alleviates uncertainties associated with spatial and temporal data matching. It is shown that observed Ka- and W-band Z DR differences are, on average, less than about 0.5 dB and do not have a pronounced trend as a function of snowfall reflectivity. The observed Z DR differences agree well with modeling results obtained using integration over nonspherical ice particle size distributions. For higher signal-to-noise ratios, KDP data derived from differential phase measurements are approximately scaled as reciprocals of corresponding radar frequencies indicating that the influence of non-Rayleigh scattering effects on this variable is rather limited. This result is also in satisfactory agreement with data obtained by modeling using realistic particle size distributions. Observed Ka- and W-band LDR differences are strongly affected by the radar hardware system polarization “leak” and are generally less than 4 dB. Smaller differences are observed for higher depolarizations, where the polarization “leak” is less pronounced. Realistic assumptions about particle canting and the system polarization isolation lead to modeling results that satisfactorily agree with observational dual-frequency LDR data.

42 ENGINEERING↗

Intensity Biased PSP Measurement

The current pressure sensitive paint (PSP) technique assumes a linear relationship (Stern-Volmer Equation) between intensity ratio (I(sub o)/I) and pressure ratio (P/P(sub o)) over a wide range of pressures (vacuum to ambient or higher). Although this may be valid for some PSPs, in most PSPs the relationship is nonlinear, particularly at low pressures (less than 0.2 psia when the oxygen level is low). This non-linearity can be attributed to variations in the oxygen quenching (de-activation) rates (which otherwise is assumed constant) at these pressures. Other studies suggest that some paints also have non-linear calibrations at high pressures; because of heterogeneous (non-uniform) oxygen diffusion and quenching. Moreover, pressure sensitive paints require correction for the output intensity due to light intensity variation, paint coating variation, model dynamics, wind-off reference pressure variation, and temperature sensitivity. Therefore to minimize the measurement uncertainties due to these causes, an insitu intensity correction method was developed. A non-oxygen quenched paint (which provides a constant intensity at all pressures, called non-pressure sensitive paint, NPSP) was used for the reference intensity (I(sub NPSP) with respect to which all the PSP intensities (I) were measured. The results of this study show that in order to fully reap the benefits of this technique, a totally oxygen impermeable NPSP must be available.

Subramanian, Chelakara S.↗

Introductory guide to backgrounds in XPS spectra and their impact on determining peak intensities

Photoelectron and Auger peaks are central to most of the important uses of x-ray photoelectron spectroscopy (XPS), and thus, they receive the most attention in many types of analysis. Quantitative chemical analysis using XPS requires the assessment of the intensities of the photoemission peaks of the elements detected. Determination of peak intensities requires separation of the photoelectron peak signals from the background on which the peaks rest. For the determination of peak area intensities, the background is subtracted from overall signal intensity. The spectral background is also critical when peak fitting is used to determine intensities of overlapping peaks, and the model of background used in this process can impact the results. In addition to the impact on quantitative analysis, information about the depth distribution of elements in the near surface region can often be obtained by visual inspection of the background and quantified using appropriate modeling. Here, this introductory guide provides some basic information about backgrounds in the XPS analysis, describes the types of background models that are commonly used, suggests some of their strengths and weaknesses, and provides examples of their use and misuse. Although the fundamental nature of some components of the background signals in XPS is not understood, indicating that none of the models in use are fully correct and the area is subject to active research, appropriate good practices have been established for most routine analysis. The guide describes good practices, identifies errors that frequently appear in the literature, and uses examples to demonstrate the impacts of background selections on determinations of peak intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intensity Biased PSP Measurement

The current pressure sensitive paint (PSP) technique assumes a linear relationship (Stern-Volmer Equation) between intensity ratio (I(sub 0)/I) and pressure ratio (P/P(sub 0)) over a wide range of pressures (vacuum to ambient or higher). Although this may be valid for some PSPs, in most PSPs the relationship is nonlinear, particularly at low pressures (less than 0.2 psia when the oxygen level is low). This non-linearity can be attributed to variations in the oxygen quenching (de-activation) rates (which otherwise is assumed constant) at these pressures. Other studies suggest that some paints also have non-linear calibrations at high pressures; because of heterogeneous (non-uniform) oxygen diffusion and c quenching. Moreover, pressure sensitive paints require correction for the output intensity due to light intensity variation, paint coating variation, model dynamics, wind-off reference pressure variation, and temperature sensitivity. Therefore to minimize the measurement uncertainties due to these causes, an in- situ intensity correction method was developed. A non-oxygen quenched paint (which provides a constant intensity at all pressures, called non-pressure sensitive paint, NPSP) was used for the reference intensity (I(sub NPSP)) with respect to which all the PSP intensities (I) were measured. The results of this study show that in order to fully reap the benefits of this technique, a totally oxygen impermeable NPSP must be available.

Subramanian, Chelakara S.↗

Viscous instability triggered by layered accretion in protoplanetary disks

Layered accretion is one of the inevitable ingredients in protoplanetary disks when disk turbulence is excited by magnetorotational instabilities (MRIs). In the accretion, disk surfaces where MRIs fully operate have a high value of disk accretion rate ( M-dot ), while the disk midplane where MRIs are generally quenched ends up with a low value of M-dot . Significant progress on understanding MRIs has recently been made by a number of dedicated MHD simulations, which requires improvement of the classical treatment of α in 1D disk models. To this end, we obtain a new expression of α by utilizing an empirical formula that is derived from recent MHD simulations of stratified disks with ohmic diffusion. It is interesting that this new formulation can be regarded as a general extension of the classical α. Armed with the new α, we perform a linear stability analysis of protoplanetary disks that undergo layered accretion, and we find that a viscous instability can occur around the outer edge of dead zones. Disks become stable in using the classical α. We identify that the difference arises from Σ-dependence of M-dot ; whereas Σ is uniquely determined for a given value of M-dot in the classical approach, the new approach leads to M-dot that is a multivalued function of Σ. We confirm our finding both by exploring a parameter space and by performing the 1D, viscous evolution of disks. We finally discuss other nonideal MHD effects that are not included in our analysis but may affect our results.

79 ASTRONOMY AND ASTROPHYSICS↗

Projective multiscale time-integration for electrostatic particle-in-cell methods

The simulation of problems in kinetic plasma physics is often challenging due to strongly coupled phenomena across multiple scales. In this work, we propose a wavelet-based coarse-grained numerical scheme, based on the framework of Equation-Free Projective Integration, for a kinetic plasma system modeled by the Vlasov–Poisson equations. A kinetic particle-in-cell (PIC) code is used to simulate the meso scale dynamics for short time intervals. This allows the extrapolation over long time-steps of the behavior of a coarse wavelet-based discretization of the system. To validate the approach and the underlying concepts, we perform two 1D1V numerical experiments: nonlinear propagation and steepening of an ion wave, and the expansion of a plasma slab in vacuum. The direct comparisons to resolved PIC simulations show good agreement. We show that the speedup of the projective integration scheme over the full particle scheme scales linearly with the system size, demonstrating efficiency while taking into account fully kinetic, non-Maxwellian effects. This suggests that the approach is potentially interesting for kinetic plasma problems with a large separation of scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interstellar Extinction, Polarization, and Grain Alignment in the Sh 2-185 (IC 59 and IC 63) Region

Optical and infrared continuum polarization from the interstellar medium is driven by radiative processes aligning the grains with the magnetic field. While a quantitative, predictive theory of radiative alignment torques (RATs) exists and has been extensively tested, several parameters of the theory remain to be fully constrained. In a recent paper, Medan & Andersson showed that the polarization efficiency (and therefore grain alignment efficiency) at different locations in the wall of the Local Bubble (LB) could be modeled as proportional to the integrated light intensity from the surrounding stars and OB associations. Here we probe that relationship at high radiation field intensities by studying the extinction and polarization in the two reflection nebulae IC 59 and IC 63 in the Sh 2-185 H ii region, illuminated by the B0 IV star γ Cassiopeia. We combine archival visual polarimetry with new seven-band photometry in the Vilnius system, to derive the polarization efficiency from the material. We find that the same linear relationship seen in the LB wall also applies to the Sh 2-185 region, strengthening the conclusion from the earlier study.

47 OTHER INSTRUMENTATION↗

Determination of the Transition Temperature of a Weak Ferromagnetic Thin Film by Means of an Evolution of the Method Based on the Arrott Plots

The ferromagnetic transition temperature (T{sub c}) of a weak ferromagnetic thin film Cu{sub 0.38}Ni{sub 0.62} has been estimated by using an evolution of the Belov, Goriaga, and Arrott (BGA) model. In fact, the typically used measurements of the zero-field-cooled-field-cooled magnetic moment (m) as a function of temperature (T), the magnetic hysteresis loops (m(H)) as a function of field, and the temperature dependence of the remanent magnetic moment, performed on our sample, did not allow us to individuate precisely the value of T{sub c}, mainly due to the influence of thermal effects and noise on the weak ferromagnetic response of the sample, resulting in a large uncertainty in the estimation of T{sub c}. On the other hand, also the commonly used method of determining the T{sub c} by means of a linear extrapolation of the high-field region of m{sup 2} as a function of the normalized magnetic field (H/m) at different temperatures (Arrott curves), working though for our analyzed material, could not completely fulfill the basic assumptions of the BGA model in general. For these reasons, we have considered the analysis of the derivatives of the Arrott curves, starting from observing that the T{sub c} of the sample corresponds to the temperature where the curvature of the Arrott curves at low fields inverts. In this way, the sensitivity in the determination of the Curie point from the inductive measurements of the magnetic moment is improved while the correct assumptions of the Arrott model are also fully respected.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

Developing VSC-HVDC Oscillation Damping Control Constraints in Unit Commitment

High-voltage direct current (HVDC) systems within the existing power grids will play a pivotal role in enabling high share of renewable power integration and transportation electrification. There is a need for transmission planning models to accurately capture the impact of VSC-HVDC oscillation damping control on the system planning and model the constraints in the security constrained unit commitment (SCUC) problem. In this paper, the SCUC problem in a hybrid DC-AC grid considering damping on electromechanical inter-area oscillation is studied. An electromechanical oscillation-driven transmission capability constraints are derived from system swing equation and incorporated into the SCUC. Multiple system damping control scenarios are set based on the extracted plausible bounds of damping coefficient and inertia constant to linearize the nonlinear expressions. The SCUC problem with proposed constraints is tested on the IEEE 73-bus reliability test system (RTS) case via PLEXOS. The impact of the proposed constraints on regional generation, inertia, power flow, annual system cost, and renewable energy curtailment are fully assessed.

damping control↗