Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “derived lattice”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Cluster Expansion Analysis of Atomic Order in Li-Ion Battery Cathode Material LiCo y Ni 1-y O 2

A modified cluster-expansion treatment was developed recently to analyze atomic order in LiCo y Ni 1-y O 2 , a model cathode material for Li-ion batteries. In this treatment, referred to as a “spin-atom” cluster expansion, the occupant of a lattice site is identified by its spin state as well as its atomic species. Further, Effective Cluster Interaction (ECI) coefficients are derived from a large training data set (i.e., the set of atomic arrangements for which DFT calculations are performed) which is filtered by an anomaly detection algorithm to eliminate poorly converged DFT calculations. The cluster expansion incorporates Li-Ni (LN) exchange as well as intralayer Co-Ni (CN) exchange. Monte Carlo simulations based on the cluster expansion were applied to the Ni-rich part of the phase diagram. The simulations predict a miscibility gap between y = 0.05 and y = 0.65.

25 ENERGY STORAGE↗
NASA NTRS1980

Generalized Grueneisen tensor from solid nonlinearity parameters

Anharmonic effects in solids are often described in terms of generalized Grueneisen parameters which measure the strain dependence of the lattice vibrational frequencies. The relationship between these parameters and the solid nonlinearity parameters measured directly in ultrasonic harmonic generation experiments is derived using an approach valid for normal-mode elastic wave propagation in any crystalline direction. The resulting generalized Grueneisen parameters are purely isentropic in contrast to the Brugger-Grueneisen parameters which are of a mixed thermodynamic state. Experimental data comparing the isentropic generalized Grueneisen parameters and the Brugger-Grueneisen parameters are presented.

Cantrell, J. H., Jr.↗
NASA NTRS1999

Modeling of Macroscopic/Microscopic Transport and Growth Phenomena in Zeolite Crystal Solutions Under Microgravity Conditions

Crystals grown from liquid solutions have important industrial applications. Zeolites, for instance, a class of crystalline aluminosilicate materials, form the backbone of the chemical process industry worldwide, as they are used as adsorbents and catalysts. Many of the phenomena associated with crystal growth processes are not well understood due to complex microscopic and macroscopic interactions. Microgravity could help elucidate these phenomena and allow the control of defect locations, concentration, as well as size of crystals. Microgravity in an orbiting spacecraft could help isolate the possible effects of natural convection (which affects defect formation) and minimize sedimentation. In addition, crystals will stay essentially suspended in the nutrient pool under a diffusion-limited growth condition. This is expected to promote larger crystals by allowing a longer residence time in a high-concentration nutrient field. Among other factors, the crystal size distribution depends on the nucleation rate and crystallization. These two are also related to the "gel" polymerization/depolymerization rate. Macroscopic bulk mass and flow transport and especially gravity, force the crystals down to the bottom of the reactor, thus forming a sedimentation layer. In this layer, the growth rate of the crystals slows down as crystals compete for a limited amount of nutrients. The macroscopic transport phenomena under certain conditions can, however, enhance the nutrient supply and therefore, accelerate crystal growth. Several zeolite experiments have been performed in space with mixed results. The results from our laboratory have indicated an enhancement in size of 30 to 70 percent compared to the best ground based controls, and a reduction of lattice defects in many of the space grown crystals. Such experiments are difficult to interpret, and cannot be easily used to derive empirical or other laws since many physical parameters are simultaneously involved in the process. At the same time, however, there is increased urgency to develop such an understanding in order to more accurately quantify the process. In order to better understand the results obtained from our prior space experiments, and design future experiments, a detailed fluid dynamic model simulating the crystal growth mechanism is required. This will not only add to the fundamental knowledge on the crystallization of zeolites, but also be useful in predicting the limits of size and growth of these important industrial materials. Our objective is to develop macro/microscopic theoretical and computational models to study the effect of transport phenomena in the growth of crystals grown in solutions. Our effort has concentrated so far in the development of separate macroscopic and microscopic models. The major highlights of our accomplishments are described.

