Search NASA⌕ Search

SEARCH · Search NASA

Results for “Derived”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Error analysis for semi-analytic displacement derivatives with respect to shape and sizing variables

Sensitivity analysis is fundamental to the solution of structural optimization problems. Consequently, much research has focused on the efficient computation of static displacement derivatives. As originally developed, these methods relied on analytical representations for the derivatives of the structural stiffness matrix (K) with respect to the design variables (b sub i). To extend these methods for use with complex finite element formulations and facilitate their implementation into structural optimization programs using the general finite element method analysis codes, the semi-analytic method was developed. In this method the matrix the derivative of K/the derivative b sub i is approximated by finite difference. Although it is well known that the accuracy of the semi-analytic method is dependent on the finite difference parameter, recent work has suggested that more fundamental inaccuracies exist in the method when used for shape optimization. Another study has argued qualitatively that these errors are related to nonuniform errors in the stiffness matrix derivatives. The accuracy of the semi-analytic method is investigated. A general framework was developed for the error analysis and then it is shown analytically that the errors in the method are entirely accounted for by errors in delta K/delta b sub i. Furthermore, it is demonstrated that acceptable accuracy in the derivatives can be obtained through careful selection of the finite difference parameter.

Fenyes, Peter A.↗
NASA NTRS1994

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗
NASA NTRS1995

Thermal inertias in the upper millimeters of the Martian surace derived using Phobus' shadow

The first thermal images of Phobos' shadow on the surface of Mars, in addition to simultaneous visible images, were obtained by the Phobus '88 Termoskan instrument. The best observed shadow occurence was on the flanks of Arsia Mons. For this occurence, we combined the observed decrease in visible illumination of the surface with the observed decrease in brightness temperature to calculate thermal inertias of the Martian surface. The most realistic of our three models of eclipse cooling improves upon our preliminary model by including nonisothermal initial conditions and downward atmospheric flux. Most of our derived inertias fall within the range 38 to 59 J/Sq m/S(exp 0.5)K (0.9 to 1.4 10(exp -3)Cal/Sq m/S(exp 0.5)/K), corresponding to dust-sized particles (for a homogeneous surface), consistent with previous theories of Tharsis as a currrent area of dust deposition. Viking infrared thermal mapper (IRTM) inertias are diurnally derived and are sensitive to centimeter depths, whereas the shadow-derived inertias sample the upper tenths of a millimeter of the surface. The shadow-derived inertias are lower than those derived from Viking IRTM measurements (84 to 147), however, uncertainties in both sets of derived inertias make conclusions about layering tenuous. Thus, near-surface millimeter versus centimeter layering may exist in this region, but if it does, it is likely not very significant. Both eclipse and diurnal inertias appear to increase near the eastern end of the shadow occurence. We also analyzed a shadow occurence near the crater Herschel that showed no observed cooling. This analysis was limited by cool morning temperatures and instrument sensitivity, but yielded a lower bound of 80 on eclipse inertias in that region. Based upon our results, we strongly recommend future spacecraft thermal observations of Phobus' shadow, and suggest that they will be most useful if they improve upon Termoskan's geographic and temporal coverage and its accuracy.

Betts, Bruce H.↗
NASA NTRS1995

The scheme machine: A case study in progress in design derivation at system levels

