Search NASA⌕ Search

SEARCH · Search NASA

Results for “compressibility factor”

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 181 records · Page 10

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Large Scale Tensor Factorization via Parallel Sketches

Tensor factorization methods have recently gained increased popularity. A key feature that renders tensors attractive is the ability to directly model multi-relational data. In this work, we propose ParaSketch, a parallel tensor factorization algorithm that enables massive parallelism, to deal with large tensors. The idea is to compress the large tensor into multiple small tensors, decompose each small tensor in parallel, and combine the results to reconstruct the desired latent factors. Prior art in this direction entails potentially very high complexity in the (Gaussian) compression and final combining stages. Adopting sketching matrices for compression, the proposed method enjoys a dramatic reduction in compression complexity, and features a much lighter combining step. Moreover, theoretical analysis shows that the compressed tensors inherit latent identifiability under mild conditions, hence establishing correctness of the overall approach. Numerical experiments corroborate the theory and demonstrate the effectiveness of the proposed algorithm.

block term decomposition↗

Virtual texture analysis to understand microstructure effects on deformation twinning and detwinning behavior in BCC metals

Understanding and predicting deformation twinning contributions to plastic deformation in BCC metals has been a long-standing challenge due to the interplay with dislocation slip and non-Schmid effects that render an asymmetry under tension and compression. This paper uses molecular dynamics simulations to understand the effect of unloading and grain orientation on deformation twinning in a nanocrystalline Fe (nc-Fe) system as a model BCC metal. The nc-Fe system is loaded under uniaxial stress tension and compression to understand the effect of grain orientation (Schmid effects) on deformation twinning behavior and the tension/compression asymmetry (non-Schmid effects). A new virtual texture analysis “VirTex” tool is used to understand the role of grain orientation on the nucleation of twins and their contributions to the observed stress–strain response. For certain grain orientations, the twinnability is observed to be different in tension and compression. In addition, the flow stress accommodation from twins in certain grains is observed to be different in tension and compression and different from that for the grains. Subsequent unloading leads to detwinning in the deformed microstructures, where the extent of detwinning depends on the strain from which the system is unloaded and on the morphology of the twin. Lastly, the simulations are carried out to analyze the role of the Schmid factor on the twinnability and asymmetry in tension and compression.

