Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial least squares”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

A Wind-Driven Snow Redistribution Module for Alpine3d V3.3.0: Adaptations Designed for Downscaling Ice Sheet Surface Mass Balance

Ice sheets gain mass via snow accumulation at the ice sheet surface, which is the primary component of surface mass balance. On the Antarctic ice sheet, winds redistribute snow resulting in surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 km. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to MERRA-2, a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 mm w.e.yr−1 (13%) and increases variance explained by 24%. Despite these improvements, Alpine3D still underestimates observed accumulation variability, thus providing an opportunity for future model improvement.

Eric Keenan↗
NASA NTRS1985

Parameter identification in distributed systems

This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.

Baruh, H.↗
NASA NTRS1985

Well-posedness of one-way wave equations and absorbing boundary conditions

A one-way wave equation is a partial differential which, in some approximate sense, behaves like the wave equation in one direction but permits no propagation in the opposite one. The construction of such equations can be reduced to the approximation of the square root of (1-s sup 2) on -1, 1 by a rational function r(s) = p sub m (s)/q sub n(s). Those rational functions r for which the corresponding one-way wave equation is well-posed are characterized both as a partial differential equation and as an absorbing boundary condition for the wave equation. We find that if r(s) interpolates the square root of (1-s sup 2) at sufficiently many points in (-1,1), then well-posedness is assured. It follows that absorbing boundary conditions based on Pade approximation are well-posed if and only if (m, n) lies in one of two distinct diagonals in the Pade table, the two proposed by Engquist and Majda. Analogous results also hold for one-way wave equations derived from Chebyshev or least-squares approximation.

Trefethen, L. N.↗
DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
DOE OSTI2020

Physics-Informed Deep Neural Network Method for Limited Observability State Estimation

The precise knowledge regarding the state of the power grid is important in order to ensure optimal and reliable grid operation. Specifically, knowing the state of the distribution grid becomes increasingly important as more renewable energy sources are connected directly into the distribution network, increasing the fluctuations of the injected power. In this paper, we consider the case when the distribution grid becomes partially observable, due to for example cyber attacks, and the state estimation problem is under-determined. We present a new methodology that leverages a deep neural network (DNN) to estimate the grid state. The standard DNN training method is modified to explicitly incorporate the physical information of the grid topology and line/shunt admittance. We show that our method leads to a superior accuracy of the estimation when compared to the case when no physical information is provided. Finally, we compare the performance of our method to the standard state estimation approach, which is based on the weighted least squares with pseudo-measurements, and show that our method performs significantly better with respect to the estimation accuracy.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Simultaneous quaternion estimation (QUEST) and bias determination

Tests of a new method for the simultaneous estimation of spacecraft attitude and sensor biases, based on a quaternion estimation algorithm minimizing Wahba's loss function are presented. The new method is compared with a conventional batch least-squares differential correction algorithm. The estimates are based on data from strapdown gyros and star trackers, simulated with varying levels of Gaussian noise for both inertially-fixed and Earth-pointing reference attitudes. Both algorithms solve for the spacecraft attitude and the gyro drift rate biases. They converge to the same estimates at the same rate for inertially-fixed attitude, but the new algorithm converges more slowly than the differential correction for Earth-pointing attitude. The slower convergence of the new method for non-zero attitude rates is believed to be due to the use of an inadequate approximation for a partial derivative matrix. The new method requires about twice the computational effort of the differential correction. Improving the approximation for the partial derivative matrix in the new method is expected to improve its convergence at the cost of increased computational effort.

Markley, F. Landis↗
NASA NTRS1986

Archaean megacrystic plagioclase units and the tectonic setting of greenstones

Large (up to 20 cm), equidimensional, commonly euhedral, plagioclase megacrysts of highly calcic composition (An sub 80-90) occur commonly in all Archean cratons in one or more of three distinct associations: (1) as cumulate crystal segregations of anorthosite or as megacrysts in basaltic dikes, sills, and flows in greenstone belts that vary in metamorphic grade from greenschist to granulite. Throughout 100's of thousands of square kilometers of northwestern Ontario and Manitoba the plagioclase megacrysts occur in pillowed and massive flows, sills, dikes, large inclusions in dikes, and intrusive anorthositic complexes with areas of up to a few 100 sq km and spanning a period of at least 100 m.y. in the 2.7 to 2.8 b.y. time frame; (2) as basaltic dike swarms in stable cratonic areas forming parallel to subparallel patterns over hundreds of thousands of square kilometers intruding both granitic gneisses and supracrustal belts including greenstones. These swams include the Ameralik-Saglek system at 3.1 to 3.4 b.y., the Matachewan system at 2.5 to 2.6 b.y., and the Beartooth-Bighorn system at 2.2 to 2.3 b.y.; and, (3) as anorthositic complexes associated with marbles and quartzites (Sittampundi, India and Messina, South Africa) in granulite grade terrains. Initial attempts to correlate tectonic settings of similar modern crystbearing units with their Archean counterparts were only partially successful.