The Scheme Machine is one of several design projects of the Digital Design Derivation group at Indiana University. It differs from the other projects in its focus on issues of system design and its connection to surrounding research in programming language semantics, compiler construction, and programming methodology underway at Indiana and elsewhere. The genesis of the project dates to the early 1980's, when digital design derivation research branched from the surrounding research effort in programming languages. Both branches have continued to develop in parallel, with this particular project serving as a bridge. However, by 1990 there remained little real interaction between the branches and recently we have undertaken to reintegrate them. On the software side, researchers have refined a mathematically rigorous (but not mechanized) treatment starting with the fully abstract semantic definition of Scheme and resulting in an efficient implementation consisting of a compiler and virtual machine model, the latter typically realized with a general purpose microprocessor. The derivation includes a number of sophisticated factorizations and representations and is also deep example of the underlying engineering methodology. The hardware research has created a mechanized algebra supporting the tedious and massive transformations often seen at lower levels of design. This work has progressed to the point that large scale devices, such as processors, can be derived from first-order finite state machine specifications. This is roughly where the language oriented research stops; thus, together, the two efforts establish a thread from the highest levels of abstract specification to detailed digital implementation. The Scheme Machine project challenges hardware derivation research in several ways, although the individual components of the system are of a similar scale to those we have worked with before. The machine has a custom dual-ported memory to support garbage collection. It consists of four tightly coupled processes--processor, collector, allocator, memory--with a very non-trivial synchronization relationship. Finally, there are deep issues of representation for the run-time objects of a symbolic processing language. The research centers on verification through integrated formal reasoning systems, but is also involved with modeling and prototyping environments. Since the derivation algebra is basd on an executable modeling language, there is opportunity to incorporate design animation in the design process. We are looking for ways to move smoothly and incrementally from executable specifications into hardware realization. For example, we can run the garbage collector specification, a Scheme program, directly against the physical memory prototype, and similarly, the instruction processor model against the heap implementation.

Johnson, Steven D.↗
NASA NTRS1995

Thermal inertias in the upper millimeters of the Martian surface derived using Phobos' shadow

The first thermal images of Phobos' shadow on the surface of Mars, in addition to simultaneous visible images, were obtained by the Phobos'88 Termoskan instrument. The best observed shadow occurrence was on the flanks of Arsia Mons. For this occurrence, we combined the observed decrease in visible illumination of the surface with the observed decrease in brightness temperature to calculate thermal inertias of the Martian surface. The most realistic of our three models of eclipse cooling improves upon our preliminary model by including nonisothermal initial conditions and downward atmospheric flux. Most of our derived inertias fall within the range 38 to 59 J/(sq m s(exp 1/2) K), (0.9 to 1.4 x 10(exp -3) cal/(sq cm s(exp 1/2) K)) corresponding to dust-sized particles (for a homogeneous surface), consistent with previous theories of Tharsis as a current area of dust deposition. Viking infrared thermal mapper (IRTM) inertias are diurnally derived and are sensitive to centimeter depths, whereas the shadow-derived inertias sample the upper tenths of a millimeter of the surface. The shadow-derived inertias are lower than those derived from Viking IRTM measurements (84 to 147), however, uncertainties in both sets of derived inertias make conclusions about layering tenuous. Thus, near-surface millimeter versus centimeter layering may exist in this region, but if it does, it is likely not very significant. Both eclipse and diurnal inertias appear to increase near the eastern end of the shadow occurrence. We also analyzed a shadow occurrence near the crater Herschel that showed no observed cooling. This analysis was limited by cool morning temperatures and instrument sensitivity, but yielded a lower bound of 80 on eclipse inertias in that region. Based upon our results, we strongly recommend future spacecraft thermal observations of Phobos' shadow, and suggest that they will be most useful if they improve upon Terinoskan's geographic and temporal coverage and its accuracy.

Betts, Bruce H.↗
NASA NTRS2000

Steady-State Computation of Constant Rotational Rate Dynamic Stability Derivatives

