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 163 records · Page 9

Performance Optimization for Relative-Error-Bounded Lossy Compression on Scientific Data

Scientific simulations in high-performance computing (HPC) environments generate vast volume of data, which may cause a severe I/O bottleneck at runtime and a huge burden on storage space for postanalysis. Unlike traditional data reduction schemes such as deduplication or lossless compression, not only can error-controlled lossy compression significantly reduce the data size but it also holds the promise to satisfy user demand on error control. Pointwise relative error bounds (i.e., compression errors depends on the data values) are widely used by many scientific applications with lossy compression since error control can adapt to the error bound in the dataset automatically. Pointwise relative-error-bounded compression is complicated and time consuming. In this article, we develop efficient precomputation-based mechanisms based on the SZ lossy compression framework. Our mechanisms can avoid costly logarithmic transformation and identify quantization factor values via a fast table lookup, greatly accelerating the relative-error-bounded compression with excellent compression ratios. In addition, we reduce traversing operations for Huffman decoding, significantly accelerating the decompression process in SZ. Experiments with eight well-known real-world scientific simulation datasets show that our solution can improve the compression and decompression rates (i.e., the speed) by about 40 and 80 p, respectively, in most of cases, making our designed lossy compression strategy the best-in-class solution in most cases.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

An iterative implicit diagonally-dominant factorization algorithm for solving the Navier-Stokes equations

Presented here is an algorithm for solving the multidimensional unsteady Navier-Stokes equations for compressible flows. It is based on a diagonally-dominant approximate factorization procedure. The factorization error and the timewise linearization error associated with this procedure are reduced by performing Newton-type inner iterations at each time step. The inviscid fluxes are evaluated by the fourth-order central differencing scheme amended with a numerical dissipation directly proportional to the entire dissipative part of the truncation error intrinsic to the third order biased upwind scheme. The important features of the proposed solution are elucidated by the numerical results of the convection of a vortex and the backward-facing step flows.

Chen, Shu-Cheng↗

Calculation of three-dimensional compressible turbulent boundary layers

From comparisons of high speed data with low speed closure procedures using variable mean density, there does not appear to be any appreciable influence of compressibility upon turbulent shear stress modeling in compressible turbulent boundary layers, even for extreme cases such as Mach 14 to 20 with a change in density across the layer of up to a factor of 100. Other evidence of apparent lack of compressibility caused new physics which may alter the shear stress for the compressible boundary layer cases including: (1) fluctuation Mach number was generally less than 1; (2) the shear stress distribution through the boundary layer was not a function of Mach number for zero pressure gradient flows; (3) the Morkovin hypothesis was valid up to Mach 5 (based on fluctuation data); (4) profile N power was not a function of Mach number, at least up to Mach 10; and (5) the nondimensional burst period was approximately the same as that for low speed.

Source record↗

Science-based Region-of-Interest Image Compression

As the number of currently active space missions increases, so does competition for Deep Space Network (DSN) resources. Even given unbounded DSN time, power and weight constraints onboard the spacecraft limit the maximum possible data transmission rate. These factors highlight a critical need for very effective data compression schemes. Images tend to be the most bandwidth-intensive data, so image compression methods are particularly valuable. In this paper, we describe a method for prioritizing regions in an image based on their scientific value. Using a wavelet compression method that can incorporate priority information, we ensure that the highest priority regions are transmitted with the highest fidelity.

Wagstaff, K. L.↗

Experimental observations and finite element analysis of the initiation of fiber microbuckling in notched composite laminates

An understanding was developed of the factors that determine the semi-circular edge-notched compressive strength and the associated failure mode(s) were identified of thermoplastic composite laminates with multidirectional stacking sequences. The experimental observations and the detailed literature review suggest at least four factors that affected the determination of the strain levels at which fiber microbuckling initiates and thus, partially control the composite's compression strength. The dependent variables studied are the compressive strength of a reduced gage section compression specimen and the compression strength of a compression specimen with two semi-circular edge notches (no opposite free edges) centered along the gage section. In this research, specimens containing two semi-circular edge notches (no opposite free edges) were loaded in compression at a relatively slow rate to provide more stable development of fiber microbuckling damage. The results indicate that the local constraints (free surfaces, supporting ply orientation, and resin-rich regions) significantly affect the strain level for the initiation of in-plane fiber microbuckling. Preliminary results at an elevated temperature, 77 C, showed the shear stress yield strength of the resin was reduced and consequently, the resistance to fiber microbuckling was also reduced. The finite element analysis of the perfectly straight fiber problem indicates that the free surface effect causes a 10 percent reduction in the critical buckling strain. However, the experimentally measured reduction for fibers with an initial fiber curvature, was 35 percent.

