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 613 records · Page 34

Instability Effects in F(R)-Modified Gravity and in Gravitational Baryogenesis

The review is devoted to the development and subsequent application of methods for studying the stability of solutions of higher-order differential equations arising in various models of modified gravity and baryosynthesis. Major attention is paid to F (R )-modified theories proposed to describe the accelerated expansion of the universe and to baryon asymmetry generation in the framework of the spontaneous and gravitational baryogenesis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Scalable Preconditioners for Structure Preserving Discretizations of Maxwell Equations in First Order Form

Here, we report multiple physical time-scales can arise in electromagnetic simulations when dissipative effects are introduced through boundary conditions, when currents follow external time-scales, and when material parameters vary spatially. In such scenarios, the time-scales of interest may be much slower than the fastest time-scales supported by the Maxwell equations, therefore making implicit time integration an efficient approach. The use of implicit temporal discretizations results in linear systems in which fast time-scales, which severely constrain the stability of an explicit method, can manifest as so-called stiff modes. This study proposes a new block preconditioner for structure preserving (also termed physics compatible) discretizations of the Maxwell equations in first order form. The intent of the preconditioner is to enable the efficient solution of multiple-time-scale Maxwell type systems. An additional benefit of the developed preconditioner is that it requires only a traditional multigrid method for its subsolves and compares well against alternative approaches that rely on specialized edge-based multigrid routines that may not be readily available. Lastly, results demonstrate parallel scalability at large electromagnetic wave CFL numbers on a variety of test problems.

97 MATHEMATICS AND COMPUTING↗

CD40/anti-CD40 antibody complexes which illustrate agonist and antagonist structural switches

CD40 is a 48 kDa type I transmembrane protein that is constitutively expressed on hematopoietic cells such as dendritic cells, macrophages, and B cells. Engagement of CD40 by CD40L expressed on T cells results in the production of proinflammatory cytokines, induces T helper cell function, and promotes macrophage activation. The involvement of CD40 in chronic immune activation has resulted in CD40 being proposed as a therapeutic target for a range of chronic inflammatory diseases. CD40 antagonists are currently being explored for the treatment of autoimmune diseases and several anti-CD40 agonist mAbs have entered clinical development for oncological indications. To better understand the mode of action of anti-CD40 mAbs, we have determined the x-ray crystal structures of the ABBV-323 (anti-CD40 antagonist, ravagalimab) Fab alone, ABBV-323 Fab complexed to human CD40 and FAB516 (anti-CD40 agonist) complexed to human CD40. These three crystals structures 1) identify the conformational CD40 epitope for ABBV-323 recognition 2) illustrate conformational changes which occur in the CDRs of ABBV-323 Fab upon CD40 binding and 3) develop a structural hypothesis for an agonist/antagonist switch in the LCDR1 of this proprietary class of CD40 antibodies. The structure of ABBV-323 Fab demonstrates a unique method for antagonism by stabilizing the proposed functional antiparallel dimer for CD40 receptor via novel contacts to LCDR1, namely residue position R32 which is further supported by a closely related agonist antibody FAB516 which shows only monomeric recognition and no contacts with LCDR1 due to a mutation to L32 on LCDR1. These data provide a structural basis for the full antagonist activity of ABBV-323.

59 BASIC BIOLOGICAL SCIENCES↗

Stable tuning of photorefractive microcavities using an auxiliary laser

Cavity nonlinear optics enables intriguing physical phenomena to occur at micro- or nano-scales with modest input powers. While this enhances capabilities in applications such as comb generation, frequency conversion, and quantum optics, undesired nonlinear effects including photorefraction and thermal bistability are exacerbated. In this Letter, we propose and demonstrate a highly effective method of achieving cavity stabilization using an auxiliary laser for controlling photorefraction in a z-cut periodically poled lithium niobate (LN) microcavity system. Our numerical study accurately models the photorefractive effect under high input powers, guiding future analyses and development of LN microcavity systems.

42 ENGINEERING↗

Moment tensor analysis of very shallow sources

An issue for moment tensor (MT) inversion of shallow seismic sources is that some components of the Green’s functions have vanishing amplitudes at the free surface, which can result in bias in the MT solution. The effects of the free surface on the stability of the MT method become important as we continue to investigate and improve the capabilities of regional full MT inversion for source–type identification and discrimination. It is important to understand free–surface effects on discriminating shallow explosive sources for nuclear monitoring purposes. It may also be important in natural systems that have very shallow seismicity, such as volcanic and geothermal systems. We examine the effects of the free surface on the MT via synthetic testing and apply the MT–based discrimination method to three quarry blasts from the HUMMING ALBATROSS experiment. These shallow chemical explosions at ~10 m depth and recorded up to several kilometers distance represent rather severe source–station geometry in terms of free–surface effects. We show that the method is capable of recovering a predominantly explosive source mechanism, and the combined waveform and first–motion method enables the unique discrimination of these events. Furthermore, recovering the design yield using seismic moment estimates from MT inversion remains challenging, but we can begin to put error bounds on our moment estimates using the network sensitivity solution technique.

