Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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 361 records · Page 20

Three-dimensional Navier-Stokes analysis of turbine passage heat transfer

The three-dimensional Reynolds-averaged Navier-Stokes equations are numerically solved to obtain the pressure distribution and heat transfer rates on the endwalls and the blades of two linear turbine cascades. The TRAF3D code which has recently been developed in a joint project between researchers from the University of Florence and NASA Lewis Research Center is used. The effect of turbulence is taken into account by using the eddy viscosity hypothesis and the two-layer mixing length model of Baldwin and Lomax. Predictions of surface heat transfer are made for Langston's cascade and compared with the data obtained for that cascade by Graziani. The comparison was found to be favorable. The code is also applied to a linear transonic rotor cascade to predict the pressure distributions and heat transfer rates.

Ameri, Ali A.↗

Augmentation of Rocket Propulsion: Physical Limits

Rocket propulsion is not ideal when the propellant is not ejected at a unique velocity in an inertial frame. An ideal velocity distribution requires that the exhaust velocity vary linearly with the velocity of the vehicle in an inertial frame. It also requires that the velocity distribution variance as a thermodynamic quantity be minimized. A rocket vehicle with an inert propellant is not optimal, because it does not take advantage of the propellant mass for energy storage. Nor is it logical to provide another energy storage device in order to realize variable exhaust velocity, because it would have to be partly unfilled at the beginning of the mission. Performance is enhanced by pushing on the surrounding because it increases the reaction mass and decreases the reaction jet velocity. This decreases the fraction of the energy taken away by the propellant and increases the share taken by the payload. For an optimal model with the propellant used as fuel, the augmentation realized by pushing on air is greatest for vehicles with a low initial/final mass ratio. For a typical vehicle in the Earth's atmosphere, the augmentation is seen mainly at altitudes below about 80 km. When drag is taken into account, there is a well-defined optimum size for the air intake. Pushing on air has the potential to increase the performance of rockets which pass through the atmosphere. This is apart from benefits derived from "air breathing", or using the oxygen in the atmosphere to reduce the mass of an on-board oxidizer. Because of the potential of these measures, it is vital to model these effects more carefully and explore technology that may realize their advantages.

Taylor, Charles R.↗

Geology of the Venus equatorial region from Pioneer Venus radar imaging

The data collected by Pioneer Venus have been used to map the broad radar characteristics of the equatorial region on the basis of radar brightness and texture. Seven units are defined and are used to assess the geologic character of the equatorial region. The correspondence between features in the 15 deg region of overlap between Pioneer Venus ans Venera 15/16 images is examined and used to extend units mapped in the northern high latitudes into the equatorial region. Using the properties of rms slope, reflectivity, the scattering behavior of the surface, and topographic signature, seven physiographic units are mapped on the basis of the distribution of the radar units in the equatorial region and are identified by increasing complexity as plains, dark halo plains, upland rises, upland plateaus, interhighland tectonic zones, tectonically segmented linear highlands, and tectonic junctions. The distribution of larger circular structures interpreted as coronae is also examined.

Senske, D. A.↗

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗

Operando visualization of porous metal additive manufacturing with foaming agents through high-speed x-ray imaging

Porous metals find extensive applications in soundproofing, filtration, catalysis, and energy-absorbing structures, thanks to their unique internal pore structure and high specific strength. In recent years, there has been an increasing interest in fabricating porous metals using additive manufacturing (AM), leveraging its unique advantages, including improved design freedom, spatial material control, and cost-effective small-batch production. In this study, we conducted pioneering operando visualization of AM porous metal using a laser powder bed fusion (L-PBF) setup combined with a high-speed synchrotron x-ray imaging system. Single track printing experiments using Ti6Al4V (Ti64) combined with titanium hydride (TiH 2 ) and sodium carbonate (Na 2 CO 3 ) as foaming agents, with varying mixing ratios were performed under different processing conditions. Here. the results elucidate the dynamic development of porosity formation. The average pore size is significantly influenced by the particle size of foaming agents when pore coalescence is absent. For all foaming agent content tested in the current study, the number of pores is found to be more sensitive to changes in laser power than in laser scanning speed. Increasing linear energy density (increasing laser power or reducing laser scanning speed) promotes the foaming agent activation thereby porosity formation. However, high linear energy density skews pore distribution towards the surface despite forming deeper melt pools. In addition, the impact of additional factors including foaming agent's laser absorptivity and decomposition kinetics with respect to AM time scales should be carefully considered to avoid ineffective activation of foaming agents during the AM of porous metals.