Gatsonis, Nikos A.↗
NASA NTRS2003

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗
NASA NTRS1922

Notes on the design of latticed columns subject to lateral loads

No rational column formula has yet been developed which gives results which are sufficiently precise for the design of airplane members, and consequently it is necessary to fall back upon experimental testing. In order to derive the maximum benefit from experiments, however, it is necessary that the experiments be guided by theory. It is the object here to modify existing formulae that may be obtained with a minimum number of tests. A discussion is given that is limited to the case of a simple column supported at both ends and subjected to uniformly distributed loads perpendicular to its axis and to end loads either axially or eccentrically applied. Discussed here are forces in the bracing of latticed columns, the application of theory to practical columns, loads in lattice members, and the strength of individual lattices.

Mccarthy, Charles J↗
DOE OSTI2020

Coulombic friction in metamaterials to dissipate mechanical energy

Product designs from a wide range of industries such as aerospace, automotive, biomedical, and others can benefit from new metamaterials for mechanical energy dissipation. In this study, we explore a novel new class of metamaterials with unit cells that absorb energy via sliding Coulombic friction. Remarkably, even materials such as metals and ceramics, which typically have no intrinsic reversible energy dissipation, can be architected to provide dissipation akin to elastomers. The concept is demonstrated at different scales (centimeter to micrometer), with different materials (metal and polymer), and in different operating environments (high and low temperatures), all showing substantial dissipative improvements over conventional non-contacting lattice unit cells. Further, as with other ‘programmable’ metamaterials, the degree of Coulombic absorption can be tailored for a given application. Additionally, an analytic expression is derived to allow rapid first-order optimization. This new class of Coulombic friction energy absorbers can apply broadly to many industrial sectors such as transportation (e.g. monolithic shock absorbers), biomedical (e.g. prosthetics), athletic equipment (e.g. skis, bicycles, etc.), defense (e.g. vibration tolerant structures), and energy (e.g. survivable electrical grid components).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Expanded Analogs of Three‐Dimensional Lead‐Halide Hybrid Perovskites

Abstract Replacing the Pb−X octahedral building unit of A I PbX 3 perovskites (X=halide) with a pair of edge‐sharing Pb−X octahedra affords the expanded perovskite analogs: A II Pb 2 X 6 . We report seven members of this new family of materials. In 3D hybrid perovskites, orbitals from the organic molecules do not participate in the band edges. In contrast, the more spacious inorganic sublattice of the expanded analogs accommodates larger pyrazinium‐based cations with low‐lying π* orbitals that form the conduction band, substantially decreasing the band gap of the expanded lattice. The molecular nature of the conduction band allows us to electronically dope the materials by reducing the organic molecules. By synthesizing derivatives with A II =pyridinium and ammonium, we can isolate the contributions of the pyrazinium‐based orbitals in the band gap transition of A II Pb 2 X 6 . The organic‐molecule‐based conduction band and the inorganic‐ion‐based valence band provide an unusual electronic platform with localized states for electrons and more disperse bands for holes upon optical or thermal excitation.

Umeyama, Daiki↗
DOE OSTI2021

Spontaneous imbibition in a square tube with corner films: theoretical model and numerical simulation

Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.

42 ENGINEERING↗
NASA NTRS1981

Effects of wingtip modifications on handling qualities of agricultural aircraft

The effect of wingtip modifications on the stability and control characteristics of an agricultural airplane has been studied by means of a nonplanar quasi-vortex-lattice method. The method is used to compute the changes in steady state and perturbed state lateral-directional stability and control derivatives produced by wingtip mounted winglets, vortex diffuser vanes, and tip extensions. The study shows that the combination of the excessive positive dihedral effect produced by the winglets and adverse yaw due to aileron deflection can have a detrimental effect on the roll control characteristics of the airplane. Introduction of an aileron-rudder-interconnect, and reduction of the effective dihedral by canting-in of the winglets, or addition of a lower winglet can eliminate the roll control problems.