Dynamic stability derivatives are essential to predicting the open and closed loop performance, stability, and controllability of aircraft. Computational determination of constant-rate dynamic stability derivatives (derivatives of aircraft forces and moments with respect to constant rotational rates) is currently performed indirectly with finite differencing of multiple time-accurate computational fluid dynamics solutions. Typical time-accurate solutions require excessive amounts of computational time to complete. Formulating Navier-Stokes (N-S) equations in a rotating noninertial reference frame and applying an automatic differentiation tool to the modified code has the potential for directly computing these derivatives with a single, much faster steady-state calculation. The ability to rapidly determine static and dynamic stability derivatives by computational methods can benefit multidisciplinary design methodologies and reduce dependency on wind tunnel measurements. The CFL3D thin-layer N-S computational fluid dynamics code was modified for this study to allow calculations on complex three-dimensional configurations with constant rotation rate components in all three axes. These CFL3D modifications also have direct application to rotorcraft and turbomachinery analyses. The modified CFL3D steady-state calculation is a new capability that showed excellent agreement with results calculated by a similar formulation. The application of automatic differentiation to CFL3D allows the static stability and body-axis rate derivatives to be calculated quickly and exactly.

Park, Michael A.↗
NASA NTRS2002

Non-Spherical Aerosol Phase Functions Derived from MODIS and AERONET Observations

We compare MODIS (Moderate Resolution Imaging Spectroradiometer) satellite aerosol retrievals of spectral optical thickness and size parameters over ocean with the same quantities derived from AERONET (Aerosol Robotic Network) observations made at island and coastal sites. Over much of the globe, the satellite-derived quantities agree well with the AERONET quantities. However, in regimes dominated by desert dust aerosol, the agreement is less robust. In the dust regimes, the MODIS retrievals show greater spectral dependence and report smaller particle sizes than do the AERONET derivations. We suggest that the reason for this discrepancy is the nonspherical nature of desert dust particles, which the initial MODIS algorithm is not able to handle. Using the discrepancy between MODIS and AERONET derived spectral optical thickness as an asset, instead of a detriment, we reconstruct the aerosol phase functions that the MODIS algorithm would have needed in order to match the AERONET retrievals. No assumptions of particle shape are used in the derivation of these functions and the results are empirical total column, ambient phase functions. We compare the empirically derived phase functions with phase functions calculated from spheres and spheroids, both situations in which assumptions about particle shape must be made. The resulting empirical nonspherical phase functions will be included in future updates of the MODIS algorithm.

Remer, L. A.↗
NASA NTRS2004

Employing Sensitivity Derivatives to Estimate Uncertainty Propagation in CFD

Two methods that exploit the availability of sensitivity derivatives are successfully employed to predict uncertainty propagation through Computational Fluid Dynamics (CFD) code for an inviscid airfoil problem. An approximate statistical second-moment method and a Sensitivity Derivative Enhanced Monte Carlo (SDEMC) method are successfully demonstrated on a two-dimensional problem. First- and second-order sensitivity derivatives of code output with respect to code input are obtained through an efficient incremental iterative approach. Given uncertainties in statistically independent, random, normally distributed flow parameters (input variables); these sensitivity derivatives enable one to formulate first- and second-order Taylor Series approximations for the mean and variance of CFD output quantities. Additionally, incorporation of the first-order sensitivity derivatives into the data reduction phase of a conventional Monte Carlo (MC) simulation allows for improved accuracy in determining the first moment of the CFD output. Both methods are compared to results generated using a conventional MC method. The methods that exploit the availability of sensitivity derivatives are found to be valid when considering small deviations from input mean values.

Putko, Michele M.↗
NASA NTRS1999

Determination of Stability and Control Derivatives using Computational Fluid Dynamics and Automatic Differentiation

With the recent interest in novel control effectors there is a need to determine the stability and control derivatives of new aircraft configurations early in the design process. These derivatives are central to most control law design methods and would allow the determination of closed-loop control performance of the vehicle. Early determination of the static and dynamic behavior of an aircraft may permit significant improvement in configuration weight, cost, stealth, and performance through multidisciplinary design. The classical method of determining static stability and control derivatives - constructing and testing wind tunnel models - is expensive and requires a long lead time for the resultant data. Wind tunnel tests are also limited to the preselected control effectors of the model. To overcome these shortcomings, computational fluid dynamics (CFD) solvers are augmented via automatic differentiation, to directly calculate the stability and control derivatives. The CFD forces and moments are differentiated with respect to angle of attack, angle of sideslip, and aircraft shape parameters to form these derivatives. A subset of static stability and control derivatives of a tailless aircraft concept have been computed by two differentiated inviscid CFD codes and verified for accuracy with central finite-difference approximations and favorable comparisons to a simulation database.

