Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orientational order”

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 541 records · Page 30

Structural Analysis Approach for the Defense Programs Package 3 (DPP-3)

The Pacific Northwest National Laboratory (PNNL) is the design authority for a new Type B hazardous materials transportation package designated as the Defense Programs Package 3 (DPP-3) for the U.S. Department of Energy (DOE) National Nuclear Security Administration (NNSA). The DPP-3 has been developed using similar materials and fabrication methods employed in previous U.S. Nuclear Regulatory Commission (NRC), DOE, and NNSA certified packages. The DPP-3 design criteria are derived from the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (BPVC), NNSA guidance and NRC regulatory guides in order to safely and securely transport a variety of payloads. Final regulatory approval by the NNSA will require physical testing to demonstrate that the containment vessel (CV) remains leaktight after enduring the entire regulatory testing sequence prescribed in Title 10 of the Code of Federal Regulations Part 71 (10 CFR 71). In order to gain confidence that the DPP-3 will remain leaktight after testing, the DPP-3 has been structurally analyzed using the Finite Element Analysis (FEA) software LS-DYNA. The FEA analyses serve two general purposes: first, they aid in design and development of the package, and second, they advise as to which drop orientations are expected to cause the most damage during regulatory testing. This paper will discuss how the design criteria are incorporated into analytical techniques needed to evaluate the FEA structural simulation results for 10 CFR 71 conditions to give confidence the DPP-3 testing campaign will be successful.

Sakalaukus Jr., Peter J.↗

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermal conductivity tensor of γ and ε -hexanitrohexaazaisowurtzitane as a function of pressure and temperature

Using reverse non-equilibrium molecular dynamics simulations, we have determined the dependences on temperature and pressure of the thermal conductivity tensors for the monoclinic γ and ε polymorphs of hexanitrohexaazaisowurtzitane (HNIW or CL20). In this work, a recently developed non-reactive force field [X. Bidault and S. Chaudhuri, RSC Adv. 9, 39649–39661 (2019)], designed to study polymorphism and phase transitions in CL20, is employed. The effects of temperature and pressure are investigated between 200 and 500 K and up to 0.5 GPa for γ-CL20 and 2 GPa for ε-CL20. In order to obtain the full thermal conductivity tensor, $κ_{ij}$, for the monoclinic crystals, four distinct heat propagation directions are used. We find that $κ_{ij}$ for both polymorphs is more isotropic than for other energetic molecular crystals, including α- and γ-RDX, β-HMX, and PETN, with a maximum difference of 9.8% between orientations observed at 300 K and 0 GPa for γ-CL20 and a maximum difference of 4.8% for ε-CL20. The average thermal conductivity, $\bar{κ}$, of ε-CL20 is 6.4% larger than that of γ-CL20 at 300 K and 0 GPa. Analytic linear functions of the inverse temperature and the pressure are provided, which fit the data well and can be used to predict the thermal conductivity of both polymorphs for any orientation, pressure, and temperature in and around the fitting range. Our predictions agree reasonably well with the limited available experimental data, for which the polymorph type is unknown.

36 MATERIALS SCIENCE↗

Bulk photospin effect: Calculation of electric spin susceptibility to second order in an electric field

Here, we compute the electric spin susceptibility of Bloch electrons with spin-orbit coupling to second order. We find that it is possible to generate a nonequilibrium spin polarization in the bulk of nonmagnetic inversion-symmetric materials using linearly polarized electric fields, but the process depends on interband coherence and produces heating. It may be possible to avoid heating with circular polarization in certain scenarios. The standard Edelstein effect and spin orientation effects are recovered in appropriate limits within the formalism. Finally, the electric spin susceptibility of metals has contributions proportional to spin multipole moments of the Fermi sea that dominate the low frequency spin response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evolution of magnetism in the magnetic topological semimetal NdSb x Te 2–x+δ

