Search NASA⌕ Search

SEARCH · Search NASA

Results for “Free surface velocity”

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 91 records · Page 5

Spall and subsequent recompaction of copper under shock loading

Research was undertaken to study spall and subsequent recompaction of oxygen-free high-thermal-conductivity copper using a single-stage large-bore light gas gun capable of planar impacts. Gun experiments were conducted that produced an initial spall in the target with a subsequent recompaction of the spall damage/layer by the use of layered impactors. Symmetric spall experiments at similar conditions were also conducted as a control to the recompaction experiments. Photonic Doppler velocimetry was used to obtain the velocity history of the target’s back surface; these velocity histories were analyzed, and the results were compared to numerical simulations. After recovery of selected samples, the microstructure was analyzed using optical imaging microscopy and electron backscatter diffraction. Recompaction waves were clearly observed in the time-resolved data obtained at the back surface of the targets. Two different series of experiments were performed to understand the effect of damage morphology on the stress required for spall recompaction. In the initial series, the peak pressure of the first shock was systematically increased to alter the amount of damage in the material. While in the second series, the peak pressure of the second shock was increased, keeping the magnitude of the first shock relatively constant to understand the effect of peak stress on the recompaction process. Findings clearly show that recompaction results in a perturbed band of the microstructure at the location of the expected spall plane in all cases except where the peak stress of the second shock is lower than the critical value.

, recompaction↗

Local models of two-temperature accretion disc coronae – II. Ion thermal conduction and the absence of disc evaporation

We use local stratified shearing-box simulations with magnetic field-aligned thermal conduction to study an idealized model of the coupling between a cold, radiatively efficient accretion disc, and an overlying, hot, two-temperature corona. Evaporation of a cold disc by conduction from the hot corona has been proposed as a means of mediating the soft-to-hard state transitions observed in X-ray binary systems. We model the coronal plasma in our local disc patch as an MHD fluid subject to both free-streaming ion conduction and a parametrized cooling function that captures the collisional transfer of energy from hot ions to colder, rapidly cooling leptons. In all of our models, independent of the initial net vertical magnetic flux (NF) threading the disc, we find no evidence of disc evaporation. The ion heat flux into the disc is radiated away before conduction can heat the disc’s surface layers. When an initial NF is present, steady-state temperature, density, and outflow velocities in our model coronae are unaffected by conduction. Instead of facilitating disc evaporation, thermal conduction is more likely to feed the disc with plasma condensing out of the corona, particularly in flows without NF. Our work indicates that uncertainties in the amount of NF threading the disc hold far greater influence over whether or not the disc will evaporate into a radiatively inefficient accretion flow compared to thermal conduction. We speculate that a change in net flux mediates disc truncation/evaporation.

79 ASTRONOMY AND ASTROPHYSICS↗

Few-femtosecond resolved imaging of laser-driven nanoplasma expansion

Abstract The free expansion of a planar plasma surface is a fundamental non-equilibrium process relevant for various fields but as-yet experimentally still difficult to capture. The significance of the associated spatiotemporal plasma motion ranges from astrophysics and controlled fusion to laser machining, surface high-harmonic generation, plasma mirrors, and laser-driven particle acceleration. Here, we show that x-ray coherent diffractive imaging can surpass existing approaches and enables the quantitative real-time analysis of the sudden free expansion of laser-heated nanoplasmas. For laser-ionized SiO 2 nanospheres, we resolve the formation of the emerging nearly self-similar plasma profile evolution and expose the so far inaccessible shell-wise expansion dynamics including the associated startup delay and rarefaction front velocity. Our results establish time-resolved diffractive imaging as an accurate quantitative diagnostic platform for tracing and characterizing plasma expansion and indicate the possibility to resolve various laser-driven processes including shock formation and wave-breaking phenomena with unprecedented resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface and sub-surface porosity effects on the initial free surface expansion in shocked single crystal aluminum

Defects and roughness at a material’s surface can allow for the formation of a material jet when a shockwave reaches a surface. The formation of a jet can lead to significant material damage of the shocked material, as well as the downstream material that the jet impacts. Here, we utilize planar shocks in single crystal aluminum to assess initial jetting formation for surface notches and sub-surface porosity. In the case of surface notches, we find that changing the width of the notch can cause the shock breakout response to transition from a standard jetting mechanism to an atomic atomization at high velocities. This transition occurs due to the mechanisms associated with lateral relaxations within the notch and how this interacts with the shock focusing event that induces jetting. For sub-surface pores, the pore collapse induces a re-shock, localized near the pore, that is stronger than the initial shock, but unsupported. Changing the pore dimensions and location relative to the surface influences the strength of the shock that emanates from its collapse when it reaches the surface, leading to a variety of different breakout conditions including jetting and localized material failure.

