Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock deformation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing. In order for large scale optimization to become routine, the benefits of parallel architectures should be exploited. Although the flow solver has been parallelized using compiler directives. The parallel efficiency is under 50 percent. Clearly, parallel versions of the codes will have an immediate impact on the ability to design realistic configurations on fine meshes, and this effort is currently underway.

Grossman, Bernard↗
NASA NTRS1999

Practical Aerodynamic Design Optimization Based on the Navier-Stokes Equations and a Discrete Adjoint Method

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Grossman, Bernard↗
NASA NTRS1994

Global, quasi-steady dynamics of the distant solar wind 2: Deformation of the heliospheric current sheet

A three-dimensional (3-D) magnetohydrodynamic (MHD) numerical model is used to trace global deformations of the heliospheric current sheet (HCS) caused by large-scale dynamical interactions associated with corotating solar wind flows. Configurations incorporating the tilted-dipole geometry are investigated out past 30 AU for a variety of dipole tilt angles, alpha. Inclusion of the full, 3-D interplanetary dynamics allows north-south displacements and the east-west warping of the HCS by advective corotational effects to be accurately assessed for the first time. It is found that large-scale spatial correlations between velocity and density imposed at the coronal source (i.e., the geometric arrangement whereby regions of high-speed low-density material lying adjacent to areas of slow, dense flow interact obliquely under the influence of solar rotation) result in a distinctive pattern of deformation of the HCS. For an alpha = 30 deg tilted-dipole example, it is shown that typical zonal variations in radial velocity lead to significant folding of the HCS within about 5 AU of the Sun. By 10 AU, additional sharp bends appear near the latitudinal extremes of the HCS surface, where it is overtaken by shock fronts driven by 3-D corotating interaction regions (CIRs). Moreover, the model suggests that inside about 20 AU, major plasma structures are systematically organized, about the HCS, such that the greatest concentrations of material and magnetic field (the centroids of the 3-D CIR structures) are coincident with the folded crests of the HCS (near heliographic latitudes lambda = +/- alpha). Thus, in these circumstances many of the more interesting dynamical features inconveniently lie well away from the heliographic equator. At larger heliocentric distances, where neighboring CIRs begin to interact strongly, the warping of the HCS abates dramatically and the association between folds in the HCS and major field and density concentrations is weakened and ultimately breaks down. On the basis of the model calculations, quasi-steady deformations and related phenomena should figure prominently in the interplanetary medium, except during those periods when the magnetic dipole and solar spin axes are nearly aligned (alpha is approximately less than 10 deg). Even then, they may be of significance if there are any substantial local warps or kinks in the streamer belt.

Pizzo, V. J.↗
NASA NTRS1991

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗
DOE OSTI2023

Comparing theoretical predictions of radiation-free velocities of edge dislocations to molecular dynamics simulations

Transonic defect motion is of interest for high strain-rate plastic deformation as well as for crack propagation. Ever since Eshelby's 1949 prediction [J. D. Eshelby, Proc. Phys. Soc. A 62, 307 (1949)] in the isotropic limit of a “radiation-free” transonic velocity $v$ RF = $\sqrt{2}$c T , where shock waves are absent, there has been speculation about the significance of radiation-free velocities (if they truly exist) for defect mobility. Here, in this work, we argue that they do not play any significant role in dislocation dynamics in metals, based on comparing theoretical predictions of radiation-free velocities for transonic edge dislocations with molecular dynamics simulations for two face-centered cubic metals: Ag, where theory predicts radiation-free states, and Cu, where it does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1986

Development of a numerical procedure to map a general 3-d body onto a near-circle

Conformal mapping is a classical technique utilized for solving problems in aerodynamics and hydrodynamics. Conformal mapping is utilized in the construction of grids around airfoils, engine inlets and other aircraft configurations. These shapes are transformed onto a near-circle image for which the equations of fluid motion are discretized on the mapped plane and solved numerically by utilizing the appropriate techniques. In comparison to other grid-generation techniques such as algerbraic or differential type, conformal mapping offers an analytical and accurate form even if the grid deformation is large. One of the most appealing features is that the grid can be constrained to remain orthogonal to the body after the transformation. Hence, the grid is suitable for analyzing the supersonic flow past a blunt object. The associated shock as a coordinate surface adjusts its position in the course of computation until convergence is reached. The present work applied conformal mapping to 3-D bodies with no axis of symmetry such as the Aerobraking Flight Experiment (AFE) vehicle, transforming the AFE shape onto a near-circle image. A numerical procedure and code are used to generate grids around the AFE body.

