Search NASA⌕ Search

SEARCH · Search NASA

Results for “projected gradient”

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 109 records · Page 6

A rotorcraft in-flight ice detection framework using computational aeroacoustics and Bayesian neural networks

Abstract This work develops a novel ice detection framework specifically suitable for rotorcraft using computational aeroacoustics and Bayesian neural networks. In an offline phase of the work, the acoustic signature of glaze and rime ice shapes on an oscillating wing are computed. In addition, the aerodynamic performance indicators corresponding to the ice shapes are also monitored. These performance indicators include the lift, drag, and moment coefficients. A Bayesian neural network is subsequently trained using projected Stein variational gradient descent to create a mapping from the acoustic signature generated by the iced wings to predict their performance indicators along with quantified uncertainty that is highly important for time- and safety-critical decision-making scenarios. While the training is carried out fully offline, usage of the Bayesian neural network to make predictions can be conducted rapidly online allowing for an ice detection system that can be used in real time and in-flight.

42 ENGINEERING↗

Bayesian inference of heterogeneous epidemic models: Application to COVID-19 spread accounting for long-term care facilities

Here we propose a high dimensional Bayesian inference framework for learning heterogeneous dynamics of a COVID-19 model, with a specific application to the dynamics and severity of COVID-19 inside and outside long-term care (LTC) facilities. We develop a heterogeneous compartmental model that accounts for the heterogeneity of the time-varying spread and severity of COVID-19 inside and outside LTC facilities, which is characterized by time-dependent stochastic processes and time-independent parameters in ~ 1500 dimensions after discretization. To infer these parameters, we use reported data on the number of confirmed, hospitalized, and deceased cases with suitable post-processing in both a deterministic inversion approach with appropriate regularization as a first step, followed by Bayesian inversion with proper prior distributions. To address the curse of dimensionality and the ill-posedness of the high-dimensional inference problem, we propose use of a dimension-independent projected Stein variational gradient descent method, and demonstrate the intrinsic low-dimensionality of the inverse problem. We present inference results with quantified uncertainties for both New Jersey and Texas, which experienced different epidemic phases and patterns. Moreover, we also present forecasting and validation results based on the empirical posterior samples of our inference for the future trajectory of COVID-19.

42 ENGINEERING↗

A Kaczmarz-inspired approach to accelerate the optimization of neural network wavefunctions

Neural network wavefunctions optimized using the variational Monte Carlo method have been shown to produce highly accurate results for the electronic structure of atoms and small molecules, but the high cost of optimizing such wavefunctions prevents their application to larger systems. We propose the Subsampled Projected-Increment Natural Gradient Descent (SPRING) optimizer to reduce this bottleneck. SPRING combines ideas from the recently introduced minimum-step stochastic reconfiguration optimizer (MinSR) and the classical randomized Kaczmarz method for solving linear least-squares problems. We demonstrate that SPRING outperforms both MinSR and the popular Kronecker-Factored Approximate Curvature method (KFAC) across a number of small atoms and molecules, given that the learning rates of all methods are optimally tuned. For example, on the oxygen atom, SPRING attains chemical accuracy after forty thousand training iterations, whereas both MinSR and KFAC fail to do so even after one hundred thousand iterations.

97 MATHEMATICS AND COMPUTING↗

Development of an Ultrahigh-bandwidth Phase Contrast Imaging System for detection of electron scale turbulence and Gigahertz Radio-Frequency Waves

