Search NASA⌕ Search

SEARCH · Search NASA

Results for “matching problem”

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 253 records · Page 14

A simulation-based model for understanding the time dependent x-ray drive asymmetries and error bars in indirectly driven implosions on the National Ignition Facility

Time-dependent low-mode asymmetries are believed to play a leading role in limiting the performance of current inertial confinement fusion implosions on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)]. These long wavelength modes are initiated and driven by asymmetries in the x-ray flux from the hohlraum; however, the underlying hydrodynamics of the implosion also act to modify and amplify these asymmetries. Here we present here a simulation-based model connecting the time-dependent drive asymmetry seen by the capsule to the measured inflight and hot spot symmetries. This approach is based on a Green's function analysis for which we evaluate the response of the capsule to impulses of drive asymmetry at a series of times. Our model sheds new light on the sensitivity to the drive asymmetry of an imploded capsule, giving a new tool for design. Inverting the problem and finding the drive asymmetry needed to match the experimental data allow us to tightly constrain the drive asymmetry seen by the capsule, providing an error estimate on the result. Doing so, we are able to point out when and how the complex hohlraum simulations start to deviate from what they should obtain to match the experimental data. Ultimately, we project to use this model to make some experimental recommendations to fix the time-dependent low-mode asymmetry of indirectly driven implosions and identify additional measurements to further constrain the asymmetries with a view to improving target design on the NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗

TH-A-9A-06: Inverse Planning of Gamma Knife Radiosurgery Using Natural Physical Models

Purpose: Treatment-planning systems rely on computer intensive optimization algorithms in order to provide radiation dose localization. We are investigating a new optimization paradigm based on natural physical modeling and simulations, which tend to evolve in time and find the minimum energy state. In our research, we aim to match physical models with radiation therapy inverse planning problems, where the minimum energy state coincides with the optimal solution. As a prototype study, we have modeled the inverse planning of Gamma Knife radiosurgery using the dynamic interactions between charged particles and demonstrate the potential of the paradigm. Methods: For inverse planning of Gamma Knife radiosurgery: (1) positive charges are uniformly placed on the surface of tumors and critical structures. (2) The Gamma Knife dose kernels of 4mm, 8mm and 16mm radii are modeled as geometric objects with variable charges. (3) The number of shots per each kernel radii is obtained by solving a constrained integer-linear problem. (4) The shots are placed into the tumor volume and move under electrostatic forces. The simulation is performed until internal forces are zero or maximum iterations are reached. (5) Finally, non-negative least squares (NNLS) is used to calculate the beam-on times for each shot. Results: A 3D C-shaped tumor surrounding a spherical critical structure was used for testing the new optimization paradigm. These tests showed that charges spread out evenly covering the tumor while keeping distance from the critical structure, resulting in a high quality plan. Conclusion: We have developed a new paradigm for dose optimization based on the simulation of physical models. As prototype studies, we applied electrostatic models to Gamma Knife radiosurgery and demonstrated the potential of the new paradigm. Further research and fine-tuning of the model are underway. NSF CBET-0853157.

60 APPLIED LIFE SCIENCES↗

Effect of a wind tunnel on the acoustic field from various aeroacoustic sources

The effects of the walls of an enclosed test section wind tunnel on measurements of sound fields from various sources has been studied. The acoustic field from a known source in a wind tunnel has been modeled as an infinitely long duct with constant cross section. The model was solved with a numerical panel technique in a control volume near the source, and matched to an outer analytic solution. Several sample problems were studied in a rectangular duct with and without flow. The results indicate that the presence of the duct affects the acoustic field, and that small changes in the product of duct cross dimensions and the source wave number can change the acoustic field significantly. It is also shown that, for low-frequency helicopter rotor harmonic noise, measured in typical wind tunnel rotor tests, the sound levels beyond one rotor diameter from the hub are unreliable indications of the free-field sound levels.

Mosher, M.↗

Large-scale motions on the sun - An overview

The history and present status of observations of large-scale velocity fields in the solar atmosphere are reviewed. Observations of the torsional oscillation and of mean meridional circulation suggest a connection of large-scale dynamics with the solar cycle. Significant problems must be solved before Doppler observations can match the precision of tracer measurements, particularly allowing for the effects of changes in line-profile asymmetries and for scattered light. Coordinated observations would establish the reliability of Doppler techniques, but Doppler measurements with precision of order 1 m/s made in a proper spatial-temporal window appear necessary for the identification of sub-global velocity fields varying with time-scales less than that of the solar cycle.

Bogart, Richard S.↗

Monolithic and mechanical multijunction space solar cells

