Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

transient↗
NASA NTRS2018

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

Column Density↗
DOE OSTI2022

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1960

A Graphical Method for Estimating Ion-Rocket Performance

Equations relating the critical temperature and ion current density for surface ionization of cesium on tungsten are derived for the cases of zero and finite electric fields at the ion-emitting surface. These equations are used to obtain a series of graphs that can be used to solve many problems relating to ion-rocket theoretical performance. The effect of operation at less than space-charge-limited current density and the effect of nonuniform propellant flux onto the ion-emitting surface are also treated.

Reynolds, Thaine W.↗
NASA NTRS1994

Manual for automatic generation of finite element models of spiral bevel gears in mesh

The goal of this research is to develop computer programs that generate finite element models suitable for doing 3D contact analysis of faced milled spiral bevel gears in mesh. A pinion tooth and a gear tooth are created and put in mesh. There are two programs: Points.f and Pat.f to perform the analysis. Points.f is based on the equation of meshing for spiral bevel gears. It uses machine tool settings to solve for an N x M mesh of points on the four surfaces, pinion concave and convex, and gear concave and convex. Points.f creates the file POINTS.OUT, an ASCI file containing N x M points for each surface. (N is the number of node points along the length of the tooth, and M is nodes along the height.) Pat.f reads POINTS.OUT and creates the file tl.out. Tl.out is a series of PATRAN input commands. In addition to the mesh density on the tooth face, additional user specified variables are the number of finite elements through the thickness, and the number of finite elements along the tooth full fillet. A full fillet is assumed to exist for both the pinion and gear.

Bibel, G. D.↗
NASA NTRS2018

Acoustic Radiation Pressure

The theoretical foundation of acoustic radiation pressure in plane wave beams is reexamined. It is shown from finite deformation theory and the Boltzmann-Ehrenfest Adiabatic Principle that the Brillouin stress tensor (BST) is the radiation stress in Lagrangian coordinates (not Eulerian coordinates) and that the terms in the BST are not the momentum flux density and mean excess Eulerian stress but are simply contributions to the variation in the wave oscillation period resulting from changes in path length and true wave velocity, respectively, from virtual variations in the strain. It is shown that the radiation stress in Eulerian coordinates is the mean Cauchy stress (not the momentum flux density, as commonly assumed) and that Langevin's second relation does not yield an assessment of the mean Eulerian pressure, since the enthalpy used in the traditional derivations is a function of the thermodynamic tensions - not the Eulerian pressure. It is shown that the transformation between Lagrangian and Eulerian quantities cannot be obtained from the commonly-used expansion of one of the quantities in terms of the particle displacement, since the expansion provides only the difference between the value of the quantity at two different points in Cartesian space separated by the displacement. The proper transformation is obtained only by employing the transformation coefficients of finite deformation theory, which are defined in terms of the displacement gradients. Finite deformation theory leads to the result that for laterally unconfined, plane waves the Lagrangian and Eulerian radiation pressures are equal with the value (1/4)(2K) along the direction of wave propagation, where (K) is the mean kinetic energy density, and zero in directions normal to the propagation direction. This is contrary to the Langevin result that the Lagrangian radiation pressure in the propagation direction is equal to (2K) and the BST result that the Eulerian radiation pressure in that direction is the momentum flux density.

Finite deformation theory↗
NASA NTRS2019

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

A NASA robotic refueling mission experiment on board the International Space Station (ISS) was designed to repeatedly transfer simulated cryogenic propellant between two dewars. After each metered transfer, the cryogen was to be vented overboard via sublimation or evaporation. ISS payload providers must conduct analyses to demonstrate that any planned gaseous venting generate no more than a certain level of material that may interfere with optical measurements from other payloads that may be located nearby. This requirement is expressed in terms of a maximum column number density (CND). For the experiment under consideration, fluid droplets likely would accompany any such vapor. Earlier related studies led to development of analytical expressions for column density along general paths for a number of practical cases under steady conditions, including droplet evaporation. These expressions are not appropriate for rapid evaporation of small droplets however, since the droplets cannot sustain mass flow rates consistent with steady conditions. In this work, analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where it is found the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next, the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction, but at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed in this framework.

Woronowicz, M. S.↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators

Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Superconductivity in the doped quantum spin liquid on the triangular lattice

Abstract Broad interest in quantum spin liquid (QSL) phases was triggered by the notion that they can be viewed as insulating phases with preexisting electron pairs, such that upon light doping they might automatically yield high temperature superconductivity. Yet despite intense experimental and numerical efforts, definitive evidence showing that doping QSLs leads to superconductivity has been lacking. We address the problem of a lightly doped QSL through a large-scale density-matrix renormalization group study of the t - J model on finite-circumference triangular cylinders with a small but nonzero concentration of doped holes. We provide direct evidences that doping QSL can naturally give rise to d -wave superconductivity. Specifically, we find power-law superconducting correlations with a Luttinger exponent, K s c ≈ 1, which is consistent with a strongly diverging superconducting susceptibility, $${\chi }_{sc} \,\sim\, {T}^{-(2\,-\,{K}_{sc})}$$ χ s c ~ T − ( 2 − K s c ) as the temperature T → 0. The spin–spin correlations—as in the undoped QSL state—fall exponentially which suggests that the superconducting pair-pair correlations evolve smoothly from the insulating parent state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Verification of the kinetic electron role in the microinstabilities in a negative triangularity model equilibrium

