Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Single crystal model”

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 271 records · Page 15

Magnetism of single crystalline breathing pyrochlore spinel AgInCr 4 ⁢S 8

Single crystals of Ag⁢In⁢Cr 4⁢ S 8 are grown by chemical vapor transport and crystallographic ordering of Ag/In that results in a breathing pyrochlore motif of Cr 3+ is verified by x-ray and neutron diffraction. Long-range antiferromagnetic order is observed below a Néel temperature of 𝑇 N ≈ 9.6 K. The magnetic properties are characterized using ac and dc magnetization, specific heat capacity, and single crystal neutron diffraction measurements. The specific heat data are characterized by a small lambda anomaly near 9.5 K and the estimated magnetic entropy reaches ≈ $\frac{1}{3}$ of the expected value by 3⁢𝑇 N , suggesting significant short-range order in the paramagnetic phase. Single crystal neutron diffraction evidences an incommensurate spin structure with propagation vector 𝒌 = (0,0,𝛿) and 𝛿 = 0.343 at 5 K. The minimal model that accounts for the data consists of ferromagnetic layers of Cr atoms, with magnetic moments lying in the plane of the layers and modulating in the perpendicular direction to form a helical structure propagating along 𝒌. Furthermore, this study represents a rare investigation of single crystals within the family of breathing pyrochlore materials.

Frustrated magnetism↗

Observation of the spiral spin liquid in a triangular-lattice material

The spiral spin liquid (SSL) is a highly degenerate state characterized by a continuous contour or surface in reciprocal space spanned by a spiral propagation vector. Although the SSL state has been predicted in a number of various theoretical models, very few materials are so far experimentally identified to host such a state. Via combined single-crystal wide-angle and small-angle neutron scattering, we report observation of the SSL in the quasi-two-dimensional delafossite-like AgCrSe 2 . We show that it is a very close realization of the ideal Heisenberg J 1 –J 2 –J 3 frustrated model on the triangular lattice. By supplementing our experimental results with microscopic spin-dynamics simulations, we demonstrate how such exotic magnetic states are driven by thermal fluctuations and exchange frustration.

36 MATERIALS SCIENCE↗

Phase Field Modeling of Crystallographic Corrosion Pits

The modeling of corrosion to understand and predict corrosion behavior is a topical issue. Here, a 3D phase field model is developed to simulate the pit morphology, primarily focusing on crystallographic pits. A crystallographic function is employed to incorporate different corrosion rates for different crystallographic planes. The model is benchmarked and validated against an analytical solution for a simple case. 3D crystallographic pits in a single crystal are simulated and the effect of substrate orientation on the pit morphology is studied. The crystallographic pit morphology changes significantly with the substrate orientation and these morphologies have a symmetry consistent with the substrate orientation. This first 3D phase field model of crystallographic pits will help in predicting the intricate shapes of pits, thereby, pushing the frontiers of pitting corrosion modeling.

36 MATERIALS SCIENCE↗

A hybrid electronically scanned pressure module for cryogenic environments

Pressure is one of the most important parameters measured when testing models in wind tunnels. For models tested in the cryogenic environment of the National Transonic Facility at NASA Langley Research Center, the technique of utilizing commercially available multichannel pressure modules inside the models is difficult due to the small internal volume of the models and the requirement of keeping the pressure transducer modules within an acceptable temperature range well above the -173 degrees C tunnel temperature. A prototype multichannel pressure transducer module has been designed and fabricated with stable, repeatable sensors and materials optimized for reliable performance in the cryogenic environment. The module has 16 single crystal silicon piezoresistive pressure sensors electrostatically bonded to a metalized Pyrex substrate for sensing the wind tunnel model pressures. An integral temperature sensor mounted on each silicon micromachined pressure sensor senses real-time temperature fluctuations to within 0.1 degrees C to correct for thermally induced non-random sensor drift. The data presented here are from a prototype sensor module tested in the 0.3 M cryogenic tunnel and thermal equilibrium conditions in an environmental chamber which approximates the thermal environment (-173 degrees C to +60 degrees C) of the National Transonic Facility.

