Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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 721 records · Page 40

Global navigation satellite system reflectometry (GNSS-R) algorithms for wetland observations

It is important to closely monitor the state of the world’s wetlands, as climate change and human encroachment in a rapid global urbanization trend threaten to cause large-scale wetland collapse. Because wetlands are often difficult to observe in situ, remote sensing is the only viable way to map wetland extent globally. However, current remote sensing methods suffer limitations in capturing wetland extent, and more importantly, wetland dynamics at appropriate spatial and temporal scales. GNSSReflectometry could help fill the current observation gap, as experimental data show that ground-reflected GNSS signals are very sensitive to changes in inundated areas. Furthermore, because this technique only requires a custom developed receiver and antenna system, a constellation of such instruments can potentially be launched at relatively low cost, providing global observations at sub-daily intervals. One challenge remains, however, which is quantitatively formulating the geophysical product of reflections over the land surface in various states of inundation. Here, we use a novel reflection dataset, derived from the SMAP radar receiver, to elucidate the sensitivity of reflections to small land surface features and their seasonal variations. Additionally, we quantify the dynamic range of reflections over both open and closed wetlands, and suggest an algorithm for wetland type classification.

UNKNOWN↗

Cross-code comparison of the impact of the fill tube on high yield implosions on the National Ignition Facility

Fill tubes are used to inject deuterium and tritium fuel into inertial confinement fusion capsules fielded on the National Ignition Facility. These fill tubes have been shown to have a detrimental effect on capsule performance, primarily by introducing a low-density pathway into the central fuel region that enables the jetting of ablation material into the hot spot. Due to the complexity of the highly nonlinear flow associated with the fill tube and the challenge of diagnosing the evolution of the fill tube jet late in the implosion experiments, the uncertainty in how this perturbation source evolves is great. In this work, we report on the results of a detailed code comparison performed to understand uncertainties in computational modeling of the impact of fill tubes on implosion performance. The study employed two radiation-hydrodynamics codes, HYDRA and xRAGE, which employ very different meshing strategies and hydrodynamics solvers, as well as two radiation transport methodologies, discrete ordinates and multi-group diffusion. Our results demonstrate generally good agreement between codes through most of the implosion although they indicate sensitivity to opacity averaging methods. Late in the implosion, differences arise in the distribution and amount of contaminant although these differences have a remarkably small impact on the amount of yield reduction due to the fill tube. While these results demonstrate sensitivity in fill tube modeling to algorithmic choices, the observed differences between codes are small relative to known sensitivities due to expected variations in the fill tube geometry. Finally, we have developed a methodology for performing multi-group diffusion simulations that show good agreement with the more accurate discrete ordinates method

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗

Low-Memory, Discrete Ordinates, Discontinuous Galerkin Methods for Radiative Transport

The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.

97 MATHEMATICS AND COMPUTING↗

Accurate and Efficient Parametric Model-Order Reduction for Turbulent Thermal Transport

This project produced new algorithms and software for low-cost exploration of turbulent thermal-fluids behavior under parametric variation by using reduced-order models (ROMs). The ROMs numerically solve the governing equations for fluid motion by using a small set of basis functions (typically, N=20-200 modes) to represent the solution. The base modes are computed as optimal combinations of solutions from expensive high-fidelity "anchor-point" solutions involving millions of unknowns, which are typically generated by solving the full Navier-Stokes equations on a supercomputer. The ROM solution is itself a combination of the base modes, where the basis coefficients are determined by evolving an NxN system of nonlinear equations. The overarching idea is to use the inexpensive ROM to predict solutions under conditions where the parameters differ from the anchor-point conditions. Several ingredients are required to make ROMs useful for thermal hydraulics analysis. These include: a stable and accurate ROM that is capable of reproducing the large-scale dynamics of turbulent flow, error indicators than can guide the choice of anchor points, and low-cost mechanisms for evaluating nonlinear terms in the reduced equations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A finite element solution algorithm for the Navier-Stokes equations

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the steady-state kinematics and thermodynamics of a variable viscosity, compressible multiple-species fluid. For an incompressible fluid, the motion may be transient as well. The primitive dependent variables are replaced by a vorticity-streamfunction description valid in domains spanned by rectangular, cylindrical and spherical coordinate systems. Use of derived variables provides a uniformly elliptic partial differential equation description for the Navier-Stokes system, and for which the finite element algorithm is established. Explicit non-linearity is accepted by the theory, since no psuedo-variational principles are employed, and there is no requirement for either computational mesh or solution domain closure regularity. Boundary condition constraints on the normal flux and tangential distribution of all computational variables, as well as velocity, are routinely piecewise enforceable on domain closure segments arbitrarily oriented with respect to a global reference frame.