Magnetic topological semimetals LnSbTe (Ln = Lanthanide) have attracted intensive attention because of the presence of interplay between magnetism, topological, and electron correlations depending on the choices of magnetic Ln elements. Recently, varying Sb-Te composition has been found to effectively control the electronic and magnetic states in LnSb x Te 2-x . Here, with this motivation, we report the evolution of magnetic properties with Sb-Te substitution in NdSb x Te 2-x+δ , (0 ≤ x ≤ 1). Our work reveals the interesting non-monotonic change in magnetic ordering temperature with varying composition stoichiometry. In addition, reducing the Sb content x drives the reorientation of moments from in-plane (ab-plane) to out-of-plane (c-axis) direction that results in the distinct magnetic structures for two end compounds NdTe 2 (x = 0) and NdSbTe (x = 1). Furthermore, the moment orientation in NdSb x Te 2-x+δ is also found to be strongly tunable upon application of weak magnetic field, leading to rich magnetic phases depending on the composition stoichiometry, temperature, and magnetic field. Such strong tuning of magnetism in this material establishes it as a promising platform for investigating tunable topological states and correlated topological physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Printing 2D Conjugated Polymer Monolayers and Their Distinct Electronic Properties

Abstract Recently, 2D monolayer films of conjugated polymers have gained increasing attention owing to the preeminence of 2D inorganic films that exhibit unique optoelectronic and mechanical properties compared to their bulk analogs. Despite numerous efforts, crystallization of semiconducting polymers into highly ordered 2D monolayer films still remains challenging. Herein, a dynamic‐template‐assisted meniscus‐guided coating is utilized to fabricate continuous, highly ordered 2D monolayer films of conjugated polymers over a centimeter scale with enhanced backbone π–π stacking. In contrast, monolayer films printed on solid substrates confer upon the 1D fiber networks strong alkyl side‐chain stacking at the expense of backbone packing. From single‐layers to multilayers, the polymer π‐stacks change from edge‐on to bimodal orientation as the film thickness reaches ≈20 nm. Spectroscopic and cyclic voltammetry analysis reveals an abrupt increase in J‐aggregation and absorption coefficient and a decrease in bandgap and highest occupied molecular orbital level until critical thickness, possibly arising from the straightened polymer backbone. This is corroborated by an abrupt increase in hole mobility with film thickness, reaching a maximum of 0.7 cm 2 V −1 s −1 near the critical thickness. Finally, fabrication of chemical sensors incorporating polymer films of various thicknesses is demonstrated, and an ultrahigh sensitivity of the ≈7 nm thick ultrathin film (bilayers) to 1 ppb ammonia is shown.

Kafle, Prapti↗

Molecular Dynamics Investigation of Nanoscale Hydrophobicity of Polymer Surfaces: What Makes Water Wet?

The wettability of a polymer surface–related to its hydrophobicity or tendency to repel water–can be crucial for determining its utility, such as for a coating or a purification membrane. While wettability is commonly associated with the macroscopic measurement of a contact angle between surface, water, and air, the molecular physics that underlie these macroscopic observations are not fully known, and anticipating relative behavior of different polymers is challenging. To address this gap in molecular-level understanding, we use molecular dynamics simulations to investigate and contrast interactions of water with six chemically distinct polymers: polytetrafluoroethylene, polyethylene, polyvinyl chloride, poly(methyl methacrylate), Nylon 66, and polyvinyl alcohol. We show that several prospective quantitative metrics for hydrophobicity agree well with experimental contact angles. Moreover, the behavior of water in proximity to these polymer surfaces can be distinguished with analysis of interfacial water dynamics, extent of hydrogen bonding, and molecular orientation–even when macroscopic measures of hydrophobicity are similar. The predominant factor dictating wettability is found to be the extent of hydrogen bonding between polymer and water, but the precise manifestation of hydrogen bonding and its impact on surface-water structure varies. In the absence of hydrogen bonding, other molecular interactions and polymer mechanics control hydrophobic ordering. Furthermore, these results provide new insights into how polymer chemistry specifically impacts water-polymer interactions and translates to surface hydrophobicity. Such factors may facilitate the design or processing of polymer surfaces to achieve targeted wetting behavior, and presented analyses can be useful in studying the interfacial physics of other systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling the symmetry of cadmium arsenide films by epitaxial strain

