The automation of optimization problems.
Digital simulation of boundary value problems of trajectory optimization, using variational and functional analysis and IBM-FORMAC language
FORTRAN↗
SEARCH · Search NASA
Results for “Value function”
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.
Pareto Optimal Spectrally Selective Emitters for Thermophotovoltaics via Weak Absorber Critical Coupling
Tailoring the emission spectra of a thermophotovoltaic (TPV) emitter away from that of a blackbody has the potential to minimize transmission and thermalization loss in a photovoltaic receiver. Selective TPV emitters can lead to solar energy conversion with efficiency greater than the Shockley-Queisser limit and can facilitate the generation of useful energy from waste heat. A new design is introduced to radically tune thermal emission that leverages the interplay between two resonant phenomena in simple planar nanostructures-absorption in weakly absorbing nanofilms and reflection in multilayer dielectric stacks. A virtual screening approach is employed to identify promising structures for a selective thermal emitter from a search space of millions, several of which approach the ideal values of a step-function selective thermal emitter. One of these structures is experimentally fabricated and evaluated, which includes a weakly absorbing alloy with tailored optical properties fabricated by atomic layer deposition (ALD). The versatility of the design and fabrication approach result in an emitter with excellent spectral density (0.8 W cm(-2) sr(-1)) and spectral efficiency (46.8%) at 1373 K. Future experimental challenges to a more accurate realization of the optimal structures calculated are also considered.
atomic layer deposition↗
Measurement of t -channel production of single top quarks and antiquarks in pp collisions at 13 TeV using the full ATLAS Run 2 data sample
The production of single top quarks and top antiquarks via the t-channel exchange of a virtual W boson is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV at the LHC using 140 fb -1 of ATLAS data. The total cross-sections are determined to be σ(tq) = 137$^{+8}_{-8}$ pb and σ($\bar{t}q$) = 84$^{+6}_{-5}$ pb for top-quark and top-antiquark production, respectively. The combined cross-section is found to be σ(tq + $\bar{t}q$) = 221$^{+13}_{-13}$ pb and the cross-section ratio is R t = σ(tq) / σ($\bar{t}q$) = 1.636$^{+0.036}_{-0.034}$. The predictions at next-to-next-to-leading-order in quantum chromodynamics are in good agreement with these measurements. The predicted value of R t using different sets of parton distribution functions is compared with the measured value, demonstrating the potential to further constrain the functions when using this result in global fits. The measured cross-sections are interpreted in an effective field theory approach, setting limits at the 95% confidence level on the strength of a four-quark operator and an operator coupling the third quark generation to the Higgs boson doublet: -0.37 < C$^{3,1}_{Qq}$ /Λ 2 < 0.06 and -0.87 < C$^{3}_{ΦQ}$ / Λ 2 < 1.42. The constraint |V tb | > 0.95 at the 95% confdence level is derived from the measured value of σ(tq + $\bar{t}q$), assuming that the Wtb interaction is a left-handed weak coupling and that |V tb | $\gg$ |V td |, |V ts |. In a more general approach, pairs of CKM matrix elements involving top quarks are simultaneously constrained, leading to confdence contours in the corresponding two-dimensional parameter spaces.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Observation of Λ Hyperon Local Polarization in 𝑝-Pb Collisions at $\sqrt{s_{NN}}$ = 8.16 TeV
The polarization of the Λ and $\bar{Λ}$ hyperons along the beam direction has been measured in proton-lead (𝑝-Pb) collisions at a center-of-mass energy per nucleon pair of 8.16 TeV. The data were obtained with the CMS detector at the LHC and correspond to an integrated luminosity of 186.0 ± 6.5 nb −1 . A significant azimuthal dependence of the hyperon polarization, characterized by the second-order Fourier sine coefficient 𝑃 𝑧,s2 , is observed. The 𝑃 𝑧,s2 values decrease as a function of charged particle multiplicity, but increase with transverse momentum. A hydrodynamic model that describes the observed 𝑃 𝑧,s2 values in nucleus-nucleus collisions by introducing vorticity effects does not reproduce either the sign or the magnitude of the 𝑝-Pb results. These observations pose a challenge to the current theoretical implementation of spin polarization in heavy ion collisions and offer new insights into the origin of spin polarization in hadronic collisions at LHC energies.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
On propagation of long waves in curved ducts.
Propagation of waves in curved ducts and pipes belong to the class of motion which is characterized by wave patterns totally different from those known in straight ducts or in unlimited space. The curvilinear boundaries are responsible for the appearance of a continuous standing radial wave which in turn affects the transmitted tangential waves. The purpose of this paper is to solve the problem of propagation of long acoustic waves in slightly and sharply bent ducts. This problem has been only partially analyzed by various authors. In this study two acoustic systems are considered which allow determination of the basic modes of motion and describe the transition and distortion of plane waves as they propagate down the curved channel. A detailed study of the behavior of waves in junctions between straight and curved ducts is also given. Solutions and expressions for principal modes of the wave are obtained by using the linearized equation of motion solved for its characteristic values. This original approach required a novel use of Bessel functions to determine the characteristic values of the steady and the decaying fields of motion.
Rostafinski, W.↗
Propagation of waves in elliptic ducts. A theoretical study
The cut-off frequencies for high order circumferential modes were calculated for various eccentricities of an elliptic duct section. The problem was studied with a view to the reduction of jet engine compressor noise by elliptic ducts, instead of circular ducts. The cut-off frequencies for even functions decrease with increasing eccentricity. The third order eigen frequencies are oscillatory as the eccentricity increases for odd functions. The eigen frequencies decrease for higher order odd functions inasmuch as, for higher orders, they assume the same values as those for even functions. Deformation of a circular pipe into an elliptic one of sufficiently large eccentricity produces only a small reduction in the cut-off frequency, provided the area of the pipe section is kept invariable.
Baskaran, S.↗
Use of shape-preserving interpolation methods in surface modeling
In many large-scale scientific computations, it is necessary to use surface models based on information provided at only a finite number of points (rather than determined everywhere via an analytic formula). As an example, an equation of state (EOS) table may provide values of pressure as a function of temperature and density for a particular material. These values, while known quite accurately, are typically known only on a rectangular (but generally quite nonuniform) mesh in (T,d)-space. Thus interpolation methods are necessary to completely determine the EOS surface. The most primitive EOS interpolation scheme is bilinear interpolation. This has the advantages of depending only on local information, so that changes in data remote from a mesh element have no effect on the surface over the element, and of preserving shape information, such as monotonicity. Most scientific calculations, however, require greater smoothness. Standard higher-order interpolation schemes, such as Coons patches or bicubic splines, while providing the requisite smoothness, tend to produce surfaces that are not physically reasonable. This means that the interpolant may have bumps or wiggles that are not supported by the data. The mathematical quantification of ideas such as physically reasonable and visually pleasing is examined.
Ftitsch, F. N.↗
A viscoplastic model with application to LiF-22 percent CaF2 hypereutectic salt
A viscoplastic model for class M (metal-like behavior) materials is presented. One novel feature is its use of internal variables to change the stress exponent of creep (where n is approximately = 5) to that of natural creep (where n = 3), in accordance with experimental observations. Another feature is the introduction of a coupling in the evolution equations of the kinematic and isotropic internal variables, making thermal recovery of the kinematic variable implicit. These features enable the viscoplastic model to reduce to that of steady-state creep in closed form. In addition, the hardening parameters associated with the two internal state variables (one scalar-valued, the other tensor-valued) are considered to be functions of state, instead of being taken as constant-valued. This feature enables each internal variable to represent a much wider spectrum of internal states for the material. The model is applied to a LiF-22 percent CaF2 hypereutectic salt, which is being considered as a thermal energy storage material for space-based solar dynamic power systems.
Freed, A. D.↗
Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets
Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Time correlations from steady-state expectation values
Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.
Górecki, Wojciech [INFN, Pavia]↗
Time correlations from steady-state expectation values
Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.
Górecki, Wojciech [INFN, Pavia]↗
Time Correlations from Steady-State Expectation Values
Recovering properties of correlation functions is typically challenging. On the one hand, experimentally, it requires measurements with a temporal resolution finer than the system’s dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a system parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable, and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to the experimental characterization of ultrafast systems and to the theoretical analysis of many-body models whose dynamics are hard to compute.
Górecki, Wojciech [INFN, Pavia] (ORCID:00000001991↗
Some remarks on elastic crack-tip stress fields.
It is shown that if the displacement field and stress intensity factor are known as functions of crack length for any symmetrical load system acting on a linear elastic body in plane strain, then the stress intensity factor for any other symmetrical load system whatsoever on the same body may be directly determined. The result is closely related to Bueckner's (1970) weight function, through which the stress intensity factor is expressed as a sum of work-like products between applied forces and values of the weight function at their points of application. An example of the method is given wherein the solution for a crack in a remotely uniform stress field is used to generate the expression for the stress intensity factor due to an arbitrary traction distribution on the faces of a crack. A corresponding theory is developed in an appendix for three-dimensional crack problems, although this appears to be directly useful chiefly for problems in which there is axial symmetry.
Rice, J. R.↗
Sounder design considerations in the selection of temperature sensing channels
Simulation studies on the performance of IR sounders under varying conditions of cloud cover and cloud heights are carried out for Nimbus 6. An analytic function is derived for calculating the relative response to cloud height errors for arbitrary cloud-sensing channels. Based on the values of the response function, the best choice of channels for determining cloud amounts are obtained. An algorithm is described for determining cloud heights, and the sensitivity of cloud-height sensing channels is tested. It is found that for the HIRS instrument, the most transparent channel in the 4.3-micron band is optimal for adjusting cloud heights, while the channel in the 15-micron band peaking closest to the surface is best for determining cloud amounts
Halem, M.↗
Structural optimization via a design space hierarchy
Mathematical programming techniques provide a general approach to automated structural design. An iterative method is proposed in which design is treated as a hierarchy of subproblems, one being locally constrained and the other being locally unconstrained. It is assumed that the design space is locally convex in the case of good initial designs and that the objective and constraint functions are continuous, with continuous first derivatives. A general design algorithm is outlined for finding a move direction which will decrease the value of the objective function while maintaining a feasible design. The case of one-dimensional search in a two-variable design space is discussed. Possible applications are discussed. A major feature of the proposed algorithm is its application to problems which are inherently ill-conditioned, such as design of structures for optimum geometry.
Vanderplaats, G. N.↗
Proposed Project Selection Method for Human Support Research and Technology Development (HSR&TD)
The purpose of HSR&TD is to deliver human support technologies to the Exploration Systems Mission Directorate (ESMD) that will be selected for future missions. This requires identifying promising candidate technologies and advancing them in technology readiness until they are acceptable. HSR&TD must select an may of technology development projects, guide them, and either terminate or continue them, so as to maximize the resulting number of usable advanced human support technologies. This paper proposes an effective project scoring methodology to support managing the HSR&TD project portfolio. Researchers strongly disagree as to what are the best technology project selection methods, or even if there are any proven ones. Technology development is risky and outstanding achievements are rare and unpredictable. There is no simple formula for success. Organizations that are satisfied with their project selection approach typically use a mix of financial, strategic, and scoring methods in an open, established, explicit, formal process. This approach helps to build consensus and develop management insight. It encourages better project proposals by clarifying the desired project attributes. We propose a project scoring technique based on a method previously used in a federal laboratory and supported by recent research. Projects are ranked by their perceived relevance, risk, and return - a new 3 R's. Relevance is the degree to which the project objective supports the HSR&TD goal of developing usable advanced human support technologies. Risk is the estimated probability that the project will achieve its specific objective. Return is the reduction in mission life cycle cost obtained if the project is successful. If the project objective technology performs a new function with no current cost, its return is the estimated cash value of performing the new function. The proposed project selection scoring method includes definitions of the criteria, a project evaluation questionnaire, and a scoring formula.
Jones, Harry↗
Coefficients of associated Legendre functions
The exact coefficients for the explicit forms of the associated Legendre functions Pm for integer values of m,n=0,1,2,...25 are presented in tabular form together with two cross-referenced listings of the zeroes for these functions rounded to five decimal places. The unfactored coefficients and the interger coefficients are presented in adjacent columns for each function. The greatest common factor and divisor have been removed and listed separately for the integer coefficients.
Peasley, Q. D.↗
Using Redundancy To Reduce Errors in Magnetometer Readings
A method of reducing errors in noisy magnetic-field measurements involves exploitation of redundancy in the readings of multiple magnetometers in a cluster. By "redundancy"is meant that the readings are not entirely independent of each other because the relationships among the magnetic-field components that one seeks to measure are governed by the fundamental laws of electromagnetism as expressed by Maxwell's equations. Assuming that the magnetometers are located outside a magnetic material, that the magnetic field is steady or quasi-steady, and that there are no electric currents flowing in or near the magnetometers, the applicable Maxwell 's equations are delta x B = 0 and delta(raised dot) B = 0, where B is the magnetic-flux-density vector. By suitable algebraic manipulation, these equations can be shown to impose three independent constraints on the values of the components of B at the various magnetometer positions. In general, the problem of reducing the errors in noisy measurements is one of finding a set of corrected values that minimize an error function. In the present method, the error function is formulated as (1) the sum of squares of the differences between the corrected and noisy measurement values plus (2) a sum of three terms, each comprising the product of a Lagrange multiplier and one of the three constraints. The partial derivatives of the error function with respect to the corrected magnetic-field component values and the Lagrange multipliers are set equal to zero, leading to a set of equations that can be put into matrix.vector form. The matrix can be inverted to solve for a vector that comprises the corrected magnetic-field component values and the Lagrange multipliers.
Kulikov, Igor↗