Park, Michael A.↗
NASA NTRS1953

Lateral Oscillatory Characteristics of the Republic F-91 Airplane Calculated by Using Low-Speed Experimental Static and Rotary Derivatives

The present investigation was conducted to determine, from low-speed tests in the Langley stability tunnel, the static and rotary derivatives of a 1/9-scale model of the Republic F-91 airplane and various of its components (including the effects of wing incidence) and to determine the accuracy with which the period and damping of the lateral oscillation of the airplane could be calculated by using these experimentally between flight and calculated period and damping of the lateral oscillation were made for Mach numbers from 0.4 to 0.9 at an altitude of 20,OOO feet for 0deg wing incidence and several other wing incidences. Some comparisons were made of the static and rotary derivatives of the model and derivatives estimated by available procedures. determined derivatives (corrected for Mach number effects). Comparisons The results of the investigation have indicated that the model did not have unusual aerodynamic characteristics except for a large (about -0.125) increment in the damping in yaw contributed by the fuselage. Changes in wing incidence, in general, had little effect on the static and rotary derivatives of the model. The static and rotary derivatives of the model could be estimated with good accuracy only in the low angle-of-attack range by using available procedures.

Jaquet, Byron M.↗
NASA NTRS2007

Derivation of Formulations 1 and 1A of Farassat

Formulations 1 and 1A are the solutions of the Ffowcs Williams-Hawkings (FW-H) equation with surface sources only when the surface moves at subsonic speed. Both formulations have been successfully used for helicopter rotor and propeller noise prediction for many years although we now recommend using Formulation 1A for this purpose. Formulation 1 has an observer time derivative that is taken numerically, and thus, increasing execution time on a computer and reducing the accuracy of the results. After some discussion of the Green's function of the wave equation, we derive Formulation 1 which is the basis of deriving Formulation 1A. We will then show how to take this observer time derivative analytically to get Formulation 1A. We give here the most detailed derivation of these formulations. Once you see the whole derivation, you will ask yourself why you did not do it yourself!

Farassat, F.↗
NASA NTRS2006

Estimation of Stability and Control Derivatives of an F-15

A technique for real-time estimation of stability and control derivatives (derivatives of moment coefficients with respect to control-surface deflection angles) was used to support a flight demonstration of a concept of an indirect-adaptive intelligent flight control system (IFCS). Traditionally, parameter identification, including estimation of stability and control derivatives, is done post-flight. However, for the indirect-adaptive IFCS concept, parameter identification is required during flight so that the system can modify control laws for a damaged aircraft. The flight demonstration was carried out on a highly modified F-15 airplane (see Figure 1). The main objective was to estimate the stability and control derivatives of the airplane in nearly real time. A secondary goal was to develop a system to automatically assess the quality of the results, so as to be able to tell a learning neural network which data to use. Parameter estimation was performed by use of Fourier-transform regression (FTR) a technique developed at NASA Langley Research Center. FTR is an equation- error technique that operates in the frequency domain. Data are put into the frequency domain by use of a recursive Fourier transform for a discrete frequency set. This calculation simplifies many subsequent calculations, removes biases, and automatically filters out data beyond the chosen frequency range. FTR as applied here was tailored to work with pilot inputs, which produce correlated surface positions that prevent accurate parameter estimates, by replacing half the derivatives with predicted values. FTR was also set up to work only on a recent window of data, to accommodate changes in flight condition. A system of confidence measures was developed to identify quality-parameter estimates that a learning neural network could use. This system judged the estimates primarily on the basis of their estimated variances and of the level of aircraft response. The resulting FTR system was implemented in the Simulink software system and auto-coded in the C programming language for use on the Airborne Research Test System (ARTS II) computer installed in the F-15 airplane. The Simulink model was also used in a control room that utilizes the Ring Buffered Network Bus hardware and software, making it possible to evaluate test points during flights. In-flight parameter estimation was done for piloted and automated maneuvers, primarily at three test conditions. Figure 2 shows results for pitching moment due to symmetric stabilator actuations for a series of three pitch doublet maneuvers (in a doublet maneuver, a command to change attitude in a given direction by a given amount is followed immediately by a command to change attitude in the opposite direction by the same amount). A time window of 5 seconds was used. The portions of the curves shown in red are those that passed the confidence tests. The technique showed good convergence for most derivatives for both kinds of maneuvers - typically within a few seconds. The confidence tests were marginally successful, and it would be necessary to refine them for use in an IFCS.