Baker, A. J.↗

The influence of NO and ClO variations at twilight on the interpretation of solar occultation measurements

Measurement of short-lived photochemically-produced species in the stratosphere by solar occultation is difficult because the rapid variation of such species near the terminator introduces ambiguities in interpreting the measured absorption in terms of meaningful atmospheric abundances. These variations produce tangent path concentrations that are asymmetric relative to the tangent point, as opposed to the symmetrical distribution usually assumed in most inversion algorithms. Neglect of this asymmetry may yield an inverted profile that deviates significantly from the true sunset/sunrise profile. In the present paper, the influence of this effect on solar occultation measurements of ClO and NO is examined. The results show that average inhomogeneity factors, which measure the concentration variation along the tangent path and which can be calculated from a photochemical model, can indicate which species require more careful data analysis.

Boughner, R.↗

Contrast reduction by the atmosphere and retrieval of nonuniform surface reflectance

A radiative transfer model is developed which gives the upward radiance at nadir for any 1-D Lambertian surface reflectance. This model is used to depict the atmospheric effect on the transmittance of contrast for any 1-D surface reflectance. Here by contrast we mean a general variation of the radiation field across the image. With the aid of this model an inversion algorithm is developed for retrieval of true surface reflectance from high resolution satellite data (e.g., Landsat). This inversion technique can be a useful tool for extraction of surface reflectance from satellite data in the case of a surface reflectance variable in one dimension only (e.g., seashore or near borders of big fields). A sensitivity study of the inversion procedure on the knowledge of atmospheric parameters and sensor calibration was performed. It is shown that this inversion technique is stable even in the presence of errors in the sensor calibration and the atmospheric parameters. The method was applied to Landsat data in two wavelengths. The results show reasonable dependence of the derived surface reflectance on the distance from the seashore.

Mekler, Y.↗

Amplitude calibration of spaceborne synthetic aperture radars

Problems encountered during attempts to calibrate SAR imagery, recent successful experiments conducted with SEASAT SAR data, and a proposed program for the calibration and validation of the radar imagery from the forthcoming SIR-B SAR are discussed. The SEASAT SAR data for 10 passes over Death valley, California, were processed with a modified digital correlator. The procedure included a preliminary screening of the data to check for raw data saturation, compensation of waveforms and estimation of the amplitude of the pilot tone. All data was normalized to this pilot tone signal to reduce the effects of variable gains in the data links and ground receivers. The digital correlation algorithm generated image data. Evaluation of 6 passes results in a maximum pass to pass gain variation of only 1.1 dB and a standard deviation amongst the passes of 0.35 dB.

Held, D. N.↗

Amplitude calibration of spaceborne synthetic aperture radars

Problems encountered during attempts to calibrate SAR imagery, recent successful experiments conducted with SEASAT SAR data, and a proposed program for the calibration and validation of the radar imagery from the forthcoming SIR-B SAR are discussed. The SEASAT SAR data for 10 passes over Death Valley, California, were processed with a modified digital correlator. The procedure included a preliminary screening of the data to check for raw data saturation, compensation of waveforms and estimation of the amplitude of the pilot tone. All data was normalized to this pilot tone signal to reduce the effects of variable gains in the data links and ground receivers. The digital correlation algorithm generated image data. Evaluation of 6 passes results in a maximum pass to pass gain variation of only 1.1 dB and a standard deviation amongst the passes of 0.35 dB. previously announced in STAR as N83-26215