36 MATERIALS SCIENCE↗

Hydrogen Adsorption, Reactivity, and Catalysis on Colloidal Iron Carbide Nanoparticles

Supported iron carbide particles have long served as catalysts for CO hydrogenation (the Fischer Tropsch synthesis, FTS) and continue to be attractive. Despite this, little is known about their chemistry. Reported here is a colloidal Fe x C nanoparticle (NP) model system that allows direct observation of surface hydrogen and CO, as well as quantification of the surface H. Dodecylamine-capped Fe x C NPs (DDA-Fe x C NPs) were synthesized through solution-phase carburization of Fe NPs and form stable colloids in low-polarity organic solvents. Treatment of these colloids with H 2 or D 2 produced highly hydrogenated materials, and FTIR spectra of DDA-Fe x C-D n showed that most of the D binds to carbides, with at least four distinct v(C-D) modes. The surface C-H(D) bonds were all reactive, transferring hydrogen to alkenes and other reagents in solution. Titration and ICP measurements showed a 0.17:1 ratio of added H:total Fe, or roughly. 40 H per 1.8 nm DDA-Fe x C NP. Conversely, CO was preferentially bound to surface Fe sites, with FTIR spectra showing a single broad v(CO) that shifted with CO coverage or co-adsorption of H 2 . The DDA-Fe x C NPs were active catalysts for both olefin hydrogenation and the FTS, under mild conditions and without catalyst pre-treatment. The CO hydrogenation reactions yielded a broad distribution of long-chain linear paraffins and olefins. Though quantitative comparisons with typical FTS results are not possible because of our use of sealed batch reactors and other factors, the observations of high catalytic reactivity demonstrate the relevance of this model system to iron-carbide catalysis. DFT calculations on model slab surfaces with varying iron carbide stoichiometries revealed that the thermodynamically preferred surface adsorption sites are C for H ads and Fe for CO ads . A variety of binding sites and binding energies were found for each adsorbate. We are unaware of previous studies indicating that a diverse array of C-H bonds is the primary source of reactive H on iron carbides. Experimentally, the diversity of *C-H sites was evident in reactions with H-atom donors and abstractors of different strengths, from both the reaction stoichiometries and IR spectra. The different surface–H binding energies correlate with the v(C-D) stretching frequencies. Furthermore, these insights into complex iron carbide surfaces and catalysis could assist catalyst design, and they showcase the importance of stoichiometric studies of reaction intermediates.

08 HYDROGEN↗

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beyond Fisher forecasting for cosmology

The planning and design of future experiments rely heavily on forecasting to assess the potential scientific value provided by a hypothetical set of measurements. The Fisher information matrix, due to its convenient properties and low computational cost, provides an especially useful forecasting tool. However, the Fisher matrix only provides a reasonable approximation to the true likelihood when data are nearly Gaussian distributed and observables have nearly linear dependence on the parameters of interest. Also, Fisher forecasting techniques alone cannot be used to assess their own validity. Thorough sampling of the exact or mock likelihood can definitively determine whether a Fisher forecast is valid, though such sampling is often prohibitively expensive. Here we propose a simple test, based on the Derivative Approximation for likelihoods (DALI) technique, to determine whether the Fisher matrix provides a good approximation to the exact likelihood. We show that the Fisher matrix becomes a poor approximation to the true likelihood in regions where two-dimensional slices of level surfaces of the DALI approximation to the likelihood differ from two-dimensional slices of level surfaces of the Fisher approximation to the likelihood. We demonstrate that our method accurately predicts situations in which the Fisher approximation deviates from the true likelihood for various cosmological models and several data combinations, with only a modest increase in computational cost compared to standard Fisher forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiparticle-hole excitations in nuclei near N = Z = 20: $^{41}$K

This experimental study of high-spin structure near the N = Z = 20 region focuses on 41 K and also reports three newly observed γ transitions in 41 Ca from the same reaction. High-spin states were populated using the 26 Mg( 18 O, p2nγ) 41 K and 26 Mg( 18 O, 3nγ) 41 Ca reactions at a beam energy of 50 MeV at the Florida State University John D. Fox Superconducting Linear Accelerator Laboratory, employing the FSU high-purity germanium detector array. The level scheme of 41 K was extended to 12325 keV, possibly up to J π = 25/2 − or 27/2 + , by means of 25 new γ-ray transitions, and that of 41 Ca to 9916 keV. Linear polarization and angular-distribution measurements were used to provide spin and parity information for several states in the 41 K decay scheme. The results are compared with spsdpf cross-shell shell-model calculations using the FSU interaction. Configurations involving zero or one nucleon promoted from the sd to the fp shell reproduce the energies of many known states reasonably well, while multi-particle excitations reveal a more complex interplay of single-particle motion, collective pairing, and deformation, posing an interesting challenge for future theoretical work.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗

Evaluating Space Object Conjunction Probabilities Using Characteristic Function Inversion

This report discusses an approach to computing the probability of a conjunction between two space objects in the short-term encounter scenario. A conjunction is defined here as an event where the miss distance between the objects is less than some specified value. The scenario assumptions are that the motion of the objects is linear, their positions are Gaussian distributed, and their velocities are known and constant. Under these assumptions, the squared-miss distance is shown to have the generalized chi-square distribution. An established statistical technique called characteristic function inversion is employed to evaluate the distribution and obtain conjunction probabilities. The method is closely related to a recent approach based on moment generating function inversion, and a qualitative comparison of the approaches is provided. Last, the method is tested on several benchmark test cases where it agrees with numerical integration on the cases with conjunction probabilities above 10 –12 . However, the exact probability in these cases is usually not needed and this probability can be bounded above using an independent Gaussian approximation. Overall, the report shows how to compute conjunction probabilities using a standard statistical method, though numerical integration seems to perform equally well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An entropy-based debiasing approach to quantifying experimental coverage for novel applications of interest in the nuclear community

This manuscript proposes a novel information-theoretic approach to the quantification of experimental relevance, i.e., coverage, to achieve optimal data assimilation results for nuclear engineering applications. Specifically, this work posits the need for a new metric, called coverage (q C ) of an application’s quantity of interest, i.e., eigenvalue or power peaking for an advanced reactor concept, defined herein as the theoretically maximum achievable reduction in the quantity’s uncertainty given measurements from a pool of experiments in a manner that is independent of the data assimilation procedure employed. Currently, reduction in a quantity’s uncertainty is strongly biased by the underlying assumptions of the assimilation procedure to account for the under-determined nature of such problems and the similarity criterion employed to identify relevant experiments. To address this challenge, this work has developed a coverage metric, q C , based on mutual information, which establishes a new conceptual framework for assessing coverage, one that is independent of the model parameters and responses degree of variations in both the experimental and application domains, i.e., linear vs non-linear, and their prior uncertainty distributions, i.e., Gaussian vs. non-Gaussian. The q C is an entropic measure capable of addressing coverage for general nonlinear problems with non-Gaussian uncertainties and inclusive of the measurement uncertainties from multiple experiments. Numerical experiments from manufactured analytical problems as well as a set of benchmarks from the ICSBEP handbook are employed to demonstrate its theoretical and practical performance as compared to the c k -based experiment selection methodology, commonly employed in the neutronic community. The manuscript then employs other well-known adaptations to existing data assimilation methodologies for nonlinear and non-Gaussian problems capable of achieving the coverage posited by q C .

Bayesian data assimilation↗

Supercooled Liquid Water Detection Capabilities from Ka-Band Doppler Profiling Radars: Moment-Based Algorithm Formulation and Assessment

The occurrence of supercooled liquid water in mixed-phase cloud (MPC) affects their cloud microphysical and radiative properties. The prevalence of MPCs in the mid- and high latitudes translates these effects to significant contributions to Earth’s radiative balance and hydrological cycle. The current study develops and assesses a radar-only, moment-based phase partition technique for the demarcation of supercooled liquid water volumes in arctic, MPC conditions. The study utilizes observations from the Ka band profiling radar, the collocated high spectral resolution lidar, and ambient temperature profiles from radio sounding deployments following a statistical analysis of 5.5 years of data (January 2014–May 2019) from the Atmospheric Radiation Measurement observatory at the North Slope of Alaska. The ice/liquid phase partition occurs via a per-pixel, neighborhood-dependent algorithm based on the premise that the partitioning can be deduced by examining the mean values of locally sampled probability distributions of radar-based observables and then compare those against the means of climatologically derived, per-phase probability distributions. Analyzed radar observables include linear depolarization ratio (LDR), spectral width, and vertical gradients of reflectivity factor and radial velocity corrected for vertical air motion. Results highlight that the optimal supercooled liquid water detection skill levels are realized for the radar variable combination of spectral width and reflectivity vertical gradient, suggesting that radar-based polarimetry, in the absence of full LDR spectra, is not as critical as Doppler capabilities. The cloud phase masking technique is proven particularly reliable when applied to cloud tops with an Equitable Threat Score (ETS) of 65%; the detection of embedded supercooled layers remains much more uncertain (ETS = 27%).