Smith, Mark↗
NASA NTRS2013

The Importance of Physical Models for Deriving Dust Masses and Grain Size Distributions in Supernova Ejecta. I. Radiatively Heated Dust in the Crab Nebula

Recent far-infrared (IR) observations of supernova remnants (SNRs) have revealed significantly large amounts of newly condensed dust in their ejecta, comparable to the total mass of available refractory elements. The dust masses derived from these observations assume that all the grains of a given species radiate at the same temperature, regardless of the dust heating mechanism or grain radius. In this paper, we derive the dust mass in the ejecta of the Crab Nebula, using a physical model for the heating and radiation from the dust. We adopt a power-law distribution of grain sizes and two different dust compositions (silicates and amorphous carbon), and calculate the heating rate of each dust grain by the radiation from the pulsar wind nebula. We find that the grains attain a continuous range of temperatures, depending on their size and composition. The total mass derived from the best-fit models to the observed IR spectrum is 0.019-0.13 Solar Mass, depending on the assumed grain composition. We find that the power-law size distribution of dust grains is characterized by a power-law index of 3.5-4.0 and a maximum grain size larger than 0.1 micron. The grain sizes and composition are consistent with what is expected for dust grains formed in a Type IIP supernova (SN). Our derived dust mass is at least a factor of two less than the mass reported in previous studies of the Crab Nebula that assumed more simplified two-temperature models. These models also require a larger mass of refractory elements to be locked up in dust than was likely available in the ejecta. The results of this study show that a physical model resulting in a realistic distribution of dust temperatures can constrain the dust properties and affect the derived dust masses. Our study may also have important implications for deriving grain properties and mass estimates in other SNRs and for the ultimate question of whether SNe are major sources of dust in the Galactic interstellar medium and in external galaxies.

Importance↗
NASA NTRS2020

Comparison of Vertical Surface Deformation Estimates Derived From Space-Based Gravimetry, Ground-Based GPS, and Model-Based Hydrologic Loading Over Snow-Dominated Watersheds in the United States

Spatiotemporal variability in Earth's terrestrial water storage (TWS) causes changes in surface deformation. The potential for using ground-based Global Positioning System (GPS) vertical displacement observations for estimating TWS is explored through a comparison of vertical displacements derived from space-based gravimetric retrievals, ground-based GPS, and model-based hydrologic estimates. The study presented here focuses on two snow-dominated basins in the Western United States for the years 2003–2016. Seasonal variations are observed in the vertical displacements derived from all three data sets, and the variation is coherent with the changes in hydrologic loading. Good consistency is observed between any two of the three data sets with gravimetric retrievals and hydrologic model estimates providing the highest level of agreement (i.e., all examined stations with correlation coefficient R > 0.70). Vertical displacements derived from gravimetric retrievals and ground-based GPS yielded R > 0.70 for more than 89% of the stations. In addition, it is found that both GPS-derived and space-based, gravimetry-derived vertical displacements clearly reflected the impact of climate variation (i.e., heavy precipitation during 2010–2011 winter followed by prolonged drought). Vertical displacements derived from the hydrologic model highlighted the relatively large precipitation convergence phase during late 2010 to early 2011 at some stations but not the prolonged drought that followed. The results indicate that ground-based GPS observations of vertical displacement have the capability to capture variations in TWS changes, which can be systematically merged in conjunction with Gravity Recovery and Climate Experiment (GRACE) into a land surface model to improve TWS estimates in a follow-up study.