Held, D. N.↗

On the implementation of a class of upwind schemes for system of hyperbolic conservation laws

The relative computational effort among the spatially five point numerical flux functions of Harten, van Leer, and Osher and Chakravarthy is explored. These three methods typify the design principles most often used in constructing higher than first order upwind total variation diminishing (TVD) schemes. For the scalar case the difference in operation count between any two algorithms may be very small and yet the operation count for their system counterparts might be vastly different. The situation occurs even though one starts with two different yet equivalent representations for the scalar case.

Yee, H. C.↗

Gain calibrating nonuniform image-array data using only the image data

An algorithm is developed for calibrating the spatial nonuniformity of image-array (CCD-type) detectors. Like other techniques this approach uses multiple, spatially displaced images. In circumstances where high-precision flat fields are not available by other means (i.e., sky flats) this technique is advantageous as it uses the data frames for gain calibration even when the array images extended, nonuniform, sources. Numerical experiments and direct observations with intrinsically uniform and quite nonuniform detectors show that this algorithm is useful when data frames are crowded with sources - circumstance where 'median filtering' flatfielding techniques often fail. The algorithm described is robust and efficiently uses information from multiple data frames to determine pixel gain variations, using visible and IR array observations of extended sources.

Kuhn, J. R.↗

SeGRAm - A practical and versatile tool for spacecraft trajectory optimization

An implementation of the Sequential Gradient/Restoration Algorithm, SeGRAm, is presented along with selected examples. This spacecraft trajectory optimization and simulation program uses variational calculus to solve problems of spacecraft flying under the influence of one or more gravitational bodies. It produces a series of feasible solutions to problems involving a wide range of vehicles, environments and optimization functions, until an optimal solution is found. The examples included highlight the various capabilities of the program and emphasize in particular its versatility over a wide spectrum of applications from ascent to interplanetary trajectories.

Rishikof, Brian H.↗

Prediction of geomagnetic activity

A simple analytic model for the directly driven, adiabatic solar wind-magnetosphere-ionosphere coupling is derived on the basis of frontside reconnection and force balance in the tail. Magnetosphere-ionosphere coupling occurs through Alfven waves which are driven by changes of electric fields in the magnetosphere. The model provides an algorithm which calculates the directly driven part of the auroral electrojets from an observed solar wind variation. For two days in May 1979 the cross correlation between the prediction and observation of AE is greater than 0.9. This suggests that most of the geomagnetic activity during these two days was directly driven by the solar wind.

Goertz, C. K.↗

Interannual Variability of Tropical Ocean Evaporation: A Comparison of Microwave Satellite and Assimilation Results

Remote sensing methodologies for turbulent heat fluxes over oceans depend on driving bulk formulations of fluxes with measured surface winds and estimated near surface thermodynamics from microwave sensors of the Special Sensor Microwave Imager (SSM/I) heritage. We will review recent work with a number of SSM/I-based algorithms and investigate the ability of current data sets to document global, tropical ocean-averaged evaporation changes in association with El Nino and La Nina SST changes. We show that in addition to interannual signals, latent heat flux increases over the period since late 1987 range from approx. .1 to .6 mm/ day are present; these represent trends 2 to 3 times larger than the NCEP Reanalysis. Since atmospheric storage cannot account for the difference, and since compensating evapotranspiration changes over land are highly unlikely to be this large, these evaporation estimates cannot be reconciled with ocean precipitation records such as those produced by the Global Precipitation Climatology Project, GPCP. The reasons for the disagreement include less than adequate intercalibration between SSM/I sensors providing winds and water vapor for driving the algorithms, biases due to the assumption that column integrated water vapor mirrors near surface water vapor variations, and other factors as well. The reanalyses have their own problems with spin-up during assimilation, lack of constraining input data at the ocean surface, and amplitude of synoptic transients.

Robertson, Franklin R.↗

Physical Retrieval of Surface Emissivity Spectrum from Hyperspectral Infrared Radiances