36 MATERIALS SCIENCE↗

Probing Alkylpolyglycoside Hydrogen Bonding and Its Destabilization by Sulfonate Hydrotropes with X‑ray and Vibrational Spectroscopy

It has been proposed that intersurfactant H-bonding networks can produce highly stable foams, specifically those generated from solutions containing alkylpolyglycoside surfactants. In this work, we aim to characterize the presence and destabilization of these networks by introducing a hydrotrope (sodium p-toluene sulfate) at various concentrations into an alkylpolyglycoside (Glucopon 225 DK) surfactant containing solution. Solution surface properties are probed by using aerosol velocity map imaging X-ray photoelectron spectroscopy (A-VMI-XPS) and tensiometer measurements. Bulk properties below the surface are probed using Fourier transform infrared spectroscopy (FTIR) and C-edge near-edge X-ray fine structure spectroscopy (NEXAFS). The surface measurements provide the framework to describe the destabilization of the solution H-bonding network, while the bulk solution measurements provide hints about the disruption of the hydrogen bonding network upon hydrotrope addition. The collected data support the hypothesis that the destabilization of the intersurfactant H-bonding network frees surfactant molecules from the bulk, increasing surfactant population at the air/water interface. This was quantified through an increase in peak area and width in XPS measurements as well as a decrease in surfactant critical micelle concentration. Three regimes with increasing amounts of hydrotrope addition are suggested, described as (1) a hydrotrope affecting only surfactant surface properties, (2) a hydrotrope affecting surfactant surface and bulk properties, and (3) hydrotrope-dominated surface and bulk properties. Future studies will characterize foam stability across the hydrotrope concentration regimes to better define correlations between intersurfactant H-bonding networks and foam stability.

Molecules↗

CFD simulations of electric motor end ring cooling for improved thermal management

Proper thermal management of an electric motor for vehicle applications extends its operating range. One cooling approach is to impinge Automatic Transmission Fluid (ATF) onto the rotor end ring. Increased ATF coverage correlates to enhanced heat transfer. Computational Fluid Dynamics (CFD) analytical tools provide a mechanism to assess motor thermal management prior to hardware fabrication. The complexity of the fluid flow (e.g., jet atomization, interface tracking, wall impingement) and heat transfer makes these simulations challenging. Computational costs are high when solving these flows on high-speed rotating meshes. Typically, a Volume-of Fluid (VOF) technique (i.e., two-fluid system) is used to resolve ATF dynamics within this rotating framework. Suitable numerical resolution of the relevant physics for thin films under strong inertial forces at high rotor speeds is computationally expensive, further increasing the run times. In this work, a numerical study of rotor-ring cooling by ATF is presented using a patent automated Cartesian cut-cell based method coupled with Automatic Mesh Refinement (AMR). This approach automatically creates the Cartesian mesh on-the-fly and can effectively handle complex rotating geometries by adaptively refining the mesh based on local gradients in the flow field which results in better resolution of the air-ATF interface. A Single non-inertial Reference Frame (SRF) approach is used to account for the rotating geometry and to further improve the overall computational efficiency. Quasi-steady state conditions are targeted in the analysis of the results. Important physics such as ATF jet structure, velocity detail near the air-jet interface, ATF coverage/accumulation on the ring surface, and cooling capacity are presented for a low-resolution Reynolds averaged Navier-Stokes (RANS), high-resolution RANS, and high-resolution Large-Eddy Simulation (LES) models. Computations are scaled over hundreds of cores on a supercomputer to maximize turnaround time. Each numerical approach is shown to capture the general trajectory of the oil jet prior to surface impingement. The high-resolution LES simulation, however, is superior in capturing small scale details and heat transfer between the free jet and surrounding air.

42 ENGINEERING↗

Emergence of liquid following laser melting of gold thin films