Effect of kinetic electrons on negative triangularity plasmas has been investigated and compared against the corresponding positive triangularity plasmas, using the global gyrokinetic code X-point Gyrokinetic Code with scale-separated delta-f option without Coulomb collisions. Our model magnetic equilibria have strong positive and negative triangularities and weak magnetic shear. However, unusually large ρ i > a and low density plasmas are chosen to maximize the nonlocal effect to investigate the finite ρ i effect and to be clearly away from kinetic ballooning modes. Similar conclusions to previous flux tube and global simulations have been obtained in this highly nonlocal model plasma: it is essential to include kinetic electrons in the micro-instability study of negative triangularity plasmas. Most physics findings agree with existing reports, with some disagreement. We offer a new “effective trapping fraction” concept that can add to the explanation of the growth rate difference between NT and PT plasmas, pointing to the significant variation in trapped particle fractions that have turning points in the mode growth regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Band Energy Dependence of Defect Formation in the Topological Semimetal Cd3As2

Cadmium Arsenide (Cd3As2) is a prototypical Dirac semimetal that manifests topological properties in a 3D bulk material. In defect-free Cd3As2, the Fermi level lies at a minimum in the density of states at the Dirac point, but experimentally it forms with excess electron carriers and an elevated EF, thereby masking the topological features. To computationally study the self-doping of Cd3As2, we combine density functional theory (DFT) calculations for defect formation energies with quasi-particle self-consistent GW (QSGW) electronic structure calculations. We demonstrate an innate dependence of the point defect formation energies on carrier concentrations and use the QSGW calculated density of states to extrapolate formation energies to arbitrary electron concentrations. This approach allows the quantitative modeling of thermodynamic defect equilibria in topological semimetals and is used to predict how Cd3As2 growth conditions affect the position of EF relative to the Dirac point.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2020

Covariant Density Functional Theory in Nuclear Physics and Astrophysics

How does subatomic matter organize itself? Neutron stars are cosmic laboratories uniquely poised to answer this fundamental question that lies at the heart of nuclear science. Newly commissioned rare isotope facilities, telescopes operating across the entire electromagnetic spectrum, and ever more sensitive gravitational wave detectors will probe the properties of neutron-rich matter with unprecedented precision over an enormous range of densities. A coordinated effort between observation, experiment, and theoretical research is of paramount importance for realizing the full potential of these investments. Theoretical nuclear physics provides valuable insights into the properties of neutron-rich matter in regimes that are not presently accessible to experiment or observation. In particular, nuclear density functional theory is likely the only tractable framework that can bridge the entire nuclear landscape by connecting finite nuclei to neutron stars. This compelling connection is the main scope of the present review.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

A design of an all digital FM discriminator.

This paper describes the design of an all digital FM discriminator. The system can be easily constructed using ordinary commercially available, digital IC's. In addition, the output SNR above threshold is computed. It is shown that the output noise consists of two terms: the first being the expected output gaussian noise having a parabolic power spectral density and the second being a quantization noise produced by the finite-level detection process used in the receiver.

Garodnick, J.↗
NASA NTRS1975

Theory of type 3b solar radio bursts

During the initial space-time evolution of an electron beam injected into the corona, the strong beam-plasma interaction occurs at the head of the beam, leading to the amplification of a quasi-monochromatic large-amplitude plasma wave that stabilizes by trapping the beam particles. Oscillation of the trapped particles in the wave troughs amplifies sideband electrostatic waves. The sidebands and the main wave subsequently decay to observable transverse electromagnetic waves through the parametric decay instability. This process gives rise to the elementary striation bursts. Owing to velocity dispersion in the beam and the density gradient of the corona, the entire process may repeat at a finite number of discrete plasma levels, producing chains of elementary bursts. All the properties of the type IIIb bursts are accounted for in the context of the theory.

Smith, R. A.↗
NASA NTRS1976

Theory of type IIIb solar radio bursts

During the initial space-time evolution of an electron beam injected into the corona, the strong beam-plasma interaction occurs at the head of the beam, leading to the amplification of a quasi-monochromatic large-amplitude plasma wave that stabilizes by trapping the beam particles. Oscillation of the trapped particles in the wave troughs amplifies sideband electrostatic waves. The sidebands and the main wave subsequently decay to observable transverse electromagnetic waves through the parametric decay instability. This process gives rise to the elementary striation bursts. Owing to velocity dispersion in the beam and the density gradient of the corona, the entire process may repeat at a finite number of discrete plasma levels, producing chains of elementary bursts. All the properties of the type IIIb bursts are accounted for in the context of the theory.

Smith, R. A.↗