Retrieval of temperature, moisture profiles and surface skin temperature from hyperspectral infrared (IR) radiances requires spectral information about the surface emissivity. Using constant or inaccurate surface emissivities typically results in large retrieval errors, particularly over semi-arid or arid areas where the variation in emissivity spectrum is large both spectrally and spatially. In this study, a physically based algorithm has been developed to retrieve a hyperspectral IR emissivity spectrum simultaneously with the temperature and moisture profiles, as well as the surface skin temperature. To make the solution stable and efficient, the hyperspectral emissivity spectrum is represented by eigenvectors, derived from the laboratory measured hyperspectral emissivity database, in the retrieval process. Experience with AIRS (Atmospheric InfraRed Sounder) radiances shows that a simultaneous retrieval of the emissivity spectrum and the sounding improves the surface skin temperature as well as temperature and moisture profiles, particularly in the near surface layer.

Li, Jun↗

Flight Crew Survey Responses from the Interval Management (IM) Avionics Phase 2 Flight Test

The Interval Management (IM) Avionics Phase 2 flight test used three aircraft over a nineteen day period to operationally evaluate a prototype IM avionics. Quantitative data were collected on aircraft state data and IM spacing algorithm performance, and qualitative data were collected through end-of-scenario and end-of-day flight crew surveys. The majority of the IM operations met the performance goals established for spacing accuracy at the Achieve-by Point and the Planned Termination Point, however there were operations that did not meet goals for a variety of reasons. While the positive spacing accuracy results demonstrate the prototype IM avionics can contribute to the overall air traffic goal, critical issues were also identified that need to be addressed to enhance IM performance. The first category was those issues that impacted the conduct and results of the flight test, but are not part of the IM concept or procedures. These included the design of arrival and approach procedures was not ideal to support speed as the primary control mechanism, the ground-side of the Air Traffic Management Technology Demonstration (ATD-1) integrated concept of operations was not part of the flight test, and the high workload to manually enter the information required to conduct an IM operation. The second category was issues associated with the IM spacing algorithm or flight crew procedures. These issues include the high frequency of IM speed changes and reversals (accelerations), a mismatch between the deceleration rate used by the spacing algorithm and the actual aircraft performance, and some spacing error calculations were sensitive to normal operational variations in aircraft airspeed or altitude which triggered additional IM speed changes. Once the issues in these two categories are addressed, the future IM avionics should have considerable promise supporting the goals of improving system throughput and aircraft efficiency.

Baxley, Brian T.↗

Bhutan Water Resources III: Analyzing Forest Disturbances and Climate Data in Bhutan to Create a Tool for Assisting the Himalayan Environment Rhythm Observation and Evaluation Systems (HEROES) Project

Forest disturbances from bark beetle outbreaks are a major concern in Bhutan, known to cause extensive tree mortality to pine and spruce forests. The NASA DEVELOP team partnered with the Ugyen Wangchuck Institute of Conservation and Environmental Research (UWICER), the Bhutan Foundation, and the Karuna Foundation to assess forest changes for the districts of Bumthang and Haa from 2000 to 2018. The project used preprocessed meteorological data from the Climate Hazards Center Infrared Precipitation with Station (CHIRPS) and Famine Early Warning System Network Land Data Assimilation System (FLDAS), along with Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper plus (ETM+), and Landsat 8 Operational Land Imager (OLI) to assess apparent forest disturbance occurrences and observed climate trends. Shuttle Radar Topography Mission (SRTM) was used to resolve variations in elevation and slope for mountainous regions. Using the Google Earth Engine LandTrendr (LT) code algorithm, along with Landsat data, the team developed an app called Forest Disturbances Detection Toolbox (FDDT) to assess forest changes in Bhutan. The app includes climate variables for the focus districts, along with LT variables, which allows the end users to further examine the cause of disturbances. The team compared geocoordinates for known disturbances with LT disturbance detection products. Although additional work is needed in the future to validate the project end products from the FDDT, the tool will be provided to the project partners to aid forest management efforts in Bhutan.

Tashi Choden↗