Chapman, J. J.↗

GCSS Idealized Cirrus Model Comparison Project

The GCSS Working Group on Cirrus Cloud Systems (WG2) is conducting a systematic comparison and evaluation of cirrus cloud models. This fundamental activity seeks to support the improvement of models used for climate simulation and numerical weather prediction through assessment and improvement of the "process" models underlying parametric treatments of cirrus cloud processes in large-scale models. The WG2 Idealized Cirrus Model Comparison Project is an initial comparison of cirrus cloud simulations by a variety of cloud models for a series of idealized situations with relatively simple initial conditions and forcing. The models (16) represent the state-of-the-art and include 3-dimensional large eddy simulation (LES) models, two-dimensional cloud resolving models (CRMs), and single column model (SCM) versions of GCMs. The model microphysical components are similarly varied, ranging from single-moment bulk (relative humidity) schemes to fully size-resolved (bin) treatments where ice crystal growth is explicitly calculated. Radiative processes are included in the physics package of each model. The baseline simulations include "warm" and "cold" cirrus cases where cloud top initially occurs at about -47C and -66C, respectively. All simulations are for nighttime conditions (no solar radiation) where the cloud is generated in an ice supersaturated layer, about 1 km in depth, with an ice pseudoadiabatic thermal stratification (neutral). Continuing cloud formation is forced via an imposed diabatic cooling representing a 3 cm/s uplift over a 4-hour time span followed by a 2-hour dissipation stage with no cooling. Variations of these baseline cases include no-radiation and stable-thermal-stratification cases. Preliminary results indicated the great importance of ice crystal fallout in determining even the gross cloud characteristics, such as average vertically-integrated ice water path (IWP). Significant inter-model differences were found. Ice water fall speed is directly related to the shape of the particle size distribution and the habits of the ice crystal population, whether assumed or explicitly calculated. In order to isolate the fall speed effect from that of the associated ice crystal population, simulations were also performed where ice water fall speed was set to the same constant value everywhere in each model. Values of 20 and 60 cm/s were assumed. Current results of the project will be described and implications will be drawn. In particular, this exercise is found to strongly focus the definition of issues resulting in observed inter-model differences and to suggest possible strategies for observational validation of the models. The next step in this project is to perform similar comparisons for well observed case studies with sufficient high quality data to adequately define model initiation and forcing specifications and to support quantitative validation of the results.

Starr, David OC.↗

Revisiting the Anisotropic Complex Refractive Indices of Sodium Nitrate for Interpretation of the Reflectance Spectrum of Pressed Pellets

Reflectance spectroscopy is notoriously confounding in that the measured spectral response is highly dependent upon the morphology of the sample. Fortunately, all such perturbations are neatly encoded by the complex refractive index of the sample. Herein, we seek to quantitatively model the measured infrared reflectance spectrum of a specularly flat pressed pellet sample of the birefringent compound, sodium nitrate. Single crystals of sodium nitrate were synthesized via a slow evaporation process and spectroscopically analyzed using polarization-dependent infrared single-angle reflectance spectroscopy. The anisotropic complex refractive index was measured from 7500 to 300 cm-1 (1.33 to 33.33 µm). The deduced anisotropic optical constants were found to be consistent with those previously reported. Once measured and validated, the optical constants were applied to model the pressed pellet reflectance spectrum. It was evident that an average of the anisotropic refractive indices was insufficient to account for the measured pellet reflectance. In order to account for contributions of all possible microcrystalline orientations within the pellet, the Python package PYELLI was used to calculate a basis set of orientation-dependent reflection spectra spanning the distinct ? and ? Euler rotations of the uniaxial crystal. When the population of orientations was allowed to vary freely in a spectral fit analysis, the fit-deduced orientations were tightly clustered along f = 45º, hinting at residual anisotropy in the pressed pellet sample. Conversely, an equally valid spectral fit (with marginally worse fit metric) was obtained when the population was constrained to an isotropic distribution of orientations. Subsequent non-zero cross-polarization reflectance measurements likewise suggested anisotropy in the pellet. However, both grazing incidence wide-angle x-ray scattering and scanning electron microscopy measurements revealed that the microcrystal orientations at the surface of the pressed pellet sample were isotropically distributed (and that the average crystallite size was larger than ?/10). Application of the measured complex refractive indices for modeling the reflectance spectrum of the pressed pellet, and rectification of these seemingly contradictory observations will be discussed.

