Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stabilization method”

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 829 records · Page 46

Melting behaviour of gold-platinum nanoalloy clusters by molecular dynamics simulations

The melting behavior of bimetallic gold-platinum nanoclusters is studied by applying Brownian-type isothermal molecular dynamics (MD) simulation, a program modified from the cubic coupling scheme (CCS). The process begins with the ground-state structures obtained from global minimum search algorithm and proceeds with the investigation of the effect of temperature on the thermal properties of gold-platinum nanoalloy clusters. N-body Gupta potential has been employed in order to account for the interactions between gold and platinum atoms. The ground states of the nanoalloy clusters, which are core-shell segregated, are heated until they become thermally segregated. The detailed melting mechanism of the nanoalloy clusters is studied via this approach to provide insight into the thermal stability of the nanoalloy clusters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Combined First-Principles and Experimental Investigation into the Reactivity of Codeposited Chromium–Carbon under Pressure

High-pressure synthesis in the diamond anvil cell suffers from the lack of a general approach for the control of precursor stoichiometry and homogeneity. Here, we present results from a new method we have developed that uses magnetron cosputtering to prepare stoichiometrically precise and atomically mixed amorphous films of Cr:C. Laser-heated diamond anvil cell experiments carried out on a flake of this sample at pressures between 13.5 and 24.3 GPa lead to the observation of Cr 3 C (Pnma) over the entire pressure range–in good agreement with our in-house theoretical predictions–but also reveal two other metastable phases that were not expected: a novel monoclinic chromium carbide phase and the NaCl-type CrC (Fm3̅m) phase. The unexpected stability of CrC is investigated by using first-principles methods, revealing a large stabilizing effect tied to substoichiometry at the carbon site. These results offer an important case study into the current limitations of crystal structure prediction methods with regard to phase complexity and bolster the growing need for advanced theoretical approaches that can more completely survey experimentally unexplored phase space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Imaginary time density-density correlations for two-dimensional electron gases at high density

We evaluate imaginary time density-density correlation functions for two-dimensional homogeneous electron gases of up to 42 particles in the continuum using the phaseless auxiliary field quantum Monte Carlo method. We use periodic boundary conditions and up to 300 plane waves as basis set elements. We show that such methodology, once equipped with suitable numerical stabilization techniques necessary to deal with exponentials, products, and inversions of large matrices, gives access to the calculation of imaginary time correlation functions for medium-sized systems. We discuss the numerical stabilization techniques and the computational complexity of the methodology and we present the limitations related to the size of the systems on a quantitative basis. We perform the inverse Laplace transform of the obtained density-density correlation functions, assessing the ability of the phaseless auxiliary field quantum Monte Carlo method to evaluate dynamical properties of medium-sized homogeneous fermion systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Upconverting Nanoparticles as Optical Sensors of Nano- to Micro-Newton Forces

Mechanical forces affect a myriad of processes, from bone growth to material fracture to touch-responsive robotics. While nano- to micro-Newton forces are prevalent at the microscopic scale, few methods have the nanoscopic size and signal stability to measure them in vivo or in situ. Here, we develop an optical force-sensing platform based on sub-25 nm NaYF 4 nanoparticles (NPs) doped with Yb 3+ , Er 3+ , and Mn 2+ . The lanthanides Yb 3+ and Er 3+ enable both photoluminescence and upconversion, while the energetically coupled d-metal Mn 2+ adds force tunability through its crystal field sensitivity. IN using a diamond anvil cell to exert up to 3.5 GPa pressure or ~10 μN force per particle, we track stress-induced spectral responses. The red (660 nm) to green (520, 540 nm) emission ratio varies linearly with pressure, yielding an observed color change from orange to red for α-NaYF 4 and from yellow–green to green for d-metal optimized β-NaYF 4 when illuminated in the near infrared. We record consistent readouts over multiple pressure cycles and hours of illumination. With the nanoscopic size, a dynamic range of 100 nN to 10 μN, and photostability, these nanoparticles lay the foundation for visualizing dynamic mechanical processes, such as stress propagation in materials and force signaling in organisms.