Epitaxial strains offer unique opportunities to obtain topological states in thin films and heterostructures that do not exist in their bulk counterparts. Here, we investigate the point group symmetries of coherently strained films of cadmium arsenide (Cd 3 As 2 ), a prototype three-dimensional Dirac semimetal, by convergent beam electron diffraction. We report a loss of the fourfold rotational axis and adoption of the orthorhombic mmm point group in (112)-oriented films under biaxial compressive stress. (001)-oriented Cd 3 As 2 films that are under a small biaxial tensile stress retain the fourfold rotational symmetry that protects the bulk nodes but adopt the non-centrosymmetric 4mm point group symmetry. This, in turn, suggests that (001) films adopt a different crystal structure in biaxial tension, one that differs in the arrangement of the ordered Cd vacancies that are an inherent feature of the crystal structure of Cd 3 As 2 and that are key to its nodal electronic structure. Density functional theory calculations confirm the experimental findings of the stability of the non-centrosymmetric structure under biaxial tension, whereas the centrosymmetric structure is stable under biaxial compression. The results show that bulk Cd 3 As 2 is already close to structural instability and showcase the extraordinary tunability of the topological states of Cd 3 As 2 .

36 MATERIALS SCIENCE↗

High-Entropy Alloys for Accelerator Beam Window Applications

Development of novel high-entropy alloys (HEAs) is currently underway for potential use as beam windows in future multi-megawatt target systems at Fermilab. HEAs encompass a new class of materials with a vast design space allowing for material properties to be tailored for particular applications and to potentially offer improved resistance to beam-induced radiation damage and thermal shock effects. The alloy systems being studied consist of several compositions of AlCoCrMnTiV with 4 6 component elements. These alloys are all predicted by CALPHAD simulation to have a single-phase BCC crystal structure and low density, with some compositions displaying ordered, nanoscale precipitates. This presentation will briefly discuss alloy design and synthesis before giving a detailed description of the characterization studies of these HEAs in both the pristine state and post-irradiation by low-energy heavy ions to high damage levels. Electron microscopy techniques to quantify elemental homogeneity and composition, determine grain size, shape, and orientation, and quantify lattice parameters, defect structures and precipitate phases are all being used to study alloy microstructures. Mechanical properties of the alloys at the microscale will be reported. The evolution of these properties as a function of radiation damage will also be described. Bulk thermal characteristics of these HEAs have been tested to measure specific heat capacity and coefficient of thermal expansion as a function of temperature. To determine bulk tensile properties a miniature tensile testing apparatus is under development; it s commissioning will be covered briefly. The talk will conclude with our plans for alloy down-selection.

Burleigh, A. [Fermilab]↗

The carboxylation status of osteocalcin has important consequences for its structure and dynamics