Wilhelm, Michael J.↗

Machine learning assisted derivation of minimal low-energy models for metallic magnets

Abstract We consider the problem of extracting a low-energy spin Hamiltonian from a triangular Kondo Lattice Model (KLM). The non-analytic dependence of the effective spin-spin interactions on the Kondo exchange excludes the use of perturbation theory beyond the second order. We then introduce a Machine Learning (ML) assisted protocol to extract effective two- and four-spin interactions. The resulting spin model reproduces the phase diagram of the original KLM as a function of magnetic field and single-ion anisotropy and reveals the effective four-spin interactions that stabilize the field-induced skyrmion crystal phase. Moreover, this model enables the computation of static and dynamical properties with a much lower numerical cost relative to the original KLM. A comparison of the dynamical spin structure factor in the fully polarized phase computed with both models reveals a good agreement for the magnon dispersion even though this information was not included in the training data set.

Sharma, Vikram (ORCID:0000000189105519)↗

Second order tensor finite element

The results of a research and software development effort are presented for the finite element modeling of the static and dynamic behavior of anisotropic materials, with emphasis on single crystal alloys. Various versions of two dimensional and three dimensional hybrid finite elements were implemented and compared with displacement-based elements. Both static and dynamic cases are considered. The hybrid elements developed in the project were incorporated into the SPAR finite element code. In an extension of the first phase of the project, optimization of experimental tests for anisotropic materials was addressed. In particular, the problem of calculating material properties from tensile tests and of calculating stresses from strain measurements were considered. For both cases, numerical procedures and software for the optimization of strain gauge and material axes orientation were developed.

Oden, J. Tinsley↗

Large Engine Technology (LET) Short Haul Civil Tiltrotor Contingency Power Materials Knowledge and Lifing Methodologies

This report documents the results of an experimental program conducted on two advanced metallic alloy systems (Rene' 142 directionally solidified alloy (DS) and Rene' N6 single crystal alloy) and the characterization of two distinct internal state variable inelastic constitutive models. The long term objective of the study was to develop a computational life prediction methodology that can integrate the obtained material data. A specialized test matrix for characterizing advanced unified viscoplastic models was specified and conducted. This matrix included strain controlled tensile tests with intermittent relaxtion test with 2 hr hold times, constant stress creep tests, stepped creep tests, mixed creep and plasticity tests, cyclic temperature creep tests and tests in which temperature overloads were present to simulate actual operation conditions for validation of the models. The selected internal state variable models where shown to be capable of representing the material behavior exhibited by the experimental results; however the program ended prior to final validation of the models.

Spring, Samuel D.↗

Genarris 2.0: A Random Structure Generator for Molecular Crystals