Gaohong Yin↗
NASA NTRSDate not supplied

Towards a Third-Order Accurate, Second-Derivative-Free, Shock-Capturing Finite-Volume Method for Hypersonic Flows on Tetrahedral Grids

In this paper, we report progress in the development of a third-order accurate, second-derivative-free, shock-capturing finite-volume solver for three-dimensional unstructured grids. The method is economical in the sense that the computation and storage of second derivatives are not required for third-order accuracy. It is based on point-valued numerical solutions stored at cells, gradients computed and stored at nodes, and an efficient projected-derivative formula that eliminates the need for second derivatives in a quadratic solution interpolation. The projected-derivative formula is also used to eliminate second derivatives from a high-order flux quadrature formula, so that it can be implemented conveniently in the form of a numerical flux at a face center plus a correction term. Similarly, a high-order source quadrature formula can also be implemented in the form of a cell-center point evaluation plus a similar correction term. These features make it relatively straightforward to extend an existing second-order finite-volume code to third-order. This paper reports progress of implementing the method in the NASA VULCAN-CFD code and discusses the implementation of a high-order accurate limiter for shock capturing.

Weighted Least-Squares↗
NASA NTRSDate not supplied

The Use of the BSRN Data as A Benchmark for the POWER Hourly DHI and DNI and In Validating Derived Hourly GTI

The satellite-based CERES SYN1deg hourly data is the source data of the POWER GIS solar data that covers 2001 to near present. The SYN1deg(Ed4.1) hourly GHI agrees well with the BSRN data, but the hourly DHI and DirHI (Direct Horizontal Irradiance) are positively and negatively, respectively, biased with appreciable magnitudes. The hourly DNI, derived by dividing the DirHI by cos(SZA), or the cosine of the solar zenith angle, is therefore negatively biased. Based on the statistics of comparisons with the BSRN data, we performed bias corrections on the hourly DHI and DNI. The corrections were executed in the 3-D phase space of latitude, cos(SZA), and cloud fraction (CLFR). The isotropic model is then used to derive the hourly global tilted irradiance (GTI). For validation purpose, we applied the isotropic model to the BSRN data at the original 1-, 2-, 3- or 5-minute interval. The satellite-based hourly GTI shows good agreement with their BSRN counterpart. We also examined two monthly-mean-based methods that empirically derive monthly mean GTI and DNI from monthly mean GHI and from both monthly mean GHI and DHI. The monthly-mean-based results compare favorably with the hourly-mean-based results. The GEWEX SRB (V4-IP) provides POWER with daily mean GHI for the years before the CERES era, and the data were corrected using quantile mapping by referencing the CERES SYN1deg data. We used the Kolmogorov -Smirnov test (K-S test) and Cramer-von Mises test to examine how well the results agree with the BSRN data. We found that if we set the lower limit for the daily mean GHI to 30 W m-2, the data can pass the K-S test at 0.01 significance level and the Cramer-von Mises test at 0.001 significance level. If no lower limit is set on the daily means, the data fail both tests. The satellite-based CERES SYN1deg hourly data is the source data of the POWER GIS solar data that covers 2001 to near present. The SYN1deg(Ed4.1) hourly GHI agrees well with the BSRN data, but the hourly DHI and DirHI (Direct Horizontal Irradiance) are positively and negatively, respectively, biased with appreciable magnitudes. The hourly DNI, derived by dividing the DirHI by cos(SZA), or the cosine of the solar zenith angle, is therefore negatively biased. Based on the statistics of comparisons with the BSRN data, we performed bias corrections on the hourly DHI and DNI. The corrections were executed in the 3-D phase space of latitude, cos(SZA), and cloud fraction (CLFR). The isotropic model is then used to derive the hourly global tilted irradiance (GTI). For validation purpose, we applied the isotropic model to the BSRN data at the original 1-, 2-, 3- or 5-minute interval. The satellite-based hourly GTI shows good agreement with their BSRN counterpart. We also examined two monthly-mean-based methods that empirically derive monthly mean GTI and DNI from monthly mean GHI and from both monthly mean GHI and DHI. The monthly-mean-based results compare favorably with the hourly-mean-based results. The GEWEX SRB (V4-IP) provides POWER with daily mean GHI for the years before the CERES era, and the data were corrected using quantile mapping by referencing the CERES SYN1deg data. We used the Kolmogorov -Smirnov test (K-S test) and Cramer-von Mises test to examine how well the results agree with the BSRN data. We found that if we set the lower limit for the daily mean GHI to 30 W m-2, the data can pass the K-S test at 0.01 significance level and the Cramer-von Mises test at 0.001 significance level. If no lower limit is set on the daily means, the data fail both tests.