Van Dam, C. P.↗
NASA NTRS1981

Clay-mediated reactions of HCN oligomers - The effect of the oxidation state of the clay

Montmorillonite clays which contain Fe(III) inhibit the oligomerization of aqueous solutions of HCN. The inhibitory effect is due to the rapid oxidation of diaminomaleonitrile, a key intermediate in HCN oligomerization, by the Fe(III) incorporated into the aluminosilicate lattice of the clay. The Fe(III) oxidizes diaminomaleonitrile to diiminosuccinonitrile, a compound which is rapidly hydrolyzed to HCN and oxalic acid derivatives. Diaminomaleonitrile is not oxidized when Fe(III) in the montmorillonite is reduced with hydrazine. The oxidation state of the clay is an important variable in experiments designed to simulate clay catalysis on the primitive earth.

Ferris, J. P.↗
NASA NTRS2011

Titan's Aerosol and Stratospheric Ice Opacities Between 18 and 500 Micrometers: Vertical and Spectral Characteristics from Cassini CIRS

Vertical distributions and spectral characteristics of Titan's photochemical aerosol and stratospheric ices are determined between 20 and 560 per centimeter (500-18 micrometers) from the Cassini Composite Infrared Spectrometer (CIRS). Results are obtained for latitudes of 15 N, 15 S, and 58 S, where accurate temperature profiles can be independently determined. In addition, estimates of aerosol and ice abundances at 62 N relative to those at 15 S are derived. Aerosol abundances are comparable at the two latitudes, but stratospheric ices are approximately 3 times more abundant at 62 N than at 15 S. Generally, nitrile ice clouds (probably HCN and HC3N), as inferred from a composite emission feature at approximately 160 per centimeter, appear to be located over a narrow altitude range in the stratosphere centered at approximately 90 km. Although most abundant at high northern latitudes, these nitrile ice clouds extend down through low latitudes and into mid southern latitudes, at least as far as 58 S. There is some evidence of a second ice cloud layer at approximately 60 km altitude at 58 S associated with an emission feature at approximately 80 per centimeter. We speculate that the identify of this cloud may be due to C2H6 ice, which in the vapor phase is the most abundant hydrocarbon (next to CH4) in the stratosphere of Titan. Unlike the highly restricted range of altitudes (50-100 km) associated with organic condensate clouds, Titan's photochemical aerosol appears to be well-mixed from the surface to the top of the stratosphere near an altitude of 300 km, and the spectral shape does not appear to change between 15 N and 58 S latitude. The ratio of aerosol-to-gas scale heights range from 1.3-2.4 at about 160 km to 1.1-1.4 at 300 km, although there is considerable variability with latitude, The aerosol exhibits a very broad emission feature peaking at approximately 140 per centimeter. Due to its extreme breadth and low wavenumber, we speculate that this feature may be caused by low-energy vibrations of two-dimensional lattice structures of large molecules. Examples of such molecules include polycyclic aromatic hydrocarbons (PAHs) and nitrogenated aromatics. Finally, volume extinction coefficients N chi EPSILON derived from 15 S CIRS data at a wavelength of lambda = 62.5 micrometers are compared with those derived from the 10 S Huygens Descent Imager/Spectral Radiometer (DISR) data at 1.583 micrometers. This comparison yields volume extinction coefficient ratios N chi EPSILON (1.583 micrometers)/N chi EPSILON (62.5 micrometers) of roughly 70 and 20, respectively, for Titan's aerosol and stratospheric ices, The inferred particle cross-section ratios chi EPSILON(1.583 micrometers)/chi EPSILON (62.5 micrometers) appear to be consistent with sub-micron size aerosol particles, and effective radii of only a few microns for stratospheric ice cloud particles.

Anderson, Carrie M.↗
DOE OSTI2024

Present and future of $\mathcal{C}$osmo $\mathcal{L}$attice