Kannan, Aadhithyan [University of Connecticut, Sto↗

Modeling ablator defects as a source of mix in high-performance implosions at the National Ignition Facility

Recent indirect drive inertial confinement fusion implosions on the National Ignition Facility (NIF) [Spaeth et al., Fusion Sci. Technol. 69, 25 (2016)] have crossed the threshold of ignition. However, performance has been variable due to several factors. One of the leading sources of variability is the quality of the high-density carbon (HDC) shells used as ablators in these experiments. In particular, these shells can have a number of defects that have been found to correlate with the appearance of ablator mix into the hot spot and a degradation in nuclear yield. These defects include pits on the ablator surface, voids in the ablator bulk, high-Z debris from the Hohlraum wall that adheres to the capsule surface, and finally the inherent granular micro-structure of the crystalline HDC itself. This paper summarizes high-resolution modeling of each of these mix sources in two recent high-performance NIF implosion experiments. The simulated impact from a range of individual capsule defects is found to be broadly consistent with the trends seen in experiment, lending credence to the modeling results and the details of the mixing process that they reveal. Interestingly, modeling of the micro-structure inherent to HDC shows that this perturbation source results in considerable mixing of the deuterium–tritium fuel with ablator material during the implosion. The reduction in fuel compression from this mix results in an approximately factor of two reduction in neutron yield in current implosions and emphasizes the importance of mitigating this significant performance degradation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Location of the radio emitting regions of interplanetary shocks

Twenty interplanetary type II radio bursts are analyzed to determine the location of the type II source region relative to the interplanetary shock. The first determination of a density-distance relationship (density model) appropriate for interplanetary type II source regions is reported. To determine source location, densities in type II source regions, derived from observed type II emission frequencies, are compared to ambient solar wind densities to determine whether source regions are in the ambient, upstream solar wind or in the compressed plasma behind the shock. Densities in the ambient solar wind upstream of each shock are computed by using a simple model to extrapolate solar wind plasma densities measured at 1 AU back to the shock front. Densities in the type II source regions are found to be enhanced relative to densities in the ambient solar wind by a factor which is close to the average shock density compression ratio. The simplest interpretation of this result is that the source is located in the compressed plasma within or behind the shock. Although it is possible that the emission is produced in enhanced density regions of the upstream solar wind, it is argued that the weight of evidence favors the compressed postshock plasma as the source site.

Lengyel-Frey, D.↗

Identifying Anomalous DESI Galaxy Spectra with a Variational Autoencoder

The tens of millions of spectra being captured by the Dark Energy Spectroscopic Instrument (DESI) provide tremendous discovery potential. In this work we show how Machine Learning, in particular Variational Autoencoders (VAE), can detect anomalies in a sample of approximately 200,000 DESI spectra comprising galaxies, quasars and stars. We demonstrate that the VAE can compress the dimensionality of a spectrum by a factor of 100, while still retaining enough information to accurately reconstruct spectral features. We then detect anomalous spectra as those with high reconstruction error and those which are isolated in the VAE latent representation. The anomalies identified fall into two categories: spectra with artefacts and spectra with unique physical features. Awareness of the former can help to improve the DESI spectroscopic pipeline; whilst the latter can lead to the identification of new and unusual objects. To further curate the list of outliers, we use the Astronomaly package which employs Active Learning to provide personalised outlier recommendations for visual inspection. In this work we also explore the VAE latent space, finding that different object classes and subclasses are separated despite being unlabelled. We demonstrate the interpretability of this latent space by identifying tracks within it that correspond to various spectral characteristics. For example, we find tracks that correspond to increasing star formation and increase in broad emission lines along the Balmer series. In upcoming work we hope to apply the methods presented here to search for both systematics and astrophysically interesting objects in much larger datasets of DESI spectra.

Nicolaou, C. [University Coll. London] (ORCID:0000↗

An Incremental Tensor Train Decomposition Algorithm

We present a new algorithm for incrementally updating the tensor train decomposition of a stream of tensor data. This new algorithm, called the tensor train incremental core expansion (TT-ICE) improves upon the current state-of-the-art algorithms for compressing in tensor train format by developing a new adaptive approach that incurs significantly slower rank growth and guarantees compression accuracy. This capability is achieved by limiting the number of new vectors appended to the TT-cores of an existing accumulation tensor after each data increment. These vectors represent directions orthogonal to the span of existing cores and are limited to those needed to represent a newly arrived tensor to a target accuracy. We provide two versions of the algorithm: TT-ICE and TT-ICE accelerated with heuristics (TT-ICE*). Here, we provide a proof of correctness for TT-ICE and empirically demonstrate the performance of the algorithms in compressing large-scale video and scientific simulation datasets. Compared to existing approaches that also use rank adaptation, TT-ICE* achieves 57× higher compression and up to 95% reduction in computational time.

97 MATHEMATICS AND COMPUTING↗

Mechanisms of Mitigating Chloride-Induced Stress Corrosion Cracking of Austenitic Steels by Laser Shock Peening

This study investigates the effect of laser shock peening (LSP) on the chloride-induced stress corrosion cracking (SCC) of 304 austenitic steels. LSP can induce a high compressive residual stress to a depth of 700 µm and plastic deformation structures of dislocations, deformation twins, and stacking faults. Constant-load SCC tests in MgCl2 solution suggested that LSP can retard the crack initiation and slow the crack growth. LSP-treated subsurface layers experience ductile fracture while the central regions exhibit intergranular SCC. The LSP-induced deformation structures may impede dislocation slips, while the LSP-induced compressive residual stress can lessen the stress intensity factor of crack tips and decrease the local stress for film rupture.

Materials Science↗

Plastic intensity factors for cracked plates subjected to biaxial loading

An elastic-plastic analysis is performed for an infinite plate under plane stress conditions which contains a finite line crack and is subjected to biaxial loading. In the large scale yielding range, the amplitude of the dominant singularity at the crack tip, the plastic stress (or strain) intensity factor, is found to depend on the magnitude of loading parallel to the crack direction. The predicted value for the fracture initiation stress which is based on a plastic intensity factor fracture criterion, increases under lateral, tensile loading. Compressive loading parallel to the crack has the opposite effect on the fracture initiation stress.

Hilton, P. D.↗

Time-marching solution of incompressible Navier-Stokes equations for internal flow

Primitive variables with central differencing on a staggered grid are used in the present, factored ADI finite-difference scheme for artificial compressibility method solution of the incompressible Navier-Stokes equations, leading to a close coupling between velocity and pressure that both enhances stability and eliminates the need for artificial damping. Computational efficiency is enhanced through the use of a spatially variable, fixed Courant number-based time-step. The numerical results obtained for a driven cavity at Re of 10,000, with local cell Re as high as 100, exhibits no flow variable spatial oscillations on a 40 x 40 stretched grid solution.

Soh, W. Y.↗

Supersonic boundary-layer transition on the LaRC F-106 and the DFRF F-15 aircraft. Part 1: Transition measurements and stability analysis

For the case of the F-15 flight tests, boundary layer transition was observed up to Mach numbers of 1.2. For very limited and specific flight conditions, laminar flow existed back to about 20 percent chord on the surface clean up glove. Hot film instrumentation was effective for locating the region of transition. For the F-106 flight tests, transition on the wing or vertical tail generally occurred very near the attachment line. Transition was believed to be caused by either attachment line contamination or strong cross flow development due to the high sweep angles of the test articles. The compressibility analysis showed that cross flow N-factors were in the range of 5 to 12 at transition.

Collier, Fayette S., Jr.↗

Enhancement of the no-vent fill process

This paper presents an analytical and experimental evaluation of an enhanced techniques for no-vent fill. The method entails injecting liquid through the top of the receiver vessel, thereby increasing surface area and agitation of the ullage/liquid interface. Both of these factors promote condensation induced ullage collapse, and reduce compressive impedance to the incoming liquid. The enhanced process was analyzed by modifying the surface area algorithm of an existing tank thermodynamic code to model a downward-pointing, conical jet impringing on a steadily rising liquid surface. Transient pressure and temperature measurements from several tests with Freon-114 were input into the revised model to calculate condensation rate as a function of fill level. By expressing these rates in dimensionless form (i.e., in terms of Stanton number and Prandtl number), an empirical correlation similar to the submerged jet model of Brown and Sonin (1989) was derived. This provided a basis for developing an expression which relates top fill to bottom fill performance.

Vaughan, David A.↗

Vortex tubes in turbulent flows: Identification, representation, reconstruction

In many cases the structure of a fluid flow is well-characterized by its vortices, especially for the purpose of visualization. In this paper we present a new algorithm for identifying vortices in complex flows. The algorithm produces a skeleton line along the center of a vortex by using a two-step predictor-corrector scheme. The vorticity vector field serves as the predictor and the pressure gradient (in the perpendicular plane) serves as the corrector. We describe an economical description of the vortex tube's cross-section: a 5-term truncated Fourier series is generally sufficient, and it compresses the representation of the flow by a factor of 4000 or more. We reconstruct the vortex tubes as generalized cylinders, providing a polygonal mesh suitable for display on a graphics workstation. We show how the reconstructed geometry of vortex tubes can be enhanced to help visualize helical motion in a static image.

Banks, David C.↗

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗

Strategies for on-chip digital data compression for X-ray pixel detectors

Here, the continued desire for X-ray pixel detectors with higher frame rates will stress the ability of application-specific integrated circuit (ASIC) designers to provide sufficient off-chip bandwidth to reach continuous frame rates in the 1 MHz regime. To move from the current 10 kHz to the 1 MHz frame rate regime, ASIC designers will continue to pack as many power-hungry high-speed transceivers at the periphery of the ASIC as possible. In this paper, however, we present new strategies to make the most efficient use of the off-chip bandwidth by utilizing data compression schemes for X-ray photon-counting and charge-integrating pixel detectors. In particular, we describe a novel in-pixel compression scheme that converts from analog to digital converter units to encoded photon counts near the photon Poisson noise level and achieves a compression ratio of >1.5x independent of the dataset. In addition, we describe a simple yet efficient zero-suppression compression scheme called "zeromask" (ZM) located at the ASIC's edge before streaming data off the ASIC chip. ZM achieves average compression ratios of >4x, >7x, and >8x for high-energy X-ray diffraction, ptychography, and X-ray photon correlation spectroscopy datasets, respectively. We present the conceptual designs, register-transfer level block diagrams, and the physical ASIC implementation of these compression schemes in 65 nm CMOS. When combined, these two digital compression schemes could increase the effective off-chip bandwidth by a factor of 6-12x.

47 OTHER INSTRUMENTATION↗

Compression through decomposition into browse and residual images

Economical archival and retrieval of image data is becoming increasingly important considering the unprecedented data volumes expected from the Earth Observing System (EOS) instruments. For cost effective browsing the image data (possibly from remote site), and retrieving the original image data from the data archive, we suggest an integrated image browse and data archive system employing incremental transmission. We produce our browse image data with the JPEG/DCT lossy compression approach. Image residual data is then obtained by taking the pixel by pixel differences between the original data and the browse image data. We then code the residual data with a form of variable length coding called diagonal coding. In our experiments, the JPEG/DCT is used at different quality factors (Q) to generate the browse and residual data. The algorithm has been tested on band 4 of two Thematic mapper (TM) data sets. The best overall compression ratios (of about 1.7) were obtained when a quality factor of Q=50 was used to produce browse data at a compression ratio of 10 to 11. At this quality factor the browse image data has virtually no visible distortions for the images tested.

Novik, Dmitry A.↗

Review - Axial compressor stall phenomena

Stall in compressors can be associated with the initiation of several types of fluid dynamic instabilities. These instabilities and the different phenomena, surge and rotating stall, which result from them, are discussed in this paper. Assessment is made of the various methods of predicting the onset of compressor and/or compression system instability, such as empirical correlations, linearized stability analyses, and numerical unsteady flow calculation procedures. Factors which affect the compressor stall point, in particular inlet flow distortion, are reviewed, and the techniques which are used to predict the loss in stall margin due to these factors are described. The influence of rotor casing treatment (grooves) on increasing compressor flow range is examined. Compressor and compression system behavior subsequent to the onset of stall is surveyed, with particular reference to the problem of engine recovery from a stalled condition. The distinction between surge and rotating stall is emphasized because of the very different consequences on recoverability. The structure of the compressor flow field during rotating stall is examined, and the prediction of compressor performance in rotating stall, including stall/unstall hysteresis, is described.

Greitzer, E. M.↗

Data volume reduction for single-look polarimetric imaging radar data

An 8-b quantization scheme to reduce the data volume for single-look complex scattering matrix data measured by polarimetric imaging radar systems is described. The scattering matrices are not symmetrized before compression, thereby retaining information about background noise and system effects. The data volume is reduced by a factor of 3.2. It is shown, with measured data, that the signal to quantization noise ratio for the compression scheme is more than 35 dB for the cross-polarized channels, and more than 45 dB for the copolarized channels.

Van Zyl, Jakob J.↗