Guynn, E. Gail↗

Skeletal model reduction with forced optimally time dependent modes

Skeletal model reduction based on local sensitivity analysis of time dependent systems is presented in which sensitivities are modeled by forced optimally time dependent (f-OTD) modes. The f-OTD factorizes the sensitivity coefficient matrix into a compressed format as the product of two skinny matrices, i.e. f-OTD modes and f-OTD coefficients. The modes create a low-dimensional, time dependent, orthonormal basis which capture the directions of the phase space associated with most dominant sensitivities. These directions highlight the instantaneous active species, and reaction paths. Evolution equations for the f-OTD modes and coefficients are derived, and the implementation of f-OTD for skeletal reduction is described. For demonstration, skeletal reduction is conducted of the constant pressure ethylene-air burning in a zero-dimensional reactor, and new reduced models are generated. The laminar flame speed, the ignition delay, and the extinction curve as predicted by the models are compared against some existing skeletal models in literature for the same detailed model. Finally, the results demonstrate the capability of f-OTD to eliminate unimportant reactions and species in a systematic, efficient and accurate manner.

42 ENGINEERING↗

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗

Thermodynamics of correlated electrons in a magnetic field

The Hofstadter–Hubbard model captures the physics of strongly correlated electrons in an applied magnetic field, which is relevant to many recent experiments on Moiré materials. Few large-scale, numerically exact simulations exists for this model. In this work, we simulate the Hubbard–Hofstadter model using the determinant quantum Monte Carlo (DQMC) algorithm. We report the field and Hubbard interaction strength dependence of charge compressibility, fermion sign, local moment, magnetic structure factor, and specific heat. The gross structure of magnetic Bloch bands and band gaps determined by the non-interacting Hofstadter spectrum is preserved in the presence of U. Incompressible regions of the phase diagram have improved fermion sign. At half filling and intermediate and larger couplings, a strong orbital magnetic field delocalizes electrons and reduces the effect of Hubbard U on thermodynamic properties of the system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Key metrics of progress in the NIF ignition implosions and future challenges on the path to higher yields

Successful demonstration of an implosion that meets the physics metrics for ignition on the National Ignition Facility (NIF) in August 2021 and demonstration of target energy gain G > 1 in December 2022 represented the culmination of more than a decade of research on the NIF. An in-depth analysis of the key ignition metrics achieved on the NIF over the past decade shows sustained progress toward identifying and mitigating degradation effects. Ignition implosions on the NIF have progressed from early experiments whose compressed fuel conditions were degraded by a factor of 10–20 below the ignition threshold with yields of only 1%–2% of the laser energy, to current experiments that significantly exceed the ignition threshold and in the best experiments achieve yields greater than twice the laser energy. It is very likely that NIF can get more than 20 MJ of yield and possibly 50–70 MJ yields when fully optimized at an ultraviolet laser energy of 2.6 MJ following the completion of the NIF Extended Yield Capability Project. Achieving these higher yields will also require improvements to capsules at larger size and possibly with alternate materials, as well as an improved understanding of hohlraum and capsule physics. Now that NIF has achieved ignition, the design of a Next Generation High Energy Density Facility that could achieve the National Nuclear Security Administration long-term goal of substantially higher yields than the NIF can proceed with a more confident physics base. We evaluate these opportunities including possible applications for inertial fusion energy and likely challenges to achieving this level of performance over the next decade.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical study of laminar and turbulent flows inside a turnaround duct with and without guide vanes

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside a turnaround duct with and without guide vanes and to investigate the effects of vanes on the flow characteristics. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. Method of segmental sweeping is developed to overcome the multiple-zone problem due to imposition of guide vanes. In the present study, several test cases have been computed. These test cases include laminar and turbulent flows inside a turnaround duct without and with two or three guide vanes. The study reveals that: (1) there exists large recirculation zones inside the duct if no vanes are present; (2) properly shaped and positional guide vanes are effective in eliminating flow separation; and (3) laminar and turbulent flows have similar flow features. But turbulent flow has less total pressure drop.