Taiping Zhang↗
DOE OSTI2026

JetGP: A derivative enhanced Gaussian process library

Derivative enhanced Gaussian Processes (DEGPs) can significantly improve surrogate model accuracy over standard Gaussian Process (GP) formulations by incorporating derivative information. However, standard implementations scale poorly with dimension, limiting their use in high dimensional engineering problems. JetGP is a Python framework that unifies existing derivative enhanced GP methodologies into a single library and extends them to support arbitrary order derivative information. The library implements four complementary formulations: standard derivative enhanced Gaussian Processes (DEGP), directional DEGP (DDEGP), generalized directional DEGP (GDDEGP), and weighted DEGP (WDEGP). By unifying these approaches in a consistent interface with robust numerical implementations, JetGP enables practitioners to balance predictive accuracy and computational efficiency for high dimensional optimization, uncertainty quantification, and sensitivity analysis in engineering design.

Derivative enhanced Gaussian process↗
DOE OSTI2026

ReMU: regional minimal updating for model-based derivative-free optimization

Derivative-free optimization (DFO) problems are optimization problems where derivative information is unavailable or extremely difficult to obtain. Model-based DFO solvers have been applied extensively in scientific computing. Powell's NEWUOA (2004) [Powell, The NEWUOA software for unconstrained optimization without derivatives, in Large-Scale Nonlinear Optimization, Nonconvex Optimization and its Applications Vol. 83, G. Di Pillo and M. Roma, eds., Springer, 2006, pp. 255–297] and Wild's POUNDerS (2014) [Wild, Solving derivative-free nonlinear least squares problems with POUNDERS, in Advances and Trends in Optimization with Engineering Applications, T. Terlaky, M.F. Anjos, and S. Ahmed, eds., SIAM, 2017, pp. 529–540] explore the numerical power of the minimal norm Hessian (MNH) model for DFO and contributed to the open discussion on building better models with fewer data to achieve faster numerical convergence. Another decade later, we propose the regional minimal updating (ReMU) models, and extend the previous models into a broader class, including the H 2 norm models [Xie and Yuan, Least H 2 norm updating of quadratic interpolation models for derivative-free trust-region algorithms, IMA J. Numer. Anal. 46 (2025), pp. 21–50]. This paper shows motivation behind ReMU models, computational details, theoretical and numerical results on particular extreme points and the barycentre of ReMU's weight coefficient region, and the associated KKT matrix error and distance. Novel metrics, such as the truncated Newton step error, are proposed to numerically understand the new models' properties. A new algorithmic strategy, based on iteratively adjusting the ReMU model type, is also proposed, and shows numerical advantages by combining and switching between the barycentric model and the classic least Frobenius norm model in an online fashion.

derivative-free trust-region methods↗