X-ray structural science is undergoing a revolution driven by the emergence of X-ray Free-electron Laser (XFEL) facilities. The structures of crystalline solids can now be studied on the picosecond time scale relevant to phonons, atomic vibrations which travel at acoustic velocities. In the work presented here, X-ray diffuse scattering is employed to characterize the time dependence of the liquid phase emerging from femtosecond laser-induced melting of polycrystalline gold thin films using an XFEL. In a previous analysis of Bragg peak profiles, we showed the supersonic disappearance of the solid phase and presented a model of pumped hot electrons carrying energy from the gold surface to scatter at internal grain boundaries. This generates melt fronts propagating relatively slowly into the crystal grains. By conversion of diffuse scattering to a partial X-ray pair distribution function, we demonstrate that it has the characteristic shape obtained by Fourier transformation of the measured F ( Q ). The diffuse signal fraction increases with a characteristic rise-time of 13 ps, roughly independent of the incident pump fluence and consequent final liquid fraction. This suggests the role of further melt-front nucleation processes beyond grain boundaries.

36 MATERIALS SCIENCE↗

Suppression of the vacuum space-charge effect in fs-photoemission by a retarding electrostatic front lens

The performance of time-resolved photoemission experiments at fs-pulsed photon sources is ultimately limited by the e–e Coulomb interaction, downgrading energy and momentum resolution. Here, we present an approach to effectively suppress space-charge artifacts in momentum microscopes and photoemission microscopes. A retarding electrostatic field generated by a special objective lens repels slow electrons, retaining the k-image of the fast photoelectrons. The suppression of space-charge effects scales with the ratio of the photoelectron velocities of fast and slow electrons. Fields in the range from -20 to 01100 V/mm for E kin = 100 eV to 4 keV direct secondaries and pump-induced slow electrons back to the sample surface. Ray tracing simulations reveal that this happens within the first 40 to 3 μm above the sample surface for E kin = 100 eV to 4 keV. An optimized front-lens design allows switching between the conventional accelerating and the new retarding mode. Time-resolved experiments at E kin = 107 eV using fs extreme ultraviolet probe pulses from the free-electron laser FLASH reveal that the width of the Fermi edge increases by just 30 meV at an incident pump fluence of 22 mJ/cm2 (retarding field -21 V/mm). For an accelerating field of +2 kV/mm and a pump fluence of only 5 mJ/cm 2 , it increases by 0.5 eV (pump wavelength 1030 nm). At the given conditions, the suppression mode permits increasing the slow-electron yield by three to four orders of magnitude. The feasibility of the method at high energies is demonstrated without a pump beam at E kin = 3830 eV using hard x rays from the storage ring PETRA III. The approach opens up a previously inaccessible regime of pump fluences for photoemission experiments.

47 OTHER INSTRUMENTATION↗

Long interior carrier lifetime in selective-area InAs nanowires on silicon

Catalyst-free, position-controlled indium arsenide (InAs) nanowires (NWs) of variable diameters were grown on Si (111) by selective-area epitaxy (SAE). Ultrafast pump-probe spectroscopy was conducted, from which carrier recombination mechanisms on the NW surface and interior were resolved and characterized. NWs grown using SAE demonstrated high optical quality, showing minority carrier lifetimes more than two-fold longer than that of the randomly-positioned (RP) NWs. The extracted SAE-InAs NW interior recombination lifetime was found to be as long as 7.2 ns , 13X longer than previous measurements on RP-NWs; and the surface recombination velocity 4154 cm · s - 1 . Transmission electron microscopy revealed a high density of stacking defects within the NWs, suggesting that interior recombination lifetime can be further increased by improving NW interior crystalline quality.

36 MATERIALS SCIENCE↗

WASP-107b’s Density Is Even Lower: A Case Study for the Physics of Planetary Gas Envelope Accretion and Orbital Migration

With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet’s low surface gravity and the star’s brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4 yr Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of the gas envelope of WASP-107b. We reveal that WASP-107b’s mass is only 1.8 Neptune masses (M {sub b} = 30.5 ± 1.7 M {sub ⊕}). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of >85% unless it is substantially inflated. The corresponding core mass of <4.6 M {sub ⊕} at 3σ is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at ≳1 au onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day a {sub b} = 0.0566 ± 0.0017 au. Beyond WASP-107b, we also detect a second, more massive planet (M{sub c}sin i=0.36±0.04M{sub J}) on a wide eccentric orbit (e {sub c} = 0.28 ± 0.07) that may have influenced the orbital migration and spin–orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system’s formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Arctic observations and numerical simulations of surface wind effects on Multi-Angle Snowflake Camera measurements