Monolithic and mechanically stacked tandem solar cells have been fabricated with encouraging AM0 efficiencies summarized as: monolithic GaAs/Ge: 19.1 percent (28 C, 4 sq cm); monolithic InP/Ga0.47In0.53As: 22.2 percent (25 C, 0.296 sq cm); monolithic AlGaAs/GaAs/InGaAs: 27.6 percent (80 C, 0.2 sq cm, 100 X); mechanically stacked GaAs/GaSb: 30.8 percent (25 C, 0.049 sq cm, 100 X); and mechanically stacked GaAs/CuInSe2: 23.1 percent (25 C, 4 sq cm). Significant improvement in tandem cell efficiencies nearing to theoretical predictions has been projected with the improvement in cell material quality and processing. Thin-film cells offer improved specific power. It is pointed out that both the monolithic and mechanically stacked cells have their own problems as to size, processing, current-voltage matching, weight, etc. More information is needed on the effect of temperature and radiation on the cell performance. Proper reference cells and full spectrum range simulators are required to measure efficiencies correctly.

Jain, Raj K.↗

Monolithic and mechanical multijunction space solar cells

High-efficiency, lightweight, radiation-resistant solar cells are essential to meet the large power requirements of future space missions. Single-junction cells are limited in efficiency. Higher cell efficiencies could be realized by developing multijunction, multibandgap solar cells. Monolithic and mechanically stacked tandem solar cells surpassing single-junction cell efficiencies have been fabricated. This article surveys the current status of monolithic and mechanically stacked multibandgap space solar cells, and outlines problems yet to be resolved. The monolithic and mechanically stacked cells each have their own problems related to size, processing, current and voltage matching, weight, and other factors. More information is needed on the effect of temperature and radiation on the cell performance. Proper reference cells and full-spectrum range simulators are also needed to measure efficiencies correctly. Cost issues are not addressed, since the two approaches are still in the developmental stage.

Jain, Raj K.↗

A Cutting Surface Algorithm for Semi-Infinite Convex Programming with an Application to Moment Robust Optimization

In this paper, we present and analyze a central cutting surface algorithm for general semi-infinite convex optimization problems and use it to develop a novel algorithm for distributionally robust optimization problems in which the uncertainty set consists of probability distributions with given bounds on their moments. Moments of arbitrary order, as well as nonpolynomial moments, can be included in the formulation. We show that this gives rise to a hierarchy of optimization problems with decreasing levels of risk-aversion, with classic robust optimization at one end of the spectrum and stochastic programming at the other. Although our primary motivation is to solve distributionally robust optimization problems with moment uncertainty, the cutting surface method for general semi-infinite convex programs is also of independent interest. The proposed method is applicable to problems with nondifferentiable semi-infinite constraints indexed by an infinite dimensional index set. Examples comparing the cutting surface algorithm to the central cutting plane algorithm of Kortanek and No demonstrate the potential of our algorithm even in the solution of traditional semi-infinite convex programming problems, whose constraints are differentiable, and are indexed by an index set of low dimension. After the rate of convergence analysis of the cutting surface algorithm, we extend the authors' moment matching scenario generation algorithm to a probabilistic algorithm that finds optimal probability distributions subject to moment constraints. The combination of this distribution optimization method and the central cutting surface algorithm yields a solution to a family of distributionally robust optimization problems that are considerably more general than the ones proposed to date.

97 MATHEMATICS AND COMPUTING↗

Optimal polarizations for achieving maximum contrast in radar images

A systematic approach is presented for obtaining the optimal polarimetric matched filter, which produces maximum contrast between two scattering classes in polarimetric radar images. The maximization procedure involves solving an eigenvalue problem when the eigenvector corresponding to the maximum contrast ratio is optimal polarimetric matched filter. It is shown how to realize the filter in terms of an equivalent transmitting and receiving polarization pair. The case where the radar transmitting polarization state is fixed is addressed. Here, the receiving polarization state is optimized so that the maximum constrast ratio is attained under this constraint.

Swartz, A. A.↗

REBUS-PC v1.5 Verification Test Summary

Motivated in part by its usage for analysis of the University of Missouri Research Reactor (MURR), the verification of the REBUS-PC software for θ-R-Z geometries is discussed and displayed in this manuscript. Of particular importance is to show that the eigenvalue, flux, power, and burnup outputs of the software given a user problem is correct. REBUS-PC contains many features and MURR only makes use of some of those aspects. As an example, the use of region specific cross section sets, interpolation of the cross section data with respect to burnup, and automatic movement of control rods during the fuel cycle is discussed and displayed. To demonstrate the accuracy, 22 verification tests based upon 7 physical geometries are discussed. Geometry descriptions are provided in Section 2 and the REBUS-PC results for each verification test are discussed in Section 3. These verification tests consist of both analytic problems with known solutions and a more complicated input consistent with usage of REBUS-PC on the MURR reactor itself called STRZ. In this verification test summary report, the REBUS-PC software is repeatedly demonstrated to reproduce, or closely match, the eigenvalue, flux, and power distributions from analytic test problem solutions, demonstrating that the implemented finite difference approximation in θ-R-Z is working as intended. Several additional REBUS-PC features are tested, such as control rod movement by detailed comparisons of the results against reference solutions obtained using an auxiliary implemented MatLab program. As expected, the REBUS-PC software produces results that are 100% consistent with the stated capabilities and thus the verification exercise is completed. These verification tests should be executed with each re-compilation of the REBUS-PC source and after any source code changes. All output results from any such new work should be rigorously compared against the stored output of the 22 verification tests cited in this report that forms the basis of the verification process.