The carboxylation status of Osteocalcin (Ocn) not only influences formation and structure in bones but also has important endocrine functions affecting energy metabolism and expenditure. In this study, the role of γ-carboxylation of the glutamate residues in the structure-dynamics-function relationship in Ocn is investigated. Three forms of Ocn, differentially carboxylated at the Glu-17, 21 and 24 residues, along with a mutated form of Ocn carrying Glu/Ala mutations, are modeled and simulated using molecular dynamics (MD) simulation in the presence of calcium ions. Characterization of the global conformational dynamics of Ocn, described in terms of the orientational variations within its 3-helical domain, highlights large structural variations in the non-carboxylated osteocalcin (nOcn). The bi-carboxylated Ocn (bOcn) and tri-carboxylated (tOcn) species, in contrast, display relatively rigid tertiary structures, with the dynamics of most regions strongly correlated. Radial distribution functions calculated for both bOcn and tOcn show long-range ordering of the calcium ion distribution around the carboxylated glutamate (γGlu) residues, likely playing an important role in promoting stability of these Ocns. Additionally, the same calcium ions are observed to coordinate with neighboring γGlu, better shielding their negative charges and in turn stabilizing these systems more than do the singly coordinating calcium ions observed in the case of nOcn. bOcn is also found to exhibit a more helical C-terminal structure, that has been shown to activate its cellular receptor GPRC6A, highlighting the allosteric role of Ocn carboxylation in modulating the stability and binding potential of the active C-terminal. Furthermore, the carboxylation status of Ocn as well and its calcium coordination appear to have a direct influence on Ocn structure and dynamics, possibly leading to the known differences in Ocn biological function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Control of liquid crystals combining surface acoustic waves, nematic flows, and microfluidic confinement

The optical properties of liquid crystals serve as the basis for display, diagnostic, and sensing technologies. Such properties are generally controlled by relying on electric fields. In this work, we investigate the effects of microfluidic flows and acoustic fields on the molecular orientation and the corresponding optical response of nematic liquid crystals. Several previously unknown structures are identified, which are rationalized in terms of a state diagram as a function of the strengths of the flow and the acoustic field. The new structures are interpreted by relying on calculations with a free energy functional expressed in terms of the tensorial order parameter, using continuum theory simulations in the Landau-de Gennes framework. Taken together, the findings presented here offer promise for the development of new systems based on combinations of sound, flow, and confinement.

74 ATOMIC AND MOLECULAR PHYSICS↗

Autoregressive long-horizon prediction of plasma edge dynamics *

Accurate modeling of scrape-off layer (SOL) and divertor-edge dynamics is vital for designing plasma-facing components in fusion devices. High-fidelity edge fluid/neutral codes such as SOLPS-ITER capture SOL physics with high accuracy, but their computational cost limits broad parameter scans and long transient studies. We present transformer-based, autoregressive surrogates for efficient prediction of 2D, time-dependent plasma edge state fields. Trained on SOLPS-ITER spatiotemporal data for the KSTAR tokamak, the surrogates forecast electron temperature, electron density, and radiated power over extended horizons. We evaluate model variants trained with increasing autoregressive horizons (1–100 steps) on short- and long-horizon prediction tasks. Longer-horizon training systematically improves rollout stability and mitigates error accumulation, enabling stable predictions over hundreds to thousands of steps and reproducing key dynamical features such as the motion of high-radiation regions. Measured end-to-end wall-clock times show the surrogate is orders of magnitude faster than SOLPS-ITER, enabling rapid parameter exploration. Prediction accuracy degrades when the surrogate enters physical regimes not represented in the training dataset, motivating future work on data enrichment and physics-informed constraints. Overall, this approach provides a fast, accurate surrogate for computationally intensive plasma edge simulations, supporting rapid scenario exploration, control-oriented studies, and progress toward real-time applications in fusion devices.

autoregressive deep learning↗

Enhanced mechanical property evaluation using innovative data analytics capability

This report focuses on efforts to improve Multiphysics Object-Oriented Simulation Environ-ment (MOOSE) for mechanical property evaluation using data analytics. These efforts involveimprovements to the stochastic tools module (STM) for stochastic simulations of MOOSE mul-tiphysics model and the development of inverse optimization capabilities. The report gives anoverview of the STM and describes recent updates to its core capabilities and theory on itsreduced-order model (ROM) schemes. Examples are also provided showing the impact of theseupdates and exhibits the usefulness of ROMs. An overview of the gradient based inverse opti-mziation algorithms are given along with examples of their application to source identification.Inverse optimization will provide a new methodology in STM for fitting model parameters toexperimental data.

