Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple shocks”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Basic Governing Equations for the Flight Regimes of Aeroassisted Orbital Transfer Vehicles

The basic governing equations for the low-density, high-enthalpy flow regimes expected in the shock layers over the heat shields of the proposed aeroassisted orbital transfer vehicles are derived by combining and extending existing theories. The conservation equations are derived from gas kinetic principles for a four-component ionized gas consisting of neutral molecules, neutral atoms, singly ionized ions, and electrons, assuming a continuum flow. The differences among translational-rotational, vibrational, and electron temperatures are accounted for, as well as chemical nonequilibrium and electric-charge separation. Expressions for convective and viscous fluxes, transport properties, and the terms representing interactions among various energy modes are explicitly given. The expressions for the rate of electron-vibration energy transfer, which violates the Landau-Teller conditions, are derived by solving the system of master equations accounting for the multiple-level transitions.

Lee, Jong-Hun↗
DOE OSTI2023

Hypervelocity impact in stellar media: heat shielding, shock fronts and ablation clouds

The Frontier Plasma Science experiments proposed by our team will model the processes occurring during spacecraft atmospheric entries and meteorite planetary collisions by inserting material targets into the Scrape off Layer (SOL) and edge plasma of a tokamak. This project will investigate important plasma-material processes during high-enthalpy atmospheric entries, including heat flux generation and material ablation. Of particular interest is to explore heat and particle flux detachment in front of the sample using concepts and techniques developed in the latest studies of divertor detachment. Exploration missions to the Solar System’s gaseous giants and hyperbolic re-entries into the Earth’s atmosphere require spacecraft that can withstand high velocity (>10 km/s), high heat fluxes (>10 MW/m2), and corresponding enthalpies. Ablative materials have been used as thermal shields to protect the spacecraft from severe heating during entry. However, developing high enthalpy ablating materials is challenging due to the lack of adequate ground testing facilities. Our team stated that the performance of candidate shielding materials could be assessed in a laboratory environment with multiple in-situ diagnostics, which was enabled by the progress in the tokamak research. Modern machines feature relatively long discharges (~10 s) with well-controlled stable plasma conditions at the edge where the heat flux and the flow speed are similar to those experienced during atmospheric entries. The experiments were carried out at the DIII-D National Fusion Facility in San Diego, CA, operated by General Atomics (GA) for the DoE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Core-collapse supernova neutrino emission and detection informed by state-of-the-art three-dimensional numerical models

ABSTRACT Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, namely Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral standing accretion shock instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate (${\gtrsim} 50 {{\ \rm per\ cent}}$), but the time-integrated signal has a relatively modest asymmetry (${\lesssim} 20 {{\ \rm per\ cent}}$). Both features are associated with the lepton-number emission self-sustained asymmetry and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavours of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of ∼20 per cent if the distance to the CCSN is ≲6 kpc.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Transonic Drag Prediction on a DLR-F6 Transport Configuration Using Unstructured Grid Solvers

A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.

Lee-Rausch, E. M.↗
DOE OSTI2023

Linear shaped-charge jet optimization using machine learning methods

Linear shaped charges are used to focus energy into rapidly creating a deep linear incision. The general design of a shaped charge involves detonating a confined mass of high explosive (HE) with a metal-lined concave cavity on one side to produce a high velocity jet for the purpose of striking and penetrating a given material target. This jetting effect occurs due to the interaction of the detonation wave with the cavity geometry, which produces an unstable fluid phenomenon known as the Richtmyer–Meshkov instability and results in the rapid growth of a long narrow jet. We apply machine learning and optimization methods to hydrodynamics simulations of linear shaped charges to improve the simulated jet characteristics. The designs that we propose and investigate in this work generally involve modifying the behavior of the detonation waves prior to interaction with the liner material. These designs include the placement of multiple detonators and the use of metal inclusions within the HE. In conclusion, we are able to produce a linear shaped-charge design with a higher penetration depth than the baseline case that we consider and accomplish this using the same amount of or less HE.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Examination of the cerium α-ε phase transition under dynamic loading with x-ray diffraction