Phinney, W. C.↗
NASA NTRS1991

Basic research for the geodynamics program

The mathematical models of space very long base interferometry (VLBI) observables suitable for least squares covariance analysis were derived and estimatability problems inherent in the space VLBI system were explored, including a detailed rank defect analysis and sensitivity analysis. An important aim is to carry out a comparative analysis of the mathematical models of the ground-based VLBI and space VLBI observables in order to describe the background in detail. Computer programs were developed in order to check the relations, assess errors, and analyze sensitivity. In order to investigate the estimatability of different geodetic and geodynamic parameters from the space VLBI observables, the mathematical models for time delay and time delay rate observables of space VLBI were analytically derived along with the partial derivatives with respect to the parameters. Rank defect analysis was carried out both by analytical and numerical testing of linear dependencies between the columns of the normal matrix thus formed. Definite conclusions were formed about the rank defects in the system.

Source record↗
NASA NTRS1974

Surface fitting three-dimensional bodies

The geometry of general three-dimensional bodies is generated from coordinates of points in several cross sections. Since these points may not be smooth, they are divided into segments and general conic sections are curve fit in a least-squares sense to each segment of a cross section. The conic sections are then blended in the longitudinal direction by fitting parametric cubic-spline curves through coordinate points which define the conic sections in the cross-sectional planes. Both the cross-sectional and longitudinal curves may be modified by specifying particular segments as straight lines and slopes at selected points. Slopes may be continuous or discontinuous and finite or infinite. After a satisfactory surface fit has been obtained, cards may be punched with the data necessary to form a geometry subroutine package for use in other computer programs. At any position on the body, coordinates, slopes and second partial derivatives are calculated. The method is applied to a blunted 70 deg delta wing, and it was found to generate the geometry very well.

Dejarnette, F. R.↗
NASA NTRS1993

Large-scale computation of incompressible viscous flow by least-squares finite element method

The least-squares finite element method (LSFEM) based on the velocity-pressure-vorticity formulation is applied to large-scale/three-dimensional steady incompressible Navier-Stokes problems. This method can accommodate equal-order interpolations and results in symmetric, positive definite algebraic system which can be solved effectively by simple iterative methods. The first-order velocity-Bernoulli function-vorticity formulation for incompressible viscous flows is also tested. For three-dimensional cases, an additional compatibility equation, i.e., the divergence of the vorticity vector should be zero, is included to make the first-order system elliptic. The simple substitution of the Newton's method is employed to linearize the partial differential equations, the LSFEM is used to obtain discretized equations, and the system of algebraic equations is solved using the Jacobi preconditioned conjugate gradient method which avoids formation of either element or global matrices (matrix-free) to achieve high efficiency. To show the validity of this scheme for large-scale computation, we give numerical results for 2D driven cavity problem at Re = 10000 with 408 x 400 bilinear elements. The flow in a 3D cavity is calculated at Re = 100, 400, and 1,000 with 50 x 50 x 50 trilinear elements. The Taylor-Goertler-like vortices are observed for Re = 1,000.

Jiang, Bo-Nan↗
DOE OSTI2020

Thermal Dehydration of Aluminum (Oxy)hydroxides on Fuel Cladding Material - 20200