We discuss the present state and planned updates of Cosmo Lattice, a cutting-edge code for lattice simulations of non-linear dynamics of scalar-gauge field theories in an expanding background. We first review the current capabilities of the code, including the simulation of interacting singlet scalars and of Abelian and non-Abelian scalar-gauge theories. We also comment on new features recently implemented, such as the simulation of gravitational waves from scalar and gauge fields. Secondly, we discuss new extensions of Cosmo Lattice that we plan to release publicly. We comment on new physics modules, which include axion-gauge interactions $\phi F\tilde{F}$, non-minimal gravitational couplings $\phi ^2R$, creation and evolution of cosmic-defect networks, and magnetohydrodynamics. Finally, we also discuss new technical features, including evolvers for non-canonical interactions, arbitrary initial conditions, simulations in 2+1 dimensions, and higher-accuracy spatial derivatives.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Site-decorated model for unconventional frustrated magnets: Ultranarrow phase crossover and two-dimensional spin reversal transition

Here, the site-decorated Ising model is introduced to advance the understanding and experimental realization of the recently discovered one-dimensional (1D) finite-temperature ultranarrow phase crossover in an external magnetic field, while mitigating the geometric complexities of traditional bond-decorated models. The unconventional frustration and physics are clarified by exactly mapping the 1D site-decorated Ising model in a magnetic field onto a zero-field bond-decorated 𝐽 1 −𝐽 2 Ising model with conventional geometrical frustration. Furthermore, although higher-dimensional Ising models in an external field remain unsolved exactly, an exact solution for a spin-reversal transition—driven by an exotic, hidden half-ice, half-fire state induced by site decoration—is derived. This transition, triggered by a slight variation in temperature or magnetic field—without changing its direction—even in the weak-field limit, offers a promising route toward energy-efficient applications such as data storage and processing. The results suggest that site decoration offers an avenue for materials and device design, particularly in systems such as mixed 𝑑−𝑓 compounds, optical lattices, and neural networks, calling for further studies with site-decorated Heisenberg models. In addition, the site-decorated model offers a rigorous test ground for artificial intelligence (AI) in science, as the analytic derivation of the present results was not only validated but also improved by a general-purpose large language model, inspiring the use of AI as scientific discoverer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

ESR Dipole Power Supply Current Ripple and Noise Specifications