Recent examination of the cerium Hugoniot with pyrometry and x-ray diffraction (XRD) has revealed a narrow solid-liquid two-phase region. It has been suggested that nonequilibrium melting may be occurring along the Ce Hugoniot, with either melt kinetics or a sluggish α-ε transition impeding the transition. In particular, the kinetics of the α-ε is unknown and the location of the phase boundary is in dispute. Static measurements suggest a nearly vertical phase boundary that intersects the Hugoniot at 6–7 GPa. This lies in direct conflict with dynamic measurements along the Hugoniot observing α-Ce through incipient melt. This work presents dynamic experiments using XRD to examine the behavior of the α-ε phase transition. The results show that the α-ε phase transition occurs through a tetragonal distortion, with the transition beginning at temperatures below the solid Hugoniot. Following the initial deviation from an ideal fcc structure, the c/a ratio is found to gradually increase with no steady value observed in the ε phase within the range of these experiments (below 17 GPa). Multiple diffraction patterns captured during the peak stress state show no significant change in c/a ratio prior to uniaxial release, upon which Ce reverts to an fcc structure. Furthermore, the results indicate that the α-ε transition occurs rapidly, both on loading and release. An examination of the c/a ratio with increasing temperatures suggests 11.5 GPa as a lower bound for the location of the α-ε-liquid triple point.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

IRIS Si iv LINE PROFILES: AN INDICATION FOR THE PLASMOID INSTABILITY DURING SMALL-SCALE MAGNETIC RECONNECTION ON THE SUN

Our understanding of the process of fast reconnection has undergone a dramatic change in the last 10 years driven, in part, by the availability of high-resolution numerical simulations that have consistently demonstrated the break-up of current sheets into magnetic islands, with reconnection rates that become independent of Lundquist number, challenging the belief that fast magnetic reconnection in flares proceeds via the Petschek mechanism which invokes pairs of slow-mode shocks connected to a compact diffusion region. The reconnection sites are too small to be resolved with images, but these reconnection mechanisms, Petschek and the plasmoid instability, have reconnection sites with very different density and velocity structures and so can be distinguished by high-resolution line-profile observations. Using IRIS spectroscopic observations we obtain a survey of typical line profiles produced by small-scale events thought to be reconnection sites on the Sun. Slit-jaw images are used to investigate the plasma heating and re-configuration at the sites. A sample of 15 events from 2 active regions is presented. The line profiles are complex with bright cores and broad wings extending to over 300 km s{sup −1}. The profiles can be reproduced with the multiple magnetic islands and acceleration sites that characterize the plasmoid instability but not by bi-directional jets that characterize the Petschek mechanism. This result suggests that if these small-scale events are reconnection sites, then fast reconnection proceeds via the plasmoid instability, rather than the Petschek mechanism during small-scale reconnection on the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SMALL-SCALE MAGNETIC ISLANDS IN THE SOLAR WIND AND THEIR ROLE IN PARTICLE ACCELERATION. I. DYNAMICS OF MAGNETIC ISLANDS NEAR THE HELIOSPHERIC CURRENT SHEET

Increases of ion fluxes in the keV–MeV range are sometimes observed near the heliospheric current sheet (HCS) during periods when other sources are absent. These resemble solar energetic particle events, but the events are weaker and apparently local. Conventional explanations based on either shock acceleration of charged particles or particle acceleration due to magnetic reconnection at interplanetary current sheets (CSs) are not persuasive. We suggest instead that recurrent magnetic reconnection occurs at the HCS and smaller CSs in the solar wind, a consequence of which is particle energization by the dynamically evolving secondary CSs and magnetic islands. The effectiveness of the trapping and acceleration process associated with magnetic islands depends in part on the topology of the HCS. We show that the HCS possesses ripples superimposed on the large-scale flat or wavy structure. We conjecture that the ripples can efficiently confine plasma and provide tokamak-like conditions that are favorable for the appearance of small-scale magnetic islands that merge and/or contract. Particles trapped in the vicinity of merging islands and experiencing multiple small-scale reconnection events are accelerated by the induced electric field and experience first-order Fermi acceleration in contracting magnetic islands according to the transport theory of Zank et al. We present multi-spacecraft observations of magnetic island merging and particle energization in the absence of other sources, providing support for theory and simulations that show particle energization by reconnection related processes of magnetic island merging and contraction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