58 GEOSCIENCES↗

Project on Nuclear Gaming: Investigating the Impact of Nuclear Weapons Capabilities on Deterrence: May 2-3, 2018 Workshop Summary

In the future, it may be possible for states to construct nuclear weapons with alternative effects that are optimized for specific military purposes. These nuclear weapon concepts might employ features such as enhanced radiation outputs (ER), electromagnetic pulse (EMP), or other approaches such as high-precision, low-yield systems (HPLY) that reduce collateral damage, but improve lethality to targets. Some have speculated that such weapons could have the effect of lowering the threshold of use for nuclear weapons by reducing the negative incentives to their use or by providing new military applications for nuclear weapons. A nuclear arsenal that includes weapons with alternate effect regimes may also have an impact on the calculus of escalation in a conflict scenario. This project of the University of California, Berkeley (UCB), Sandia National Laboratories (SNL), and Lawrence Livermore National Laboratory (LLNL) investigates the potential impacts of alternate nuclear weapons effects regimes on deterrence and strategic stability using experimental gaming methods. The purpose of this document is to capture the research questions, assumptions, scenario characteristics, and other details that provide the framework used to address this challenge presented during the first workshop discussion hosted by the Center for Global Security Research at LLNL with subject matter experts from government, academia, and the private sector. The report also summarizes feedback from participants. Topics discussed during the workshop included alternative effects of nuclear weapons, existing approaches for modeling and measuring conflict, and the state of the art in experimental game design.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

One-Step Spark Plasma Erosion Processing of Carbon-Coated Sn-Si Nanoparticles for Lithium-Ion Battery Anodes

High density portable energy storage is desirable owing to the energy requirements of portable electronics and electric vehicles. The Li-ion battery’s high energy density could be even further improved through the utilization of alternative materials (instead of carbon) for the anode, such as Sn or Si. Nonetheless, the large volume expansion upon lithiation, up to ~300% for Li 22 Si 5 , causes pulverization and rapid capacity degradation during cycling. Sn also forms a Li 22 Sn 5 compound with the equivalent stoichiometric Li capacity but with enhanced ductility. Nano-sized Si and Sn have demonstrated distinctive nanoscale properties, facilitating the retention of higher capacities, particularly when coated with carbon, which improves mechanical stability. To date, the methods of synthesizing coated Si, Sn, or Si-Sn alloyed nanoparticles are complicated, costly, and not readily scalable to meet the demands of cost-effective manufacturing. Spark plasma erosion in a hydrocarbon dielectric has been explored as a one-step process to produce Sn-Si alloy nanoparticles coated with a thin carbon film, offering a scalable and cost-effective processing route. The resulting Sn-Si particles exhibited a bi-modal size distribution at ~5 nm and ~500 nm and were carbon-coated, as intended, from the hydrocarbon dielectric breakdown. The spark-eroded nanoparticles were thoroughly characterized using TEM/EDS, XPS, AES, SSNMR, and TGA, and their improved electrochemical performance was assessed through half-cell experiments.

25 ENERGY STORAGE↗

Fission Product Diffusion in Silicon Carbide: As Revealed by Computer Simulations and Experimental Measurements

The diffusion and release of fission products through silicon carbide in the tri-structural isotropic (TRISO) fuel particles remain unsolved for decades. The underlying mechanism is quite challenging to be determined. To help unveil the mysterious story, the current work applies molecular dynamics method to show the stability of silver, palladium, ruthenium and iodine as an interstitial and their atomic diffusion along coincident site lattice (CSL) boundary, especially ?3 grain boundary (GB). The major finding presents a much faster diffusion along GB than in bulk for all elements considered. The reasonably close estimate to experiments and simulations where available has confirmed the important role of grain boundary diffusion of Ag and Pd in SiC. However, the discrepancy addressed in Ag with measurements from fuel studies suggest a more complicated mechanism, which might be in correlation with high energy grain boundaries or the presence of crack. The subsequent characterization of Ru and I distribution in SiC-PyC-SiC diffusion couples, which have been ion irradiated at 900?C to 10 dpa and 20 dpa, has performed by secondary ion mass spectrometry (SIMS) analysis. The experiment measurements correlate well with the grain boundary diffusion by simulation, which provide further evidences that the grain boundary diffusion cannot be neglected once the fission products are accessible at the grain boundary.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