97 MATHEMATICS AND COMPUTING↗

An outflow acoustic boundary condition for internal duct flows

A boundary condition for the linear acoustic equation has been developed that allows the acoustic pressure waves to propagate out of the computational domain boundary, just as they would propagate in an infinitely long duct. The problem is divided into two domains: numerical and analytical. The boundary condition provides a matching of the two domains. Examples show this method works well in an acoustic panel program for a model problem (simple source in a rectangular duct with several propagating modes present). The boundary condition is described so that it can be used with various duct various duct geometries and numerical methods.

Mosher, M.↗

Time-dependent mixed convection heat transfer from a sphere in a micro-gravity environment

A fundamental problem of interest for crystal growth in micro-gravity applications involves the mixed convection heat transfer from a sphere in a uniform flow of fluid at a differing temperature. Under the combined influence of the imposed free stream as well as an induced buoyancy force due to thermal expansion of the fluid, the heat transfer from the sphere will be different from that of either the pure forced convection flow or the pure free convection flow. For the present study, the method of matched asymptotic expansions is applied to the laminar flow problem of an impulsively heated, impulsively started sphere in an originally quiescent fluid. Time series expansions are developed for the dependent variables by acknowledging the existence of two distinct regions: one, an inner region, near the sphere, in which viscous effects are significant; and two, an outer region in which the fluid may be treated as inviscid. The time series expansions are developed in terms of the Reynolds number and Richardson number (Buoyancy Parameter), and the relevant heat transfer and drag coefficients are calculated and plotted.

Hommel, Mark J.↗

Comparison of Fault Detection Algorithms for Real-time Diagnosis in Large-Scale System

In this paper, we present a review of different real-time capable algorithms to detect and isolate component failures in large-scale systems in the presence of inaccurate test results. A sequence of imperfect test results (as a row vector of I's and O's) are available to the algorithms. In this case, the problem is to recover the uncorrupted test result vector and match it to one of the rows in the test dictionary, which in turn will isolate the faults. In order to recover the uncorrupted test result vector, one needs the accuracy of each test. That is, its detection and false alarm probabilities are required. In this problem, their true values are not known and, therefore, have to be estimated online. Other major aspects in this problem are the large-scale nature and the real-time capability requirement. Test dictionaries of sizes up to 1000 x 1000 are to be handled. That is, results from 1000 tests measuring the state of 1000 components are available. However, at any time, only 10-20% of the test results are available. Then, the objective becomes the real-time fault diagnosis using incomplete and inaccurate test results with online estimation of test accuracies. It should also be noted that the test accuracies can vary with time --- one needs a mechanism to update them after processing each test result vector. Using Qualtech's TEAMS-RT (system simulation and real-time diagnosis tool), we test the performances of 1) TEAMSAT's built-in diagnosis algorithm, 2) Hamming distance based diagnosis, 3) Maximum Likelihood based diagnosis, and 4) HidderMarkov Model based diagnosis.

Kirubarajan, Thiagalingam↗

The optimal polarizations for achieving maximum contrast in radar images

There is considerable interest in determining the optimal polarizations that maximize contrast between two scattering classes in polarimetric radar images. A systematic approach is presented for obtaining the optimal polarimetric matched filter, i.e., that filter which produces maximum contrast between two scattering classes. The maximization procedure involves solving an eigenvalue problem where the eigenvector corresponding to the maximum contrast ratio is an optimal polarimetric matched filter. To exhibit the physical significance of this filter, it is transformed into its associated transmitting and receiving polarization states, written in terms of horizontal and vertical vector components. For the special case where the transmitting polarization is fixed, the receiving polarization which maximizes the contrast ratio is also obtained. Polarimetric filtering is then applies to synthetic aperture radar images obtained from the Jet Propulsion Laboratory. It is shown, both numerically and through the use of radar imagery, that maximum image contrast can be realized when data is processed with the optimal polarimeter matched filter.

Swartz, A. A.↗