STATISTICALLY DETERMINED DISPERSION RELATIONS OF MAGNETIC FIELD FLUCTUATIONS IN THE TERRESTRIAL FORESHOCK

We obtain dispersion relations of magnetic field fluctuations for two crossings of the terrestrial foreshock by Cluster spacecraft. These crossings cover plasma conditions that differ significantly in their plasma β and in the density of the reflected ion beam, but not in the properties of the encountered ion population, both showing shell-like distribution function. Dispersion relations are reconstructed using two-point instantaneous wave number estimations from pairs of Cluster spacecraft. The accessible range of wave vectors, limited by the available spacecraft separations, extends to ≈2 × 10{sup 4} km. Results show multiple branches of dispersion relations, associated with different powers of magnetic field fluctuations. We find that sunward propagating fast magnetosonic waves and beam resonant modes are dominant for the high plasma β interval with a dense beam, while the dispersions of the interval with low beam density include Alfvén and fast magnetosonic modes propagating sunward and anti-sunward.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A method for assessing economic, environmental, and reliability tradeoffs of interregional transmission connecting ERCOT (the Texas grid) to the eastern and western grids

Reliable development of the power grid is an evolving concern for humanity due to extreme weather that frequently threatens power sector infrastructure. The state of Texas is a uniquely structured testbed for grid planners to study when looking for solutions to development, innovation, and overcoming such challenges. Because of its size and islanded structure, Texas is small enough to model, but big enough to matter. Texas is a global leader in energy production, energy consumption, and maintains an unusually diverse fuel mix. In addition, the state has experienced winter freezes, heat waves, wind storms, droughts and floods that have threatened power sector infrastructure or caused recent blackouts and calls for demand side conservation. One of the most devastating of these events was the North American winter storm, dubbed “Winter Storm Uri” by the Weather Channel, that froze the region in February 2021 and led to an extended power outage event that put the majority of Texan residents in darkness for days. While preparing to avoid such outage events in the future, various tools have been proposed to improve grid reliability, including energy efficiency, demand response, and distributed energy resources. An additional option would be to develop interregional transmission that connects the Texas grid to other national grids. To assess the merits of this idea, we developed a novel, universally-applicable and internationally-relevant framework to study how the Texas grid would evolve alongside access to various interregional ties. This method allows us to stress the synthetic grid structure and analyze how it would respond to the shock of a simulated winter storm event. Our method leverages open-source modeling tools, such as PowerGenome, pyGRETA, and GenX to synthesize unique zonal grid data, construct a consolidated network of model regions, and simulate different developmental pathways of capacity expansion and operational dispatch. We demonstrate our method with an analysis connecting the Electric Reliability Council of Texas (ERCOT), the grid that serves most of Texas, the Western Electricity Coordinating Council (WECC), the grid that serves the western half of the contiguous U.S., and the Eastern Interconnect, the grid that serves the eastern half of the contiguous U.S. Our results indicate that the cost-optimal capacity of interregional transmission connecting the ERCOT grid to other grids lies between 9–13 GW assuming baseline conditions. Building this amount of connecting capacity in one or multiple directions lowers the costs and emissions of development and operation by up to $16 billion and 257 million metric tonnes (MMT) respectively. Additionally, our results show that the interregional connections between ERCOT and other national grids reduce the amount of total load shed required through mild winter storm events. However, our results also show that there is a threshold of very extreme winter storm conditions, spanning multiple service areas, above which the connections exacerbate resource adequacy problems. Therefore, the results indicate that the connections need to be carefully planned alongside the rest of the grid infrastructure to avoid over-reliance on specific resources or technology options.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Progenitor, precursor, and evolution of the dusty remnant of the stellar merger M31-LRN-2015

