Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave prediction”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Diverse Manifestations of Electron-Phonon Coupling in a Kagome Superconductor

Recent angle-resolved photoemission spectroscopy (ARPES) experiments on the kagome metal CsV_{3}Sb_{5} revealed distinct multimodal dispersion kinks and nodeless superconducting gaps across multiple electron bands. The prominent photoemission kinks suggest a definitive coupling between electrons and certain collective modes, yet the precise nature of this interaction and its connection to superconductivity remain to be established. Here, employing the state-of-the-art ab initio many-body perturbation theory computation, we present direct evidence that electron-phonon (e-ph) coupling induces the multimodal photoemission kinks in CsV_{3}Sb_{5}, and profoundly, drives the nodeless s-wave superconductivity, showcasing the diverse manifestations of the e-ph coupling. Our calculations well capture the experimentally measured kinks and their fine structures, and reveal that vibrations from different atomic species dictate the multimodal behavior. Results from anisotropic GW-Eliashberg equations predict a phonon-mediated superconductivity with nodeless s-wave gaps, in excellent agreement with various ARPES and scanning tunneling spectroscopy measurements. Despite the universal origin of the e-ph coupling, the contributions of several characteristic phonon vibrations vary in different phenomena, highlighting a versatile role of e-ph coupling in shaping the low-energy excitations of kagome metals.

Electron-phonon coupling↗
NASA NTRS2011

Imaging the Moon's Core with Seismology

Constraining the structure of the lunar core is necessary to improve our understanding of the present-day thermal structure of the interior and the history of a lunar dynamo, as well as the origin and thermal and compositional evolution of the Moon. We analyze Apollo deep moonquake seismograms using terrestrial array processing methods to search for the presence of reflected and converted energy from the lunar core. Although moonquake fault parameters are not constrained, we first explore a suite of theoretical focal spheres to verify that fault planes exist that can produce favorable core reflection amplitudes relative to direct up-going energy at the Apollo stations. Beginning with stacks of event seismograms from the known distribution of deep moonquake clusters, we apply a polarization filter to account for the effects of seismic scattering that (a) partitions energy away from expected components of ground motion, and (b) obscures all but the main P- and S-wave arrivals. The filtered traces are then shifted to the predicted arrival time of a core phase (e.g. PcP) and stacked to enhance subtle arrivals associated with the Moon s core. This combination of filtering and array processing is well suited for detecting deep lunar seismic reflections, since we do not expect scattered wave energy from near surface (or deeper) structure recorded at varying epicentral distances and stations from varying moonquakes at varying depths to stack coherently. Our results indicate the presence of a solid inner and fluid outer core, overlain by a partial-melt-containing boundary layer (Table 1). These layers are consistently observed among stacks from four classes of reflections: P-to-P, S-to-P, P-to-S, and S-to-S, and are consistent with current indirect geophysical estimates of core and deep mantle properties, including mass, moment of inertia, lunar laser ranging, and electromagnetic induction. Future refinements are expected following the successful launch of the GRAIL lunar orbiter and SELENE 2 lunar lander missions.

Weber, Renee C.↗
NASA NTRS2012

Experimental Validation of a Forward Looking Interferometer for Detection of Clear Air Turbulence due to Mountain Waves

The Forward-Looking Interferometer (FLI) is an airborne sensor concept for detection and estimation of potential atmospheric hazards to aircraft. The FLI concept is based on high-resolution Infrared Fourier Transform Spectrometry technologies that have been developed for satellite remote sensing. The FLI is being evaluated for its potential to address multiple hazards, during all phases of flight, including clear air turbulence, volcanic ash, wake vortices, low slant range visibility, dry wind shear, and icing. In addition, the FLI is being evaluated for its potential to detect hazardous runway conditions during landing, such as wet or icy asphalt or concrete. The validation of model-based instrument and hazard simulation results is accomplished by comparing predicted performance against empirical data. In the mountain lee wave data collected in the previous FLI project, the data showed a damped, periodic mountain wave structure. The wave data itself will be of use in forecast and nowcast turbulence products such as the Graphical Turbulence Guidance and Graphical Turbulence Guidance Nowcast products. Determining how turbulence hazard estimates can be derived from FLI measurements will require further investigation.