All SiC/SiC cladding irradiation in LWR environments

The Muroran Institute of Technology has been in charge of developing SiC/SiC claddings by the DEMO-NITE method, and to evaluate stability in reactor water environment, in high temperature steam environment and in neutron exposure environment. Those issues have been studied during five years and carried out in 3 experiments (Scarlet-1, Scarlet-2 and Inspire-1) irradiated in the Halden Boiling Water Reactor (HBWR) in Norway. The main conclusions follow. The corrosion of SiC/SiC is accelerated by the neutron irradiation. The corrosion of SiC/SiC significantly depends on the condition of water chemistry, and the corrosion of SiC/SiC is drastically decreased under hydrogen water chemistry. Under boiling water condition, the corrosion rate of SiC/SiC is slightly increased. The significant degradation at a joint was confirmed by the galvanic corrosion caused by the joining of dissimilar materials. As expected, there is no degradation of full SiC/SiC segment prepared by diffusion bonding. Microstructural change of CVD-SiC coat on SiC/SiC tubes was not observed. The effectiveness of CVD (Chemical Vapor Deposition) SiC coat as environmental barrier coating is confirmed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Geometric multigrid for an implicit-time immersed boundary method

The immersed boundary (IB) method is an approach to fluid-structure interaction that uses Lagrangian variables to describe the deformations and resulting forces of the structure and Eulerian variables to describe the motion and forces of the fluid. Explicit time stepping schemes for the IB method require solvers only for Eulerian equations, for which fast Cartesian grid solution methods are available. Such methods are relatively straightforward to develop and are widely used in practice but often require very small time steps to maintain stability. Implicit-time IB methods permit the stable use of large time steps, but efficient implementations of such methods require significantly more complex solvers that effectively treat both Lagrangian and Eulerian variables simultaneously. Moreover, several different approaches to solving the coupled Lagrangian-Eulerian equations have been proposed, but a complete understanding of this problem is still emerging. This paper presents a geometric multigrid method for an implicit-time discretization of the IB equations. This multigrid scheme uses a generalization of box relaxation that is shown to handle problems in which the physical stiffness of the structure is very large. Numerical examples are provided to illustrate the effectiveness and efficiency of the algorithms described herein. Finally, these tests show that using multigrid as a preconditioner for a Krylov method yields improvements in both robustness and efficiency as compared to using multigrid as a solver. They also demonstrate that with a time step 100–1000 times larger than that permitted by an explicit IB method, the multigrid-preconditioned implicit IB method is approximately 50–200 times more efficient than the explicit method.

97 MATHEMATICS AND COMPUTING↗

Stability analysis and discretization of A –$Φ$ time domain integral equations for multiscale electromagnetics

The growth of applications at the intersection between electromagnetic and quantum physics is necessitating the creation of novel computational electromagnetic solvers. Work in this paper presents a new set of time domain integral equations (TDIEs) formulated directly in terms of the magnetic vector and electric scalar potentials that can be used to meet many of the requirements of this emerging area. Stability for this new set of TDIEs is achieved by leveraging an existing rigorous functional framework that can be used to determine suitable discretization approaches to yield stable results in practice. The basics of this functional framework are reviewed before it is shown in detail how it may be applied in developing the TDIEs of this work. Numerical results are presented which validate the claims of stability and accuracy of this method over a wide range of frequencies where traditional methods would fail.

97 MATHEMATICS AND COMPUTING↗

Optimal Power Hardware-in-the-Loop Interfacing: Applying Modern Control For Design and Verification of High-Accuracy Interfaces

Here, we develop a novel approach to designing power hardware-in-the-loop (PHIL) simulation interfaces that maximizes simulation bandwidth and accuracy by leveraging a modern control framework that explicitly considers objectives on accuracy and can automatically synthesize an optimal controller that meets these objectives while stabilizing the closed-loop system. The method developed improves upon common approaches to PHIL interface design that typically involve multiple design steps for manual compensation and stabilization design that can result in interfaces that are stable but have suboptimal bandwidth and accuracy. The approach developed is general and can be applied to most PHIL system configurations. The modeling framework allows for the inclusion of scaling factors and hardware-under-test current injection models that are common to practical PHIL simulations and have significant impacts on stability. We also present practical methods and metrics for verifying the absolute accuracy of PHIL simulation results without relying on relative comparisons to potentially inaccurate models or previous simulation results. We demonstrate the accuracy evaluation method and show the improved performance achievable when using the optimal PHIL interface design approach in an experimental case study involving a 100-kVA battery inverter.

42 ENGINEERING↗