ABSTRACT M31-LRN-2015 is a likely stellar merger discovered in the Andromeda Galaxy in 2015. We present new optical to mid-infrared photometry and optical spectroscopy for this event. Archival data show that the source started to brighten ∼2 yr before the nova event. During this precursor phase, the source brightened by ∼3 mag. The light curve at 6 and 1.5 months before the main outburst may show periodicity, with periods of 16 ± 0.3 and 28.1 ± 1.4 d, respectively. This complex emission may be explained by runaway mass-loss from the system after the binary undergoes Roche lobe overflow, leading the system to coalesce in tens of orbital periods. While the progenitor spectral energy distribution shows no evidence of pre-existing warm dust in the system, the remnant forms an optically thick dust shell at approximately four months after the outburst peak. The optical depth of the shell increases dramatically after 1.5 yr, suggesting the existence of shocks that enhance the dust formation process. We propose that the merger remnant is likely an inflated giant obscured by a cooling shell of gas with mass ∼0.2 M⊙ ejected at the onset of the common envelope phase.

(ISM:) dust, extinction↗
NASA NTRS2017

Ames Engineering Directorate

The Ames Engineering Directorate is the principal engineering organization supporting aerospace systems and spaceflight projects at NASA's Ames Research Center in California's Silicon Valley. The Directorate supports all phases of engineering and project management for flight and mission projects-from R&D to Close-out-by leveraging the capabilities of multiple divisions and facilities.The Mission Design Center (MDC) has full end-to-end mission design capability with sophisticated analysis and simulation tools in a collaborative concurrent design environment. Services include concept maturity level (CML) maturation, spacecraft design and trades, scientific instruments selection, feasibility assessments, and proposal support and partnerships. The Engineering Systems Division provides robust project management support as well as systems engineering, mechanical and electrical analysis and design, technical authority and project integration support to a variety of programs and projects across NASA centers. The Applied Manufacturing Division turns abstract ideas into tangible hardware for aeronautics, spaceflight and science applications, specializing in fabrication methods and management of complex fabrication projects. The Engineering Evaluation Lab (EEL) provides full satellite or payload environmental testing services including vibration, temperature, humidity, immersion, pressure/altitude, vacuum, high G centrifuge, shock impact testing and the Flight Processing Center (FPC), which includes cleanrooms, bonded stores and flight preparation resources. The Multi-Mission Operations Center (MMOC) is composed of the facilities, networks, IT equipment, software and support services needed by flight projects to effectively and efficiently perform all mission functions, including planning, scheduling, command, telemetry processing and science analysis.

Small Satellites↗
DOE OSTI2021

Collective optical Thomson scattering in pulsed-power driven high energy density physics experiments (invited)

Optical collective Thomson scattering is used to diagnose magnetized high energy density physics experiments at the Magpie pulsed-power generator at Imperial College London. The system uses an amplified pulse from the 2nd harmonic of a Nd:YAG laser (3 J, 8 ns, 532 nm) to probe a wide diversity of hightemperature plasma objects; with densities in the range of 10 17 -10 19 cm -3 and temperatures between 10 eV and a few keV. The scattered light is collected from 100 μm-scale volumes within the plasmas, which are imaged onto optical fiber arrays. Multiple collection systems observe these volumes from different directions, providing simultaneous probing with different scattering K-vectors (and different associated α-parameters, typically in the range 0.5 – 3) allowing independent measurements of separate velocity components of the bulk plasma flow. The fiber arrays are coupled to an imaging spectrometer with a gated ICCD. The spectrometer is configured to view the ion-acoustic waves (IAWs) of the collective Thomson scattered spectrum. Fits to the spectra with the theoretical spectral density function S(K,ω) yield measurements of the local plasma temperatures and velocities. Fitting is constrained by independent measurements of the electron density from laser interferometry, and the corresponding spectra for different scattering vectors. Furthermore, this TS diagnostic has been successfully implemented on a wide range of experiments, revealing temperature and flow velocity transitions across magnetized shocks, inside rotating plasma jets and imploding wire arrays, as well as providing direct measurements of drift velocities inside a magnetic reconnection current sheet.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Metagenomic insights into the taxonomy, function, and dysbiosis of prokaryotic communities in octocorals

