Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discrete velocity”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Heat transfer in a discrete-velocity gas

Heat transfer of a discrete-velocity gas composed of identical hard sphere molecules between two parallel plates is studied numerically. A finite-difference method recently developed by the authors is used for solving the discrete Boltzmann equations. The temperature distributions between the plates of temperature ratios 2.0 and 4.0, and the heat transfer of temperature ratio 4.0 are presented and are compared with results of other methods.

Inamuro, Takaji↗
NASA NTRS2007

Theory of the Lattice Boltzmann Equation: Symmetry properties of Discrete Velocity Sets

In the lattice Boltzmann equation, continuous particle velocity space is replaced by a finite dimensional discrete set. The number of linearly independent velocity moments in a lattice Boltzmann model cannot exceed the number of discrete velocities. Thus, finite dimensionality introduces linear dependencies among the moments that do not exist in the exact continuous theory. Given a discrete velocity set, it is important to know to exactly what order moments are free of these dependencies. Elementary group theory is applied to the solution of this problem. It is found that by decomposing the velocity set into subsets that transform among themselves under an appropriate symmetry group, it becomes relatively straightforward to assess the behavior of moments in the theory. The construction of some standard two- and three-dimensional models is reviewed from this viewpoint, and procedures for constructing some new higher dimensional models are suggested.

Rubinstein, Robert↗
NASA NTRS1982

Mass loss at discrete velocities in Be stars: Evidence for non-radiatively driven wind

The temporal behavior of the high velocity, relatively narrow absorption components in C IV, Si III, IV observed in the pole-one Be stars psi Ori, 66 Oph, and HR 4009 and what these observations contribute toward our understanding of the Be phenomenon are discussed. Ground based and X-ray data on these stars are compared as well as their UV differences with similarly active Be stars with larger projected rotational velocities.

Peters, G. P.↗
NASA NTRS2000

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre↗
NASA NTRS1972

Numerical methods for solving the Boltzmann equation (a review)

The methods are reviewed which are utilized in principal attempts to obtain the numerical solution or modeling of the Boltzmann equation over a broad range of Knudsen numbers. The primary methods considered are the Monte Carlo and the discrete velocities methods. The conculsions drawn from the analysis include the following: (1) The Monte Carlo methods are not well suited in the area of small Knudsen numbers. (2) Among the Monte Carlo methods, the Bird method appears to be the most attractive, since it is more directly related to the Boltzmann equation. (3) The deterministic methods, which include the discrete ordinate technique, offer great possibilities but require exceedingly large computer times. (4) The use of approximating equations in combination with the discrete velocities method will possibly improve computation time and reduce the required memory volume.

Limar, Y. F.↗
NASA NTRS1997

X-Ray Flare Characteristics in the B2e Star Lambda Eridani (ROSAT)

We document the results of a simultaneous wavelength monitoring on the B2e star (lambda) Eri. This campaign was carried out from ground stations and with the ROSAT, ASCA, IUE, and Voyager 2 space platforms during a week in February-March 1995; a smaller follow-up was conducted in September 1995. During the first of these intervals (lambda) Eri exhibited extraordinary wind and disk-ejection activity. The ROSAT/HRI X-ray light curves showed no large flares such as the one the ROSAT/PSCA observed in 1991. However, possible low level fluctuations in the February-March ROSAT data occurred at the same time as unusual activity in H(alpha) He I (lambda)6678, He II (lambda)1640, and the C IV doublet. For example, the hydrogen and helium lines exhibited an emission in the blue half of their profiles, probably lasting several hours. The C IV lines showed a strong high-velocity Discrete Absorption Component (DAC) accompanied by unusually strong absorption at lower velocities. The helium line activity suggests that a mass ejection occurred at the base of the wind while the strong C III (Voyager) and C IV (IUE) lines implies that shock interactions occurred in the wind flow. It is not clear that the X-ray elevations are directly related to the strong C IV absorptions because the former changed on a much more rapid timescale than absorptions in the C IV lines. Within hours of the mild X-ray flux variations found by ROSAT on February 28, the Voyager UVS observed a "ringing" that decayed over three 3-hr. cycles. The amplitude of these fluctuations was strong (50%) at (lambda)(lambda)950-1100, decreased rapidly with wavelength, and faded to nondetection longward of (lambda)1300. Various considerations indicate that these continuum variations were not due to an instrumental pathology in the UVS. Rather, they appear to be due to a time-dependent flux deficit in the (lambda)(lambda)950-1250 region. We outline a scenario in which a dense plasma structure over the star's surface is heated and cooled quasi-periodically to produce such flux changes. Observations of new examples of this phenomenon are badly needed. Amateur astronomers can make a significant contribution to its understanding by searching for ringing in light curves of Be stars during their outburst phases. Finally we draw attention to an increase in the emission of the H(alpha) line that occurred at about the time the FUV ringing started. This increased emission hints that approximately 50,000K plasma near the star's surface can influence the circumstellar disc at approximately 12R. by its increased Lyman continuum flux.

Smith, Myron A.↗
NASA NTRS1994

Internal motions of HII regions and giant HII regions

We report new echelle observations of the kinematics of 30 HII regions in the Large Magellanic Clouds (LMC), including the 30 Doradus giant HII region. All of the HII regions possess supersonic velocity dispersions, which can be attributed to a combination of turbulent motions and discrete velocity splitting produced by stellar winds and/or embedded supernova remnants (SNRs). The core of 30 Dor is unique, with a complex velocity structure that parallels its chaotic optical morphology. We use our calibrated echelle data to measure the physical properties and energetic requirements of these velocity structures. The most spectacular structures in 30 Dor are several fast expanding shells, which appear to be produced at least partially by SNRs.

