Search NASA⌕ Search

SEARCH · Search NASA

Results for “MONOTONE METHOD”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A O(1/k{sup 3/2}) hybrid proximal extragradient primal–dual interior point method for nonlinear monotone mixed complementarity problems

We present a Newton-type hybrid proximal extragradient primal–dual interior point method for solving smooth monotone mixed complementarity problems. Dual variables for the nonnegativity constraints are introduced. The ergodic complexity of the method is O(1/k{sup 3/2}). The method performs two types of iterations: underelaxed hybrid proximal extragradient iterations and short-steps primal–dual interior point iterations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

An Unconditionally Monotone C 2 Quartic Spline Method with Nonoscillation Derivatives

Here, a one-dimensional monotone interpolation method based on interface reconstruction with partial volumes in the slope-space utilizing the Hermite cubic-spline, is proposed. The new method is only quartic, however is C 2 and unconditionally monotone. A set of control points is employed to constrain the curvature of the interpolation function and to eliminate possible nonphysical oscillations in the slope space. An extension of this method in two-dimensions is also discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Evaluation of kriging based surrogate models constructed from mesoscale computations of shock interaction with particles

Macro-scale computations of shocked particulate flows require closure laws that model the exchange of momentum/energy between the fluid and particle phases. Closure laws are constructed in this work in the form of surrogate models derived from highly resolved mesoscale computations of shock-particle interactions. The mesoscale computations are performed to calculate the drag force on a cluster of particles for different values of Mach Number and particle volume fraction. Two Kriging-based methods, viz. the Dynamic Kriging Method (DKG) and the Modified Bayesian Kriging Method (MBKG) are evaluated for their ability to construct surrogate models with sparse data; i.e. using the least number of mesoscale simulations. It is shown that if the input data is noise-free, the DKG method converges monotonically; convergence is less robust in the presence of noise. The MBKG method converges monotonically even with noisy input data and is therefore more suitable for surrogate model construction from numerical experiments. This work is the first step towards a full multiscale modeling of interaction of shocked particle laden flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

On the Measurements of Numerical Viscosity and Resistivity in Eulerian MHD Codes

We propose a simple ansatz for estimating the value of the numerical resistivity and the numerical viscosity of any Eulerian MHD code. We test this ansatz with the help of simulations of the propagation of (magneto)sonic waves, Alfvén waves, and the tearing mode (TM) instability using the MHD code Aenus. By comparing the simulation results with analytical solutions of the resistive-viscous MHD equations and an empirical ansatz for the growth rate of TMs, we measure the numerical viscosity and resistivity of Aenus. The comparison shows that the fast magnetosonic speed and wavelength are the characteristic velocity and length, respectively, of the aforementioned (relatively simple) systems. We also determine the dependence of the numerical viscosity and resistivity on the time integration method, the spatial reconstruction scheme and (to a lesser extent) the Riemann solver employed in the simulations. From the measured results, we infer the numerical resolution (as a function of the spatial reconstruction method) required to properly resolve the growth and saturation level of the magnetic field amplified by the magnetorotational instability in the post-collapsed core of massive stars. Our results show that it is most advantageous to resort to ultra-high-order methods (e.g., the ninth-order monotonicity-preserving method) to tackle this problem properly, in particular, in three-dimensional simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Globally convergent algorithms for finding zeros of duplomonotone mappings

We introduce a new class of mappings, called duplomonotone, which is strictly broader than the class of monotone mappings.We study some of themain properties of duplomonotone functions and provide various examples, including nonlinear duplomonotone functions arising from the study of systems of biochemical reactions. Finally,we present three variations of a derivative-free line search algorithm for finding zeros of systems of duplomonotone equations, and we prove their linear convergence to a zero of the function.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Mean Field Type Control with Congestion (II): An Augmented Lagrangian Method

This work deals with a numerical method for solving a mean-field type control problem with congestion. It is the continuation of an article by the same authors, in which suitably defined weak solutions of the system of partial differential equations arising from the model were discussed and existence and uniqueness were proved. Here, the focus is put on numerical methods: a monotone finite difference scheme is proposed and shown to have a variational interpretation. Then an Alternating Direction Method of Multipliers for solving the variational problem is addressed. It is based on an augmented Lagrangian. Two kinds of boundary conditions are considered: periodic conditions and more realistic boundary conditions associated to state constrained problems. Various test cases and numerical results are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Nonlinear Integral Equations with Potential-Type Kernels on a Segment

We study various classes of nonlinear equations containing operators of potential type (Riesz potential). By the method of monotone operators in the Lebesgue spaces of real-valued functions L{sub p}(a, b) we prove global theorems on the existence, uniqueness, estimates, and methods of construction of their solutions. We present applications that illustrate the results obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

CFD simulation of simultaneous monotonic cooling and surface heat transfer coefficient

The monotonic heating regime method for determination of thermal diffusivity is based on the analysis of an unsteady-state (stabilised) thermal process characterised by an independence of the space-time temperature distribution on initial conditions. At the first kind of the monotonic regime a sample of simple geometry is heated / cooled at constant ambient temperature. The determination of thermal diffusivity requires the determination rate of a temperature change and simultaneous determination of the first eigenvalue. According to a characteristic equation the first eigenvalue is a function of the Biot number defined by a surface heat transfer coefficient and thermal conductivity of an analysed material. Knowing the surface heat transfer coefficient and the first eigenvalue the thermal conductivity can be determined. The surface heat transport coefficient during the monotonic regime can be determined by the continuous measurement of long-wave radiation heat flow and the photoelectric measurement of the air refractive index gradient in a boundary layer. CFD simulation of the cooling process was carried out to analyse local convective and radiative heat transfer coefficients more in detail. Influence of ambient air flow was analysed. The obtained eigenvalues and corresponding surface heat transfer coefficient values enable to determine thermal conductivity of the analysed specimen together with its thermal diffusivity during a monotonic heating regime.

