Interplay of martensitic phase transformation and plastic slip in polycrystals

Inspired by key experimental and analytical results regarding Shape Memory Alloys (SMAs), we propose a modelling framework to explore the interplay between martensitic phase transformations and plastic slip in polycrystalline materials, with an eye towards computational efficiency. The resulting framework uses a convexified potential for the internal energy density to capture the stored energy associated with transformation at the meso-scale, and introduces kinetic potentials to govern the evolution of transformation and plastic slip. The framework is novel in the way it treats plasticity on par with transformation.

We implement the framework in the setting of anti-plane shear, using a staggered implicit/explict update: we first use a Fast-Fourier Transform (FFT) solver based on an Augmented Lagrangian formulation to implicitly solve for the full-field displacements of a simulated polycrystal, then explicitly update the volume fraction of martensite and plastic slip using their respective stick-slip type kinetic laws. We observe that, even in this simple setting with an idealized material comprising four martensitic variants and four slip systems, the model recovers a rich variety of SMA type behaviors. We use this model to gain insight into the isothermal behavior of stress-stabilized martensite, looking at the effects of the relative plastic yield strength, the memory of deformation history under non-proportional loading, and several others.

We extend the framework to the generalized 3-D setting, for which the convexified potential is a lower bound on the actual internal energy, and show that the fully implicit discrete time formulation of the framework is governed by a variational principle for mechanical equilibrium. We further propose an extension of the method to finite deformations via an exponential mapping. We implement the generalized framework using an existing Optimal Transport Mesh-free (OTM) solver. We then model the $\alpha$--$\gamma$ and $\alpha$--$\varepsilon$ transformations in pure iron, with an initial attempt in the latter to account for twinning in the parent phase. We demonstrate the scalability of the framework to large scale computing by simulating Taylor impact experiments, observing nearly linear (ideal) speed-up through 256 MPI tasks. Finally, we present preliminary results of a simulated Split-Hopkinson Pressure Bar (SHPB) experiment using the $\alpha$--$\varepsilon$ model.

Slender-body hypervelocity boundary-layer instability

With novel application of optical techniques, the slender-body hypervelocity boundary-layer instability is characterized in the previously unexplored regime where thermo-chemical effects are important. Narrowband disturbances (500-3000~kHz) are measured in boundary layers with edge velocities of up to 5~km/s at two points along the generator of a 5 degree half angle cone. Experimental amplification factor spectra are presented. Linear stability and PSE analysis is performed, with fair prediction of the frequency content of the disturbances; however, the analysis over-predicts the amplification of disturbances. The results of this work have two key implications: 1) the acoustic instability is present and may be studied in a large-scale hypervelocity reflected-shock tunnel, and 2) the new data set provides a new basis on which the instability can be studied.

The instability of fluids with time dependent heating

The stability of a fluid having a non-uniform temperature stratification is examined analytically for the response of infinitesimal disturbances. The growth rates of disturbances have been established for a semi-infinite fluid for Rayleigh numbers of 10³, 10⁴, and 10⁵ and for Prandtl numbers of 7.0 and 0.7.

The critical Rayleigh number for a semi-infinite fluid, based on the effective fluid depth, is found to be 32, while it is shown that for a finite fluid layer the critical Rayleigh number depends on the rate of heating. The minimum critical Rayleigh number, based on the depth of a fluid layer, is found to be 1340.

The stability of a finite fluid layer is examined for two special forms of heating. The first is constant flux heating, while in the second, the temperature of the lower surface is increased uniformly in time. In both cases, it is shown that for moderate rates of heating the critical Rayleigh number is reduced, over the value for very slow heating, while for very rapid heating the critical Rayleigh number is greatly increased. These results agree with published experimental observations.

The question of steady, non-cellular convection is given qualitative consideration. It is concluded that, although the motion may originate from infinitesimal disturbances during non-uniform heating, the final flow field is intrinsically non-linear.

Supernovae explosions induced by pair production instability

Stars with a core mass greater than about 30 M_⊙ become dynamically unstable due to electron-positron pair production when their central temperature reaches 1.5-2.0 x 10^{9 0}K. The collapse and subsequent explosion of stars with core masses of 45, 52, and 60 M_⊙ is calculated. The range of the final velocity of expansion (3,400 – 8,500 km/sec) and of the mass ejected (1 – 40 M_⊙) is comparable to that observed for type II supernovae.

An implicit scheme of hydrodynamic difference equations (stable for large time steps) used for the calculation of the evolution is described.