The aluminum cladding of research-reactor fuel undergoes general corrosion with resulting formation of adherent aluminum (oxy)hydroxide films during in-reactor and post-discharge exposure to water under various conditions and temperatures. These (oxy)hydroxides contain chemically-bound water that poses challenges for extended dry storage due to the risk of thermal or radiolytic decomposition releasing free water and/or hydrogen and oxygen gases. This study comprised laboratory experiments of thermal drying behavior of aluminum (oxy)hydroxides to identify approaches for reducing bound water on fuel cladding prior to sealed dry storage. Effective drying strategies for adherent (oxy)hydroxides on fuel cladding will improve the safety of dry storage by mitigating potential avenues for additional corrosion and/or generation of flammable gases inside the storage canister. This work is part of a broader investigation to address knowledge gaps and technical data needs for dry storage of aluminum-clad spent nuclear fuel (ASNF), which included investigation of (oxy)hydroxide formation on aluminum alloy substrates immersed in water, characterization of service-grown films from ASNF, and radiolytic yield of hydrogen from (oxy)hydroxide powders and films. The current work comprises experimental thermal dehydration of aluminum trihydroxides (bayerite and gibbsite) characteristic of low-temperature (<80 deg. C) corrosion in water and aluminum oxyhydroxide (boehmite) characteristic of high temperature (>80 deg. C) corrosion in water. First, commercially produced (oxy)hydroxide powders (boehmite and gibbsite) were tested via thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC). The effort aimed to identify or confirm key temperature ranges/thresholds for the thermal decomposition reactions, as well as impacts of ramp rate and hold times, for drying of isolated, high-surface-area (oxy)hydroxides. The information gleaned from powder tests was used to guide subsequent drying tests on adherent (oxy)hydroxide films grown on aluminum alloy substrates. TGA was used to analyze small samples of the aluminum coupons with adherent (oxy)hydroxide films. Specimens were characterized both pre- and post-drying using X-ray diffraction (XRD) to determine the film composition and scanning electron microscopy (SEM) to determine its morphology. Drying of adherent (oxy)hydroxide films was anticipated to pose additional complications relative to drying of (oxy)hydroxide powders. Anticipated challenges include achieving sufficient drying of the (oxy)hydroxide at low-enough temperature to avoid melting or other undesirable phase changes in the aluminum substrate and overcoming transport limitations of water through the thickness of a dense film or through a tortuous pore structure to reach the outer surface of the film. Inn addition, dehydration of (oxy)hydroxides significantly changes the film density, which may lead to substantial alterations in the morphology, including potential cracking and spalling of the film. TGA/DSC tests of gibbsite powders resulted in successful conversion to boehmite or alumina, depending on the maximum temperature reached. At low ramp rates (≤5 deg. C/min), the conversion to boehmite occurred around 300 deg. C (210-340 deg. C). XRD confirmed that boehmite was the only phase detected after tests reaching 450 deg. C. The transition to alumina occurred around 510 deg. C (470-550 deg. C), with XRD detecting only alumina after tests reaching at least 600 deg. C. Boehmite powders dehydrated to alumina at about 400 deg. C (330-460 deg. C) for coarse (77-μm particle diameter) powder and about 490 deg. C (420-520 deg. C) for fine (0.7-μm particle diameter) powder. For all powders tested, the maximum percent mass loss after drying to high temperature slightly exceeded the theoretical mass loss for complete dehydration of the stoichiometric (oxy)hydroxide to alumina, which is likely attributable to physisorbed water in the powder. Drying of a thick (∼8.6 μm), predominantly bayerite adherent film displayed a dramatic change in film morphology, with the initially continuous trihydroxide film cracking into sections on the order of 50 μm square and partially delaminating. The outermost layer of the film completely faked of in some regions. The layer exposed under the spalled oxide also displayed prominent cracking, with spacing on the order of 10 μm, but this surface layer appeared to remain completely adhered to the aluminum substrate. TGA showed mass losses per unit surface area up to ∼1.0 mg/cm{sup 2}. XRD characterization of the remaining oxide was inconclusive, with no crystalline phases detected. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2003

A New Fitting Procedure for the Blended He I 1083 nm Multiplet

The He I 1083 nm multiplet is a powerful tool for observing the outer solar atmosphere but is difficult to analyze because the lines are weak, highly variable, and spectrally blended, both internally and with other neighboring solar and telluric lines. After separation from nearby spectral features, two components of the He I multiplet are resolved. Fitting these lines with two unconstrained Gaussian profiles always gives a ratio of major to minor component of less than half the value which would be expected for optically thin lines. One possibility for explaining the discrepancy between the weakness of the lines and the ratio of the spectral components is to assume that the line formation regions are concentrated in laterally unresolved, optically thick structures with small filling factor. However, we present here a least-squares fitting technique using cubic splines with fixed breakpoints with the constraint that the blend is the sum of three identically shaped profiles shifted in wavelength according to the atomic structure of the blend and weighted by the corresponding statistics weights, in agreement with optically thin line formation. The basis functions for the fitting procedure have no built-in spectral symmetry or shape. The resulting underlying profiles tend to be asymmetric with excess absorption to the blue, consistent with formation by "hot" and "cool" spatial elements within the observed volume, with the hotter regions having differential motion toward the observer. The results agree well with NASA/XSO Spectromagnetograph observations in quiet sun and coronal holes. Partial funding of this research was provided through the NASA Sun-Earth Connection SR&T program.

