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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Spectral methods for the Euler equations

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collocation methods are put in historical context. The basic concepts of both Fourier and Chebyshev spectral collocation methods are provided. Filtering strategies for both shock-fitting and shock-capturing approaches are also presented. Fourier shock capturing techniques are evaluated using a one-dimensional, periodic astrophysical 'nozzle' problem. Examples of shock-fitting approaches include a shock/acoustic wave interaction, shock/vortex interaction, and the classical blunt body problem. While the shock capturing spectral method does not yet show a clear advantage over second-order finite differences, equivalent accuracy can be obtained using shock fitting with far fewer grid points.

Hussaini, M. Y.↗
DOE OSTI2023

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The simulation of 3D hypervelocity spallation using a hydrocode PAGOSA with FLIP+MPM

n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.

42 ENGINEERING↗
NASA NTRS1989

Review of image correlation systems - Hybrid and optical

The operating principles, design, and development histories of three image correlation systems, the Portable Correlator Unit (PCU), the Miniature Correlator Unit (MCU), and the Coherence Interferometer (CI) are reviewed. The PCU and MCU are compact versions of the classical Vander Lught filter optical correlator, essentially parallel optical processors by light valve; the CI is a white-light correlator which functions entirely by interfering light waves, with electronic detection only at the final output, and thus represents a parallel optical processor by rotational-shear interferometer. Details of the design and construction are discussed and illustrated with drawings and photographs. Also considered are applications of the MCU to NASA tasks such as rendezvous and docking on the Mars Rover and Sample Return mission, autonomous robotic exploration of planetary surfaces, and compression and pattern-recognition of remote-sensing data.

Breckinridge, James B.↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Finite-volume and integral-equation techniques for transonic and supersonic vortex-dominated flows

Two computational techniques are developed to calculate the compressible vortex-dominated flows. The first technique is a finite-volume Euler Solver which uses four-Stage Runge-Kutta time stepping with second- and fourth-order dissipation terms. The technique is applied to supersonic conical and three-dimensional flows about sharp- and round-edged delta wings. Attached and separated-flow solutions have been obtained depending on the values of damping coefficients. The second technique is an integral-equation solver of the full potential equation which uses a volume-integral term in addition to the classical surface-integral terms. The technique is applied to transonic three-dimensional flows about sharp-edged delta wings. A hybrid technique which combines the finite-volume and the integral-equation solvers is also presented.

Kandil, O. A.↗
DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Simulating a pulsed-power-driven plasma with ideal MHD

We describe a simple practical numerical method for simulating plasma driven within a vacuum chamber by a pulsed power generator. Typically, in this type of simulation, the vacuum region adjacent to the plasma is approximated as a highly resistive, light fluid; this involves computationally expensive solvers describing the diffusion of the magnetic field through this fluid. Instead, we provide a recipe for coupling pulsed power generators to the magnetohydrodynamics (MHD) domain by approximating the perfectly insulating vacuum as a light, perfectly conducting, inviscid MHD fluid and discuss the applicability of this counter-intuitive technique. This much more affordable ideal MHD representation is particularly useful in situations where a plasma exhibits interesting three-dimensional phenomena, either due to the design of the experiment or due to developing instabilities. We verified that this coupling recipe works by modeling an exactly solvable flux compression generator as well as a self-similar Noh-like solution and demonstrated convergence to the theoretical solution. We also showed examples of simulating complex three-dimensional pulsed power devices with this technique. We release our code implementation to the public.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shock propagation in aerogel and TPP foams for inertial fusion energy target design

Achieving practical inertial fusion energy (IFE) requires the development of target designs with well-characterized microstructure and compression response. We measured shock dynamics in low-density (17.5–500 mg/cm 3 ) aerogel and two-photon polymerization (TPP) foams using x-ray phase contrast imaging (XPCI) methods and the Velocity Interferometer System for Any Reflector. By analyzing shock front evolution, we examined how target type and density influence shock propagation and energy dissipation. Talbot-XPCI shows that aerogels support a smooth, bowed shock front due to their homogeneous nanometer-scale pore network. In contrast, TPP foams exhibit irregular, stepwise propagation driven by interactions with their periodic micrometer-scale lattice. Shock velocity follows a power-law relation: aerogels deviate from classical ρ −1/2 scaling due to pore-collapse dissipation, while TPP foams follow the trend with larger uncertainties from density variations. Comparisons with xRAGE simulations reveal systematic underestimation of shock speeds. These results provide the first experimental constraints on shock propagation in TPP foams over a wide density range and highlight the influence of internal structure on anisotropic shock behavior. Our findings support improved benchmarking of EOS and hydrodynamic models and inform the design of foam architectures that promote implosion symmetry in IFE capsules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
NASA NTRS1994