Genarris is an open source Python package for generating random molecular crystal structures with physical constraints for seeding crystal structure prediction algorithms and training machine learning models. Here we present a new version of the code, containing several major improvements. A MPI-based parallelization scheme has been implemented, which facilitates the seamless sequential execution of user-defined workflows. A new method for estimating the unit cell volume based on the single molecule structure has been developed using a machine-learned model trained on experimental structures. A new algorithm has been implemented for generating crystal structures with molecules occupying special Wyckoff positions. A new hierarchical structure check procedure has been developed to detect unphysical close contacts efficiently and accurately. New intermolecular distance settings have been implemented for strong hydrogen bonds. To demonstrate these new features, we study two specific cases: benzene and glycine. Genarris finds the experimental structures of the two polymorphs of benzene and the three polymorphs of glycine. Program summary Program Title: Genarris 2.0 Program Files doi: http://dx.doi.org/10.17632/grx6mz4pjn.1 Licensing provisions: BSD-3 Clause Programming language: Python, C External routines/libraries: Spglib, ASE, pymatgen, SciPy, mpi4py, scikit-learn, PyTorch, FHI-aims. Nature of problem: Molecular crystal structure prediction. Solution method: Genarris 2.0 generates molecular crystal structures over the 230 space groups, on general and special Wyckoff positions, using physical constraints. Down-sampling of the generated structures may be performed subsequently, based on molecular crystal packing descriptors and an unsupervised machine learning algorithm. Lastly, ab initio structure relaxation may be performed for the final pool. Depending on the user-defined workflow implemented, Genarris may be used to generate diverse molecular crystal datasets to seed evolutionary algorithms or to train machine learning algorithms or as a standalone crystal structure prediction method. Restrictions: For crystal structure generation, the molecule of interest must be semi-rigid with no bond rotational degrees of freedom. Unusual features: Genarris 2.0 is a highly distributed program, making use of MPI for Python parallelization. The user has the ability to design and implement workflows by executing a user-defined list of procedures. Genarris 2.0 offers new features including a machine learning model for estimating the molecular volume in the solid state from the single molecule structure, structure generation in special Wyckoff positions of space groups, hierarchical structure checks including rigorous treatment of non-orthogonal structures, and clustering and down-selection workflows combining first principles simulations with machine learning. (C) 2020 Elsevier B.V. All rights reserved.

Crystal structure prediction↗

An Efficient Single Frequency Ho:YLF Laser for IPDA Lidar Applications

