Search NASA⌕ Search

SEARCH · Search NASA

Results for “standard 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 199 records · Page 11

Analytical and Computational Aspects of Collaborative Optimization

Bilevel problem formulations have received considerable attention as an approach to multidisciplinary optimization in engineering. We examine the analytical and computational properties of one such approach, collaborative optimization. The resulting system-level optimization problems suffer from inherent computational difficulties due to the bilevel nature of the method. Most notably, it is impossible to characterize and hence identify solutions of the system-level problems because the standard first-order conditions for solutions of constrained optimization problems do not hold. The analytical features of the system-level problem make it difficult to apply conventional nonlinear programming algorithms. Simple examples illustrate the analysis and the algorithmic consequences for optimization methods. We conclude with additional observations on the practical implications of the analytical and computational properties of collaborative optimization.

Alexandrov, Natalia M.↗

Nyx: A MASSIVELY PARALLEL AMR CODE FOR COMPUTATIONAL COSMOLOGY

We present a new N-body and gas dynamics code, called Nyx, for large-scale cosmological simulations. Nyx follows the temporal evolution of a system of discrete dark matter particles gravitationally coupled to an inviscid ideal fluid in an expanding universe. The gas is advanced in an Eulerian framework with block-structured adaptive mesh refinement; a particle-mesh scheme using the same grid hierarchy is used to solve for self-gravity and advance the particles. Computational results demonstrating the validation of Nyx on standard cosmological test problems, and the scaling behavior of Nyx to 50,000 cores, are presented.

79 ASTRONOMY AND ASTROPHYSICS↗

Standard Model–axion–seesaw–Higgs portal inflation. Five problems of particle physics and cosmology solved in one stroke

We present a minimal extension of the Standard Model (SM) providing a consistent picture of particle physics from the electroweak scale to the Planck scale and of cosmology from inflation until today. Three right-handed neutrinos N {sub i} , a new color triplet Q and a complex SM-singlet scalar σ, whose vacuum expectation value v {sub σ} ∼ 10{sup 11} GeV breaks lepton number and a Peccei-Quinn symmetry simultaneously, are added to the SM. At low energies, the model reduces to the SM, augmented by seesaw generated neutrino masses and mixing, plus the axion. The latter solves the strong CP problem and accounts for the cold dark matter in the Universe. The inflaton is comprised by a mixture of σ and the SM Higgs, and reheating of the Universe after inflation proceeds via the Higgs portal. Baryogenesis occurs via thermal leptogenesis. Thus, five fundamental problems of particle physics and cosmology are solved at one stroke in this unified Standard Model—axion—seesaw—Higgs portal inflation (SMASH) model. It can be probed decisively by upcoming cosmic microwave background and axion dark matter experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluating the effects of real power losses in optimal power flow based storage integration