Jones, H. P.↗
DOE OSTI2020

Lift & Learn: Physics-informed machine learning for large-scale nonlinear dynamical systems

In this work, we present Lift & Learn, a physics-informed method for learning low-dimensional models for large-scale dynamical systems. The method exploits knowledge of a system’s governing equations to identify a coordinate transformation in which the system dynamics have quadratic structure. This transformation is called a lifting map because it often adds auxiliary variables to the system state. The lifting map is applied to data obtained by evaluating a model for the original nonlinear system. This lifted data is projected onto its leading principal components, and low-dimensional linear and quadratic matrix operators are fit to the lifted reduced data using a least-squares operator inference procedure. Analysis of our method shows that the Lift & Learn models are able to capture the system physics in the lifted coordinates at least as accurately as traditional intrusive model reduction approaches. This preservation of system physics makes the Lift & Learn models robust to changes in inputs. Numerical experiments on the FitzHugh–Nagumo neuron activation model and the compressible Euler equations demonstrate the generalizability of our model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Dynamic calibration of differential equations using machine learning, with application to turbulence models

We present a methodology for calibration of parametric ordinary and partial differential equation models, using off-the-shelf software for back-propagation in Neural Networks (NN). As a prototypical example, we consider calibration of a Reynolds-averaged Navier-Stokes (RANS) turbulence closure model, against ground truth data from direct numerical simulations (DNS) of two different turbulent flows. Numerical time integration is represented as a custom NN, where only the RANS model parameters are trainable. A loss function is defined to quantify the mismatch between the NN prediction and the ground truth over a predefined, finite time integration window. This loss function is then minimized using a gradient descent method utilizing the back-propagation algorithm. Furthermore, this dynamic approach to training is to be contrasted with a static approach, wherein a least square regression estimate for parameters is obtained in the limit of an infinitesimal time integration window. In a first test of static and dynamic approaches against ground truth data generated by the model, the former proves to be significantly faster and more accurate than the latter at recovering the parameters. When both calibration approaches are tested against DNS data, for which it is known that the model cannot achieve a perfect fit, the static approach yields a good prediction only for short times, while the dynamic approach results in physical and stable predictions over the entire integration window. After optimization of the dynamic approach for time step, spatial resolution, stability, and physics-based constraints, we obtain a 50% improvement of outcomes over those obtained from the existing, manually calibrated set of parameters, demonstrating the merits of this systematic and automated procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Extended Galerkin Neural Network Approximation of Singular Variational Problems with Error Control

We present extended Galerkin neural networks, a variational framework for approximating general boundary value problems (BVPs) with error control. The main contributions of this work are (1) a rigorous theory guiding the construction of new weighted least squares variational formulations suitable for use in neural network approximation of general BVPs, and (2) an “extended” feedforward network architecture which incorporates and is even capable of learning singular solution structures, thus greatly improving approximability of singular solutions. Furthermore, numerical results are presented for several problems, including steady Stokes flow around reentrant corners and in convex corners with Moffatt eddies in order to demonstrate efficacy of the method.

a posteriori error estimate↗
NASA NTRS2016

Carbon Monoxide Accumulation in the Extravehicular Mobility Unit