A highly efficient, versatile, single frequency 2-micron pulsed laser can be used in a pulsed Differential Absorption Lidar (DIAL) / Integrated Path Differential Absorption (IPDA) instrument to make precise, high-resolution measurements to investigate sources and sinks of CO2. For a direct detection IPDA lidar, the desired 2 m Ho:YLF laser should generate 30-40 mJ pulses at the repetition rate of 100 to 200 Hz, with short pulse length (<100 ns) and better than 2% wall plug efficiency. A Tm fiber laser in-band pumped Ho:YLF laser has been developed to meet this technical challenge. This Ho:YLF laser is designed in a four mirror ring resonator with bow tie configuration, which helps to obtain high beam quality. It is end-pumped by a 40 W linearly polarized Tm fiber laser at 1.94μm. The resonator length is 1.10 meters with output coupler reflectivity at 45%. The laser crystal size is 3 x 3 x 60 mm (w, h, l) with a doping concentration of 0.5% Holmium. The laser beam and pump beam are mode-matched in the active medium. Thus, the pump and laser beams have the same confocal parameters. Mode-matching is also helpful for operating the laser in a single transverse mode. The laser beam waist is slightly less than 0.5 mm at the center of the laser crystal. Based on quasi-four level modeling, pump absorption and saturation depend on laser intensity. Laser amplification and saturation also depend on the pump intensity in the crystal. The laser is injection seeded to obtain the single frequency required by an IPDA lidar measurement. The seed beam is entered into the resonator through an output coupler. The laser is mounted on a water cooled optical bench for stable and reliable operation. The size of the optical bench is 22.16 x 9.20 x 1.25 inches. It is stiffened so that the laser can be operated in any orientation of the optical bench. This packaged Ho:YLF laser is designed for either mobile trailer or airborne platform operation. The engineering prototype Ho:YLF laser has been fully characterized to demonstrate laser performance. Figure 1 shows the laser output power as a function of pump power at different pulse repetition rates from 100 Hz to 333 Hz. The threshold of the laser is less than 14 W. The slope efficiencies are 28%, 40%, 41% and 43% for pulse repetition rates of 100, 200, 250 and 333 Hz, respectively. Maximum power increases with the pulse repetition rate. Output power of 4.2 W, 6 W, 6.7 W, and 7.7 W is achieved for pulse repetition rates of 100, 200, 250 and 333Hz, respectively. This represents the optical conversion efficiency of 16.7%, 22.4%, 23.7%, and 26.5% at these various pulse repletion rates. It is the most efficient and compact Ho:YLF laser demonstrated in the high pulse energy (>20mJ) and moderate pulse repetition rate (100's Hz) operation range. As shown in Figure 1, the maximum pulse energy at 100 Hz is 42 mJ. This is limited due to optical damage. The laser stability is characterized and found to be very stable. A relative pulse energy standard deviation of 2% was measured. The beam quality of the Ho:YLF was measured by a Spiricon infrared laser beam camera. Figure 2 shows the beam profile image of the laser. Both the X-profile of the beam (horizontal direction) and the Y-profile of the beam (vertical direction) are well fitted by a Gaussian profile. The qualitative beam quality measurement shows excellent beam quality in both axis. The M-square value for the laser beam is measured at 1.06 and 1.09 for the x and y axis respectively.-

Yu, J.↗

Dose-dependent structural and electron-density features in the lytic polysaccharide monooxygenase NcAA9D

Structural studies of copper-containing lytic polysaccharide monooxygenases (LPMOs) by X-ray crystallography are often complicated by radiation damage. In this study, we analyze a series of 36 X-ray crystal structures of NcAA9D, a Neurospora crassa AA9-family LPMO, determined from data collected at cryogenic temperature from a single crystal to investigate the progressive effects of radiation damage at the active site of this enzyme. We report new insights into the dose-dependence of active-site geometry in LPMOs and utilize the unique pre-bound dioxygen site of NcAA9D to analyze the impact of X-ray dose on the electron density of this species. It is well established that photoreduction of the LPMO active-site copper(II) leads to expulsion of its water ligands. We further characterize this displacement and the corresponding electron-density smearing, a phenomenon that can lead to the erroneous modeling of copper-bound dioxygen species. These findings suggest that radiation-dose series collected from a single crystal provide invaluable data to support unambiguous assignment of radiation-sensitive intermediates at the active site of LPMOs and other radiation-sensitive redox enzymes.

Miller, Samuel [ORNL] (ORCID:0009000459491817)↗

Finite Element Modeling of a Cylindrical Contact Using Hertzian Assumptions

The turbine blades in the high-pressure fuel turbopump/alternate turbopump (HPFTP/AT) are subjected to hot gases rapidly flowing around them. This flow excites vibrations in the blades. Naturally, one has to worry about resonance, so a damping device was added to dissipate some energy from the system. The foundation is now laid for a very complex problem. The damper is in contact with the blade, so now there are contact stresses (both normal and tangential) to contend with. Since these stresses can be very high, it is not all that difficult to yield the material. Friction is another non-linearity and the blade is made out of a Nickel-based single-crystal superalloy that is orthotropic. A few approaches exist to solve such a problem and computer models, using contact elements, have been built with friction, plasticity, etc. These models are quite cumbersome and require many hours to solve just one load case and material orientation. A simpler approach is required. Ideally, the model should be simplified so the analysis can be conducted faster. When working with contact problems determining the contact patch and the stresses in the material are the main concerns. Closed-form solutions for non-conforming bodies, developed by Hertz, made out of isotropic materials are readily available. More involved solutions for 3-D cases using different materials are also available. The question is this: can Hertzian1 solutions be applied, or superimposed, to more complicated problems-like those involving anisotropic materials? That is the point of the investigation here. If these results agree with the more complicated computer models, then the analytical solutions can be used in lieu of the numerical solutions that take a very long time to process. As time goes on, the analytical solution will eventually have to include things like friction and plasticity. The models in this report use no contact elements and are essentially an applied load problem using Hertzian assumptions to determine the contact patch dimensions.

Knudsen, Erik↗

Bayesian calibration of a physics-based crystal plasticity and damage model

In this work, we present a model parameter calibration procedure for a physics-based crystal plasticity model. The calibration process utilizes a powerful statistics-based Bayesian calibration method. Calibration of the crystal plasticity parameters makes use of experimentally-measured data, i.e. compressive stress–strain response, from $\langle$100$\rangle$ and $\langle$123$\rangle$ single crystal copper dynamically loaded via Kolsky bar tests. The calibration of damage parameters is achieved using experimentally-measured free-surface velocity history data from plate impact test on $\langle$100$\rangle$ and $\langle$110$\rangle$ single crystal copper, which generates shock compression followed by dynamic tensile failure. A validation assessment is then carried out by comparing the calibrated model predictions and experimental measurements of the dynamic tensile damage generated in an impacted bicrystal copper plate. Lastly, a model-informed rationale for the experimentally-observed dependence of the spatial distribution of ductile damage (porosity) on crystallography is provided.

42 ENGINEERING↗

Vortex phases and glassy dynamics in the highly anisotropic superconductor HgBa 2 CuO 4+δ

We present an extensive study of vortex dynamics in a high-quality single crystal of HgBa 2 CuO 4+ δ , a highly anisotropic superconductor that is a model system for studying the effects of anisotropy. From magnetization M measurements over a wide range of temperatures T and fields H , we construct a detailed vortex phase diagram. We find that the temperature-dependent vortex penetration field H p ( T ), second magnetization peak H smp ( T ), and irreversibility field H irr ( T ) all decay exponentially at low temperatures and exhibit an abrupt change in behavior at high temperatures T / T c >~0.5. By measuring the rates of thermally activated vortex motion (creep) S ( T , H ) = | d ln M ( T , H )/ d ln t |, we reveal glassy behavior involving collective creep of bundles of 2D pancake vortices as well as temperature- and time-tuned crossovers from elastic (collective) dynamics to plastic flow. Based on the creep results, we show that the second magnetization peak coincides with the elastic-to-plastic crossover at low T , yet the mechanism changes at higher temperatures.

36 MATERIALS SCIENCE↗

Growth of zinc selenide single crystals by physical vapor transport in microgravity

The goals of this research were the optimization of growth parameters for large (20 mm diameter and length) zinc selenide single crystals with low structural defect density, and the development of a 3-D numerical model for the transport rates to be expected in physical vapor transport under a given set of thermal and geometrical boundary conditions, in order to provide guidance for an advantageous conduct of the growth experiments. In the crystal growth studies, it was decided to exclusively apply the Effusive Ampoule PVT technique (EAPVT) to the growth of ZnSe. In this technique, the accumulation of transport-limiting gaseous components at the growing crystal is suppressed by continuous effusion to vacuum of part of the vapor contents. This is achieved through calibrated leaks in one of the ground joints of the ampoule. Regarding the PVT transport rates, a 3-D spectral code was modified. After introduction of the proper boundary conditions and subroutines for the composition-dependent transport properties, the code reproduced the experimentally determined transport rates for the two cases with strongest convective flux contributions to within the experimental and numerical error.

Rosenberger, Franz↗

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗

Beyond single tetrahedron physics of the breathing pyrochlore compound Ba 3 ⁢Yb 2 ⁢Zn 5 ⁢O 11

Recently, a new class of quantum magnets, the so-called breathing pyrochlore spin systems, have attracted much attention due to their potential to host exotic emergent phenomena. Here, we present magnetometry, heat capacity, thermal conductivity, muon-spin relaxation, and polarized inelastic neutron scattering measurements performed on high-quality single crystal samples of the breathing pyrochlore compound Ba 3 ⁢Yb 2 ⁢Zn 5 ⁢O 11 . We interpret these results using a simplified toy model and provide insight into the low-energy physics of this system beyond the single tetrahedron physics proposed previously.

36 MATERIALS SCIENCE↗