54 ENVIRONMENTAL SCIENCES↗

Stochastic Averaging of Radiative Transfer Coefficients for Relativistic Electrons

Abstract Synchrotron emissivities, absorptivities, and Faraday rotation and conversion coefficients are needed in modeling a variety of astrophysical sources, including Event Horizon Telescope (EHT) sources. We develop a method for estimating transfer coefficients that exploits their linear dependence on the electron distribution function, decomposing the distribution function into a sum of parts each of whose emissivity can be calculated easily. We refer to this procedure as stochastic averaging and apply it in two contexts. First, we use it to estimate the emissivity of an isotropicκdistribution function with a high-energy cutoff. The resulting coefficients can be evaluated efficiently enough to be used directly in ray-tracing calculations, and we provide an example calculation. Second, we use stochastic averaging to assess the effect of subgrid turbulence on the volume-averaged emissivity and along the way provide a prescription for a turbulent emissivity. We find that for parameters appropriate to EHT sources turbulence reduces the emissivity slightly. In the infrared, turbulence can dramatically increase the emissivity.

Astronomy & Astrophysics↗

GENERAL INSTABILITY OF ORTHOGONALLY STIFFENED CYLINDRICAL SHELLS

Earlier research at the National Aeronautics Research Institute (N.L.R.), Amsterdam, which forms the basis of recent work is reviewed. This early work refers to 2 schemes: the orthotropic shell and, in view of buckling modes where the half wave length is of the order of the ring distance, the shell with continuously distributed stringers and discrete rings. Linear theory is considered to be adequate for these structures, where the imperfections are small in comparison to the height of the ring sections. Recent developments account for pressure difference in addition to axial compression, for the correct stiffness matrix of skin panels in the post-buckling stage and for stringer bending due to hoop stresses in the skin, which are of importance as has been shown by the investigation of the post-buckling behaviour. Numerical data for the stiffness matrix of skin panels have been established. Numerical evaluation of the stability equation has not been performed as yet.

SHELL STABILITY↗

Comparison of Experimental and Theoretical Zero-Lift Wave-Drag Results for Various Wing-Body-Tail Combinations at Mach Numbers up to 1.9

Comparisons are made of experimental and theoretical zero-lift wave drag for several nose shapes, wing-body combinations, and models of current airplanes at Mach numbers up to 1.0. The experimental data were obtained from tests in the Ames 6- by6-foot supersonic wind tunnel and at the NACA Wallops Island facility. The theoretical drag was found by use of linear theory utilizing model area distributions. The agreement between theoretical and experimental zero-lift wave-drag coefficients was generally very good, especially for a fuselage or for fuselage-wing combinations that were vertically symmetrical. For other models that had rapid changes in body shape and/or were not vertically symmetrical, the agreement of theory with experiment ranged from fair to poor, depending on the severity of the change in shape.

DRAG COEFFICIENT↗

Analysis of the Flow About Delta Wings with Leading Edge Separation at Supersonic Speeds

A research program was conducted to develop an improved theoretical flow model for the flow about sharp edge delta wings with leading-edge separation at supersonic speeds. The flow model incorporates a representation of the secondary separation region which occurs just inboard of the leading edge on such wings and is based on a slender-wing theory whereby the full three-dimensional problem is reduced to a quasi two-dimensional problem in the cross-flow plane. The secondary separation region was modeled by a surface distribution of singularities or a linearized type of cavity representation. The primary vortex and separation were modeled by a concentrated vortex and cut in the cross-flow potential which represents its feeding sheet. The cross-flow solutions for the cavity model were obtained, but these solutions have physical significance only in a very restricted range of angle of attack. The reasons for the failure of the flow model are discussed. The analysis is presented so that other interested researchers may critically review the work.

Nenni, J. P.↗

Study of a light curve of Beta Persei at 3428 A

Photometric light-curve data obtained during Copernicus satellite observations of Beta Persei in the Balmer continuum at 3428 A are analyzed and reduced. A small but significant fluctuation with a period of very nearly 0.069167 day in the photometric counts is tentatively attributed to a small perturbation in spacecraft attitude introduced each time the satellite passes from daylight into darkness. The light curve is intensity-rectified, phase-rectified, and solved on the assumption that the limb darkening follows the linear cosine law. The resulting distributions of luminosity among the three components of Beta Persei are plotted, and it is shown that the effective temperature of Algol B may be higher than 5000 K.

Chen, K.-Y.↗