Magnetotransport and magnetic textures in Ho/FeCoGd/ β -W multilayers
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.
We present elucidation of homoepitaxial growth mechanisms on vicinal non-polar surfaces of GaN that is highly important for gaining an understanding of and control thin film surface morphology and properties. Using first-principles calculations, we study the step-flow growth in m-plane GaN based on atomic row nucleation and kink propagation kinetics. Ga–N dimer adsorption onto the m-plane is energetically more favorable than that of Ga and N isolated adatoms. Therefore, we have treated the dimers as the dominant growth species attached to the step edges. By calculating the free energies of sequentially attached Ga–N dimers, we have elucidated that the a-step edge kink growth proceeds by parallel attachment rather than by across the step edge approach. We found a series of favorable configurations of kink propagation and calculated the free energy and nucleation barriers for kink evolution on five types of step edges (a, +c, -c, +a + c, and -a - c). By changing the chemical potential μGa and the excess chemical potential Δμ, the growth velocities at the five types of edges are controlled by the corresponding kink pair nucleation barrier E* in their free energy profiles. To explore the kink-flow growth instability observed at different Ga/N flux ratios, calculations of kink pairs on the incompact -c and +c-step edges are further performed to study their formation energies. Variations of these step edge morphologies with a tuned chemical environment are consistent with previous experimental observations, including stable diagonal ±a ± c-direction steps. In conclusion, our work provides a first-principles approach to explore step growth and surface morphology of the vicinal m-plane GaN, which is applicable to analyze and control the step-flow growth of other binary thin films.
Alpha channeling uses waves to extract hot ash from a fusion plasma, transferring energy from the ash to the wave. It has been proposed that this process could create a radial electric field, efficiently driving E × B rotation. However, existing theories ignore the nonresonant particles, which play a critical role in enforcing momentum conservation in quasilinear theory. Because cross field charge transport and momentum conservation are fundamentally linked, this non-consistency throws the rotation drive into question. This paper has two main goals. First, we provide a pedantic and cohesive introduction to the recently developed simple, general, self-consistent quasilinear theory for electrostatic waves that explains the torques which allow for current drive parallel to the magnetic field, and charge extraction across it; a theory that has largely resolved the question of rotation drive by alpha channeling. We show how the theory reveals a fundamental difference between the reaction of nonresonant particles to plane waves that grow in time vs steady-state waves that have a nonuniform spatial structure, allowing rotation drive in the latter case while precluding it in the former, and we review the local and global conservation laws that lead to this result. Second, we provide two new results in support of the theory. First, we provide a novel two-particle Hamiltonian model that rigorously establishes the relationship between charge transport and momentum conservation. Second, we compare the new quasilinear theory to the oscillation-center theories of ponderomotive forces, showing how the latter often obscure the time-dependent nonresonant recoil, but ultimately lead to similar results.
The spin-flop phase of a honeycomb antiferromagnet of Mn 0.84 Mg 0.16 TiO 3 has been studied under magnetic field along normal direction to the honeycomb plane by neutron scattering. With increasing a magnetic field from 0 to 5 T, the inelastic scattering intensity increased by 1.5 times at (1 0 0), where the magnetic structure was found to be collinear-antiferromagnet with the spin parallel to [1 1 0] direction. The spin direction is consistent with a 1.5 times increase of the inelastic intensity. Furthermore, the negligible energy difference between K and K' points suggests small values for both Dzyaloshinskii–Moriya interaction and bond-dependent symmetric anisotropic exchange interaction in this material.
We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching ~100,000 electrons) using Kohn-Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn-Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn-Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes---both in CPU-times and wall-times, and on both non-periodic and periodic systems---at systems sizes beyond a few thousand electrons, with over 5-10 fold speedups in systems with more than 10,000 electrons. In conclusion, the benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated on up to 192,000 MPI tasks.
Magnetic reconnection in a quasi-parallel bow shock is investigated with two-dimensional local particle-in-cell simulations. In the shock transition and downstream regions, large amplitude magnetic fluctuations exist, and abundant current sheets form. In some current sheets, reconnection occurs, and ion-scale and electron-scale magnetic islands are generated. In electron-scale island regions, only electron outflow jets are observed, producing a quadrupolar out-of-plane magnetic field pattern, while in ion-scale islands, both ions and electrons are involved and energized in reconnection. Finally, normalized reconnection rates are obtained to be between around 0.1 to 0.2, and particle acceleration signatures are seen in distribution functions.
Long schlieren imaging systems, where the parallel light test section is longer than the focal length of the focusing schlieren optics, have a limited region in the test section in which any occluding object, such as a wind tunnel model, is in focus in the schlieren image. Corrector lenses are introduced here to alter the location in the test section where image focus is achieved. Corrector lenses with focal lengths varying from −1000 to −100mm were studied. Lens-type inline and z -type mirror schlieren systems were experimentally tested with multiple collecting optic focal lengths to characterize the changes in focal position. An automated image processing method was used to determine the plane of best focus from image sequences. The introduction of the corrector lens was observed to move the location of best focus within the schlieren system test section while also causing a decrease in the focal sharpness of the images and altering the magnification. The thin lens equation was found to provide a good estimate of the focal location change in the schlieren imaging systems with the addition of the corrector lens.
Diffusive transport processes in magnetized plasmas are highly anisotropic, with fast parallel transport along the magnetic field lines sometimes faster than perpendicular transport by orders of magnitude. This constitutes a major challenge for describing non-grid-aligned magnetic structures in Eulerian (grid-based) simulations. Here, the present paper describes and validates a new method for parallel diffusion in magnetized plasmas based on the anti-symmetry representation [Halpern and Waltz, Phys. Plasmas 25, 060703 (2018)]. In the anti-symmetry formalism, diffusion manifests as a flow operator involving the logarithmic derivative of the transported quantity. Qualitative plane wave analysis shows that the new operator naturally yields better discrete spectral resolution compared to its conventional counterpart. Numerical simulations comparing the new method against existing finite difference methods are carried out, showing significant improvement. In particular, we find that combining anti-symmetry with finite differences in diagonally staggered grids essentially eliminates the so-called “artificial numerical diffusion” that affects conventional finite difference and finite volume methods.
The templated growth of the basic porphyrin unit, free-base porphine (2H-P), is characterized by means of X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine-structure (NEXAFS) spectroscopy measurements and density functional theory (DFT). The DFT simulations allow the deconvolution of the complex XPS and NEXAFS signatures into contributions originating from five inequivalent carbon atoms, which can be grouped into C–N and C–C bonded species. Polarization-dependent NEXAFS measurements reveal an intriguing organizational behavior: On both Cu(111) and Ag(111), for coverages up to one monolayer, the molecules adsorb undeformed and parallel to the respective metal surface. Upon increasing the coverage, however, the orientation of the molecules in the thin films depends on the growth conditions. Multilayers deposited at low temperatures exhibit a similar average tilting angle (30° relative to the surface plane) on both substrates. Conversely, for multilayers grown at room temperature a markedly different scenario exists. On Cu(111) the film thickness is self-limited to a coverage of approximately two layers, while on Ag(111) multilayers can be grown easily and, in contrast to the bulk 2H-P crystal, the molecules are oriented perpendicular to the surface. This difference in molecular orientation results in a modified line-shape of the C 1s XPS signatures, which depends on the incident photon energy and is explained by comparison with depth-resolved DFT calculations. Simulations of ionization energies for differently stacked molecules show no indication for a packing-induced modification of the multilayer XP spectra, thus indicating that the comparison of single molecule calculations to multilayer data is justified.
Electrochemical deposition technique has been used to deposit cobalt nanowires into the nano sized channels of Anodized Aluminium Oxide (AAO) templates. CoCl{sub 2}М‡6H2O salt solution was used, which was buffered with H{sub 3}BO{sub 3} and acidified by dilute H{sub 2}SO{sub 4} to increase the plating life and control pH of the solution. Thin film of copper around 150 nm thick on one side of AAO template coated by e-beam evaporation system served as cathode to create electrical contact. FESEM analysis shows that the as-deposited nanowires are highly aligned, parallel to one another and have high aspect ratio with a reasonably high pore-filing factor. The TEM results show that electrodeposited cobalt nanowires are crystalline in nature. The Hysteresis loop shows the magnetization properties for in and out of plane configuration. The in plane saturation magnetization (Ms) is lower than out of plane configuration because of the easy axis of magnetization is perpendicular to nanowire axis. These magnetic nanowires could be utilized for applications such as spintronic devices, high density magnetic storage, and magnetic sensor applications.
Using the recently developed Bond Order Potential for Cr, we investigate the mobility of isolated 1/2[111] screw dislocations under stress at 0 K. These studies provide the dependence of the critical resolved shear stress (i.e., the Peierls stress) on the orientation of the maximum resolved shear stress plane and on the shear stress perpendicular to the slip direction. Similar studies are made for loading in tension and compression to show that the onset of yielding in Cr single crystals can be described using a combination of the shear stresses parallel and perpendicular to the slip direction. Based on these results, we develop a single crystal yield criterion for Cr that includes all non-Schmid contributions to the Peierls stress. This criterion is used to predict the shape of the yield surface and compare it with the regular Tresca's hexagon arising from the Schmid law. For all directions of uniaxial load, the criterion is used to quantify the character and magnitude of the tension-compression asymmetry.
SUMMARY The seismic anisotropy of the Earth's solid inner core has been the topic of much research. It could be explained by the crystallographic preferred orientation (CPO) developing during convection. The likely phase is hexagonal close-packed iron (hcp), alloyed with nickel and some lighter elements. Here we use high energy synchrotron X-rays to study CPO in Fe-9wt%Si, uniaxially compressed in a diamond anvil cell in radial geometry. The experiments reveal that strong preferred orientation forms in the low-pressure body-centred cubic (bcc) phase that appears to be softer than pure iron. CPO is attributed to dominant {110}<111> slip. The onset of the bcc→hcp transition occurs at a pressure of ≈15 GPa, and the alloy remains in a two phase bcc + hcp state up to 40 GPa. The hcp phase forms first with a distinct {11$\bar{2}$0} maximum perpendicular to compression. Modelling shows that this is a transformation texture, which can be described by Burgers orientation relationship with variant selection. Experimental results suggest that bcc grains oriented with <100> parallel to compression transform into hcp first. The CPO of the hcp changes only slowly during further pressure and deviatoric stress increase at ambient temperature. After heating to 1600 K, a change in the hcp CPO is observed with alignment of (0001) planes perpendicular to compression that can be interpreted as dominant (0001)<11$\bar{2}$0> slip, combined with {10$\bar{1}$2}<$\bar{1}$011> mechanical twinning, which is similar to the deformation modes suggested previously for pure hcp iron at inner core conditions.
Code to organize, prepare, and segment multi-scan, larger-than-memory, volumetric dataset. Deep learning segmentation uses pytorch-lightning for code organization, logging and metrics, and multi-GPU parallelization of training and prediction. Also uses moani for volumetric data augmentation and model implementation (currently U-Net++). Training can be done slice-by-slice with 2D models trained on 2D patches either aligned with or perpendicular to the electrode plane, or with 3D models on 3D patches. Sparse labels can be used to minimize required hand-labeling -- less than 0.1% of the entire dataset was labeled. Also includes code for data preparation (alignment/cropping/patching).
BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.
The impact of the finite grid size in SOLPS-ITER edge plasma simulations is assessed for JET H-mode discharges with a metal wall. For a semi-horizontal divertor configuration it is shown that the separatrix density is at least 30% higher when a narrow scrape-off layer (SOL) grid width is chosen in SOLPS-ITER compared to the case for which the SOL grid width is maximised. The density increase is caused by kinetic neutrals being not confined inside the divertor region because of the reduced extent of the plasma grid. In this case, an enhanced level of reflections of energetic neutrals at the low-field side (LFS) metal divertor wall is observed. This leads to a shift of the ionisation source further upstream which must be accounted for as a numerical artefact. An overestimate in the cooling at the divertor entrance is observed in this case, identified by a reduced heat flux decay parameters λ q div . Otherwise and further upstream the mid-plane heat decay length λ q parameter is not affected by any change in divertor dissipation. This confirms the assumptions made for the ITER divertor design studies, i.e. that λ q upstream is essentially set by the assumptions for the ratio radial to parallel heat conductivity. It is also shown that even for attached conditions the decay length relations λ ne > λ Te > λ q hold in the near-SOL upstream. Thus for interpretative edge plasma simulations one must take the (experimental) value of λ ne into account, rather than λ q , as the former actually defines the required minimum upstream SOL grid extent.
Fe 4 Si 2 Sn 7 O 16 hosts an undistorted kagome lattice of Fe 2+ (3d 6 , S = 2) ions. We present results of bulk magnetization and Sn nuclear magnetic resonance (NMR) measurements on an oriented Fe 4 Si 2 Sn 7 O 16 powder sample oriented in geometries parallel (∥) and perpendicular (⊥) to the external applied magnetic field used for orienting the powder (Bori). The bulk susceptibility χ shows a broad peak at T N ~ 3 K associated with antiferromagnetic ordering. NMR spectra indicate the presence of planar anisotropy in the kagome planes. From an analysis of the static NMR shift (K) and dynamic spin-lattice relaxation rate (1/T 1 ) we conclude the presence of dominant magnetic fluctuations in the kagome planes. For the ∥ orientation, K scales linearly with the bulk susceptibility for temperatures down to ~ 4 K, while in the ⊥ orientation K starts to deviate strongly below T ~ 30 K. We associate this deviation with the onset of spin-tilting towards the kagome planes. These correlations are also reflected in the 1/T 1 data for the ∥ orientation, which starts to decrease below T ~ 30 K. Here, in this correlated regime, T N < T < ~ 30 K, we discuss the formation of positive chiral spin correlations in the kagome planes.
: As the rapid growth of marine renewable energy and off-shore wind energy, there have been concerns that the noises generated from construction and operation of the devices may interfere marine animals’ communication. In this research, a underwater sound model is developed to simulate sound prorogation generated by marine-hydrokinetic energy (MHK) devices or offshore wind (OSW) energy platforms. Finite volume and finite difference methods are developed to solve the 3D Helmholtz equation of sound propagation in the coastal environment. For finite volume method, the grid system consists of triangular grids in horizontal plane and sigma-layers in vertical dimension. A 3D sparse matrix solver with complex coefficients is formed for solving the resulting acoustic pressure field. The Complex Shifted Laplacian Preconditioner (CSLP) method is applied to efficiently solve the matrix system iteratively with MPI parallelization using a high performance cluster. The sound model is then coupled with the Finite Volume Community Ocean Model (FVCOM) for simulating sound propagation generated by human activities in a range-dependent setting, such as offshore wind energy platform constructions and tidal stream turbines. As a proof of concept, initial validation of the finite difference solver is presented for two coastal wedge problems. Validation of finite volume method will be reported separately.