36 MATERIALS SCIENCE↗

Design and implementation of grid multi-scroll fractional-order chaotic attractors

This paper proposes a novel approach for generating multi-scroll chaotic attractors in multi-directions for fractional-order (FO) systems. The stair nonlinear function series and the saturated nonlinear function are combined to extend equilibrium points with index 2 in a new FO linear system. With the help of stability theory of FO systems, stability of its equilibrium points is analyzed, and the chaotic behaviors are validated through phase portraits, Lyapunov exponents, and Poincaré section. Choosing the order 0.96 as an example, a circuit for generating 2-D grid multiscroll chaotic attractors is designed, and 2-D 9 × 9 grid FO attractors are observed at most. Numerical simulations and circuit experimental results show that the method is feasible and the designed circuit is correct.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Influence of Modification on Crystal Lattice Stability of Austenite in Stainless Steel

Using the methods of electron diffraction microscopy and X-ray diffraction analysis, the influence of alloying of the austenitic steel, Grade 110H13, with chromium and vanadium, as well as high-melting, ultrafine-grained TiO{sub 2}, ZrO{sub 2} powders and Na{sub 3}AlF{sub 6} cryolite on its structural-phase state and microstructure is investigated. It is shown that the matrix of non-modified steel is completely austenitic and consists of an iron-based solid solution and the interstitial (C, N, O and other) and substitutional (Cr, V and other) atoms simultaneously. Alloying with chromium and vanadium changes neither its phase composition nor defect structure, while alloy modification results in qualitatively new structural features: γ → ε-transformation, high-intensity microtwinning, defect structure changes, and a sharp increase in the scalar dislocation density. The features of the deformation-induced microtwinning and ε-martensite plates identified in the modified steel promote revealing additional microtwin systems in the matrix γ-phase, which result in structural changes making it possible to classify it as a γ′-phase. It is found out that an introduction of modifying additions gives rise to the following sequence of structural-phase transformations: γ→γ′→(γ′ +ε). The experimental data obtained demonstrate that as a result of modification the crystal lattice transits into a low-stability state. This transition is accompanied by marked structural-phase changes consisting in the formation of several microtwin systems and γ → ε-transformation. These structural-phase changes in the modified steel are due to the crystal-lattice transition into the low-stability state, followed by new structural-phase alterations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Early-stage human islet amyloid polypeptide aggregation: Mechanisms behind dimer formation

Amyloid aggregates of human islet amyloid polypeptide (hIAPP or human amylin) have long been implicated in the development of type II diabetes. While hIAPP is known to aggregate into amyloid fibrils, it is the early-stage prefibrillar species that have been proposed to be cytotoxic. A detailed picture of the early-stage aggregation process and relevant intermediates would be valuable in the development of effective therapeutics. Here, we use atomistic molecular dynamics simulations with a combination of enhanced sampling methods to examine the formation of the hIAPP dimer in water. Bias-exchange metadynamics calculations reveal relative conformational stabilities of the hIAPP dimer. Finite temperature string method calculations identify pathways for dimer formation, along with relevant free energy barriers and intermediate structures. We show that the initial stages of dimerization involve crossing a substantial free energy barrier to form an intermediate structure exhibiting transient beta-sheet character, before proceeding to form an entropically stabilized dimer structure.

59 BASIC BIOLOGICAL SCIENCES↗

Stability and response of trapped solitary wave solutions of coupled nonlinear Schrödinger equations in an external, $\mathcal{PT}$- and supersymmetric potential

