Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “classical compression”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A description of structured waves in shock compressed particulate composites

Dynamic compression of composite materials is of scientific interest because the mechanical mismatch between internal phases challenges continuum theories. Typical assumptions about steady wave propagation and quasi-instantaneous state changes require reexamination along with the need for time-dependent models. To that end, data and models are presented here for the shock compression of an idealized particulate composite. To serve as a generic representative of this material class, a polymer matrix was filled with tungsten particles, ranging from 1 to 50 vol. %. This creates a simple microstructure containing randomly scattered particles with an extreme impedance mismatch to the binding matrix. These materials were parallel plate impact loaded by Al flyers traveling at 1.8–5.0 km/s. Velocimetry provided records of the equilibrium state and the compression wave structure for each case with trends quantified by an empirical fit. The same quantities were also studied as a function of the wave's propagation distance. A homogenized viscoelastic model then made it possible to progress from cataloging trends to making predictions. Starting from a Mie–Gruneisen equation of state, additional time varying terms were added to capture the transient response. After calibration, accurate predictions of the steady wave structure were possible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

Analytical characterization of AOTV perigee aerothermodynamic regime

In preliminary design studies of AOTV (Aeroassisted Orbital Transfer Vehicle) configurations, it is desirable to have a simple analytical method for characterizing the aerothermodynamic regime at skip trajectory perigee as a function of vehicle L/D and m/CDA (ballistic coefficient). The present study derives an approximate perigee solution from the generalized equations of motion. For a prescribed entry velocity vector, the solution determines a Zeta function (proportional to density divided by m/CDA) at perigee as a parametric function of the exit velocity vector and (constant) L/D. The perigee Zeta function then determines perigee density (or altitude) as a parametric function of m/CDA. The solution allows the following classic aerothermodynamic parameters to be determined at perigee as parametric functions of m/CDA, L/D, and the exit velocity vector: Reynolds Number (viscous effects), Mach Number (compressibility effects), Knudsen Number (rarefaction effects), sphere and disk bow shock standoff distance, Damkohler Number (relaxation effects), viscous correlation parameter (viscid-inviscid interactions), and Stanton Number (convective heat transfer). Results of the analysis are given for low L/D and mid L/D AOTV configurations on return from geosynchronous and L5 orbits. It is concluded the method successfully provides preliminary estimates of the aerothermodynamic parameters through the use of simple algebraic equations and plots.

Desautel, D.↗
DOE OSTI2023

Quantum JPEG

The joint photographic expert group algorithm compresses a digital image by filtering its high spatial-frequency components. Similarly, we introduce a quantum algorithm that uses the quantum Fourier transform to discard the high spatial-frequency qubits of an image, downsampling it to a lower resolution. This allows one to capture, compress, and send images even with limited quantum resources for storage and communication. We show under which conditions this protocol is advantageous with respect to its classical counterpart.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Fluid-structure finite-element vibrational analysis

A fluid finite element has been developed for a quasi-compressible fluid. Both kinetic and potential energy are expressed as functions of nodal displacements. Thus, the formulation is similar to that used for structural elements, with the only differences being that the fluid can possess gravitational potential, and the constitutive equations for fluid contain no shear coefficients. Using this approach, structural and fluid elements can be used interchangeably in existing efficient sparse-matrix structural computer programs such as SPAR. The theoretical development of the element formulations and the relationships of the local and global coordinates are shown. Solutions of fluid slosh, liquid compressibility, and coupled fluid-shell oscillation problems which were completed using a temporary digital computer program are shown. The frequency correlation of the solutions with classical theory is excellent.

Feng, G. C.↗
NASA NTRS1994

Structural design criteria for filament-wound composite shells

Advanced composite cylinders, manufactured by filament winding, provide a cost effective solution to many present structural applications; however, the compressive performance of filament-wound cylinders is lower than comparable shells fabricated from unidirectional tape. The objective of this study was to determine the cause of this reduction in thin filament-wound cylinders by relating the manufacturing procedures to the quality of the cylinder and to its compressive performance. The experiments on cylinder buckling were complemented by eigenvalue buckling analysis using a detailed geometric model in a finite element analysis. The applicability of classical buckling analyses was also investigated as a design tool.

Hahn, H. T.↗
DOE OSTI2022

Measurement of shock roughness due to phase plate speckle imprinting relevant for x-ray diffraction experiments on 3rd and 4th generation light sources

