SEARCH · Engineering Papers
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.
In this paper the problem of synthesizing dual reflector antennas for both amplitude and phase control of the final aperture distribution is discussed. An approximate procedure for the offset synthesis problem is presented and applications of the procedure to the shaping of beam waveguides and reflectors for high-gain antenna systems are illustrated.
Nuclear Computational Low-Energy Initiative (NUCLEI) collaboration uses Density Functional Theory (DFT) simulations to predict the structure and binding energies of nuclei over a wide range of proton (Z) and neutron (N) numbers. The DFT simulations utilize a particular parameterization of a Skyrme energy density functional called UNEDF1 which depends on 12 free parameters that must be fit to data (M Kortelainen et al 2014). Fitting involves comparing (e.g.) predicted binding energies of nuclei to experimentally measured values. We use only binding energies as observables, but DFT with UNEDF1 will predict structure (shape) observables as well. In this work, assessing the capability of approximate GP emulators to balance emulator accuracy with computational speed to facilitate improved UNEDF1 calibration. Sparse GPs are straightforward to train and accurate. Calibration is not straightforward with MCMC (using MH or HMC/NUTS). We produced reusable software for continuing and building on this work as well as accessing and using Darwin cluster compute resources
NASA’s Orbital Debris Program Office (ODPO) and Hypervelocity Impact Technology (HVIT) team have coordinated to better understand the risks to upper stages and spacecraft from non-spherical orbital debris. It is well understood that fragmentation (collision or explosion) events in orbit produce fragments of various materials, sizes, and shapes. To further characterize these parameters, the ODPO is developing the next-generation Orbital Debris Engineering Model (ORDEM) version 4.0 to include orbital debris shape distributions. Ground-based assets, such as radar and optical sensors, can provide size estimates and some insight into material based on radar return or optical filter photometry/spectroscopy, respectively. Characterizing an object’s shape requires more laboratory analyses to infer how shape affects these measurements. More importantly, in addition to size and material/density, the shape of fragments in orbit will alter the ballistic limit equations used in orbital debris risk assessments with NASA’s Bumper Code. The ODPO plans to release ORDEM 4.0 in the coming years. Performing ground-based laboratory impact tests on high-fidelity spacecraft mockups provides the means to directly measure size, mass, material/density, and shape of fragments, all key parameters needed to characterize real-world break up events. The DebriSat test, the results of which are provided, showcases the details of this type of experiment. The goal of this collaborative research between the ODPO and the HVIT team is to include a shape parameter in the environmental and breakup models used to assess risk for various space structures. This paper examines ground-based laboratory impact tests and the associated fragment shape categories. Provided these defined shapes, the approach is simplified by assuming a right circular cylinder (RCC) approximation with varying length-to-diameter ratios. Highlights of impact tests conducted by the HVIT team using non-spherical projectiles based on the RCC approximation are presented. Hydrocode simulations have also been performed to expand on the complexity of variations with non-spherical projectiles. Lastly, ray-tracing simulations of various RCCs of known material are provided to support the ongoing research on optical reflectance distributions with known shapes and to highlight how this may modify the current optical size estimation model. The status and plan forward are outlined for NASA's orbital debris shape effect investigation using a multidisciplinary approach by the ODPO and the HVIT team.
NASA’s Orbital Debris Program Office (ODPO) and Hypervelocity Impact Technology (HVIT) team have coordinated to better understand the risks to upper stages and spacecraft from non-spherical orbital debris. It is well understood that fragmentation (collision or explosion) events in orbit produce fragments of various materials, sizes, and shapes. To further characterize these parameters, the ODPO is developing the next-generation Orbital Debris Engineering Model (ORDEM) version 4.0 to include orbital debris shape distributions. Ground-based assets, such as radar and optical sensors, can provide size estimates and some insight into material based on radar return or optical filter photometry/spectroscopy, respectively. Characterizing an object’s shape requires more laboratory analyses to infer how shape affects these measurements. More importantly, in addition to size and material/density, the shape of fragments in orbit will alter the ballistic limit equations used in orbital debris risk assessments with NASA’s Bumper Code. The ODPO plans to release ORDEM 4.0 in the coming years. Performing ground-based laboratory impact tests on high-fidelity spacecraft mockups provides the means to directly measure size, mass, material/density, and shape of fragments, all key parameters needed to characterize real-world break up events. The DebriSat test, the results of which are provided, showcases the details of this type of experiment. The goal of this collaborative research between the ODPO and the HVIT team is to include a shape parameter in the environmental and breakup models used to assess risk for various space structures. This paper examines ground-based laboratory impact tests and the associated fragment shape categories. Provided these defined shapes, the approach is simplified by assuming a right circular cylinder (RCC) approximation with varying length-to-diameter ratios. Highlights of impact tests conducted by the HVIT team using non-spherical projectiles based on the RCC approximation are presented. Hydrocode simulations have also been performed to expand on the complexity of variations with non-spherical projectiles. Lastly, ray-tracing simulations of various RCCs of known material are provided to support the ongoing research on optical reflectance distributions with known shapes and to highlight how this may modify the current optical size estimation model. The status and plan forward are outlined for NASA's orbital debris shape effect investigation using a multidisciplinary approach by the ODPO and the HVIT team.
Physical band models that use simplifying assumptions about the structure of the CO2 15-micron band are described and compared with the accurate line-by-line model of Fels and Schwarzkopf (1981). The methods have application to the models of the atmospheres of Venus, earth, and Mars. The narrow-band methods include a general random model and techniques employing the Malkmus (1967) line-strength distribution, and the accuracy of these models at low pressure is found to be sensitive to the line-strength distribution and the line shape. When combined with a simple Voigt parameterization, the more efficient approximate methods represented by an exponential wideband model are found to be almost as accurate as the expensive narrow-band models. Empirical logarithmic wideband models used in climate-modeling studies are found to be very efficient, but of limited range of validity.
A method is presented to calculate analytically the sensitivity derivatives of wing static aeroelastic characteristics with respect to wing shape parameters. The wing aerodynamic response under fixed total load is predicted with Weissinger's L-method; its structural response is obtained with Giles' equivalent plate method. The characteristics of interest include the spanwise distribution of lift, trim angle of attack, rolling and pitching moments, wind induced drag, as well as the divergence dynamic pressure. The shape parameters considered are the wing area, aspect ratio, taper ratio, sweep angle, and tip twist angle. Results of sensitivity studies indicate that: (1) approximations based on analytical sensitivity derivatives can be used over wide ranges of variations of the shape parameters considered, and (2) the analytical calculation of sensitivity derivatives is significantly less expensive than the conventional finite-difference alternative.
A method is presented to calculate analytically the sensitivity derivatives of wing static aeroelastic characteristics with respect to wing shape parameters. The wing aerodynamic response under fixed total load is predicted using Weissinger's L-method; its structural response is obtained with Giles' equivalent plate method. The characteristics of interest in this study include the spanwise distribution of lift, trim angle of attack, rolling and pitching moments, wing induced drag as well as the divergence dynamic pressure. The shape parameters considered are the wing area, aspect ratio, taper ratio, sweep angle and tip twist angle. Results of sensitivity studies indicate that (1) approximations based on analytical sensitivity derivatives can be used for wide ranges of variations of the shape parameters considered and (2) the analytical calculation of sensitivity derivatives is significantly less expensive than the conventional finite difference alternative.
This article explains results developed from the following research: 'The Stability of Plates and Shells beyond the Elastic Limit.' A significant improvement is found in the derivation of the relations between the stress factors and the strains resulting from the instability of plates and shells. In a strict analysis, the problem reduces to the solution of two simultaneous nonlinear partial differential equations of the fourth order in the deflection and stress function, and in the approximate analysis to a single linear equation of the Bryan type. Solutions are given for the special cases of a rectangular plate buckling into a cylindrical form, and of an arbitrarily shaped plate under uniform compression. These solutions indicate that the accuracy obtained by the approximate method is satisfactory.
We have analyzed a large set of Rossi X-ray Timing Explorer/Proportional Counter Array (RXTE/PCA) scanning and slewing observations performed between April 1996 and March 1999. We obtained the 3-20 keV spectrum of the cosmic X-ray background (CXB) by subtracting Earth-occulted observations from observations of the X-ray sky at high galactic latitude and far away from sources. The sky coverage is approximately 22.6 x 10(exp 3) square degrees. The PCA spectrum of CXB in 3-20 keV energy band is adequately approximated by a single power law with photon index GAMMA approximately 1.4 and normalization at 1 keV approximately 9.5 phot/s/square centimeter/keV/sr. Instrumental background uncertainty precludes accurate RXTE/PCA measurements of the spectrum of cosmic X-ray background at energies above 15 keV and therefore we cannot detect the high energy cutoff observed by the High Energy Astronomical Observatory (HEAO)-1 A2 experiment. Deep observations of the 6 high latitude points used to model the PCA background provide a coarse measure of the spatial variation of the CXB. The CXB variations are consistent with a fixed spectral shape and variable normalization characterized by a fractional rms amplitude of approximately 7% on angular scales of approximately 1 square deg.
Air-to-refrigerant heat exchangers (HXs) are fundamental components in HVAC&R systems, and considerable research has been dedicated designing continually smaller, lighter, and more efficient HX designs. In recent years, researchers have leveraged advancements in Computational Fluid Dynamics (CFD), Finite Element Analysis (FEA), and optimization algorithms to consider primary tube shape and topology optimization to design highly compact, high performance HXs for a multitude of applications. In this research, we present a computationally efficient, comprehensive, multi-scale, and multi-physics analysis and optimization method for air-to-refrigerant HXs featuring automated CFD and FEA simulations and approximation-assisted optimization. This methodology was utilized to design HXs with shape-optimized, non-round tubes which outperform current state-of-the-art tube-fin HXs without compromising structural integrity. The optimal HXs were shown to deliver similar thermal performance to the baseline HXs while also achieving more than 20% reductions in airside pressure drop and core envelope volume and more than 30% reduction in internal volume. Comprehensive experimental validation of the optimization methodology was conducted through the testing of two prototypes in a standardized wind tunnel facility under multiple operating conditions. For prototype #1 under dry evaporator conditions, the predicted heat load agreed within ± 10% of measured values and the predicted airside pressure drop agreed within ± 30%, while for dehumidifying conditions, the predicted sensible and latent heat loads agreed within ± 10% and ± 20% of the measured values, respectively. For prototype #2, the predicted condenser heat load agreed within ± 3.0% of measured values, and the predicted airside pressure drop agreed within ± 27%. The acceptable agreement between simulation and experimental results for the present application highlights the flexibility of the novel optimization methodology to design next generation HXs with improved performance and reduced volume, weight, and environmental impact.
Radio-occultation and thermal-infrared measurements are complementary investigations for sounding planetary atmospheres. The vertical resolution afforded by radio occultations is typically approximately 1 km or better, whereas that from infrared sounding is often comparable to a scale height. On the other hand, an instrument like CIRS can easily generate global maps of temperature and composition, whereas occultation soundings are usually distributed more sparsely. The starting point for radio-occultation inversions is determining the residual Doppler-shifted frequency, that is the shift in frequency from what it would be in the absence of the atmosphere. Hence the positions and relative velocities of the spacecraft, target atmosphere, and DSN receiving station must be known to high accuracy. It is not surprising that the inversions can be susceptible to sources of systematic errors. Stratospheric temperature profiles on Titan retrieved from Cassini radio occultations were found to be very susceptible to errors in the reconstructed spacecraft velocities (approximately equal to 1 mm/s). Here the ability to adjust the spacecraft ephemeris so that the profiles matched those retrieved from CIRS limb sounding proved to be critical in mitigating this error. A similar procedure can be used for Saturn, although the sensitivity of its retrieved profiles to this type of error seems to be smaller. One issue that has appeared in inverting the Cassini occultations by Saturn is the uncertainty in its equatorial bulge, that is, the shape in its iso-density surfaces at low latitudes. Typically one approximates that surface as a geopotential surface by assuming a barotropic atmosphere. However, the recent controversy in the equatorial winds, i.e., whether they changed between the Voyager (1981) era and later (after 1996) epochs of Cassini and some Hubble observations, has made it difficult to know the exact shape of the surface, and it leads to uncertainties in the retrieved temperature profiles of one to a few kelvins. This propagates into errors in the retrieved helium abundance, which makes use of thermal-infrared spectra and synthetic spectra computed with retrieved radio-occultation temperature profiles. The highest abundances are retrieved with the faster Voyager-era winds, but even these abundances are somewhat smaller than those retrieved from the thermal-infrared data alone (albeit with larger formal errors). The helium abundance determination is most sensitive to temperatures in the upper troposphere. Further progress may include matching the radio-occultation profiles with those from CIRS limb sounding in the upper stratosphere.
A molded magnetic article and fabrication method are provided. Particles of ferromagnetic material embedded in a polymer binder are molded under heat and pressure into a geometric shape. Each particle is an oblate spheroid having a radius-to-thickness aspect ratio approximately in the range of 15-30. Each oblate spheroid has flattened poles that are substantially in perpendicular alignment to a direction of the molding pressure throughout the geometric shape.
A molded magnetic article and fabrication method are provided. Particles of ferromagnetic material embedded in a polymer binder are molded under heat and pressure into a geometric shape. Each particle is an oblate spheroid having a radius-to-thickness aspect ratio approximately in the range of 15-30. Each oblate spheroid has flattened poles that are substantially in perpendicular alignment to a direction of the molding pressure throughout the geometric shape.
Dark matter scattering off a nucleus has a small probability of inducing an observable ionization through the inelastic excitation of an electron, called the Migdal effect. We use an effective field theory to extend the computation of the Migdal effect in semiconductors to regions of small momentum transfer to the nucleus, where the final state of the nucleus is no longer well described by a plane wave. Our analytical result can be fully quantified by the measurable dynamic structure factor of the semiconductor, which accounts for the vibrational degrees of freedom (phonons) in a crystal. We show that, due to the sum rules obeyed by the structure factor, the inclusive Migdal rate and the shape of the electron recoil spectrum is well captured by approximating the nuclei in the crystal as free ions; however, the exclusive differential rate with respect to energy depositions to the crystal depends on the phonon dynamics encoded in the dynamic structure function of the specific material. Our results now allow the Migdal effect in semiconductors to be evaluated even for the lightest dark matter candidates (m χ ≳ 1 MeV) that can kinematically excite electrons.
Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.
The availability of actinide standard materials for use in nuclear safeguard applications is critical, as is thorough characterization thereof. Although accurate trace element compositions and isotopic considerations are paramount for deployment of reference standards, structural characterization is also essential towards accurately describing the chemical form and potential matrix effects in candidate materials. Here, to this end, samples of NpO 2 were synthesized via a direct denitration (DD) method and probed with powder X-ray diffraction (PXRD), Raman spectroscopy, and scanning electron microscopy (SEM) for structural and morphological characterization and comparison with NpO 2 materials produced via modified direct denitration (MDD). PXRD confirmed the bulk identity of NpO 2 , and no additional phases were identified using this method. Analysis of Raman data collected using a 532 nm excitation wavelength indicates that samples are mostly phase pure; however, some variability in spectral features is observed. Analysis of additional spectroscopic data collected with a 785 nm excitation wavelength revealed variability in the relative intensity of spectral features. Raman spectroscopy indicates that the sample is primarily NpO 2 ; however, additional signals indicate possible structural disorder, oxidized species, or potential contributions from other Np phases. To further investigate the possibility of additional phase contributions within the sample of NpO 2 , Raman spectroscopic mapping was employed to examine the homogeneity of the sample produced via DD. From this analysis, we determined that despite variability in the intensity of Raman-active vibrational modes, consistent spectra are obtained throughout the area of the sample investigated. SEM images show aggregates with variable sizes and shapes, with rounded, primary particles possessing an average diameter of approximately 100 nm. Comparison of the results of these multimodal analyses to the literature indicates that the crystal chemical, spectroscopic, and microstructural properties of NpO 2 vary based on synthesis method, even if X-ray diffraction data indicate that the bulk phase is NpO 2 .
Pixel tracking detectors at upcoming collider experiments will see unprecedented charged-particle densities. Real-time data reduction on the detector will enable higher granularity and faster readout, possibly enabling the use of the pixel detector in high-rate online event selection, such as the ATLAS or CMS first-level trigger systems. This data reduction can be accomplished with a neural network (NN) in the readout chip bonded with the sensor that recognizes and rejects tracks with low transverse momentum (p T ) based on the geometrical shape of the charge deposition (“cluster”). To design viable detectors for deployment, the dependence of the NN as a function of the sensor geometry, external magnetic field, irradiation, and noise must be understood. In this paper, we present first studies of the efficiency and data reduction for planar pixel sensors exploring these parameters. For the CMS HL-LHC sensor geometry, we obtain a signal efficiency of (91.9 ± 0.7)% and a data reduction of (29.7 ± 1.0)%. A smaller sensor pitch in the bending direction improves the p T discrimination, but a larger pitch can be partially compensated with detector thickness. Any accumulated radiation damage also changes the cluster shape, reducing the signal efficiency compared to the baseline by approximately 30–60% in absolute terms, but nearly all of the performance can be recovered through retraining of the network and updating the weights. Finally, the impact of noise was investigated, and retraining the network on noise-injected datasets was found to maintain performance within 6% of the baseline network trained and evaluated on noiseless data. •ASIC-compatible track-momentum classifier is robust in realistic detector conditions.•About 90% signal efficiency and 30% data reduction per layer for CMS HL-LHC geometry.•Single-layer signal efficiency increases for smaller pixel pitch or thicker sensors.•Performance with noise or after radiation damage mostly recovered by retraining.
Abstract Addressing human anatomical and physiological variability is a crucial component of human health risk assessment of chemicals. Experts have recommended probabilistic chemical risk assessment paradigms in which distributional adjustment factors are used to account for various sources of uncertainty and variability, including variability in the pharmacokinetic behavior of a given substance in different humans. In practice, convenient assumptions about the distribution forms of adjustment factors and human equivalent doses (HEDs) are often used. Parameters such as tissue volumes and blood flows are likewise often assumed to be lognormally or normally distributed without evaluating empirical data for consistency with these forms. In this work, we performed dosimetric extrapolations using physiologically based pharmacokinetic (PBPK) models for dichloromethane (DCM) and chloroform that incorporate uncertainty and variability to determine if the HEDs associated with such extrapolations are approximately lognormal and how they depend on the underlying distribution shapes chosen to represent model parameters. We accounted for uncertainty and variability in PBPK model parameters by randomly drawing their values from a variety of distribution types. We then performed reverse dosimetry to calculate HEDs based on animal points of departure for each set of sampled parameters. Corresponding samples of HEDs were tested to determine the impact of input parameter distributions on their central tendencies, extreme percentiles, and degree of conformance to lognormality. This work demonstrates that the measurable attributes of human variability should be considered more carefully and that generalized assumptions about parameter distribution shapes may lead to inaccurate estimates of extreme percentiles of HEDs.