In this work, we present trapped solitary wave solutions of a coupled nonlinear Schrödinger (NLS) system in 1 + 1 dimensions in the presence of an external, supersymmetric and complex $\mathcal{PT}$-symmetric potential. The Schrödinger system this work focuses on possesses exact solutions whose existence, stability, and spatio-temporal dynamics are investigated by means of analytical and numerical methods. Two different variational approximations are considered where the stability and dynamics of the solitary waves are explored in terms of eight and twelve time-dependent collective coordinates (CCs). We find regions of stability for specific potential choices as well as analytic expressions for the small oscillation frequencies in the CC approximation. Our findings are further supported by performing systematic numerical simulations of the NLS system.

97 MATHEMATICS AND COMPUTING↗

Thermodynamic and mechanical stabilities of α- and β-Ta{sub 4}AlC{sub 3} via first-principles investigations

Recently, it has been predicted that Ta{sub 4}AlC{sub 3} is likely to exhibit an (α–β) polymorphic transformation at temperatures above 1873 K. However, recent X-ray diffraction and transmission electron microscopy experiments suggest on the other hand that the α phase remains stable up to temperatures close to the limit of experimental capabilities and no transition has yet been observed. While the matter has already been settled experimentally, in this work, we re-investigate the phase stability problem in Ta{sub 4}AlC{sub 3} by using first-principles methods. The study was carried out by considering both thermodynamic and mechanical stabilities of the Ta{sub 4}AlC{sub 3} polymorphs. Particularly, finite-temperature Gibbs free energies and elastic properties of the polymorphs were calculated using density functional theory. Calculation results reveal that the α phase continue to be stable even at temperatures exceeding 1875 K, which is in agreement with experimental results reported in literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Scalable implicit incompressible resistive MHD with stabilized FE and fully-coupled Newton–Krylov-AMG

Here, we discuss that the computational solution of the governing balance equations for mass, momentum, heat transfer and magnetic induction for resistive magnetohydrodynamics (MHD) systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena, as well as the significant range of time- and length-scales that the interactions of these physical mechanisms produce. This paper explores the development of a scalable, fully-implicit stabilized unstructured finite element (FE) capability for 3D incompressible resistive MHD. The discussion considers the development of a stabilized FE formulation in context of the variational multiscale (VMS) method, and describes the scalable implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton–Krylov methods, which are preconditioned using fully-coupled algebraic multilevel preconditioners. These preconditioners are shown to enable a robust, scalable and efficient solution approach for the large-scale sparse linear systems generated by the Newton linearization. Verification results demonstrate the expected order-of-accuracy for the stabilized FE discretization. The approach is tested on a variety of prototype problems, that include MHD duct flows, an unstable hydromagnetic Kelvin–Helmholtz shear layer, and a 3D island coalescence problem used to model magnetic reconnection. Initial results that explore the scaling of the solution methods are also presented on up to 128K processors for problems with up to 1.8B unknowns on a CrayXK7.

97 MATHEMATICS AND COMPUTING↗

Electrochemical Stability of Li 10 GeP 2 S 12 and Li 7 La 3 Zr 2 O 12 Solid Electrolytes

The electrochemical stability window of solid electrolyte is overestimated by the conventional experimental method using a Li/electrolyte/inert metal semiblocking electrode because of the limited contact area between solid electrolyte and inert metal. Since the battery is cycled in the overestimated stability window, the decomposition of the solid electrolyte at the interfaces occurs but has been ignored as a cause for high interfacial resistances in previous studies, limiting the performance improvement of the bulk-type solid-state battery despite the decades of research efforts. Thus, there is an urgent need to identify the intrinsic stability window of the solid electrolyte. The thermodynamic electrochemical stability window of solid electrolytes is calculated using first principles computation methods, and an experimental method is developed to measure the intrinsic electrochemical stability window of solid electrolytes using a Li/electrolyte/electrolyte-carbon cell. The most promising solid electrolytes, Li10GeP2S12 and cubic Li-garnet Li7La3Zr2O12, are chosen as the model materials for sulfide and oxide solid electrolytes, respectively. The results provide valuable insights to address the most challenging problems of the interfacial stability and resistance in high-performance solid-state batteries.

25 ENERGY STORAGE↗