SU-D-BRC-01: An Automatic Beam Model Commissioning Method for Monte Carlo Simulations in Pencil-Beam Scanning Proton Therapy

Purpose: Monte Carlo (MC) simulation is typically regarded as the most accurate dose calculation method for proton therapy. Yet for real clinical cases, the overall accuracy also depends on that of the MC beam model. Commissioning a beam model to faithfully represent a real beam requires finely tuning a set of model parameters, which could be tedious given the large number of pencil beams to commmission. This abstract reports an automatic beam-model commissioning method for pencil-beam scanning proton therapy via an optimization approach. Methods: We modeled a real pencil beam with energy and spatial spread following Gaussian distributions. Mean energy, and energy and spatial spread are model parameters. To commission against a real beam, we first performed MC simulations to calculate dose distributions of a set of ideal (monoenergetic, zero-size) pencil beams. Dose distribution for a real pencil beam is hence linear superposition of doses for those ideal pencil beams with weights in the Gaussian form. We formulated the commissioning task as an optimization problem, such that the calculated central axis depth dose and lateral profiles at several depths match corresponding measurements. An iterative algorithm combining conjugate gradient method and parameter fitting was employed to solve the optimization problem. We validated our method in simulation studies. Results: We calculated dose distributions for three real pencil beams with nominal energies 83, 147 and 199 MeV using realistic beam parameters. These data were regarded as measurements and used for commission. After commissioning, average difference in energy and beam spread between determined values and ground truth were 4.6% and 0.2%. With the commissioned model, we recomputed dose. Mean dose differences from measurements were 0.64%, 0.20% and 0.25%. Conclusion: The developed automatic MC beam-model commissioning method for pencil-beam scanning proton therapy can determine beam model parameters with satisfactory accuracy.

60 APPLIED LIFE SCIENCES↗

Adjoint-Based Inversion of Geodetic Data for Sources of Deformation and Strain

An adjoint-based formulation leads to a particularly efficient approach for inverting geodetic measurements for the source of the deformation. Specifically, the quantities necessary to iteratively improve the fit to the observations can be computed with just three forward calculations, one to obtain the current residuals, another to solve the adjoint problem, and a third to compute the step length. An inversion algorithm utilizing the adjoint-based gradient is applied to a set of Interferometric Synthetic Aperture Radar (InSAR) data gathered between 2016 and 2018 over the Tulare Basin in California's Central Valley. Because the measured deformation is due to groundwater withdrawal, a penalty function is included in the inversion to avoid placing aquifer volume change in locations that are far from any documented wells. The solution of the inverse problem provides estimates of aquifer compaction that provide a match to the observed range changes while honoring the well data. The solution indicates an average aquifer volume loss of 2.17 km 3 /year over the two year period from January 2016 to January 2018, encompassing one drought year (2016) and one wet year (2017). Finally, this magnitude of lost volume is compatible with the 3.1 km 3 /year decrease in water volume for the entire Central Valley, estimated from GRACE satellite gravity data.

58 GEOSCIENCES↗

Real-Time Multi-Vehicle Multi-Camera Tracking with Graph-Based Tracklet Features

An essential application in intelligent transportation systems is multi-target multi-camera tracking (MTMCT), where the target’s activity is tracked from different cameras. Although the tracking-by-detection scheme is the primary paradigm in MTMCT, the object association information from the video frames is lost. This is mainly because the multi-camera multi-object matching uses the information from the video frames separately. To solve this problem and leverage this association information, we propose an MTMCT framework, where features are built in the form of a graph and a graph similarity algorithm is used to match multi-camera objects. In this paper, we focus on the real-time scenario, where only the past images are used to match an object. Our method achieves an IDF1 score (the ratio of the number of correctly identified objects to the number of ground truth and average objects) of 0.75 with a rate of 14 frames per second (fps).

Engineering↗

Surface matching for correlation of virtual models: Theory and application

Virtual reality can enable a robot user to off line generate and test in a virtual environment a sequence of operations to be executed by the robot in an assembly cell. Virtual models of objects are to be correlated to the real entities they represent by means of a suitable transformation. A solution to the correlation problem, which is basically a problem of 3-dimensional adjusting, has been found exploiting the surface matching theory. An iterative algorithm has been developed, which matches the geometric surface representing the shape of the virtual model of an object, with a set of points measured on the surface in the real world. A peculiar feature of the algorithm is to work also if there is no one-to-one correspondence between the measured points and those representing the surface model. Furthermore the problem of avoiding convergence to local minima is solved, by defining a starting point of states ensuring convergence to the global minimum. The developed algorithm has been tested by simulation. Finally, this paper proposes a specific application, i.e., correlating a robot cell, equipped for biomedical use with its virtual representation.

Caracciolo, Roberto↗