Introduction: Life support technology in large closed systems like submarines and space stations catalyzes carbon monoxide (CO) to carbon dioxide, which is easily removed. However, in a small system like the Extravehicular Mobility Unit (EMU), spacesuit, CO from exogenous (contaminated oxygen (O (sub 2) supply) and endogenous (human metabolism) sources will accumulate in the free suit volume. The free volume becomes a sink for CO that is rebreathed by the astronaut. The accumulation through time depends on many variables: the amount absorbed by the astronaut, the amount produced by the astronaut (between 0.28 and 0.34 moles per hour per kilogram)[1], the amount that enters the suit from contaminated O (sub 2), the amount removed through suit leak, the free volume of the suit, and the O (sub 2) partial pressure[2], just to list a few. Contamination of the EMU O (sub 2) supply with no greater than 1 part per million CO was the motivation for empirical measurements from CO pulse oximetry (SpCO) as well as mathematical modeling of the EMU as a rebreather for CO. Methods: We developed a first-order differential mixing equation as well as an iterative method to compute CO accumulation in the EMU. Pre-post measurements of SpCO (Rad-57, Masimo Corporation) from EMU ground training and on-orbit extravehicular activities (EVAs) were collected. Results: Initial modeling without consideration of the astronaut as a sink but only the source of CO showed that after 8 hours breathing 100 percent O (sub 2) with a 10 milliliter per minute (760 millimeters Hg at 21 degrees Centigrade standard) suit leak, an endogenous production rate of 0.23 moles per hour per kilogram for a 70 kilogram person with 42 liters (1.5 cubic feet) free suit volume resulted in a peak CO partial pressure (pCO) of 0.047 millimeters Hg at 4.3 pounds per square inch absolute (222 millimeters Hg). Preliminary results based on a 2008 model[3] with consideration of the astronaut as a sink and source of CO suggests that most of the rebreathed CO stays bound to hemoglobin and myoglobin and; therefore, pCO only increased to 0.002 millimeters Hg in the EMU. Hemoglobin saturation after 8 hours was an insignificant 0.4 percent compared to about 4 percent for cigarette smokers in the general population. This preliminary modeling result supplements 11 pre-post index finger SpCO measurements from EMU ground training (mean 0.5 percent versus 1 percent, probability equal to 0.41 from paired t-test) and 10 on-orbit pre-post EVAs (mean 1.5 percent versus 1.1 percent, probability equal to 0.17 from paired t-test) that showed no consistent increase, at least no increase outside the accuracy of the oximeter (1 percent display resolution with plus or minus 3 percent Standard Deviation). Discussion: Simulations continue, but a preliminary conclusion is that rebreathing endogenous CO accumulating in the EMU is not a serious medical issue during EVA. The absence of CO poisoning signs or symptoms following hundreds of EVAs is also good empirical evidence that corroborates the limited SpCO measurements and preliminary modeling results.

Conkin, J.↗
NASA NTRS2013

The NuSTAR Extragalactic Survey: A First Sensitive Look at the High-Energy Cosmic X-Ray Background Population

We report on the first 10 identifications of sources serendipitously detected by the Nuclear Spectroscopic Telescope Array (NuSTAR) to provide the first sensitive census of the cosmic X-ray background source population at approximately greater than 10 keV. We find that these NuSTAR-detected sources are approximately 100 times fainter than those previously detected at approximately greater than 10 keV and have a broad range in redshift and luminosity (z = 0.020-2.923 and L(sub 10-40 keV) approximately equals 4 × 10(exp 41) - 5 × 10(exp 45) erg per second; the median redshift and luminosity are z approximately equal to 0.7 and L(sub 10-40 keV) approximately equal to 3 × 10(exp 44) erg per second, respectively. We characterize these sources on the basis of broad-band approximately equal to 0.5 - 32 keV spectroscopy, optical spectroscopy, and broad-band ultraviolet-to-mid-infrared spectral energy distribution analyses. We find that the dominant source population is quasars with L(sub 10-40 keV) greater than 10(exp 44) erg per second, of which approximately 50% are obscured with N(sub H) approximately greater than 10(exp 22) per square centimeters. However, none of the 10 NuSTAR sources are Compton thick (N(sub H) approximately greater than 10(exp 24) per square centimeters) and we place a 90% confidence upper limit on the fraction of Compton-thick quasars (L(sub 10-40 keV) greater than 10(exp 44) erg per second) selected at approximately greater than 10 keV of approximately less than 33% over the redshift range z = 0.5 - 1.1. We jointly fitted the rest-frame approximately equal to 10-40 keV data for all of the non-beamed sources with L(sub 10-40 keV) greater than 10(exp 43) erg per second to constrain the average strength of reflection; we find R less than 1.4 for gamma = 1.8, broadly consistent with that found for local active galactic nuclei (AGNs) observed at approximately greater than 10 keV. We also constrain the host-galaxy masses and find a median stellar mass of approximately 10(exp 11) solar mass, a factor approximately 5 times higher than the median stellar mass of nearby high-energy selected AGNs, which may be at least partially driven by the order of magnitude higher X-ray luminosities of the NuSTAR sources. Within the low source-statistic limitations of our study, our results suggest that the overall properties of the NuSTAR sources are broadly similar to those of nearby high-energy selected AGNs but scaled up in luminosity and mass.

Spectroscopic↗