Background. In octocorals (Cnidaria Octocorallia), the functional relationship between host health and its symbiotic consortium has yet to be determined. Here, we employed comparative metagenomics to uncover the distinct functional and phylogenetic features of the microbiomes of healthy Eunicella gazella, Eunicella verrucosa, and Leptogorgia sarmentosa tissues, in contrast with the microbiomes found in seawater and sediments. We further explored how the octocoral microbiome shifts to a pathobiome state in E. gazella. Results. Multivariate analyses based on 16S rRNA genes, Clusters of Orthologous Groups of proteins (COGs), Protein families (Pfams), and secondary metabolite-biosynthetic gene clusters annotated from 20 Illumina-sequenced metagenomes each revealed separate clustering of the prokaryotic communities of healthy tissue samples of the three octocoral species from those of necrotic E. gazella tissue and surrounding environments. While the healthy octocoral microbiome was distinguished by so-far uncultivated Endozoicomonadaceae, Oceanospirillales, and Alteromonadales phylotypes in all host species, a pronounced increase of Flavobacteriaceae and Alphaproteobacteria, originating from seawater, was observed in necrotic E. gazella tissue. Increased abundances of eukaryotic-like proteins, exonucleases, restriction endonucleases, CRISPR/Cas proteins, and genes encoding for heat-shock proteins, inorganic ion transport, and iron storage distinguished the prokaryotic communities of healthy octocoral tissue regardless of the host species. An increase of arginase and nitric oxide reductase genes, observed in necrotic E. gazella tissues, suggests the existence of a mechanism for suppression of nitrite oxide production by which octocoral pathogens may overcome the host’s immune system. Conclusions. This is the first study to employ primer-less, shotgun metagenome sequencing to unveil the taxonomic, functional, and secondary metabolism features of prokaryotic communities in octocorals. Our analyses reveal that the octocoral microbiome is distinct from those of the environmental surroundings, is host genus (but not species) specific, and undergoes large, complex structural changes in the transition to the dysbiotic state. Host-symbiont recognition, abiotic-stress response, micronutrient acquisition, and an antiviral defense arsenal comprising multiple restriction endonucleases, CRISPR/Cas systems, and phage lysogenization regulators are signatures of prokaryotic communities in octocorals. We argue that these features collectively contribute to the stabilization of symbiosis in the octocoral holobiont and constitute beneficial traits that can guide future studies on coral reef conservation and microbiome therapy.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2006

Heliosheath Space Environment Interactions with Icy Bodies in the Outermost Solar System

The Voyager 1 and 2 spacecraft are exploring the space environment of the outermost solar system at the same time that earth-based astronomy continues to discover new icy bodies, one larger than Pluto, in the transitional region outward from the Classical Kuiper Belt to the Inner Oort Cloud. Some of the Scattered Disk Objects in this region periodically pass through the heliosheath, entered by Voyager 1 in Dec. 2004 and later expected to be reached by Voyager 2, and out even beyond the heliopause into the Very Local Interstellar Medium. The less energetic heliosheath ions, important for implantation and sputtering processes, are abundant near and beyond the termination shock inner boundary, but the source region of the more penetrating anomalous cosmic ray component has not yet been found. Advantageous for modeling of icy body interactions, the measured heliosheath flux spectra are relatively more stable within this new regime of isotropic compressional magnetic turbulence than in the upstream heliospheric environment. The deepest interactions and resultant radiation-induced chemistry arise from the inwardly diffusing component of the galactic cosmic ray ions with significant intensity modulation also arising in the heliosheath beyond Voyager 1. Surface gardening by high-velocity impacts of smaller bodies (e.g., fragments of previous KBO collisions) and dust is a further space weathering process setting the time scales for long term exposure of different regolith layers to the ion irradiation. Sputtering and ionization of impact ejecta grains may provide a substantial feedback of pickup ions for multiple cycles of heliosheath acceleration and icy body interaction. Thus the space weathering interactions are potentially of interest not only for effects on sensible surface composition of the icy bodies but also for evolution of the heliosheath plasma energetic ion, and neutral emission environment.

Cooper, John F.↗
DOE OSTI2022

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration

Turbulent Richtmyer–Meshkov instability (RMI) is investigated through a series of high resolution three-dimensional simulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regimes combined. Overall, the growth rate exponent is determined as θ=0.292±0.009, in good agreement with prior studies; however, the exponent is decaying slowly in time. Also, θ is shown to be relatively insensitive to the choice of mixing layer width measurements. The asymptotic integral molecular mixing measures Θ=0.792±0.014, Ξ=0.800±0.014, and Ψ=0.782±0.013 are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49% and 66% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single θ. In conclusion, quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