Chu, You-Hua↗
NASA NTRS1993

Goddard high-resolution spectrograph observations of narrow discrete stellar wind absorption features in the ultraviolet spectrum of the O7.5 III star Xi Persei

We report the observation of transient narrow absorption components (NACs) in the stellar wind of the O giant Xi Per. Two sets of GHRS observations of the Si IV ultraviolet resonance doublet have been obtained. These features are extremely weak, with column densities of approximately 10 exp 12/sq cm and optical depths of order 0.1. The features are narrow, less than 30 km/s, and seem to occur in groups. If the NACs are due to the 1393 A component, they represent previously undetected low-velocity discrete absorption components at V(rad) below -600 km/s. If they are high-velocity features on the 1402 A doublet component, they may represent the decay phase of the discrete absorption components at the terminal velocity. In either case, they are a new aspect of the NAC phenomenon that could not have been detected with previous ultraviolet spectrographs.

Shore, Steven N.↗
NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗
NASA NTRS1988

The detection of high-velocity outflows from M8E-IR

A high-resolution (0.059/cm) M band (4.6 micron) spectrum of the embedded young stellar object M8E-IR is presented and discussed. The spectrum shows strong absorption to large blueshifts in the rotational lines of the fundamental vibrational band, v = 1-0, of CO. The absorption is interpreted as being due to gas near to, and flowing from, the central object. The outflowing gas is warm (95-330 K) and consists of discrete velocity components with the very high velocities of 90, 130, 150, and 160 km/s. On the basis of a simple model, it is estimated that the observed outflows are less than 100 yr old.

Mitchell, George F.↗
DOE OSTI2024

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Structured ionized winds shooting out from a quasar at relativistic speeds

Evidence indicates that supermassive black holes (SMBHs) exist at the centres of most galaxies. Their mass correlates with the galactic bulge mass, suggesting a coevolution with their host galaxies, most likely through powerful winds. X-ray observations have detected highly ionized winds outflowing at sub-relativistic speeds from the accretion disks around SMBHs. However, the limited spectral resolution of present X-ray instruments has left the physical structure and location of the winds poorly understood, hindering accurate estimates of their kinetic power. Here, in this study, the first X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the luminous quasar PDS 456 is reported. The high-resolution spectrometer Resolve aboard XRISM enabled the discovery of five discrete velocity components outflowing at 20–30% of the speed of light. This demonstrates that the wind structure is highly inhomogeneous, which probably consists of up to a million clumps. The mass outflow rate is estimated to be 60–300 solar masses per year, with the wind kinetic power exceeding the Eddington luminosity limit. Compared with the galaxy-scale outflows, the kinetic power is more than three orders of magnitude larger, whereas the momentum flux is ten times larger. These estimates disfavour both energy-driven and momentum-driven outflow models. This suggests that such wind activity occurs in less than 10% of the quasar phase and/or that its energy/momentum is not efficiently transferred to the galaxy-scale outflows owing to the clumpiness of the wind and the interstellar medium.

Audard, Marc [Univ. of Geneva, Versoix (Switzerlan↗
NASA NTRS1989

A multidomain spectral collocation method for the Stokes problem

A multidomain spectral collocation scheme is proposed for the approximation of the two-dimensional Stokes problem. It is shown that the discrete velocity vector field is exactly divergence-free and we prove error estimates both for the velocity and the pressure.

Landriani, G. Sacchi↗
NASA NTRS1994

Spectropolarimetric variability in OB stars

Data from the Wisconsin spectropolarimetric monitoring program provide information about the polarimetric variability of numerous types of stars. Examining the data that exist thus far, we find, not surprisingly, that many OB stars, especially OB supergiants and Be stars, are polarimetric variables. Be stars show the most extreme variable as a class. Results from the Be star pi Aqr are presented which illustrate the dramatic polarimetric variability that such stars can exhibit. Contemporaneous IUE data show that the development of low velocity discrete absorption components in the UV wind lines may be associated with the onset of polarimetric outbursts. These preliminary findings add further support to the picture that winds from hot stars are non-homogeneous and variable.

Bjorkman, K. S.↗
NASA NTRS2010

Navier-Stokes Dynamics by a Discrete Boltzmann Model

This work investigates the possibility of particle-based algorithms for the Navier-Stokes equations and higher order continuum approximations of the Boltzmann equation; such algorithms would generalize the well-known Pullin scheme for the Euler equations. One such method is proposed in the context of a discrete velocity model of the Boltzmann equation. Preliminary results on shock structure are consistent with the expectation that the shock should be much broader than the near discontinuity predicted by the Pullin scheme, yet narrower than the prediction of the Boltzmann equation. We discuss the extension of this essentially deterministic method to a stochastic particle method that, like DSMC, samples the distribution function rather than resolving it completely.

Rubinstein, Robet↗
DOE OSTI2026

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species f s on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where f s is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6–60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

Discontinuous Galerkin↗
NASA NTRS1982

Flute mode waves near the lower hybrid frequency excited by ion rings in velocity space

Discrete emissions at the lower hybrid frequency are often seen on the S3-3 satellite. Simultaneous observation of perpendicularly heated ions suggests that these ions may provide the free energy necessary to drive the instability. Studies of the dispersion relation for flute modes excited by warm ion rings in velocity space show that waves are excited with real frequencies near the lower hybrid frequency and with growth rates ranging from about 0.01 to 1 times the ion cyclotron frequency. Numerical results are therefore consistent with the possibility that the observed ions are the free energy source for the observed waves.

Cattell, C.↗