Hommel, M. J.↗
NASA NTRS1999

An Adaptive Unstructured Grid Method by Grid Subdivision, Local Remeshing, and Grid Movement

An unstructured grid adaptation technique has been developed and successfully applied to several three dimensional inviscid flow test cases. The approach is based on a combination of grid subdivision, local remeshing, and grid movement. For solution adaptive grids, the surface triangulation is locally refined by grid subdivision, and the tetrahedral grid in the field is partially remeshed at locations of dominant flow features. A grid redistribution strategy is employed for geometric adaptation of volume grids to moving or deforming surfaces. The method is automatic and fast and is designed for modular coupling with different solvers. Several steady state test cases with different inviscid flow features were tested for grid/solution adaptation. In all cases, the dominant flow features, such as shocks and vortices, were accurately and efficiently predicted with the present approach. A new and robust method of moving tetrahedral "viscous" grids is also presented and demonstrated on a three-dimensional example.

Pirzadeh, Shahyar Z.↗
DOE OSTI2024

New dynamic diamond anvil cell for time-resolved radial x-ray diffraction

The dynamic diamond anvil cell (dDAC) is a recently developed experimental platform that has shown promise for studying the behavior of materials at strain rates ranging from intermediate to quasi-static and shock compression regimes. Combining dDAC with time-resolved x-ray diffraction (XRD) in the radial geometry (i.e., with incident x-rays perpendicular to the axis of compression) enables the study of material properties such as strength, texture evolution, and deformation mechanisms. This work describes a radial XRD dDAC setup at beamline P02.2 (Extreme Conditions Beamline) at DESY’s PETRA III synchrotron. Time-resolved radial XRD data are collected for titanium, zirconium, and zircon samples, demonstrating the ability to study the strength and texture of materials at compression rates above 300 GPa/s. In addition, the simultaneous optical imaging of the DAC sample chamber is demonstrated. The ability to conduct simultaneous radial XRD and optical imaging provides the opportunity to characterize plastic strain and deviatoric strain rates in the DAC at intermediate rates, exploring the strength and deformation mechanisms of materials in this regime.

47 OTHER INSTRUMENTATION↗
NASA NTRS1994

The formation of small grains in shocks in the ISM