Schaffner, Philip R.↗
NASA NTRS1990

Strong Langmuir turbulence generated by electron beams - Electric-field distributions and electron scattering

Strong turbulence and transit-time scattering theory are applied here to calculate the statistical distribution of intense Langmuir fields and the consequent beam scattering in plasma turbulence driven by an electron beam. The experimentally observed electric-field distributions are compared with predictions of strong-turbulence theory, concentrating on the wave levels, the Gaussian tail of the high-field distribution observed in one experiment, the arrest scale of collapse, and the fractional volume occupied by the highest fields. The Guassian form of the tail is confirmed, and the results imply that the collapse is arrested at a scale where the peak electrostatic energy density is of the same order as the thermal energy density. The theory of transit-time interactions is generalized to include relativistic particle dynamics and is applied to predict the scattering of the beam electrons in energy and angle as they pass through strong Langmuir turbulence. The results support the validity of the recently developed scaling theory of strong turbulence.

Robinson, P. A.↗
DOE OSTI2024

A-B Transition in Superfluid $^3$He and Cosmological Phase Transitions

First-order phase transitions in the very early universe are a prediction of many extensions of the Standard Model of particle physics and could provide the departure from equilibrium needed for a dynamical explanation of the baryon asymmetry of the Universe. They could also produce gravitational waves of a frequency observable by future space-based detectors such as the Laser Interferometer Space Antenna. All calculations of the gravitational wave power spectrum rely on a relativistic version of the classical nucleation theory of Cahn-Hilliard and Langer, due to Coleman and Linde. The high purity and precise control of pressure and temperature achievable in the laboratory made the first-order A to B transition of superfluid $^3$He ideal for test of classical nucleation theory. As Leggett and others have noted, the theory fails dramatically. The lifetime of the metastable A phase is measurable, typically of order minutes to hours, far faster than classical nucleation theory predicts. If the nucleation of B phase from the supercooled A phase is due to a new, rapid intrinsic mechanism that would have implications for first-order cosmological phase transitions as well as predictions for gravitational wave production in the early universe. Here we discuss studies of the A-B phase transition dynamics in $^3$He, both experimental and theoretical, and show how the computational technology for cosmological phase transition can be used to simulate the dynamics of the A-B transition, support the experimental investigations of the A-B transition in the QUEST-DMC collaboration with the goal of identifying and quantifying the mechanism(s) responsible for nucleation of stable phases in ultra-pure metastable quantum phases.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Spin Splitting Energy of Transition Metals: A New, More Affordable Wave Function Benchmark Method and Its Use to Test Density Functional Theory

Accurately predicting the spin splitting energy of chemical species is important for understanding their reactivity and magnetic properties, but it is very challenging, especially for molecules containing transition metals. One impediment to progress is the scarcity of accurate benchmark data. Here we report a set of calculations designed to yield reliable benchmarks for simple transition-metal complexes that can be used to test density functional methods that are affordable for large systems of more practical interest. Various wave function methods are tested against experiment for Fe 2+ , Fe 3+ , and Co 3+ , including CASSCF, CASPT2, CASPT3, MRCISD, MRCISD+Q, ACPF, AQCC, CCSD(T), and CASPT2/CCSD(T) and also a new method called CASPT2.5, which is performed by taking the average of the CASPT2 and CASPT3 energies. Furthermore, we find that MRCISD+Q, ACPF, and AQCC require smaller active spaces for good accuracy than are required by CASPT2 and CASPT3, and this aspect may be important for calculations on larger molecules; here we find that CASPT2.5 extrapolated to a complete basis set is the most suitable method—in terms of computational cost and in terms of accuracy on monatomic systems—and therefore we chose this method for molecular benchmarks. Then Kohn–Sham density functional calculations with 60 exchange-correlation functionals are tested for FeF 2 , FeCl 2 , and CoF 2 . We find that MN15-L, M06-SX, and revM06 have very good agreement with CASPT2.5 benchmarks in terms of both the spin splitting energy and the optimized geometry for each spin state. In addition, we recommend def2-TZVP as the most suitable basis set to perform density functional calculations for molecular spin splitting energies; extra polarization functions in the basis set do not help to increase the accuracy of the spin splitting energy in KS calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Hydrodynamic Stability of Liquid-Propellant Combustion: Landau's Problem Revisited