97 MATHEMATICS AND COMPUTING↗

Deep-Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile

Polyacrylonitrile (PAN) is an important commercial polymer, bearing atactic stereochemistry resulting from nonselective radical polymerization. As such, an accurate, fundamental understanding of governing interactions among PAN molecular units is indispensable for advancing the design principles of final products at reduced processability costs. While ab initio molecular dynamics (AIMD) simulations can provide the necessary accuracy for treating key interactions in polar polymers, such as dipole–dipole interactions and hydrogen bonding, and analyzing their influence on the molecular orientation, their implementation is limited to small molecules only. Herein, we show that the neural network interatomic potentials (NNIPs) that are trained on the small-scale AIMD data (acquired for oligomers) can be efficiently employed to examine the structures and properties at large scales (polymers). NNIP provides critical insight into intra- and interchain hydrogen-bonding and dipolar correlations and accurately predicts the amorphous bulk PAN structure validated by modeling the experimental X-ray structure factor. Furthermore, the NNIP-predicted PAN properties, such as density and elastic modulus, are in good agreement with their experimental values. Overall, the trend in the elastic modulus is found to correlate strongly with the PAN structural orientations encoded in the Hermans orientation factor. In conclusion, this study enables the ability to predict the structure–property relations for PAN and analogues with sustainable ab initio accuracy across scales.

36 MATERIALS SCIENCE↗

A control-oriented combustion model framework for compression ignition engines operating on low-reactivity fuel

This work focuses on zero-dimensional modeling of the heat release rate in a compression ignition engine operating on gasoline-like fuels. Due to the properties of gasoline, such as high volatility and longer ignition delay than diesel, the injection strategies can vary significantly from the operation with conventional diesel fuel. Different injection strategies are commonly used to achieve varying degrees of in-cylinder stratification in order to shape the combustion event and maximize efficiency. The proposed zero-dimensional combustion model was developed to account for the different stages in combustion caused by the fuel stratification. As the ignition delay model is an integral part of the entire combustion process and significantly affects the prediction accuracy, special attention has been paid to local phenomena influencing ignition delay. A one-dimensional spray model by Musculus and Kattke was employed in conjunction with a Lagrangian tracking approach in order to estimate the local air–fuel ratio within the spray tip, as a proxy for reactivity. The local air–fuel ratio, in-cylinder temperature and pressure were used in an integral fashion to estimate the ignition delay. Heat release rates were modeled using first-order non-linear differential equations. The proposed combustion model was validated against experimental data of a heavy-duty compression ignition engine with up to three injection events at mostly 1038 r/min and 14 bar brake mean effective pressure. Further validation of the model was carried out at other engine loads and speeds. Model prediction errors in CA50 of less than 1 °CA across all conditions were found. Modeling results of other combustion metrics such as combustion duration and indicated mean effective pressure are also highly satisfactory. In addition, the model has been shown to be capable of estimating the ringing intensity for most conditions.

Pamminger, Michael↗

Nondestructive measurements of residual 235 U mass of Israeli Research Reactor-1 fuel using the Advanced Experimental Fuel Counter