Lin, S.-J.↗

A Navier-Stokes flow simulation of the Space Shuttle Main Engine Hot Gas Manifold

Incompressible viscous flow inside the turnaround duct, the fuel bowl, the transfer duct and the racetrack of the Space Shuttle Main Engine (SSME) Hot Gas Manifold (HGM) has been computed using the method of pseudo-compressibility together with an implicit, approximate-factorization algorithm. A multiple-zone method is used to make solution of flows in complex geometries easy. A model which predicts the pressure loading for the shield and the injector post arrangement without solving the complex flow field in the main injector region is proposed. The computed results show good qualitative agreement with experimental data.

Yang, Ruey-Jen↗

Numerical simulation of flow path in the oxidizer side hot gas manifold of the Space Shuttle main engine

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside the oxidizer side Hot Gas Manifold of the Space Shuttle Main Engine. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is developed to overcome the complexity of the geometry. In the present study, the laminar and turbulent flows in the oxidizer side Hot Gas Manifold have been computed. The study reveals that: (1) there exists large recirculation zones inside the bowl if no vanes are present; (2) strong secondary flows are observed in the transfer tube; and (3) properly shaped and positioned guide vanes are effective in eliminating flow separation.

Lin, S. J.↗

Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks with Gyrotropic Pressure

Shock jump conditions, i.e., the specification of the downstream parameters of the gas in terms of the upstream parameters, are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. This is done by combining the continuity of particle number flux and the electromagnetic boundary conditions at the shock with the magnetohydrodynamic conservation laws derived from the stress-energy tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. Examples assuming isotropic pressure illustrate how the shock compression ratio depends on the shock speed and obliquity. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically with- out additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., approximately 20%) produce significant changes in the shock compression ratio, r , at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks; such models include the modification of the fluid velocity profile due to the contribution of energetic particles to the momentum and energy fluxes.

Double, Glen P.↗

Current concepts in the pathophysiology, evaluation, and diagnosis of compartment syndrome

This article reviews present knowledge of the pathophysiology and diagnosis of acute compartment syndromes. Recent results using compression of legs in normal volunteers provide objective data concerning local pressure thresholds for neuromuscular dysfunction in the anterior compartment. Results with this model indicate that a progression of neuromuscular deficits occurs when IMP increases to within 35 to 40 mm Hg of diastolic blood pressure. These findings provide useful information on the diagnosis and compression thresholds for acute compartment syndromes. Time factors are also important, however, and usually are incompletely known in most cases of acute compartment syndrome. Although the slit catheter is a very good technique for monitoring IMP during rest, these catheters and their associated extracorporeal transducer systems are not ideal. Recently developed miniature transducer-tipped catheters and, perhaps, future development of noninvasive techniques may provide accurate recordings of IMP in patients with acute compartment syndromes.

NASA Discipline Cardiopulmonary↗

Mechanics of Preloaded Bolt Tensile Loading With Focus on Load Introduction Factor