Hydrodynamic, or Landau, instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. As its name suggests, it stems from hydrodynamic effects connected with thermal expansion across the reaction region. In the context of liquid-propellant combustion, the classical models that originally predicted this phenomenon have been extended to include the important effects that arise from a dynamic dependence of the burning rate on the local pressure and temperature fields. Thus, the onset of Landau instability has now been shown to occur for sufficiently small negative values of the pressure sensitivity of the burning rate, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. It has also been shown that the onset of instability occurs for decreasing values of the disturbance wave number as the gravitational-acceleration parameter decreases. Consequently, in an appropriate weak-gravity limit, Landau instability becomes a long-wave phenomena associated with the formation of large cells on the liquid-propellant surface. Additionally, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. This instability occurs for sufficiently large negative values of the pressure sensitivity, and is enhanced by increasing values of the burning-rate temperature sensitivity. It is further shown that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the (disturbance-wave number, pressure-sensitivity) plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗
NASA NTRS1997

Validation of an Accurate Three-Dimensional Helical Slow-Wave Circuit Model

The helical slow-wave circuit embodies a helical coil of rectangular tape supported in a metal barrel by dielectric support rods. Although the helix slow-wave circuit remains the mainstay of the traveling-wave tube (TWT) industry because of its exceptionally wide bandwidth, a full helical circuit, without significant dimensional approximations, has not been successfully modeled until now. Numerous attempts have been made to analyze the helical slow-wave circuit so that the performance could be accurately predicted without actually building it, but because of its complex geometry, many geometrical approximations became necessary rendering the previous models inaccurate. In the course of this research it has been demonstrated that using the simulation code, MAFIA, the helical structure can be modeled with actual tape width and thickness, dielectric support rod geometry and materials. To demonstrate the accuracy of the MAFIA model, the cold-test parameters including dispersion, on-axis interaction impedance and attenuation have been calculated for several helical TWT slow-wave circuits with a variety of support rod geometries including rectangular and T-shaped rods, as well as various support rod materials including isotropic, anisotropic and partially metal coated dielectrics. Compared with experimentally measured results, the agreement is excellent. With the accuracy of the MAFIA helical model validated, the code was used to investigate several conventional geometric approximations in an attempt to obtain the most computationally efficient model. Several simplifications were made to a standard model including replacing the helical tape with filaments, and replacing rectangular support rods with shapes conforming to the cylindrical coordinate system with effective permittivity. The approximate models are compared with the standard model in terms of cold-test characteristics and computational time. The model was also used to determine the sensitivity of various circuit parameters including typical manufacturing dimensional tolerances and support rod permittivity. By varying the circuit parameters of an accurate model using MAFIA, these sensitivities can be computed for manufacturing concerns, and design optimization previous to fabrication, thus eliminating the need for costly experimental iterations. Several variations were made to a standard helical circuit using MAFIA to investigate the effect that variations on helical tape and support rod width, metallized loading height and support rod permittivity, have on TWT cold-test characteristics.

Kory, Carol L.↗
NASA NTRS1990

Laser driven launch vehicles for continuous access to space