Acceleration of GPU-based Krylov solvers via data transfer reduction

Krylov subspace iterative solvers are often the method of choice when solving large sparse linear systems. At the same time, hardware accelerators such as graphics processing units continue to offer significant floating point performance gains for matrix and vector computations through easy-to-use libraries of computational kernels. However, as these libraries are usually composed of a well optimized but limited set of linear algebra operations, applications that use them often fail to reduce certain data communications, and hence fail to leverage the full potential of the accelerator. In this study, we target the acceleration of Krylov subspace iterative methods for graphics processing units, and in particular the Biconjugate Gradient Stabilized solver that significant improvement can be achieved by reformulating the method to reduce data-communications through application-specific kernels instead of using the generic BLAS kernels, e.g. as provided by NVIDIA’s cuBLAS library, and by designing a graphics processing unit specific sparse matrix-vector product kernel that is able to more efficiently use the graphics processing unit’s computing power. Furthermore, we derive a model estimating the performance improvement, and use experimental data to validate the expected runtime savings. Finally, considering that the derived implementation achieves significantly higher performance, we assert that similar optimizations addressing algorithm structure, as well as sparse matrix-vector, are crucial for the subsequent development of high-performance graphics processing units accelerated Krylov subspace iterative methods.

97 MATHEMATICS AND COMPUTING↗

GEMPIC: geometric electromagnetic particle-in-cell methods

We present a novel framework for finite element particle-in-cell methods based on the discretization of the underlying Hamiltonian structure of the Vlasov–Maxwell system. We derive a semi-discrete Poisson bracket, which retains the defining properties of a bracket, anti-symmetry and the Jacobi identity, as well as conservation of its Casimir invariants, implying that the semi-discrete system is still a Hamiltonian system. In order to obtain a fully discrete Poisson integrator, the semi-discrete bracket is used in conjunction with Hamiltonian splitting methods for integration in time. Techniques from finite element exterior calculus ensure conservation of the divergence of the magnetic field and Gauss’ law as well as stability of the field solver. The resulting methods are gauge invariant, feature exact charge conservation and show excellent long-time energy and momentum behaviour. Due to the generality of our framework, these conservation properties are guaranteed independently of a particular choice of the finite element basis, as long as the corresponding finite element spaces satisfy certain compatibility conditions.

Physics↗

Thermodynamical stability of the Bardeen black hole

We analyze the stability of the regular magnetic Bardeen black hole both thermodynamically and dynamically. For the thermodynamical analysis we consider a microcanonical ensemble and apply the turning point method. This method allows to decide a change in stability (or instability) of a system, requiring only the assumption of smoothness of the area functional. The dynamical stability is asserted using criteria based on the signs of the Lagrangian and its derivatives. It turns out from our analysis that the Bardeen black hole is both thermodynamically and dynamically stable.

36 MATERIALS SCIENCE↗

An efficient parallel implementation of explicit multirate Runge–Kutta schemes for discontinuous Galerkin computations

Although explicit time integration schemes require small computational efforts per time step, their efficiency is severely restricted by their stability limits. Indeed, the multi-scale nature of some physical processes combined with highly unstructured meshes can lead some elements to impose a severely small stable time step for a global problem. Multirate methods offer a way to increase the global efficiency by gathering grid cells in appropriate groups under local stability conditions. These methods are well suited to the discontinuous Galerkin framework. The parallelization of the multirate strategy is challenging because grid cells have different workloads. The computational cost is different for each sub-time step depending on the elements involved and a classical partitioning strategy is not adequate any more. In this paper, we propose a solution that makes use of multi-constraint mesh partitioning. It tends to minimize the inter-processor communications, while ensuring that the workload is almost equally shared by every computer core at every stage of the algorithm. Particular attention is given to the simplicity of the parallel multirate algorithm while minimizing computational and communication overheads. Our implementation makes use of the MeTiS library for mesh partitioning and the Message Passing Interface for inter-processor communication. Performance analyses for two and three-dimensional practical applications confirm that multirate methods preserve important computational advantages of explicit methods up to a significant number of processors.