Idealized textile composites for experimental/analytical correlation

Textile composites are fiber reinforced materials produced by weaving, braiding, knitting, or stitching. These materials offer possible reductions in manufacturing costs compared to conventional laminated composites. Thus, they are attractive candidate materials for aircraft structures. To date, numerous experimental studies have been performed to characterize the mechanical performance of specific textile architectures. Since many materials and architectures are of interest, there is a need for analytical models to predict the mechanical properties of a specific textile composite material. Models of varying sophistication have been proposed based on mechanics of materials, classical laminated plate theory, and the finite element method. These modeling approaches assume an idealized textile architecture and generally consider a single unit cell. Due to randomness of the textile architectures produced using conventional processing techniques, experimental data obtained has been of limited use for verifying the accuracy of these analytical approaches. This research is focused on fabricating woven textile composites with highly aligned and accurately placed fiber tows that closely represent the idealized architectures assumed in analytical models. These idealized textile composites have been fabricated with three types of layer nesting configurations: stacked, diagonal, and split-span. Compression testing results have identified strength variations as a function of nesting. Moire interferometry experiments are being used to determine localized deformations for detailed correlation with model predictions.

Adams, Daniel O.↗
DOE OSTI2020

Equations of State in Computational Physics: What, Why, and How

Summary: 1. Equations of state furnish thermodynamic properties of a material. 2. Usually involve P, p (or V), and T, or P, r, and E 3. Necessary ingredient to solve compressible flow problems. 4. Varieties: i. Analytical expressions (ideal gas, Grüneisen) – fast, convenient for verification and test problems ii. Graphs and charts – good for hand calculations iii. Tables and databases – generally most accurate, require interpolation i. Steam tables ii. NIST database iii. SESAME EOS Library 5. Constructed by using semi-empirical models from Statistical Mechanics and fitting to available data

36 MATERIALS SCIENCE↗
NASA NTRS1980

Thermal expansion of composites using Moire interferometry

An experimental technique for precise measurement of the thermal response of fiber-reinforced composite materials uses moire interferometry with fringe multiplication which yield a sensitivity of 833 nm (32.8 mu in.) per fringe. Results from the technique are compared with those obtained from electrical resistance strain gages, and also those predicted from classical lamination theory. Temperature dependent coefficients of thermal expansion for composite materials subjected to thermal cycling in the temperature range of 297 K (75 F) to 422 K (300 F) were determined for four laminate configurations (0, 90, 0/ + or - 45/90 sub s and 0/90/ + or - 45 sub s) of T300/5208 graphite epoxy, and ranged from -0.107 mu epsilon K/1 (-0.059 mu epsilon deg F/-) for the 0 laminate to 32.18 mu epsilon K/1 (17.88 mu epsilon F/1) for the 90 laminate. Moisture was found to greatly influence the thermal response of a quasi-isotropic laminate, resulting in hysteresis and residual compressive strain as the moisture content was reduced. Comparisons between moire and strain gage measurements were inconclusive with both techniques giving consistent but systematically different results. Differences of as much as 29% were observed.

Bowles, D. E.↗
DOE OSTI2020

Automated Fitting of a Semi-empirical Multiphase Equation of State for Carbon

The equation of state (EOS) of carbon is important in high explosive, geophysical, and inertial confinement fusion applications. Within the semi-empirical Sesame framework, the EOS of each phase is represented by a sum of cold, vibrational, and thermal electronic Helmholtz free energy contributions. Each phase has ~5-10 independent parameters that are adjusted to reproduce single-phase data (e.g., thermal expansion, isothermal compression) as well as experimentally- and computationally- derived phase boundaries. Manual calibration of the full multiphase EOS is arduous. We present our progress in development of automated EOS parameterization based on minimization of an objective function. Here, this function encodes deviation of model EOS results from experimental/computational benchmarks. Optimization is implemented as a combination of global (particle swarm) and local optimization techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Scale Invariance of the Homentropic Inviscid Euler Equations with Application to the Noh Problem