The availability of megawatt laser systems in the next century will make laser launch systems from ground to orbit feasible and useful. Systems studies indicate launch capabilities of 1 ton payload per gigawatt laser power. Recent research in ground to orbit laser propulsion has emphasized laser supported detonation wave thrusters driven by repetitively pulsed infrared lasers. In this propulsion concept each laser repetition cycle consists of two pulses. A lower energy first pulse is used to vaporize a small amount of solid propellant and then after a brief expansion period, a second and higher energy laser pulse is used to drive a detonation wave through the expanded vapor. The results are reported of numerical studies comparing the detonation wave properties of various candidate propellants, and the simulation of thruster performance under realistic conditions. Experimental measurements designed to test the theoretical predictions are also presented. Measurements are discussed of radiance and opacity in absorption waves, and mass loss and momentum transfer. These data are interpreted in terms of specific impulse and energy conversion efficiency.

Rollins, C. J.↗
NASA NTRS2010

The Prediction of Broadband Shock-Associated Noise from Dualstream and Rectangular Jets Using RANS CFD

Supersonic jets operating off-design produce broadband shock-associated noise. Broadband shock-associated noise is characterized by multiple broadband peaks in the far-field and is often the dominant source of noise towards the sideline and upstream direction relative to the jet axis. It is due to large scale coherent turbulence structures in the jet shear layers interacting with the shock cell structure. A broadband shock-associated noise model recently developed by the authors predicts this noise component from solutions to the Reynolds averaged Navier-Stokes equations using a two-equation turbulence model. The broadband shock-associated noise model is applied to dualstream and rectangular nozzles operating supersonically, heated, and off-design. The dualstream jet broadband shock-associated noise predictions are conducted for cases when the core jet is supersonic and the fan jet is subsonic, the core jet is subsonic and the fan jet is supersonic, and when both jet streams operate supersonically. Rectangular jet predictions are shown for a convergent-divergent nozzle operating both over- and under-expanded for cold and heated conditions. The original model implementation has been heavily modified to make accurate predictions for the dualstream jets. It is also argued that for over-expanded jets the oblique shock wave attached to the nozzle lip contributes little to broadband shock-associated noise. All predictions are compared with experiments.

Miller, Steven A.↗
DOE OSTI2023

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗
NASA NTRS2010

Assessment of the Breakup of the Antarctic Polar Vortex in Two New Chemistry-Climate Models

Successful simulation of the breakup of the Antarctic polar vortex depends on the representation of tropospheric stationary waves at Southern Hemisphere middle latitudes. This paper assesses the vortex breakup in two new chemistry-climate models (CCMs). The stratospheric version of the UK Chemistry and Aerosols model is able to reproduce the observed timing of the vortex breakup. Version 2 of the Goddard Earth Observing System (GEOS V2) model is typical of CCMs in that the Antarctic polar vortex breaks up too late; at 10 hPa, the mean transition to easterlies at 60 S is delayed by 12-13 days as compared with the ERA-40 and National Centers for Environmental Prediction reanalyses. The two models' skill in simulating planetary wave driving during the October-November period accounts for differences in their simulation of the vortex breakup, with GEOS V2 unable to simulate the magnitude and tilt of geopotential height anomalies in the troposphere and thus underestimating the wave driving. In the GEOS V2 CCM the delayed breakup of the Antarctic vortex biases polar temperatures and trace gas distributions in the upper stratosphere in November and December.

Hurwitz, M. M.↗
NASA NTRS1986

Analysis of rainfall over northern Peru during El Nino: A PCDS application

In an examination of GOES satellite data during the 1982 through 1983 El Nino period, the appearance of lee wave cloud patterns was revealed. A correlation was hypothesized relating an anomalous easterly flow across the Andes with the appearance of these wave patterns and with the subsequent onset of intense rainfall. The cloud patterns are belived to be associated with the El Nino period and could be viewed as precursors to significant changes in weather patterns. The ultimate goal of the researchers will be the ability to predict occurrences of rainstorms associated with the appearance of lee waves and related cloud patterns as harbingers of destruction caused by flooding, huaycos, and other catastrophic consequences of heavy and abnormal rainfall. Rainfall data from about 70 stations in northern Peru from 1980 through 1984 were formatted to be utilized within the Pilot Climate Data System (PCDS). This time period includes the 1982 through 1983 El Nino period. As an example of the approach, a well-pronounced lee wave pattern was shown from a GOES satellite image of April 4, 1983. The ground truth data were then displayed via the PCDS to graphically demonstrate the increase in intensity and areal distribution of rainfall in the northern Peruvian area in the next 4 to 5 days.