The study of waves and turbulence is vital to the development of future reactor-grade plasma devices developed in the quest for fusion energy. These fluctuations are responsible for moving heat and particles across the magnetic field, and a predictive understanding of them is needed to achieve the density and temperature required to sustain a plasma fusion reaction. While many techniques have been developed for measuring waves and fluctuations, every measurement method has limitations. There are relatively few techniques for measuring very high frequency fluctuations, such as radio frequency waves injected to heat the plasma, unstable waves driven by suprathermal particles, or short wavelength electrostatic waves driven by electron temperature or density gradients. The present project builds upon the proven phase contrast imaging (PCI) technique to extend the response of the diagnostic by orders of magnitude in frequency and almost a factor of ten in spatial resolution. PCI provides a measurement of electron density based on small angle scattering of a CO-2 laser beam by using optical techniques to render a phase shift as an intensity change on a detector. Due to the telecommunications revolution, technological development by manufacturers has focused on components in the near -infrared, so that high- quality lasers and detectors at 1.55 µm are readily available. Shifting PCI design to a new, shorter wavelength has numerous advantages and challenges. Similar detector performance is available with room-temperature arrays with GHz bandwidth, while the detector arrays for 10.6 µm required liquid nitrogen cooling and were therefore limited to a bandwidth of about 1 MHz. Shifting to a shorter wavelength reduces the angle at which the laser beam scatters off of plasma waves, allowing more such scattered components to pass through the aperture of the vacuum vessel port and be collected by the PCI, which increases the spatial resolution. Concomitant with these benefits, various questions of performance arise. At shorter laser wave- length, the sensitivity to mirror and lens quality is increased, the contribution of the laser to the overall system noise is increased, and the sensitivity to vibrations is increased. The custom optical components at the heart of the PCI technique were required to be properly scaled for the shorter wavelength, so fabrication technologies needed to be explored. This project was designed to show that a low noise, high response PCI system at 1.55 µm could be constructed and operated, and then to quantify potential issues to allow extrapolation to a full-size production system providing physics measurements on a large plasma device. The first stage, producing the custom optical component called a Phase Plate, was successfully achieved using two methods. First, an easily reproduced masking and coating technique was able to produce good phase plates with the required parameters. Second, a nanofabrication technique was found to produce extremely high quality phase plates at a competitive price. The PCI constructed with the new phase plates and 1.55 µm laser was found to provide excellent wavelength measurements with the theoretically expected response. The sensitivity to optical surface quality was found to be in line with previous measurements at 10.6 µm. The observed signal-to-noise ratio was similar to the theoretically expected value. The effect of vibrations on PCI was studied with the first measurement of the effect of beam motion on PCI response and comparison to theory, allowing for a quantitative prediction of the effect of vibrations on a production PCI system and the requirements for improved beam stabilization. PCI is an extremely cost-effective method to provide a low noise, absolutely calibrated measurement of plasma fluctuations across a wide spatial scale. This project has shown that a 1.55 µm PCI using modern techniques and components is less expensive than the 10.6 µm PCI of a few years ago, with the largest savings in phase plate fabrication and the infrared detector array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pattern Search Algorithms for Bound Constrained Minimization

We present a convergence theory for pattern search methods for solving bound constrained nonlinear programs. The analysis relies on the abstract structure of pattern search methods and an understanding of how the pattern interacts with the bound constraints. This analysis makes it possible to develop pattern search methods for bound constrained problems while only slightly restricting the flexibility present in pattern search methods for unconstrained problems. We prove global convergence despite the fact that pattern search methods do not have explicit information concerning the gradient and its projection onto the feasible region and consequently are unable to enforce explicitly a notion of sufficient feasible decrease.

Lewis, Robert Michael↗

Excellent performance of 650 MHz single-cell niobium cavity after electropolishing

Electropolishing process and cathodes have undergone modification and optimization for both low- and high-beta 650 MHz five-cell niobium cavities for PIP-II. Cavities treated with these modified electropolishing conditions exhibited smooth surfaces and good performance in baseline tests. Nonetheless, due to administrative constraints on project cavities, maximum gradient performance testing was not conducted. This paper presents a study conducted on a single-cell 650 MHz cavity utilizing the optimized electropolishing conditions, highlighting the maximum performance attained for this specific cavity. The cavity tested at 2 K in a vertical cryostat reached a superior accelerating field gradient of 53.3 MV/m at Q0 of 1.6x1010, which is the highest gradient attained for this type of large-sized cavities.

43 PARTICLE ACCELERATORS↗

A multivariate variational objective analysis-assimilation method. Part 1: Development of the basic model

The variational method of undetermined multipliers is used to derive a multivariate model for objective analysis. The model is intended for the assimilation of 3-D fields of rawinsonde height, temperature and wind, and mean level temperature observed by satellite into a dynamically consistent data set. Relative measurement errors are taken into account. The dynamic equations are the two nonlinear horizontal momentum equations, the hydrostatic equation, and an integrated continuity equation. The model Euler-Lagrange equations are eleven linear and/or nonlinear partial differential and/or algebraic equations. A cyclical solution sequence is described. Other model features include a nonlinear terrain-following vertical coordinate that eliminates truncation error in the pressure gradient terms of the horizontal momentum equations and easily accommodates satellite observed mean layer temperatures in the middle and upper troposphere. A projection of the pressure gradient onto equivalent pressure surfaces removes most of the adverse impacts of the lower coordinate surface on the variational adjustment.