We investigate the inviscid compressible flow (Euler) equations constrained by an ”isentropic” equation of state (EOS), whose functional form in pressure is an arbitrary function of density alone. Under the aforementioned condition, we interrogate using symmetry methods the scale-invariance of the homentropic inviscid Euler equations. We find that under general conditions, we can reduce the inviscid Euler equations into a system of two coupled ordinary differential equations. To exemplify the utility of these results, we formulate two example scale-invariant, self-similar solutions. The first example includes a shock-free expanding bubble scenario, featuring a modified Tait EOS. The second example features the classical Noh problem, coupled to an arbitrary isentropic EOS. In this case, in order to satisfy the conditions set forth in the classical Noh problem, we find that the solution for the flow is given by a transcendental algebraic equation in the shocked density.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

A New Approach to Fibrous Composite Laminate Strength Prediction

A method of predicting the strength of cross-plied fibrous composite laminates is based on expressing the classical maximum-shear-stress failure criterion for ductile metals in terms of strains. Starting with such a formulation for classical isotropic materials, the derivation is extended to orthotropic materials having a longitudinal axis of symmetry, to represent the fibers in a unidirectional composite lamina. The only modification needed to represent those same fibers with properties normalized to the lamina rather than fiber is a change in axial modulus. A mirror image is added to the strain-based lamina failure criterion for fiber-dominated failures to reflect the cutoffs due to the presence of orthogonal fibers. It is found that the combined failure envelope is now identical with the well-known maximum-strain failure model in the tension-tension and compression-compression quadrants but is truncated in the shear quadrants. The successive application of this simple failure model for fibers in the 0/90 degree and +/- 45 degree orientations, in turn, is shown to be the necessary and sufficient characterization of the fiber-dominated failures of laminates made from fibers having the same tensile and compressive strengths. When one such strength is greater than the other, the failure envelope is appropriately truncated for the lesser direct strain. The shear-failure cutoffs are now based on the higher axial strain to failure since they occur at lower strains than and are usually not affected by such mechanisms as microbuckling. Premature matrix failures can also be covered by appropriately truncating the fiber failure envelope. Matrix failures are excluded from consideration for conventional fiber/polymer composites but the additional features needed for a more rigorous analysis of exotic materials are covered. The new failure envelope is compared with published biaxial test data. The theory is developed for unnotched laminates but is easily shrunk to incorporate reductions to allow for bolt holes, cutouts, reduced compressive strength after impact, and the like.

Hart-Smith, L. J.↗
NASA NTRS2007

Nonstandard Analysis and Shock Wave Jump Conditions in a One-Dimensional Compressible Gas

Nonstandard analysis is a relatively new area of mathematics in which infinitesimal numbers can be defined and manipulated rigorously like real numbers. This report presents a fairly comprehensive tutorial on nonstandard analysis for physicists and engineers with many examples applicable to generalized functions. To demonstrate the power of the subject, the problem of shock wave jump conditions is studied for a one-dimensional compressible gas. It is assumed that the shock thickness occurs on an infinitesimal interval and the jump functions in the thermodynamic and fluid dynamic parameters occur smoothly across this interval. To use conservations laws, smooth pre-distributions of the Dirac delta measure are applied whose supports are contained within the shock thickness. Furthermore, smooth pre-distributions of the Heaviside function are applied which vary from zero to one across the shock wave. It is shown that if the equations of motion are expressed in nonconservative form then the relationships between the jump functions for the flow parameters may be found unambiguously. The analysis yields the classical Rankine-Hugoniot jump conditions for an inviscid shock wave. Moreover, non-monotonic entropy jump conditions are obtained for both inviscid and viscous flows. The report shows that products of generalized functions may be defined consistently using nonstandard analysis; however, physically meaningful products of generalized functions must be determined from the physics of the problem and not the mathematical form of the governing equations.

Baty, Roy S.↗