Goldberg, R.↗
NASA NTRSDate not supplied

NASA’s Orion Crew Module Seakeeping Simulation and Test Comparisons in Ocean Wave Environments

The NASA Orion Crew Module (CM) post-landing seakeeping and uprighting capabilities are being evaluated across a variety of ocean wave conditions. This paper provides an overview of quarter-scale and full-scale test data and how it compares with the CM seakeeping dynamic simulation results evaluated using Wave Energy Convertor Simulator (WEC-Sim).This includes a discussion of the status of the dynamic modeling effort, correlation of model inputs with the wave field from the test, and a comparison of the results against the physical test data with emphasis on the full-scale results. Analysis in the frequency domain showed similar response characteristics between the model and the test and shows energy peaks corresponding to both the natural response and the forced response due to the wave environment. Time domain comparisons showed that the linear WEC-Sim model under predicted the amplitude of the pitch and roll motion of the test article, with improved matching for smaller, longer crested waves. An attempt was made to correlate the pitch/roll response to incident wave slope to account for a possible delay in the BTA’s response. The results did not show a strong relationship in the response characteristics to incident wave slope, but better correlation was seen by comparing responses in the frequency domain. Evaluation of the simulation’s ability to match test results aided in understanding the current fidelity of the CM model in WEC-Sim and identifying areas for model improvement

Abigail Lockard↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
DOE OSTI2023

Gravitational Waves from Rapid Structure Formation on Microscopic Scales before Matter-Radiation Equality

We propose a new scenario for the formation of asteroid-mass primordial black holes (PBHs). Our mechanism is based on the annihilation of the string-wall network associated with thebreaking of a U(1) global symmetry into a discrete ZN symmetry. If the potential hasmultiple local minima (N > 1) the network is stable, and the annihilation is guaranteed by a biasamong the different vacua. The collapse of the string-wall network is accompanied bycatastrogenesis, a large production of pseudo-Goldstone bosons (pGBs) — e.g. axions,ALPs, or majorons — gravitational waves, and PBHs. If pGBs rapidly decay into products thatthermalize, as predicted e.g. in the high-quality QCD axion and heavy majoron models, they do notcontribute to the dark matter population, but we show that PBHs can constitute 100% of the darkmatter. The gravitational wave background produced by catastrogenesis with heavy unstable axions,ALPs, or majorons could be visible in future interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Analysis and prediction of ocean swell using instrumented buoys

During the period 20-23 September 1990, the remnants of Supertyphoon Flo moved into the central North Pacific Ocean with sustained wind speeds of 28 m/s. The strong wind and large fetch area associated with this storm generated long-period swell that propagated to the west coast of North America. National Data Buoy Center moored-buoy stations, located in a network that ranged from the Gulf of Alaska to the California Bight, provided wave spectral estimates of the swell from this storm. The greatest dominant wave periods measured were approximately 20-25 s, and significant wave heights measured ranged from 3 to 8 m. Wave spectra from an array of three nondirectional buoys are used to find the source of the long-period swell. Directional wave spectra from a heave-pitch-roll buoy are also used to make an independent estimate of the source of the swell. The ridge-line method, using time-frequency contour plots of wave spectral energy density, is used to determine the time of swell generation, which is used with the appropriate surface pressure analysis to infer the swell generation area. The diagnosed sources of the swell are also compared with nowcasts from the Global Spectral Ocean Wave Model of the Fleet Numerical Oceanography Center. A simple method of predicting the propagation of ocean swell, by applying a simple kinematic model of wave propagation to the estimated point and time source, is demonstrated.

Mettlach, Theodore↗
DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