This study proposes a DC optimal power flow (DCOPF) with losses formulation (the `-DCOPF+S problem) and uses it to investigate the role of real power losses in OPF based grid-scale storage integration. We derive the `- DCOPF+S problem by augmenting a standard DCOPF with storage (DCOPF+S) problem to include quadratic real power loss approximations. This procedure leads to a multi-period nonconvex quadratically constrained quadratic program, which we prove can be solved to optimality using either a semidefinite or second order cone relaxation. Our approach has some important benefits over existing models. It is more computationally tractable than ACOPF with storage (ACOPF+S) formulations and the provably exact convex relaxations guarantee that an optimal solution can be attained for a feasible problem. Adding loss approximations to a DCOPF+S model leads to a more accurate representation of locational marginal prices, which have been shown to be critical to determining optimal storage dispatch and siting in prior ACOPF+S based studies. Case studies demonstrate the improved accuracy of the `-DCOPF+S model over a DCOPF+S model and the computational advantages over an ACOPF+S formulation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Quantum Machine Learning with SQUID

In this work we present the Scaled QUantum IDentifier (SQUID), an open-source framework for exploring hybrid Quantum-Classical algorithms for classification problems. The classical infrastructure is based on PyTorch and we provide a standardized design to implement a variety of quantum models with the capability of back-propagation for efficient training. We present the structure of our framework and provide examples of using SQUID in a standard binary classification problem from the popular MNIST dataset. In particular, we highlight the implications for scalability for gradient-based optimization of quantum models on the choice of output for variational quantum models.

97 MATHEMATICS AND COMPUTING↗

CCSDS Spacecraft Monitor and Control Mission Operations Interoperability Prototype

We are entering a new era in space exploration. Reduced operating budgets require innovative solutions to leverage existing systems to implement the capabilities of future missions. Custom solutions to fulfill mission objectives are no longer viable. Can NASA adopt international standards to reduce costs and increase interoperability with other space agencies? Can legacy systems be leveraged in a service oriented architecture (SOA) to further reduce operations costs? The Operations Technology Facility (OTF) at the Johnson Space Center (JSC) is collaborating with Deutsches Zentrum fur Luft- und Raumfahrt (DLR) to answer these very questions. The Mission Operations and Information Management Services Area (MOIMS) Spacecraft Monitor and Control (SM&C) Working Group within the Consultative Committee for Space Data Systems (CCSDS) is developing the Mission Operations standards to address this problem space. The set of proposed standards presents a service oriented architecture to increase the level of interoperability among space agencies. The OTF and DLR are developing independent implementations of the standards as part of an interoperability prototype. This prototype will address three key components: validation of the SM&C Mission Operations protocol, exploration of the Object Management Group (OMG) Data Distribution Service (DDS), and the incorporation of legacy systems in a SOA. The OTF will implement the service providers described in the SM&C Mission Operation standards to create a portal for interaction with a spacecraft simulator. DLR will implement the service consumers to perform the monitor and control of the spacecraft. The specifications insulate the applications from the underlying transport layer. We will gain experience with a DDS transport layer as we delegate responsibility to the middleware and explore transport bridges to connect disparate middleware products. A SOA facilitates the reuse of software components. The prototype will leverage the capabilities of existing legacy systems. Various custom applications and middleware solutions will be combined into one system providing the illusion of a set of homogenous services. This paper will document our journey as we implement the interoperability prototype. The team consists of software engineers with experience on the current command, telemetry and messaging systems that support the International Space Station (ISS) and Space Shuttle programs. Emphasis will be on the objectives, results and potential cost saving benefits.

Lucord, Steve↗

An improved simulated annealing algorithm for standard cell placement

Simulated annealing is a general purpose Monte Carlo optimization technique that was applied to the problem of placing standard logic cells in a VLSI ship so that the total interconnection wire length is minimized. An improved standard cell placement algorithm that takes advantage of the performance enhancements that appear to come from parallelizing the uniprocessor simulated annealing algorithm is presented. An outline of this algorithm is given.

Jones, Mark↗

A method for nonlinear optimization with discrete design variables

A numerical method is presented for the solution of nonlinear discrete optimization problems. The applicability of discrete optimization to engineering design is discussed, and several standard structural optimization problems are solved using discrete design variables. The method uses approximation techniques to create subproblems suitable for linear mixed-integer programming methods. The method employs existing software for continuous optimization and integer programming.

Olsen, Gregory R.↗

SOLAR MODELS WITH ACCRETION. I. APPLICATION TO THE SOLAR ABUNDANCE PROBLEM

We generate new standard solar models using newly analyzed nuclear fusion cross sections and present results for helioseismic quantities and solar neutrino fluxes. The status of the solar abundance problem is discussed. We investigate whether nonstandard solar models with accretion from the protoplanetary disk might alleviate this problem. We examine a broad range of models, analyzing metal-enriched and metal-depleted accretion and three scenarios for the timing of accretion. Only partial solutions are found. For metal-rich accreted material (Z{sub ac} {approx}> 0.018) there exist combinations of accreted mass and metallicity that bring the depth of the convective zone into agreement with the helioseismic value. For the surface helium abundance, the helioseismic value is reproduced if metal-poor or metal-free accretion is assumed (Z{sub ac} {approx}< 0.09). In both cases a few percent of the solar mass must be accreted. Precise values depend on when accretion takes place. We do not find a simultaneous solution to both problems but speculate that changing the hydrogen-to-helium mass ratio in the accreted material may lead to more satisfactory solutions. We also show that, with current data, solar neutrinos are already a very competitive source of information about the solar core and can help constraining possible accretion histories. Even without helioseismic constraints, solar neutrinos rule out the possibility that more than 0.02 M{sub Sun} from the protoplanetary disk were accreted after the Sun settled on the main sequence. Finally, we discuss how measurements of neutrinos from the CN cycle could shed light on the interaction between the early Sun and its protoplanetary disk.

79 ASTRONOMY AND ASTROPHYSICS↗

An archetype hydrogen atmosphere problem

Populations for the first three bound states and the continuum of hydrogen are determined for an isothermal hydrostatic atmosphere at 20,000 K. The atmosphere is treated as optically thin in the Balmer and Paschen continua and illuminated by continuum radiation at these wavelengths with prescribed radiation temperatures. The atmosphere is optically thick in the 2-1, 3-1, 3-2 and c-1 transitions. Three stages of approximation are treated: (1) radiative detailed balance in the 2-1, 3-1 and 3-2 transitions, (2) radiative detailed balance in the 3-1 and 3-2 transitions, and (3) all transitions out of detailed balance. The solution of this problem is nontrivial and presents sufficient difficulty to have caused the failure of at least one rather standard technique. The problem is thus a good archetype against which new methods or new implementations of old methods may be tested.

Athay, R. G.↗

HamPerf: A Hamiltonian-Oriented Approach to Quantum Benchmarking

Quantum computing technologies are undergoing rapid development. The different qubit modalities being considered for quantum computing each have their strengths and weaknesses, making it challenging to compare their performance relative to each other and the state-of-the-art in classical high-performance computing. To better understand the utility of a given quantum processor and to assess when and how it will be able to advance the frontiers of computational science, researchers need a robust approach to quantum benchmarking. A variety of approaches have been proposed, many of which characterize the presence of noise in current quantum devices. These efforts include component-level performance metrics, such as randomized benchmarking and gate set tomography; high-level application-dependent metrics; and devicelevel metrics, such as the Quantum Volume. However, it remains unclear how low-level metrics, such as fidelities and decoherence times, and global device metrics, such as Quantum Volume, relate to the computational utility and practical limitations of quantum processors to solve useful problems. In this paper, we describe our Hamiltonian-oriented approach to quantum benchmarking called HamPerf. Where previous application-dependent approaches specify a suite of benchmarking circuits inspired by applications, we place the problem Hamiltonian at the center. Our strategy allows us to probe the computational performance of a quantum processor on standardized and relevant problem sets, agnostic of the algorithms and hardware used to solve them; it also provides fundamental insights into how device characteristics correlate with computational utility.

Butko, Anastasiia↗

A Singlet Extension of the Minimal Supersymmetric Standard Model: Towards a More Natural Solution to the Little Hierarchy Problem

In this work, I present a generalization of the Next-to-Minimal Supersymmetric Standard Model (NMSSM), with an explicit μ-term and a supersymmetric mass for the singlet superfield, as a route to alleviating the little hierarchy problem of the Minimal Supersymmetric Standard Model (MSSM). I analyze two limiting cases of the model, characterized by the size of the supersymmetric mass for the singlet superfield. The small and large limits of this mass parameter are studied, and I find that I can generate masses for the lightest neutral Higgs boson up to 140 GeV with top squarks below the TeV scale, all couplings perturbative up to the gauge unification scale, and with no need to fine tune parameters in the scalar potential. This model, which I call the S-MSSM is also embedded in a gauge-mediated supersymmetry breaking scheme. I find that even with a minimal embedding of the S-MSSM into a gauge mediated scheme, the mass for the lightest Higgs boson can easily be above 114 GeV, while keeping the top squarks below the TeV scale. Furthermore, I also study the forward-backward asymmetry in the t¯t system within the framework of the S-MSSM. For this purpose, non-renormalizable couplings between the first and third generation of quarks to scalars are introduced. The two limiting cases of the S-MSSM, characterized by the size of the supersymmetric mass for the singlet superfield is analyzed, and I find that in the region of small singlet supersymmetric mass a large asymmetry can be obtained while being consistent with constraints arising from flavor physics, quark masses and top quark decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spectral Solutions of Self-adjoint Elliptic Problems with Immersed Interfaces

This paper describes a spectral representation of solutions of self-adjoint elliptic problems with immersed interfaces. The interface is assumed to be a simple non-self-intersecting closed curve that obeys some weak regularity conditions. The problem is decomposed into two problems, one with zero interface data and the other with zero exterior boundary data. The problem with zero interface data is solved by standard spectral methods. The problem with non-zero interface data is solved by introducing an interface space H{sub {Gamma}}({Omega}) and constructing an orthonormal basis of this space. This basis is constructed using a special class of orthogonal eigenfunctions analogously to the methods used for standard trace spaces by Auchmuty (SIAM J. Math. Anal. 38, 894-915, 2006). Analytical and numerical approximations of these eigenfunctions are described and some simulations are presented.

97 MATHEMATICS AND COMPUTING↗

A highly parallel multigrid-like method for the solution of the Euler equations

We consider a highly parallel multigrid-like method for the solution of the two dimensional steady Euler equations. The new method, introduced as filtering multigrid, is similar to a standard multigrid scheme in that convergence on the finest grid is accelerated by iterations on coarser grids. In the filtering method, however, additional fine grid subproblems are processed concurrently with coarse grid computations to further accelerate convergence. These additional problems are obtained by splitting the residual into a smooth and an oscillatory component. The smooth component is then used to form a coarse grid problem (similar to standard multigrid) while the oscillatory component is used for a fine grid subproblem. The primary advantage in the filtering approach is that fewer iterations are required and that most of the additional work per iteration can be performed in parallel with the standard coarse grid computations. We generalize the filtering algorithm to a version suitable for nonlinear problems. We emphasize that this generalization is conceptually straight-forward and relatively easy to implement. In particular, no explicit linearization (e.g., formation of Jacobians) needs to be performed (similar to the FAS multigrid approach). We illustrate the nonlinear version by applying it to the Euler equations, and presenting numerical results. Finally, a performance evaluation is made based on execution time models and convergence information obtained from numerical experiments.

Tuminaro, Ray S.↗

Development of tearing instability in a current sheet forming by sheared incompressible flow

Sweet–Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvénic incompressible flows which narrow the sheet. Standard tearing theory is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal magnetohydrodynamics, a coordinate transformation reveals that the flow compresses and stretches perturbations. Here, a multiple scale formulation describes how linear tearing mode theory can be applied to an equilibrium changing under flow, showing that the flow affects the separable exponential growth only implicitly, by making the standard scalings time dependent. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Tearing Instability in a Current Sheet Forming by Sheared Incompressible Flow

Sweet-Parker current sheets in high Lundquist number plasmas are unstable to tearing, suggesting they will not form in physical systems. Understanding magnetic reconnection thus requires study of the stability of a current sheet as it forms. Formation can occur due to sheared, sub-Alfvenic incompressible flows which narrow the sheet. Standard tearing theory (Furth et al. 1963; Coppi et al. 1976; Rutherford 1973) is not immediately applicable to such forming sheets for two reasons: first, because the flow introduces terms not present in the standard calculation; second, because the changing equilibrium introduces time dependence to terms which are constant in the standard calculation, complicating the formulation of an eigenvalue problem. This paper adapts standard tearing mode analysis to confront these challenges. In an initial phase when any perturbations are primarily governed by ideal MHD, a coordinate transformation reveals that the flow compresses and stretches perturbations. A multiple scale formulation describes how linear tearing mode theory (Furth et al. 1963; Coppi et al. 1976) can be applied to an equilibrium changing under flow, showing that tearing eigenfunctions grow slowly due to the flow and undergo exponential growth at a rate given by standard scalings with time dependence added. In the nonlinear Rutherford stage, the coordinate transformation shows that standard theory can be adapted by adding to the stationary rates time dependence and an additional term due to the strengthening equilibrium magnetic field. Overall, this understanding supports the use of flow-free scalings with slight modifications to study tearing in a forming sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Problem solving and decisionmaking: An integration

An attempt was made to redress a critical fault of decisionmaking and problem solving research-a lack of a standard method to classify problem or decision states or conditions. A basic model was identified and expanded to indicate a possible taxonomy of conditions which may be used in reviewing previous research or for systematically pursuing new research designs. A generalization of the basic conditions was then made to indicate that the conditions are essentially the same for both concepts, problem solving and decisionmaking.

Dieterly, D. L.↗

The new inflationary universe, 1984

The present status of the new inflationary theory in cosmology is discussed. The standard scenario of the very early universe, the problems of this scenario, and the basics of the new inflationary theory are reviewed. The ways in which this theory solves the problems connected with the standard scenario are described.

Guth, Alan H.↗