97 MATHEMATICS AND COMPUTING↗

De Sitter decays to infinity

Bubbles of nothing are a class of vacuum decay processes present in some theories with compactified extra dimensions. We investigate the existence and properties of bubbles of nothing in models where the scalar pseudomoduli controlling the size of the extra dimensions are stabilized at positive vacuum energy, which is a necessary feature of any realistic model. We map the construction of bubbles of nothing to a four-dimensional Coleman-De Luccia problem and establish necessary conditions on the asymptotic behavior of the scalar potential for the existence of suitable solutions. We perform detailed analyses in the context of five-dimensional theories with metastable dS 4 x S 1 vacua, using analytic approximations and numerical methods to calculate the decay rate. We find that bubbles of nothing sometimes exist in potentials with no ordinary Coleman-De Luccia decay process, and that in the examples we study, when both processes exist, the bubble of nothing decay rate is typically faster. Our methods can be generalized to other stabilizing potentials and internal manifolds.

classical theories of gravity↗

Preparation of meta-stable phases of barium titanate by Sol-hydrothermal method

Two low-cost chemical methods of sol–gel and the hydrothermal process have been strategically combined to fabricate barium titanate (BaTiO{sub 3}) nanopowders. This method was tested for various synthesis temperatures (100 °C to 250 °C) employing barium dichloride (BaCl{sub 2}) and titanium tetrachloride (TiCl{sub 4}) as precursors and sodium hydroxide (NaOH) as mineralizer for synthesis of BaTiO{sub 3} nanopowders. The as-prepared BaTiO{sub 3} powders were investigated for structural characteristics using x-ray diffraction (XRD), scanning electron microscopy (SEM) and transmission electron microscopy (TEM). The overall analysis indicates that the hydrothermal conditions create a gentle environment to promote the formation of crystalline phase directly from amorphous phase at the very low processing temperatures investigated. XRD analysis showed phase transitions from cubic - tetragonal - orthorhombic - rhombohedral with increasing synthesis temperature and calculated grain sizes were 34 – 38 nm (using the Scherrer formula). SEM and TEM analysis verified that the BaTiO{sub 3} nanopowders synthesized by this method were spherical in shape and about 114 - 170 nm in size. The particle distribution in both SEM and TEM shows that as the reaction temperature increases from 100 °C to 250 °C, the particles agglomerate. Selective area electron diffraction (SAED) shows that the particles are crystalline in nature. The study shows that choosing suitable precursor and optimizing pressure and temperature; different meta-stable (ferroelectric) phases of undoped BaTiO{sub 3} nanopowders can be stabilized by the sol-hydrothermal method.

36 MATERIALS SCIENCE↗

Identification of multiple eigenmode growth rates in DIII-D and EAST tokamak plasmas

Here, the successful application of three-dimensional (3D) magnetohydrodynamic (MHD) spectroscopy in the stable DIII-D and EAST plasmas enables to directly extract the multi-mode plasma transfer function in tokamak experiments. The transfer function not only reveals the contribution of each dominant MHD eigenmode in the multi-model plasma response, but also quantifies the corresponding eigenvalue which is the critical stability index of MHD mode. The method performs the active detection of stable plasma by utilizing the upper and lower rows of internal coils to scan the frequency and poloidal spectrum of the applied 3D field, and reconstructing the multi-mode transfer functions model through the least square fitting of the plasma response measured by 3D-field magnetic sensors distributed at different poloidal locations. The results point to the potential development of an advanced strategy for tracking the plasma stability based on the extracted eigenvalues of stable modes. The improved understanding of the dominant eigenmodes' behavior in multi-mode plasma response through transfer function can also help to optimize the applied 3D fields for the purposes, such as the type-I edge localized mode suppression and the core stability control in future fusion reactors.

3D MHD spectroscopy↗