This note presents key findings for the ESR main magnet dipole power supplies (PS), where we find the current ripple specification to be close to or beyond the state-of-the-art. These specifications originate from beam-beam considerations, with the requirement to limit the ripple-induced hadron emittance growth to below 10%/hour. Beam dynamics that drive this PS ripple specification arise from the beam motions at the Interaction Point (IP). The frequency of the motions can be separated into "low", compared to the betatron frequency, and "high", i.e. around the betatron frequency and harmonics. In terms of the driving frequency, "low" implies ƒ<<ƒ 0 v x,y and "fast" means ƒ≈{ƒ 0 v x,y , ƒ 0 (1- v x,y , etc.}, where ƒ 0 =1/ T 0 =78.2 kHz is the revolution frequency, and v x,y are the fractional parts of the betatron tunes. Frequencies higher than ƒ 0 /2 will be folded back due to the particles sampling the field once per turn. To provide flexibility for future lattice adjustments and working point variations, we do not consider the tunes as fixed. Instead, we assume a certain margin and allow them to potentially fall within the range of 0.1< v x,y <0.5. In other words, the high-frequency region spans approximately from 8kHz to 40 kHz. Consequently,, we definte the dipole PS ripple in two distinct frequency ranges: the low-frequency range of [1-8000] Hz and the high-frequency range of [8-40] kHz. For the physics effects we analyzed in this note, there is no distinction between the ripple (which can be approximately reproduced in the frequency domain) and the random noise if both have some power within the frequency bandwidth of interest. Therefore, while we will use the term "ripple" for short, it should always be understood that we are referring to "ripple plus noise." Except for the lower end of the low-frequency range, the impact of the rippling PS current on the beam will be considerably reduced due to the eddy currents induced in the walls of the vacuum chamber. We will account for this effect in the PS ripple specifications to follow. The remaining sections of this note are structure as follows: Section II outlines the beam-beam physics requirements for the positional stability of the beam at the IP. Section III describes the anticipated shielding effect of the eddy currents induced in the vacuum chamber. Section IV derives the ripple requirement for the low-frequency range by propagating the closed orbit ripple resulting from the rippling dipoles to the IP (relevant lattice simulation results are summarized in the Appendix). In Section V, we present the analytical criterion for the ripple in the high-frequency range by considering resonant oscillations of the electron beam around a stable closed orbit near the betatron frequency. Finally, Section VI provides a summary of our findings and discusses related work.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

A combined ensemble-volume average homogenization method for lattice structures with defects under dynamic and static loading

In the study of lattices structures, both experiments and numerical simulations are often conducted with small samples. Using combined ensemble and volume averaging, this work introduces a method to extract a macroscopic constitutive response of a lattice material from numerical simulations performed in periodic domains. The domain size needed to obtain statistically accurate results is investigated. Similar to molecular dynamics, the concept of the virial stress is introduced after homogenized equations are derived using the ensemble averaging method. Under static conditions, the virial stress is shown to agree with the volume averaged solid stress. Using the homogenization method, constitutive relations for this stress can be obtained from systems with uniform strains. Application of such obtained constitutive relations to more general cases results in an error proportional to the square of the ratio between the lattice length scale and the macroscopic length scale. Taking advantage of this property, numerical simulations are performed in systems with a uniform gradient of the average velocity. The volume average method is then used to accelerate convergence when studying lattices with defects. To avoid the artificial numerical time scale from the size of a representative volume element divided by the wave speed, a numerical scheme is developed to enforce a spatially uniform velocity gradient within the computational domain while allowing fluctuations of the velocity or displacement to develop naturally. To account for probability distribution of lattice defects, the stress is calculated as the ensemble-volume averaged value. For dynamic systems, energy dissipation properties are also studied.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Phonons and phase symmetries in bulk CrCl 3 from scattering measurements and theory

Phonon-derived behaviors are important indicators of novel phenomena in transition metal trihalides, including spin liquid behavior, two-dimensional magnetism, and spin-lattice coupling. However, phonons and their dependence on spin structure and excitations have not been adequately explored. In this work, we probe and critically examine the vibrational properties of the prototype ferromagnetic honeycomb lattice material CrCl 3 using inelastic neutron scattering and density functional theory. We demonstrate that magnetic and van der Waals interactions are essential to describing the structure and phonons in CrCl 3 ; however, the specific spin configuration is unimportant. Further, this provides context for understanding thermal transport measurements as governed by dynamical spin-lattice couplings. More importantly, we introduce an efficient dynamic method that exploits translational symmetries in large conventional unit cells that generates insights into phonon dispersions, interactions, and measured spectra in terms of quantum phase interference conditions. This work opens new avenues for understanding phonons in layered magnets and more generally in conventional cell geometries of a variety of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Qubit Lattice Algorithms based on the Schrodinger-Dirac representation of Maxwell Equations and their Extensions

It is well known that Maxwell equations can be expressed in a unitary Schrodinger-Dirac representation for homogeneous media. However, difficulties arise when considering inhomoge- neous media. A Dyson map points to a unitary field qubit basis, but the standard qubit lattice algorithm of interleaved unitary collision-stream operators must be augmented by some sparse non-unitary potential operators that recover the derivatives on the refractive indices. The effect of the steepness of these derivatives on two dimensional scattering is examined with simulations showing quite complex wavefronts emitted due to transmissions/reflections within the dielectric objects. Maxwell equations are extended to handle dissipation using Kraus operators. Then, our theoretical algorithms are extended to these open quantum systems. A quantum circuit diagram is presented as well as estimates on the required number of quantum gates for implementation on a quantum computer.

Vahala, George↗