In 2018, a measurement campaign took place with participants from Los Alamos National Laboratory (LANL), the Nuclear Research Centre-Negev (NRCN) and Soreq Nuclear Research Center (SNRC) at the Israeli Research Reactor-1 (IRR-1) in which 14 of the reactor’s used fuel assemblies (FAs) with varied amount of depletion were measured with the nondestructive assay instrument Advanced Experimental Fuel Counter (AEFC). Designed for safeguards purposes, the AEFC measures both neutrons emitted from the FA (passive neutrons) and fission neutrons induced by an external neutron source (in this experiment, 252 Cf). Signals recorded with the AEFC include total neutron count rates (Singles), time-correlated neutron count rates (Doubles), and total gamma-ray count rates. The 235 U content of the FAs was previously assessed by two independent methods: (1) measurement of the transparency of the FA to low-energy gamma rays from an activated rhenium source (rhenium gamma transmission, or the RGT method) and (2) calculation of the 30-year burnup history of the core using detailed three-dimensional Monte-Carlo core depletion calculations. The results from the FAs that had been measured via the RGT method were used to construct the calibration curves, which translate the AEFC count rates to 235 U mass. Then, the calibration was evaluated using AEFC measurements of six additional FAs that were not measured via the RGT method. From the results, it was determined the Doubles calibration curve was more reliable than that of the Singles and follows a simple second-order polynomial fit for the whole range of residual 235 U mass content, albeit with larger statistical uncertainty. Detailed uncertainties quantification was conducted for both the AEFC Singles and Doubles. This includes the analysis of statistical uncertainties, calibration uncertainty, and random uncertainties due to the sensitivity of the AEFC to several sources of uncertainty, namely the FA position, FA orientation, interrogation source position, and ambient pool temperature. Finally, an overall total uncertainty of 6 g of 235 U is estimated for the Singles and Doubles, which is mainly due to calibration uncertainty (for the Singles) and statistical uncertainty (for the Doubles), and which constitutes 3%–6% of the 235 U total mass in the FAs, depending on their level of depletion.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Creep Behavior of a High-$\gamma^{\prime}$ Ni-Based Superalloy Fabricated via Electron Beam Melting

Additive manufacturing enables the fabrication of complex engineering components previously inaccessible through traditional processes. Nickel-base superalloys with large γ' volume fraction are typically considered non-weldable and therefore exhibit a propensity for cracking during the fusion process. These crack-prone materials, however, are of great importance in gas turbine engines due to their excellent high temperature creep resistance. In this study we investigate the creep behavior of IN738LC produced by the electron beam melting process. Here, we find that with appropriate post-build heat treatment the creep response of material oriented in the build direction exhibits deformation and rupture behavior comparable to that of conventionally cast IN738 & IN738LC. In the transverse direction properties fall below the expected cast behavior, however, we argue this is likely due to differences in grain scale and crystallographic texture. It may be possible to coarsen the grain morphology with appropriate process-parameter optimization in order to reduce the severity of intergranular fracture in the transverse direction. These results illustrate that high temperature properties exhibited by additively manufactured IN738LC are suitable for the hot section of gas turbine engines.

36 MATERIALS SCIENCE↗

Three-dimensional modeling and simulations of single-crystal and bi-crystal titanium for high-strain-rate loading conditions

High purity single crystal titanium (Ti) under shock wave loading is modeled under both one-dimensional and three-dimensional cylindrical conditions. Cylinder sizes of 10 μm and 20 μm radius are both considered in order to assess influence of boundary conditions. A thermodynamically consistent single crystal model for application to shock conditions is presented. The model accounts for the coupled non-linear elastic, dislocation slip, deformation twinning, and structural phase transformation response of the titanium material. Plate impact experiment results using a copper flyer are used to compare against the simulations for crystals oriented in [0001] and [10$\overline1$1] crystallographic directions. The one-dimensional and threedimensional simulations of the two differently oriented single crystals indicate differences between the one-dimensional and three-dimensional representation, especially for the [10$\overline1$1] oriented single crystal. This orientation breaks the relative orientation symmetry between the crystal and cylinder which otherwise exists for the [0001] oriented single crystal. A significant amount of heterogeneity in the field response of the[10$\overline1$1] oriented simulation was demonstrated due to the highly coupled nature of the deformation. A bi-crystal model composed of both the [0001] and [10$\overline1$1] orientations with the boundary between the two along the axis of the cylinder is also considered for a cylinder model size of 10 μm. The results indicate a strong interaction between the two grains that affects the ω phase volume fraction achieved relative to the single crystal calculations.

36 MATERIALS SCIENCE↗