Carbonaceous and silicate grains swept up, and betatron accelerated, by supernova-generated shock waves in the interstellar medium are exposed to grain destructive processing. The degree of grain destruction is determined by the differential gas-grain and grain-grain velocities, which lead to sputtering of the grain surface and grain core disruption (deformation, vaporization and shattering), respectively. The threshold pressure for grain shattering in grain-grain collisions (100 k bar) is considerably lower than that for vaporization (approximately 5 M bar). Therefore, collisions between grains shatter large grains into smaller fragments (i.e., small grains and PAH's). Using a new algorithms for the destructive processes, it was possible to model the formation fo small grain fragments in grain-grain collisions in the warm phase of the interstellar medium. It was found that in one cycle through the warm medium (approximately 3 x 10(sup 6) years) of order 1-2% of the total grain mass is shattered into particles with radii of less than 50 A.

Jones, Anthony P.↗
DOE OSTI2024

High pressure suppression of plasticity due to an overabundance of shear embryo formation

Abstract High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is known to occur in a variety of these materials, the governing dynamics of the mechanisms are not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands, “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high-strain rate application materials such as those for high explosives and hypersonic systems.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Experimental Hypersonic Aerodynamic Characteristics of the Space Shuttle Orbiter for a Range of Damage Scenarios

Aerodynamic tests in support of the Columbia accident investigation were conducted in two hypersonic wind tunnels at the NASA Langley Research Center, the 20-Inch Mach 6 Air Tunnel and the 20-Inch CF4 Tunnel. The primary purpose of these tests was to measure the forces and moments generated by a variety of outer mold line alterations (damage scenarios) using 0.0075-scale models of the Space Shuttle Orbiter. Simultaneously acquired global heat transfer mappings were obtained for a majority of the configurations tested. Test parametrics included angles of attack from 38 to 42 deg, unit Reynolds numbers from 0.3 x 10(exp 6) to 3.0 x 10(exp 6) per foot, and normal shock density ratios of 5 (Mach 6 air) and 12 (CF4). The damage scenarios evaluated included asymmetric boundary layer transition, gouges in the windward surface thermal protection system tiles, wing leading edge damage (partially and fully missing reinforced carbon-carbon (RCC) panels), deformation of the wing windward surface, and main landing gear and/or door deployment. The measured aerodynamic increments for the damage scenarios examined were generally small in magnitude, as were the flight-derived values during most of the entry prior to loss of communication. A progressive damage scenario is presented that qualitatively matches the flight observations for the STS-107 entry.

Brauckmann, Gregory J.↗
NASA NTRS1993

Stratigraphy and sedimentology of the K/T boundary deposit in Haiti

The K/T boundary sequence is exposed in uplifted carbonate sediments of the southwest peninsula of Haiti. It is found at 15 localities within the Beloc formation, a sequence of limestone and marls interpreted as a monoclinal nappe structure thrust to the north. This tectonic deformation has affected the K/T boundary deposit to varying degrees. In some cases the less competent K/T deposit has acted as a slip plane leading to extensive shearing of the boundary layer, as well as duplication of the section. The presence of glassy tektites, shocked quartz, and an Ir anomaly directly link the deposit to a bolide impact. Stratigraphic and sedimentological features of the tripartite sequence indicate that it was formed by deposition from ballistic fallout of coarse tektites, emplacement of particle gravity flows and fine grained fallout of widely dispersed impact ejecta.

Carey, S.↗
NASA NTRS2002

Flutter Analysis of a Transonic Fan

This paper describes the calculation of flutter stability characteristics for a transonic forward swept fan configuration using a viscous aeroelastic analysis program. Unsteady Navier-Stokes equations are solved on a dynamically deforming, body fitted, grid to obtain the aeroelastic characteristics using the energy exchange method. The non-zero inter-blade phase angle is modeled using phase-lagged boundary conditions. Results obtained show good correlation with measurements. It is found that the location of shock and variation of shock strength strongly influenced stability. Also, outboard stations primarily contributed to stability characteristics. Results demonstrate that changes in blade shape impact the calculated aerodynamic damping, indicating importance of using accurate blade operating shape under centrifugal and steady aerodynamic loading for flutter prediction. It was found that the calculated aerodynamic damping was relatively insensitive to variation in natural frequency.

Srivastava, R.↗
DOE OSTI2021

Dynamic flow stress of pure polycrystalline aluminum: Pressure-shear plate impact experiments and extension of dislocation-based modeling to large strains

The dynamic thermo-mechanical behavior of pure aluminum has attracted renewed interest lately due to experimental observations of an anomalous increase in Hugoniot Elastic Limit (HEL) at incipient plasticity and elevated temperatures in polycrystalline pure metals. In context of current dislocation-mediated plasticity models for metals, this increase in dynamic strength is indicative of a transition in the rate controlling mechanism for dislocation glide from being thermally assisted to phonon-drag restricted due to increase in phonon viscosity at elevated temperatures. Though these studies have helped to shed light on these important mechanisms operative in FCC metals at incipient plasticity, the extent to which they contribute to flow stress, particularly at larger strains, remains unclear. Here, we address these questions through a combined experimental and modeling effort focused on investigating the evolution of dynamic flow stress in polycrystalline aluminum using experimental data gathered from a series of combined pressure-and-shear plate impact (PSPI) experiments designed to reveal the flow stress of pure aluminum at strain rates ~ 10 5 /s, plastic strains of up to 40% and temperatures ranging from room to 866 K. In all cases, the flow stress of aluminum, as inferred from the measured transverse particle velocity histories at the free surface of an elastic tungsten carbide target plate, reveals saturation with increasing plastic strains at stress levels that decrease with increasing test temperatures. Numerical simulations are performed to correlate the experimentally observed temperature and strain rate dependence of flow stress at small and large plastic strains using the Austin-McDowell dislocation-mediated plasticity model parametrized to normal plate impact experiments conducted in an earlier study by Zaretsky and Kanel (2012). Extensions to the model are made to better represent dynamic behavior of pure aluminum at larger plastic strains as observed in elevated temperature split Hopkinson Pressure bar (SPHB) experiments of Samanta (1971) and Lindholm and Yeakly (1965), and the combined pressure-and-shear plate-impact experiments conducted in the present study. The main theoretical extension to the Austin-McDowell model made in this paper is the introduction of a new rate- and temperature- dependent dynamic recovery function, which can potentially allow for an effective reduction in the rate of accumulation of dislocations at large plastic strains. The numerical predictions of the revised and re-calibrated plasticity model are brought into agreement with the experimental observations, correlating sufficiently well with dynamic yield stress at incipient plasticity and the flow stress levels at larger plastic strains, plastic strain rates in the range 10 3 – 10 6 /s, and elevated temperatures up to near melt. The model, in concert with the experimental measurements, suggests that at incipient plasticity the rate governing mechanism for plastic flow is phonon-drag restricted dislocation glide, whereas, at higher magnitudes of plastic strain, it transitions to stress-assisted thermally activated glide. The transition between these two rate governing mechanisms is controlled by the evolution of dislocations throughout the deformation process.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Rapid activation of non-oriented mechanophores via shock loading and spallation

Mechanophores (MPs), stimuli-responsive molecules that respond chromatically to mechanochemical reactions, are important for understanding the coupling between mechanics and chemistry as well as in engineering applications. However, the atomic-level understanding of their activation originates from gas phase studies or under simple linear elongation forces directly on molecules or polymer chains containing MPs. The effect of many-body distortions, pervasive in condensed-phase applications, is not understood. Therefore, we performed large-scale molecular dynamics (MD) simulations of a poly(methyl methacrylate)-spiropyran copolymer under dynamic mechanical loading and studied the activation of the MP under various conditions from dynamical compression to tension during unloading. Detailed analysis of the all-atom MD trajectories shows that the MP blocks experience significant many-body intramolecular distortion that can significantly decrease the activation barrier as compared with when deformation rates are slow relative to molecular relaxation time scales. Here we find that the reactivity of MPs under material compression states is governed by many-body effects of intramolecular torsions, whereas under tension, the reactions are governed by tensile stresses.

36 MATERIALS SCIENCE↗
DOE OSTI2023

In-Situ X-Ray Imaging High Strain Rate Compression of Laminate Al-Graphene Composite and Mechanical Property Characterization

Owing to the opaque nature of the laminated structures, traditional high-speed optical camera cannot be used to detect the dynamic process of sub-surface deformation. In this article, we report a study of using high speed X-ray imaging to study the high strain rate deformation in laminated Al structures. We used a Kolsky bar apparatus to apply dynamic compression and a high-speed synchrotron X-ray phase contrast imaging (PCI) setup to conduct the in situ X-ray imaging study. The in situ X-ray imaging captures the shock wave propagation in the laminated structures. After shock compression, we characterized the microstructures by using transmission electron microscopy (TEM), which demonstrates an increase of dislocation density. In conclusion, the micro-pillar compression tests show that the yield strength at 0.2% offset of laminated Al-graphene composite has a significant increase of 67%, from 30 to 50 MPa, compared to laminate Al after shock loading.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Laser-driven ramp-compression experiments on the national ignition facility

This report details the analyses and related uncertainties in measuring longitudinal-stress–density paths in indirect laser-driven ramp equation-of-state (EOS) experiments. Experiments were conducted at the National Ignition Facility (NIF) located at the Lawrence Livermore National Laboratory. The NIF can deliver up to 2 MJ of laser energy over 30 ns and provide the necessary laser power and control to ramp compress materials to TPa pressures (1 TPa = 10 × 10 6 atmospheres). These data provide low-temperature solid-state EOS data relevant to the extreme conditions found in the deep interiors of giant planets. In these experiments, multi-stepped samples with thicknesses in the range of 40–120 µm experience an initial shock compression followed by a time-dependent ramp compression to peak pressure. Interface velocity measurements from each thickness combine to place a constraint on the Lagrangian sound speed as a function of particle velocity, which in turn allows for the determination of a continuous stress–density path to high levels of compressibility. In this report, we present a detailed description of the experimental techniques and measurement uncertainties and describe how these uncertainties combine to place a final uncertainty in both stress and density. Here, we address the effects of time-dependent deformation and the sensitivity of ramp EOS techniques to the onset of phase transformations.

47 OTHER INSTRUMENTATION↗