Laser-shock compression experiments at 3rd and 4th generation light sources generally employ phase plates, which are inserted into the beamline to achieve a repeatable intensity distribution at the focal plane. Here, the laser intensity profile is characterized by a high-contrast, high-frequency laser speckle. Without sufficient smoothing, these laser non-uniformities can translate to a significant pressure distribution within the sample layer and can affect data interpretation in x-ray diffraction experiments. Here, we use a combination of one- and two-dimensional velocity interferometry to directly measure the extent to which spatial frequencies within the laser focal spot intensity pattern are smoothed out during propagation within the laser plasma and a polyimide ablator. We find that the use of thicker polyimide layers results in spatially smoother shock fronts, with the greatest degree of smoothing associated with the highest spatial frequencies. Focal spots with the smallest initial speckle separation produce the most rapid smoothing. Laser systems that employ smoothing by spectral dispersion techniques to rapidly modulate the focal plane intensity distribution are shown to be the most effective ones in producing a spatially smooth shock front. We show that a simple transport model combined with the known polyimide Hugoniot adequately describes the extent of shock smoothness as a function of polyimide thickness. Finally, our results provide a description of spatial structure smoothing across a shock front, which can be used to design targets on x-ray free electron laser facilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantum Monte Carlo determination of the principal Hugoniot of deuterium

Here, we present coupled electron-ion Monte Carlo results for the principal Hugoniot of deuterium together with an accurate study of the initial reference state of shock-wave experiments. We discuss the influence of nuclear quantum effects, thermal electronic excitations, and the convergence of the potential energy surface by wave-function optimization within variational Monte Carlo and projection quantum Monte Carlo methods. Compared to a previous study, our calculations also include low pressure-temperature (P,T) conditions resulting in close agreement with experimental data, while our revised results at higher (P,T) conditions still predict a more compressible Hugoniot than experimentally observed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Novel method to control explosive shock sensitivity: A mesoscale study to understand the effect of thermally expandable microsphere (TEM) inclusions in high explosives (HE) microstructure

When the void content and/or void structure of a high explosive (HE) is altered by some means (i.e., bulk heating or mechanical damage), the shock initiation behavior of the material changes. The ability to precisely predict the change in shock sensitivity after an HE has undergone microstructural changes is a crucial capability in multi-scale reactive flow models. Here, we utilize thermally expandable microspheres (TEMs) as a dopant in a polymer bonded explosive (PBX) matrix to alter the shock initiation properties in a controlled fashion. Here, using a mesoscale modeling approach, we evaluated how a single TEM (before and after thermal expansion) behaves under shock compression, as well as how the matrix PBX in the direct vicinity of the TEM is affected. We first examined the effect of an unexpanded TEM in the explosive matrix and found that its presence does not significantly perturb the bulk flow and by extension will not affect bulk sensitivity. Next, we examined the effect of an expanded TEM and found that its presence significantly perturbs the flow via hydrodynamic jetting, which causes a secondary shock wave with a strength that exceeds that of the incident wave. Finally, we showed that this secondary shock interacts with the downstream porosity to ignite a larger fraction of the overall pore volume, commensurate with the secondary shock strength and the affected volume, increasing the global (bulk) shock sensitivity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

HTR-1.3 solver: Predicting electrified combustion using the hypersonic task-based research solver

Here this manuscript presents an updated open-source version of the Hypersonics Task-based Research (HTR) solver. The solver, whose main features are presented in Di Renzo et al. (2020) and Di Renzo & Pirozzoli (2021), is designed for direct numerical simulation of reacting flows at high Reynolds numbers. This new version extends the applications of the HTR solver to turbulent combustion in the presence of external electric fields. In particular, a new distributed Poisson solver compatible with heterogeneous architectures has been incorporated in the algorithm to compute the electric potential distribution in bi-periodic configurations. The drift fluxes of the electrically charged species are now included in the transport equations using a targeted essentially non-oscillatory scheme. A verification of these new features of the solver is provided using one-dimensional burner stabilized flames, whereas a three dimensional turbulent flame is utilized to discuss the scalability of the proposed numerical tool.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Hybrid algorithm for the time-dependent Hartree–Fock method using the Yang–Baxter equation on quantum computers *

Abstract The time-dependent Hartree–Fock (TDHF) method is an approach to simulate the mean field dynamics of electrons within the assumption that the electrons move independently in their self-consistent average field and within the space of single Slater determinants. One of the major advantages of performing time dynamics within Hartree–Fock theory is the free fermionic nature of the problem, which makes TDHF classically simulatable in polynomial time. Here, we present a hybrid TDHF implementation for quantum computers. This quantum circuit grows with time; but with our recent work on circuit compression via the Yang–Baxter equation (YBE), the resulting circuit is constant depth. This study provides a new way to simulate TDHF with the aid of a quantum device as well as provides a new direction for the application of YBE symmetry in quantum chemistry simulations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Rotational relaxation of molecular H2 at moderate temperatures

A coupled rotation-vibration-dissociation model is used to compute the rotational relaxation times for molecular hydrogen as a function of final temperatures (500-5000 K) in a hypothetical scenario of sudden compression. The theoretical model is based on a master equation solver. The bound-bound and bound-free transition rates have been computed using a quasi-classical trajectory method. A review of the available experimental data on the rotational relaxation of hydrogen is presented, with a critical overview of the method of measurements and data reduction, including the sources of errors. The experimental data are compared with the computed results.

Sharma, S. P.↗
DOE OSTI2022

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Understanding the evolution of liquid and solid microjets from grooved Sn and Cu samples using radiography