Achtemeier, Gary L.↗

Unprecedented Drought Challenges for Texas Water Resources in a Changing Climate: What Do Researchers and Stakeholders Need to Know?

Long‐range water planning is complicated by factors that are rapidly changing in the 21st century, including climate, population, and water use. Here, we analyze climate factors and drought projections for Texas as an example of a diverse society straddling an aridity gradient to examine how the projections can best serve water stakeholder needs. We find that climate models are robust in projecting drying of summer‐season soil moisture and decreasing reservoir supplies for both the eastern and western portions of Texas during the 21st century. Further, projections indicate drier conditions during the latter half of the 21st century than even the most arid centuries of the last 1,000 years that included megadroughts. To illustrate how accounting for drought non‐stationarity may increase water resiliency, we consider generalized case studies involving four key stakeholder groups: agricultural producers, large surface water suppliers, small groundwater management districts, and regional water planning districts. We also examine an example of customized climate information being used as input to long‐range water planning. We find that while stakeholders value the quantitative capability of climate model outputs, more specific climate‐related information better supports resilience planning across multiple stakeholder groups. New suites of tools could provide necessary capacity for both short and long‐term, stakeholder‐specific adaptive planning.

drought↗

Experimental and Computational Studies of Crystal Nucleation in Composition Gradients (Final Report)

The major goals of this project were to develop and validate a predictive nucleation model that incorporates composition gradients and is applicable to a large range of metallic systems which form both stochiometric and non-stochiometric compounds. Computationally we planned to expand the extent of thermodynamics-based theories and lay the groundwork to improve the general understanding and predictive capabilities of the role that gradients play in phase formation, glass formability and stability, and nucleation and growth events. To address these goals, we used Molecular Dynamics (MD) simulations and both isothermal and isochronal nanocalorimetric experiments on amorphous phases with a controlled composition gradient with the hopes of validating and improving the classical nucleation model and its use in solid solutions. While the computations were successful and identified ways to improve the classical nucleation model, the in situ nanocalorimetry studies proved very challenging due to unexpected difficulties in fabricating effective calorimeters and samples. Thus, we could not experimentally validate our predicted influence of composition gradients on nucleation. Nonetheless, important insights were gained and modifications to the classical nucleation theory were suggested.

36 MATERIALS SCIENCE↗

Intercomparison of methods of coupling between convection and large-scale circulation: 2. Comparison over nonuniform surface conditions

As part of an international intercomparison project, the weak temperature gradient (WTG) and damped gravity wave (DGW) methods are used to parameterize large-scale dynamics in a set of cloud-resolving models (CRMs) and single column models (SCMs). The WTG or DGW method is implemented using a configuration that couples a model to a reference state defined with profiles obtained from the same model in radiative-convective equilibrium. We investigated the sensitivity of each model to changes in SST, given a fixed reference state. We performed a systematic comparison of the WTG and DGW methods in different models, and a systematic comparison of the behavior of those models using the WTG method and the DGW method. The sensitivity to the SST depends on both the large-scale parameterization method and the choice of the cloud model. In general, SCMs display a wider range of behaviors than CRMs. All CRMs using either the WTG or DGW method show an increase of precipitation with SST, while SCMs show sensitivities which are not always monotonic. CRMs using either the WTG or DGW method show a similar relationship between mean precipitation rate and column-relative humidity, while SCMs exhibit a much wider range of behaviors. DGW simulations produce large-scale velocity profiles which are smoother and less top-heavy compared to those produced by the WTG simulations. Lastly, these large-scale parameterization methods provide a useful tool to identify the impact of parameterization differences on model behavior in the presence of two-way feedback between convection and the large-scale circulation.

54 ENVIRONMENTAL SCIENCES↗

SUPG stabilization for the nonconforming virtual element method for advection–diffusion–reaction equations