Efficient and robust stability analysis in the internal energy, volume, and moles (UVN) space

Hydrodynamic simulation of fluids in vessels is often based on the internal energy (U), volume (V), and the number of moles of all components ( re N = N 1 , N 2 , ..., N n ) as the working variables to determine the thermodynamic state by finding the global maximum of entropy. The procedure of finding the equilibrium state may be divided into two steps: stability analysis and phase-split computations. Stability analysis is performed first, and the phase-split computation proceeds if the system is unstable. Most of the past work on UVN space has been based on phase-split computations, while the stability analysis is barely investigated. Furthermore, previous studies for the stability analysis use only Newton's method to solve the non-linear algebraic equations. We present an efficient and robust approach for stability analysis where UVN is specified. The successive substitution iteration (SSI) is used to provide good initial guesses for Newton's method. The proposed approach results in a reduced number of unknowns and does not require a large number of iterations to achieve convergence in Newton's method. Overall, our proposed formulation is compatible with different equations of state and is applicable to both pure component and multi-components. The robustness and efficiency of the algorithm for stability analysis are demonstrated in various examples.

02 PETROLEUM↗

Robust H{sub ∞} control for singular time-delay systems with saturating actuators via static output feedback

The problem of H{sub ∞} robust exponential stability and static output feedback stabilization of uncertain singular time-delay systems with saturating actuators are investigated. First, the delay-dependent condition is given to guarantee that the singular time-delay systems are not only exponential admissible, but also satisfies H{sub ∞} performance γ. Moreover, the estimation of stability region and the design method of H{sub ∞} static output feedback controller are given by solving a convex optimization problem. Some simulation examples are provided to demonstrate the merit of the proposed design methods.

97 MATHEMATICS AND COMPUTING↗

A Strong Stability Preserving Analysis for Explicit Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions

High-order strong stability preserving (SSP) time discretizations are usually needed to ensure the nonlinear (and sometimes non-inner-product) strong stability properties of spatial discretizations specially designed for the solution of hyperbolic PDEs. Multi-derivative time-stepping methods have recently been increasingly used for evolving hyperbolic PDEs, and the strong stability properties of these methods are of interest. In our prior work we explored time discretizations that preserve the strong stability properties of spatial discretizations coupled with forward Euler and a second-derivative formulation. Yet, many spatial discretizations do not satisfy strong stability properties when coupled with this second-derivative formulation, but rather with a more natural Taylor series formulation. Here, we demonstrate sufficient conditions for an explicit two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and Taylor series formulation. We call these strong stability preserving Taylor series (SSP-TS) methods. We also prove that the maximal order of SSP-TS methods is p=6, and define an optimization procedure that allows us to find such SSP methods. Several types of these methods are presented and their efficiency compared. Finally, these methods are tested on several PDEs to demonstrate the benefit of SSP-TS methods, the need for the SSP property, and the sharpness of the SSP time-step in many cases.

97 MATHEMATICS AND COMPUTING↗

A second-order ensemble method based on a blended backward differentiation formula timestepping scheme for time-dependent Navier-Stokes equations

Here, we present a second-order ensemble method based on a blended three-step backward differentiation formula (BDF) timestepping scheme to compute an ensemble of Navier–Stokes equations. Compared with the only existing second-order ensemble method that combines the two-step BDF timestepping scheme and a special explicit second-order Adams–Bashforth treatment of the advection term, this method is more accurate with nominal increase in computational cost. We give comprehensive stability and error analysis for the method. Numerical examples are also provided to verify theoretical results and demonstrate the improved accuracy of the method.

97 MATHEMATICS AND COMPUTING↗

Nonlinear asymmetric tearing mode evolution in cylindrical geometry

The growth of a tearing mode is described by reduced MHD equations. For a cylindrical equilibrium, tearing mode growth is governed by the modified Rutherford equation, i.e., the nonlinear Δ'(w). For a low beta plasma without external heating, Δ'(w) can be approximately described by two terms, Δ' ql (w), Δ'A(w). In this work, we present a simple method to calculate the quasilinear stability index Δ'ql rigorously, for poloidal mode number m ≥ 2. Δ' ql is derived by solving the outer equation through the Frobenius method. Δ'ql is composed of four terms proportional to: constant Δ' 0 , w, wlnw, and w2. Δ' A is proportional to the asymmetry of island that is roughly proportional to w. The sum of Δ' ql and Δ' A is consistent with the more accurate expression calculated perturbatively. The reduced MHD equations are also solved numerically through a 3D MHD code M3D-C1. The analytical expression of the perturbed helical flux and the saturated island width agree with the simulation results. Lastly, it is also confirmed by the simulation that the Δ' A has to be considered in calculating island saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