Experiments were performed on grooved Sn and Cu samples to study the temporal evolution of microjets. Jets were generated by the impact of gun-launched flyer plates against the back of grooved targets made from either Cu or Sn (groove depth of ~250 μm). The Hugoniot states in the various Sn targets encompassed conditions where solid phases are maintained throughout (7 and 16 GPa) and also conditions where melting occurs upon the release of compression (25 and 34 GPa); the transition occurs near a Hugoniot pressure of 23 GPa. Cu targets at 27 and 56 GPa provide comparisons in which the jets move at similar speeds but remain solid. In all cases, the spatial distribution of mass within the microjets was measured using high-speed synchrotron radiography. The result is a time history of the jet thickness profile from which quantities like total jet mass and jet velocity can be derived. In both the solid and liquid states, we generally observe that an increase in the shock strength leads to an increase in jet mass. However, this trend breaks down for Hugoniot states near the transition from continuously solid to melted-on-release. This is evidenced by the observation that there was no difference in the rate of mass flow in Sn jets at 16 and 25 GPa, while similar pressure jumps on either side of this range caused substantial changes in the jet mass. This contrasts with the behavior of smaller polishing defects that were present on the same samples (~1 μm deep). From these, no ejecta mass was detected below the melt boundary, but obvious microjets were generated once melting occurred. This indicates that crossing the bulk melt-on-release threshold can alternately promote or inhibit the flow of mass into microjets based on the amplitude of the initial perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

Numerical analysis and parametric studies of the buckling of composite orthotropic compression and shear panels

A computer program is presented which was developed for the combined compression and shear of stiffened variable thickness orthotropic composite panels on discrete springs: boundary conditions are general and include elastic boundary restraints. Buckling solutions are obtained by using a newly developed trigonometric finite difference procedure which improves the solution convergence rate over conventional finite difference methods. The classical general shear buckling results which exist only for simply supported panels over a limited range of orthotropic properties, were extended to the complete range of these properties for simply supported panels and, in addition, to the complete range of orthotropic properties for clamped panels. The program was also applied to parametric studies which examine the effect of filament orientation upon the buckling of graphite-epoxy panels. These studies included an examination of the filament orientations which yield maximum shear or compressive buckling strength for panels having all four edges simply supported or clamped over a wide range of aspect ratios. Panels with such orientations had higher buckling loads than comparable, equal weight, thin skinned aluminum panels. Also included among the parameter studies were examinations of combined axial compression and shear buckling and examinations of panels with rotational elastic edge restraints.

Housner, J. M.↗
DOE OSTI2026

Mechanistic Insights into Defect-Mediated Crystallization Revealed by Lattice Strain Evolution

Structural defects and lattice strain are intrinsic to many crystalline materials, yet their roles in controlling chemical reaction mechanisms and directing crystallization pathways remain poorly understood. Here, in this study, we revealed the three-dimensional evolution of strain and dislocation defects at the nanoscale during the growth of heterogeneously nucleated barite (BaSO 4 ) and calcite (CaCO 3 ) crystals by using coherent X-ray scattering, electron microscopy, and molecular simulations. Unlike barite, which formed with minimal internal strain, calcite developed dislocation defects and exhibited spatially varying strain that increased during growth. During growth in Sr-rich solutions, calcite likely incorporates Sr 2+ into the defects, which further modulates the local lattice structure and increases both the compressive and tensile strain. These findings suggest that calcite crystallization was likely dominated by attachment of precursor phases, which gave rise to defect-enriched domain structures not predicted by classical growth models. By linking defect formation to ion incorporation and growth dynamics, this work provides fundamental insight into how lattice-level strain heterogeneity governs the chemical reactivity of ionic crystals.

Bragg coherent diffractive imaging↗
DOE OSTI2022

Direct cosmological inference from three-dimensional correlations of the Lyman α forest

ABSTRACT When performing cosmological inference, standard analyses of the Lyman α (Ly α) three-dimensional correlation functions only consider the information carried by the distinct peak produced by baryon acoustic oscillations (BAO). In this work, we address whether this compression is sufficient to capture all the relevant cosmological information carried by these functions. We do this by performing a direct fit to the full shape, including all physical scales without compression, of synthetic Ly α autocorrelation functions and cross-correlations with quasars at effective redshift zeff = 2.3, assuming a DESI-like survey, and providing a comparison to the classic method applied to the same data set. Our approach leads to a $3.5{{\ \rm per\ cent}}$ constraint on the matter density ΩM, which is about three to four times better than what BAO alone can probe. The growth term fσ8(zeff) is constrained to the $10{{\ \rm per\ cent}}$ level, and the spectral index ns to $\sim 3-4{{\ \rm per\ cent}}$. We demonstrate that the extra information resulting from our ‘direct fit’ approach, except for the ns constraint, can be traced back to the Alcock–Paczynski effect and redshift space distortion information.

79 ASTRONOMY AND ASTROPHYSICS↗