We present here the design, convergence analysis and numerical investigations of the nonconforming virtual element method with Streamline Upwind/Petrov–Galerkin (VEM-SUPG) stabilization for the numerical resolution of convection–diffusion–reaction problems in the convective-dominated regime. According to the virtual discretization approach, the bilinear form is split as the sum of a consistency and a stability term. The consistency term is given by substituting the functions of the virtual space and their gradients with their polynomial projection in each term of the bilinear form (including the SUPG stabilization term). Polynomial projections can be computed exactly from the degrees of freedom. The stability term is also built from the degrees of freedom by ensuring the correct scalability properties with respect to the mesh size and the equation coefficients. The nonconforming formulation relaxes the continuity conditions at cell interfaces and a weaker regularity condition is considered involving polynomial moments of the solution jumps at cell interface. Optimal convergence properties of the method are proved in a suitable norm, which includes contribution from the advective stabilization terms. Experimental results confirm the theoretical convergence rates.

97 MATHEMATICS AND COMPUTING↗

Inferring morphology and strength of magnetic fields from proton radiographs

Proton radiography is an important diagnostic method for laser plasma experiments and is particularly important in the analysis of magnetized plasmas. The theory of radiographic image analysis has heretofore only permitted somewhat limited analysis of the radiographs of such plasmas. We furnish here a theory that remedies this deficiency. We show that to linear order in magnetic field gradients, proton radiographs are projection images of the MHD current along the proton trajectories. We demonstrate that in the linear regime (i.e., the small image contrast regime), the full structure of the projected perpendicular magnetic field can be reconstructed by solving a steady-state inhomogeneous 2-dimensional diffusion equation sourced by the radiograph fluence contrast data. Here, we explore the validity of the scheme with increasing image contrast, as well as limitations of the inversion method due to the Poisson noise, discretization errors, radiograph edge effects, and obstruction by laser target structures. We also provide a separate analysis that is suited to the inference of isotropic-homogeneous magnetic turbulence spectra. Finally, we discuss extension of these results to the nonlinear regime (i.e., the order unity image contrast regime).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Western Sector In-Situ Chemical Oxidation Project: Additional Injection Pilot Testing Supplemental Results (U)

Trichloroethylene (TCE) and tetrachloroethylene (PCE) are persistent at elevated concentrations in the Western Sector of the M-Area Hazardous Waste Management Facility (HWMF) groundwater plume. Migration of these chlorinated ethenes as Dense Non-Aqueous Phase Liquids (DNAPLs) occurred from the M-Area Settling Basin along geologic strata resulting in a high concentration groundwater plume in Western Sector. Pump-and-treat has been active for over 30 years at the M-Area HWMF to hydraulically contain the primary sources of the plume in addition to other remediation strategies throughout the contaminated zone, like soil vapor extraction, dynamic underground stripping, groundwater recirculation wells, and bioremediation. In-situ chemical oxidation (ISCO) has been identified as a viable treatment method for removal of these volatile organic compounds (VOCs) due to the complex lithology and high concentrations in this part of the aquifer. In addition to ISCO, a new recovery well (i.e., RWM018) was also installed in Western Sector to contain the high concentration VOC plume. The ISCO projects have utilized the gradient induced by RWM018 to help plan the injection and monitoring locations. The combination of ISCO and the additional pump and treat in Western Sector with the installation of RWM018 will help to control the plume and maximize remediation efforts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

ON THE NATURE OF X-RAY SURFACE BRIGHTNESS FLUCTUATIONS IN M87

X-ray images of galaxy clusters and gas-rich elliptical galaxies show a wealth of small-scale features that reflect fluctuations in density and/or temperature of the intracluster medium. In this paper we study these fluctuations in M87/Virgo to establish whether sound waves/shocks, bubbles, or uplifted cold gas dominate the structure. We exploit the strong dependence of the emissivity on density and temperature in different energy bands to distinguish between these processes. Using simulations we demonstrate that our analysis recovers the leading type of fluctuation even in the presence of projection effects and temperature gradients. We confirm the isobaric nature of cool filaments of gas entrained by buoyantly rising bubbles, extending to 7′ to the east and southwest, and the adiabatic nature of the weak shocks at 40″ and 3′ from the center. For features of ∼5–10 kpc, we show that the central 4′ × 4′ region is dominated by cool structures in pressure equilibrium with the ambient hotter gas while up to 30% of the variance in this region can be ascribed to adiabatic fluctuations. The remaining part of the central 14′ × 14′ region, excluding the arms and shocks described above, is dominated by apparently isothermal fluctuations (bubbles) with a possible admixture (at the level of ∼30%) of adiabatic (sound waves) and by isobaric structures. Larger features, of about 30 kpc, show a stronger contribution from isobaric fluctuations. The results broadly agree with a model based on feedback from an active galactic nucleus mediated by bubbles of relativistic plasma.