Abstract. Ground-based measurements of frozen precipitation are heavily influenced by interactions of surface winds with gauge-shield geometry. The Multi-Angle Snowflake Camera (MASC), which photographs hydrometeors in free-fall from three different angles while simultaneously measuring their fall speed, has been used in the field at multiple midlatitude and polar locations both with and without wind shielding. Here, we present an analysis of Arctic field observations – with and without a Belfort double Alter shield – and compare the results to computational fluid dynamics (CFD) simulations of the airflow and corresponding particle trajectories around the unshielded MASC. MASC-measured fall speeds compare well with Ka-band Atmospheric Radiation Measurement (ARM) Zenith Radar (KAZR) mean Doppler velocities only when winds are light (≤5ms-1) and the MASC is shielded. MASC-measured fall speeds that do not match KAZR-measured velocities tend to fall below a threshold value that increases approximately linearly with wind speed but is generally <0.5ms-1. For those events with wind speeds ≤1.5ms-1, hydrometeors fall with an orientation angle mode of 12∘ from the horizontal plane, and large, low-density aggregates are as much as 5 times more likely to be observed. Simulations in the absence of a wind shield show a separation of flow at the upstream side of the instrument, with an upward velocity component just above the aperture, which decreases the mean particle fall speed by 55 % (74 %) for a wind speed of 5 m s−1 (10 m s−1). We conclude that accurate MASC observations of the microphysical, orientation, and fall speed characteristics of snow particles require shielding by a double wind fence and restriction of analysis to events where winds are light (≤5ms-1). Hydrometeors do not generally fall in still air, so adjustments to these properties' distributions within natural turbulence remain to be determined.

54 ENVIRONMENTAL SCIENCES↗

Observed volatilization behavior of silicon carbide in flowing hydrogen above 2000 K

In this work, the intrinsic compatibility of silicon carbide (SiC) and hydrogen (H 2 ) at high temperatures (2000-2473 K) and pressure near one atmosphere was evaluated through a combination of thermodynamic calculations and hot hydrogen exposure testing. Thermodynamic calculations predict the decomposition of SiC in a hydrogen environment to form free silicon (Si) and free carbon (C). Free Si is predicted to vaporize from the surface as a volatile species, while free C may interact with H 2 to form the hydrocarbons CH 4 (T < 2100 K) or C 2 H 2 (T > 2100 K). Coupons of high purity chemical vapor deposition (CVD) β-SiC were exposed to slowly flowing hydrogen at temperatures ranging between 2000 and 2473 K. SiC experienced active attack as the result of H 2 exposure, exhibiting linear weight loss kinetics and an Arrhenius dependence of weight loss on exposure temperature. The linear volatilization constant was experimentally evaluated to correspond with an activation energy of 370 ± 18 kJ/mol. Due to the dependence of observed corrosion rates on gas velocity, corrosion of SiC in flowing H 2 was determined to be governed by external mass transfer of volatile Si species through the boundary layer. Experimentally derived mass losses were in good agreement with mass losses predicted by a boundary layer limited gas diffusion model.

36 MATERIALS SCIENCE↗

Hydrokinetic Power Conversion Using Vortex-Induced Oscillation with Cubic Restoring Force

A cubic-spring restoring function with high-deformation stiffening is introduced to passively improve the harnessed marine hydrokinetic power by using flow-induced oscillations/vibrations (FIO/V) of a cylinder. In these FIO/V experiments, a smooth, rigid, single-cylinder on elastic end-supports is tested at Reynolds numbers ranging from 24,000 < Re < 120,000. The parameters of the tested current energy converter (CEC) are cubic stiffness and linear damping. Using the second generation of digital virtual spring-damping (Vck) controller developed by the Marine Renewable Energy Laboratory (MRELab), the cubic modeling of the oscillator stiffness is tested. Experimental results show the influence of the parameter variation on the amplitude, frequency, energy conversion, energy efficiency, and power of the converter. All experiments are conducted in the low turbulence-free surface water (LTFSW) channel of the MRELab of the University of Michigan. The main conclusions are: (1) The nonlinearity in the cubic oscillator is an effective way to extend the vortex-induced vibration (VIV) upper branch, which results in higher harnessing power and efficiency compared to the linear stiffness cylinder converter. (2) Compared to the linear converter, the overall power increase is substantial. The nonlinear power optimum, occurring at the end of the VIV upper branch, is 63% higher than its linear counterpart. (3) The cubic stiffness converter with low harnessing damping achieves consistently good performance in all the VIV regions because of the hardening restoring force, especially at higher flow velocity.

13 HYDRO ENERGY↗

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗

Hydrokinetic energy harvesting from slow currents using flow-induced oscillations