The bolt tensile and joint separation loads are directly influenced by the locations at which the external loads enter the clamped members of a preloaded bolted joint (PBJ) and the associated load-paths through the joint. This physical load introduction mechanism affecting the bolt tensile loading is typically represented in the bolt tensile load equation, in part, by a load introduction factor (LIF), which was shown by H.M. Lee of Marshall Spaceflight Center to be a natural product of the bolt tensile load equation using a linear spring stiffness model. This LIF, being a function of load-path stiffness, has subsequently been denoted as the stiffness-based LIF (SBLIF), providing a framework to calculate the LIF using whatever load-path stiffness approximations are appropriate. Expanding upon the work of Lee, it is shown that the SBLIF and the joint stiffness factor are functions of the stiffnesses of the same load-paths and regions within a PBJ, and thus they should not be treated as independent variables. Mathematical expressions for the SBLIF are presented. Comparisons are shown between the analytically calculated SBLIF, the analytically calculated geometric LIF (GLIF), which is a simple clamped-member thickness ratio, the experimentally derived LIF, and the LIF determined by finite element analysis (FEA). Using experiment and FEA as a benchmark, the SBLIF, using traditional load-path stiffness approximations, enables a more accurate prediction of bolt tensile loading than the GLIF, although it can be unconservative near joint separation. The GLIF generally attributes more of the externally applied tensile load to the bolt than does the SBLIF, potentially resulting in heavier and/or more costly bolted joints. Mathematical relationships between the SBLIF and the GLIF are developed. Supplemental material is provided in the appendixes where the historical practice of using the joint compressive stiffness in place of the joint tensile stiffness is evaluated. The appendixes include step-by-step examples demonstrating the calculation of the SBLIF using traditional stiffness approximations and conclude with the development of the joint diagram in terms of the SBLIF, culminating into formulas for the key features of a joint diagram, which is useful for programming.

Load Path↗

Rapid-Rate Compression Testing of Sheet Materials at High Temperatures

This Report describes the test equipment that was developed and the procedures that were used to evaluate structural sheet-material compression properties at preselected constant strain rates and/or loads. Electrical self-resistance was used to achieve a rapid heating rate of 200 F/sec. Four materials were tested at maximum temperatures which ranged from 600 F for the aluminum alloy to 2000 F for the Ni-Cr-Co iron-base alloy. Tests at 0.1, 0.001, and 0.00001 in./in./sec showed that strain rate has a major effect on the measured strength, especially at the high temperatures. The tests, under conditions of constant temperature and constant compression stress, showed that creep deformation can be a critical factor even when the time involved is on the order of a few seconds or less. The theoretical and practical aspects of rapid-rate compression testing are presented, and suggestions are made regarding possible modifications of the equipment which would improve the over-all capabilities.

Bernett, E. C.↗

Turbulence models and Reynolds analogy for two-dimensional supersonic compression ramp flow

Results of the application of turbulence models and the Reynolds analogy to the Navier-Stokes computations of Mach 2.9 two-dimensional compression ramp flows are presented. The Baldwin-Lomax eddy viscosity model and the kappa-epsilon turbulence transport equations for the turbulent momentum flux modeling in the Navier-Stokes equations are studied. The Reynolds analogy for the turbulent heat flux modeling in the energy equation was also studied. The Navier-Stokes equations and the energy equation were numerically solved for the flow properties. The Reynolds shear stress, the skin friction factor, and the surface heat transfer rate were calculated and compared with their measurements. It was concluded that with a hybrid kappa-epsilon turbulence model for turbulence modeling, the present computations predicted the skin friction factors of the 8 deg and 16 deg compression ramp flows and with the turbulent Prandtl number Pr(sub t) = 0.93 and the ratio of the turbulent thermal and momentum transport coefficients mu(sub q)/mu(sub t) = 2/Prt, the present computations also predicted the surface heat transfer rates beneath the boundary layer flow of the 16 compression ramp.

Wang, Chi R.↗

Renormalization of states and quasiparticles in many-body downfolding

We explore the principles of many-body Hamiltonian complexity reduction via downfolding on an effective low-dimensional representation. We show that the renormalization factor provides a unique measure of the quality of the compression as it directly represents the projection between the approximate stationary state of the many-body Hamiltonian and the full many-body wavefunction. Hence, the renormalization factor is a measure of fidelity between the effective (reduced-rank) description and the full many-body treatment for arbitrary (i.e., ground and excited) states. When the entire problem is mapped on a system of interacting quasiparticles [Romanova et al., npj Comput. Mater. 9, 126 (2023)], the effective Hamiltonians can faithfully reproduce the physics only when a clear energy scale separation exists between the subsystems and their environment. We also demonstrate that it is necessary to include quasiparticle renormalization at distinct energy scales, capturing the distinct interaction between subsystems and their surrounding environments. Numerical results from simple, exactly solvable models highlight the limitations and strengths of this approach, particularly for ground and low-lying excited states. This work lays the groundwork for applying dynamical downfolding techniques to problems concerned with (quantum) interfaces.

Green-functions technique↗