79 ASTRONOMY AND ASTROPHYSICS↗

A Numerical Method for High Energy X-Ray Source Shape Reconstruction

Quantifying the source shape of a pulsed power X-ray imaging system is essential in determining a lower bound for the overall resolution of the system. A method of reconstructing the source shape involves solving an inverse problem with an operator matrix of size n x n, where n is the number of pixels in the image. Building and inverting such a large matrix quickly becomes intractable as n increases due to both memory and computational constraints, even in well-posed problems. The ill-posed nature of this problem adds further complications, making direct inversion unstable. Taking advantage of the structure of the operator matrix and utilizing the Projected Newton and Conjugate Gradient algorithms, we are able to perform reconstructions of images that are orders of magnitude larger than previously possible.

97 MATHEMATICS AND COMPUTING↗

MO-FG-CAMPUS-JeP3-03: Detection of Unpredictable Patient Movement During SBRT Using a Single KV Projection of An On-Board CBCT System: Simulation Study

Purpose: An unpredictable movement of a patient can occur during SBRT even when immobilization devices are applied. In the SBRT treatments using a conventional linear accelerator detection of such movements relies heavily on human interaction and monitoring. This study aims to detect such positional abnormalities in real-time by assessing intra-fractional gantry mounted kV projection images of a patient’s spine. Methods: We propose a self-CBCT image based spine tracking method consisting of the following steps: (1)Acquire a pre-treatment CBCT image; (2)Transform the CBCT volume according to the couch correction; (3)Acquire kV projections during treatment beam delivery; (4)Simultaneously with each acquisition generate a DRR from the CBCT volume based-on the current projection geometry; (5)Perform an intensity gradient-based 2D registration between spine ROI images of the projection and the DRR images; (6)Report an alarm if the detected 2D displacement is beyond a threshold value. To demonstrate the feasibility, retrospective simulations were performed on 1,896 projections from nine CBCT sessions of three patients who received lung SBRT. The unpredictable movements were simulated by applying random rotations and translations to the reference CBCT prior to each DRR generation. As the ground truth, the 3D translations and/or rotations causing >3 mm displacement of the midpoint of the thoracic spine were regarded as abnormal. In the measurements, different threshold values of 2D displacement were tested to investigate sensitivity and specificity of the proposed method. Results: A linear relationship between the ground truth 3D displacement and the detected 2D displacement was observed (R{sup 2} = 0.44). When the 2D displacement threshold was set to 3.6 mm the overall sensitivity and specificity were 77.7±5.7% and 77.9±3.5% respectively. Conclusion: In this simulation study, it was demonstrated that intrafractional kV projections from an on-board CBCT system have a potential to detect unpredictable patient movement during SBRT. This research is funded by Interfractional Imaging Research Grant from Elekta.

60 APPLIED LIFE SCIENCES↗

Salton Sea solar pond project

The feasibility of constructing salt gradient solar ponds within the Salton Sea is being studied. These ponds would serve a dual purpose: (1) become a depository for unwanted salt and (2) supply thermal energy for driving turbine electric power systems. Under present circumstances, the rise in salinity is expected to eliminate fish life and create other unfavorable conditions. The proposed concept would have a power generation potential of 600 MWe.

French, R. L.↗

ECP Milestone Report WBS 2.3.4.13 ECP/VTK-m FY18Q1 [MS-18/01-03] Multiblock / Gradients / Release STDA05-5.

The FY18Q1 milestone of the ECP/VTK-m project includes the implementation of a multiblock data set, the completion of a gradients filtering operation, and the release of version 1.1 of the VTK-m software. With the completion of this milestone, the new multiblock data set allows us to iteratively schedule algorithms on composite data structures such as assemblies or hierarchies like AMR. The new gradient algorithms approximate derivatives of fields in 3D structures with finite differences. Finally, the release of VTK-m version 1.1 tags a stable release of the software that can more easily be incorporated into external projects.

97 MATHEMATICS AND COMPUTING↗