42 ENGINEERING↗
DOE OSTI2023

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗
DOE OSTI2016

Solution of monotone complementarity and general convex programming problems using a modified potential reduction interior point method

We present a homogeneous algorithm equipped with a modified potential function for the monotone complementarity problem. We show that this potential function is reduced by at least a constant amount if a scaled Lipschitz condition (SLC) is satisfied. A practical algorithm based on this potential function is implemented in a software package named iOptimize. The implementation in iOptimize maintains global linear and polynomial time convergence properties, while achieving practical performance. It either successfully solves the problem, or concludes that the SLC is not satisfied. When compared with the mature software package MOSEK (barrier solver version 6.0.0.106), iOptimize solves convex quadratic programming problems, convex quadratically constrained quadratic programming problems, and general convex programming problems in fewer iterations. Moreover, several problems for which MOSEK fails are solved to optimality. In addition, we also find that iOptimize detects infeasibility more reliably than the general nonlinear solvers Ipopt (version 3.9.2) and Knitro (version 8.0).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Pan-Arctic seasonal cycles and long-term trends of aerosol properties from 10 observatories

Even though the Arctic is remote, aerosol properties observed there are strongly influenced by anthropogenic emissions from outside the Arctic. This is particularly true for the so-called Arctic haze season (January through April). In summer (June through September), when atmospheric transport patterns change, and precipitation is more frequent, local Arctic sources, i.e., natural sources of aerosols and precursors, play an important role. Over the last few decades, significant reductions in anthropogenic emissions have taken place. At the same time a large body of literature shows evidence that the Arctic is undergoing fundamental environmental changes due to climate forcing, leading to enhanced emissions by natural processes that may impact aerosol properties. In this study, we analyze 9 aerosol chemical species and 4 particle optical properties from 10 Arctic observatories (Alert, Kevo, Pallas, Summit, Thule, Tiksi, Barrow/Utgiavik, Villum, and Gruvebadet and Zeppelin Observatory – both at Ny-Ålesund Research Station) to understand changes in anthropogenic and natural aerosol contributions. Variables include equivalent black carbon, particulate sulfate, nitrate, ammonium, methanesulfonic acid, sodium, iron, calcium and potassium, as well as scattering and absorption coefficients, single scattering albedo and scattering Ångström exponent. First, annual cycles are investigated, which despite anthropogenic emission reductions still show the Arctic haze phenomenon. Second, long-term trends are studied using the Mann–Kendall Theil–Sen slope method. We find in total 41 significant trends over full station records, i.e., spanning more than a decade, compared to 26 significant decadal trends. The majority of significantly declining trends is from anthropogenic tracers and occurred during the haze period, driven by emission changes between 1990 and 2000. For the summer period, no uniform picture of trends has emerged. Twenty-six percent of trends, i.e., 19 out of 73, are significant, and of those 5 are positive and 14 are negative. Negative trends include not only anthropogenic tracers such as equivalent black carbon at Kevo, but also natural indicators such as methanesulfonic acid and non-sea-salt calcium at Alert. Positive trends are observed for sulfate at Gruvebadet. No clear evidence of a significant change in the natural aerosol contribution can be observed yet. However, testing the sensitivity of the Mann–Kendall Theil–Sen method, we find that monotonic changes of around 5 % yr –1 in an aerosol property are needed to detect a significant trend within one decade. This highlights that long-term efforts well beyond a decade are needed to capture smaller changes. It is particularly important to understand the ongoing natural changes in the Arctic, where interannual variability can be high, such as with forest fire emissions and their influence on the aerosol population. To investigate the climate-change-induced influence on the aerosol population and the resulting climate feedback, long-term observations of tracers more specific to natural sources are needed, as well as of particle microphysical properties such as size distributions, which can be used to identify changes in particle populations which are not well captured by mass-oriented methods such as bulk chemical composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Recovering the QNEC from the ANEC

We study the relative entropy in QFT comparing the vacuum state to a special family of purifications determined by an input state and constructed using relative modular flow. We use this to prove a conjecture by Wall that relates the shape derivative of relative entropy to a variational expression over the averaged null energy (ANE) of possible purifications. Furthermore, this variational expression can be used to easily prove the quantum null energy condition (QNEC). We formulate Wall’s conjecture as a theorem pertaining to operator algebras satisfying the properties of a half-sided modular inclusion, with the additional assumption that the input state has finite averaged null energy. We also give a new derivation of the strong superadditivity property of relative entropy in this con- text. We speculate about possible connections to the recent methods used to strengthen monotonicity of relative entropy with recovery maps.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗
DOE OSTI2022

A Quantized State Integrator With Second Order Errors Over Monotonic Segments

We introduce a novel quantized state integration method that is analogous to the explicit, second order Runge-Kutta technique. A feature of this method is that, for a single differential equation over an interval where the solution is monotonic, it exhibits an error proportional to the square of the integration quantum. We offer a theoretical proof of this behavior and demonstrate it with a numerical example. At the same time, an extension of the method to a system with input causes the errors to become proportional to the integration quantum. The latter observation fits established theory; the former appears to be a new result that might offer new insights into the construction of quantized state integration methods.

Mahawattege, Rasika↗