For fast evolution the turbulence caused by convective instability does not produce the zero entropy gradient and perfect mixing found for slower evolution. A dynamical model of the convection is derived from the equations of motion and then incorporated into the difference equations.

Experimental and theoretical studies of the electrothermal instability in a nonequilibrium plasma

Experimental and theoretical studies have been made of the electrothermal waves occurring in a nonequilibrium MHD plasma. These waves are caused by an instability that occurs when a plasma having a dependence of conductivity on current density is subjected to crossed electric and magnetic fields. Theoretically, these waves were studied by developing and solving the equations of a steady, one-dimensional nonuniformity in electron density. From these nonlinear equations, predictions of the maximum amplitude and of the half width of steady waves could be obtained. Experimentally, the waves were studied in a nonequilibrium discharge produced in a potassium-seeded argon plasma at 2000°K and 1 atm. pressure. The behavior of such a discharge with four different configurations of electrodes was determined from photographs, photomultiplier measurements, and voltage probes. These four configurations were chosen to produce steady waves, to check the stability of steady waves, and to observe the manifestation of the waves in a MHD generator or accelerator configuration.

Steady, one-dimensional waves were found to exist in a number of situations, and where they existed, their characteristics agreed with the predictions of the steady theory. Some extensions of this theory were necessary, however, to describe the transient phenomena occurring in the inlet region of a discharge transverse to the gas flow. It was also found that in a discharge away from the stabilizing effect of the electrodes, steady waves became unstable for large Hall parameters. Methods of prediction of the effective electrical conductivity and Hall parameter of a plasma with nonuniformities caused by the electrothermal waves were also studied. Using these methods and the values of amplitude predicted by the steady theory, it was found that the measured decrease in transverse conductivity of a MHD device, 50 per cent at a Hall parameter of 5, could be accounted for in terms of the electrothermal instability.

A Continuum Model for Slip-Twinning Interactions in Magnesium and Magnesium Alloys

Due to their high specific strength and low density, magnesium and magnesium-based alloys have gained great technological importance in recent years. However, their underlying hexagonal crystal structure furnishes Mg and its alloys with a complex mechanical behavior because of their comparably smaller number of energetically favorable slip systems. Besides the commonly studied slip mechanism, another way to accomplish general deformation is through the additional mechanism of deformation-induced twinning. The main aim of this thesis research is to develop an efficient continuum model to understand and ultimately predict the material response resulting from the interaction between these two mechanisms.

The constitutive model we present is based on variational constitutive updates of plastic slips and twin volume fractions and accounts for the related lattice reorientation mechanisms. The model is applied to single- and polycrystalline pure magnesium. We outline the finite-deformation plasticity model combining basal, pyramidal, and prismatic dislocation activity as well as a convexification based approach for deformation twinning. A comparison with experimental data from single-crystal tension-compression experiments validates the model and serves for parameter identification. The extension to polycrystals via both Taylor-type modeling and finite element simulations shows a characteristic stress-strain response that agrees well with experimental observations for polycrystalline magnesium. The presented continuum model does not aim to represent the full details of individual twin-dislocation interactions, yet it is sufficiently efficient to allow for finite element simulations while qualitatively capturing the underlying microstructural deformation mechanisms.

Dynamic response of hysteretic systems with application to a system containing limited slip

A general class of single degree of freedom systems possessing rate-independent hysteresis is defined. The hysteretic behavior in a system belonging to this class is depicted as a sequence of single-valued functions; at any given time, the current function is determined by some set of mathematical rules concerning the entire previous response of the system. Existence and uniqueness of solutions are established and boundedness of solutions is examined.

An asymptotic solution procedure is used to derive an approximation to the response of viscously damped systems with a small hysteretic nonlinearity and trigonometric excitation. Two properties of the hysteresis loops associated with any given system completely determine this approximation to the response: the area enclosed by each loop, and the average of the ascending and descending branches of each loop.

The approximation, supplemented by numerical calculations, is applied to investigate the steady-state response of a system with limited slip. Such features as disconnected response curves and jumps in response exist for a certain range of system parameters for any finite amount of slip.

To further understand the response of this system, solutions of the initial-value problem are examined. The boundedness of solutions is investigated first. Then the relationship between initial conditions and resulting steady-state solution is examined when multiple steady-state solutions exist. Using the approximate analysis and numerical calculations, it is found that significant regions of initial conditions in the initial condition plane lead to the different asymptotically stable steady-state solutions.