To harness marine hydrokinetic energy from slow flows, which constitute the majority of currents, tides, and rivers, new Passive Turbulence Control (PTC), consisting of large turbulence stimulators, is tested experimentally on circular cylinders on springs. Here, this study experimentally investigates the effect of PTC on the onset of Flow-Induced Oscillations (FIO) and particularly the relative onset of Vortex-Induced Vibrations (VIV) and galloping. Experiments are conducted in the Low Turbulence Free Surface Water Channel, University of Michigan. Fixed are: mass ratio m* = 1.48, aspect ratio l/D = 10.29, and total damping ratio ζ = 0.04. Parameters are: cylinder diameter D, spring stiffness K, PTC location and height, and flow speed U$\in$[0.36 m/s-1.45 m/s]. Placing the leading edge of PTC at 40–60° induces high amplitude FIO while placement at 10–20° suppresses FIO. As PTC height increases, VIV and galloping initiate earlier and exhibit higher amplitude with a steeper slope. Lower spring stiffness initiates VIV earlier by reducing the oscillator natural frequency in water. Even though large PTC maintained its effectiveness in initiating galloping early, it has no effect on the earlier initiation of VIV, which starts at a nearly fixed reduced velocity. Lower spring stiffness and large PTC enable power generation at low current speed (0.2 m/s).

16 TIDAL AND WAVE POWER↗

NUMERICAL CALCULATION OF LOSSES OF TRAPPED VORTICES UNDER STRONG RF MEISSNER CURRENT AND DC SUPERHEATING FIELD IN TYPE II SUPERCONDUCTORS

Research on the vortex dynamics and enhancing of superheating field in superconductors has attracted much attention in accelerator physics community to develop next-generation high-performance accelerator cavities. However, the extreme dynamics of curvilinear elastic vortices driven by very strong currents close to the depairing limit or superheating field of a superconductor with a nanostructured surface has not been well understood. We calculated the superheating field Hsh and critical momentum kc characterizing the wavelength of the instability ?m of the Meissner state to flux penetration by solving numerically the Ginzburg?Landau equations. A bulk superconductor, superconductor with the inhomogeneous surface disorder (S-S), and multilayered surface (S-I-S) have been thoroughly investigated in this work. Our result showed that S-S and S-I-S structures can enhance the superheating field well above their clean limit. In this work extensive numerical simulation of the power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was performed. Our simulations take into account the nonlinear vortex line tension, vortex mass, Bardeen?Stephen viscous vortex drag applicable at low fields, and nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient ?(v) at high fields and pinning force. We showed that the LO decrease of ?(v) with the vortex velocity v could radically change the field dependence of the surface resistance Ri(H) caused by trapped vortices. At low frequencies Ri(H) exhibits a conventional increase with H. However, as frequency increases, the surface resistance becomes a nonmonotonic function of H which decreases with H at higher fields irrespective of the pinning distribution. Overheating can mask the descending field dependence of Ri(H) as frequency increases. Our numerical simulations also show that the LO effect can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along with the vortex when a vortex is pinned strongly to one end. Nonlinear losses of trapped vortices in thick films under high-amplitude RF fields as functions of frequency, mean free path and pinning characteristics have been calculated.

Pathirana, Walive↗

Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions

This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.

54 ENVIRONMENTAL SCIENCES↗

Mitigation of spall fracture by evolving porosity

Ductile materials subject to impact loading conditions can undergo spall fracture when an incoming compressive stress wave reflects off interfaces and free surfaces as a tensile stress wave. Experimental observations suggest that in ductile materials, spall fracture is driven by the evolution of porosity. However, the presence of initial porosity in ductile materials also introduces plastic compressibility, which can attenuate the incoming compressive stress wave and, as a result, reduce the amplitude of the reflected tensile stress wave. This, in turn, can mitigate or delay spall fracture. In this work, we report on finite deformation finite element calculations that analyze the response of porous ductile materials subjected to impact loading conditions. Two sets of calculations are carried out, in the first set the material contains initial porosity values ranging from 0% to 5% while in the second set the material also undergoes stress-controlled porosity nucleation. Both sets of calculations are carried out for a wide range of imposed impact velocities. Here, our results show that porosity in ductile materials can, under certain circumstances, mitigate spall fracture by attenuating stress waves. Results correlating the effects of impact velocity, initial porosity, and porosity nucleation on spall fracture are presented and the underlying mechanisms are discussed.

36 MATERIALS SCIENCE↗