Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

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 379 records · Page 21

A Limiting Case of Constant Counterion Electrochemical Potentials in the Membrane for Examining Ion Transfer at Ion-Exchange Membranes and Patches

We present that ion passage through ion-exchange membranes is vital in electrodialysis desalination, batteries and fuel cells, and water splitting. Simplified models of ion transport through such membranes frequently assume complete exclusion of coions (ions with the same sign of charge as the fixed charge in the membrane) from the membrane. However, a second assumption of constant counterion electrochemical potentials across the membrane leads to simple analytical expressions for ion fluxes and transmembrane potentials. Moreover, linear corrections to account for a small membrane electrical resistance yield analytical expressions with a wider applicability. For bi-ionic potential measurements and current-induced concentration polarization at low salt concentrations, these analytical solutions match the fluxes and potentials obtained numerically without the limiting assumptions. This gives confidence in both the limiting assumptions (under appropriate conditions) and the numerical solutions. At low ion concentrations, the analytical solutions may enable rapid characterization of membrane coatings or boundary layers in solution, and such boundary layers are important in many applications of ion-exchange membranes. In fact, the assumption of complete coion exclusion is sometimes more limiting than the constraint of constant electrochemical potentials of counterions across the membrane. Remarkably, this limiting case readily yields the ion accumulation and depletion regions above “ion-exchange patches” that reside beneath a solution with an applied electric field. Finally, such regions are important for sample pre-concentration in microfluidic devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analysis of plasma instabilities and verification of the BOUT code for the Large Plasma Device

The properties of linear instabilities in the Large Plasma Device [W. Gekelman et al., Rev. Sci. Instrum. 62, 2875 (1991)] are studied both through analytic calculations and solving numerically a system of linearized collisional plasma fluid equations using the three-dimensional fluid code BOUT[M. Umansky et al., Contrib. Plasma Phys. 180, 887 (2009)], which has been successfully modified to treat cylindrical geometry. Instability drive from plasma pressure gradients and flows is considered, focusing on resistive drift waves and the Kelvin-Helmholtz and rotational interchange instabilities. A general linear dispersion relation for partially ionized collisional plasmas including these modes is derived and analyzed. For Large Plasma Device relevant profiles including strongly driven flows, it is found that all three modes can have comparable growth rates and frequencies. Detailed comparison with solutions of the analytic dispersion relation demonstrates that BOUT accurately reproduces all characteristics of linear modes in this system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Stochastic Resonance and First Arrival Time for Excitable Systems

We study the noise induced thermally activated barrier crossing of Brownian particles that hop in a piecewise linear potential. Using the exact analytic solutions and via numerical simulations not only we explore the dependence for the first passage time of a single particle but also we calculate the first arrival time for one particle out of N particles. The first arrival time decreases as the number of particles increases as expected. We then explore the thermally activated barrier crossing rate of the system in the presence of time varying signal. The dependence of signal to noise ratio SNR as well as the power amplification (η) on model parameters is explored. η and SNR depict a pronounced peak at particular noise strength. In the presence of N particles, η is considerably amplified as N steps up showing the weak periodic signal plays a vital role in controlling the noise induced dynamics of the system. Moreover, for the sake of generality, the viscous friction γ is considered to decrease exponentially when the temperature T of the medium increases (γ=Be{sup −AT}) as proposed originally by Reynolds (Philos Trans R Soc Lond 177:157, 1886).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Connecting the irreversible capacity loss in Li-ion batteries with the electronic insulating properties of solid electrolyte interphase (SEI) components.

The formation and continuous growth of a solid electrolyte interphase (SEI) layer are responsible for the irreversible capacity loss of batteries in the initial and subsequent cycles, respectively. In this article, the electron tunneling barriers from Li metal through three insulating SEI components, namely Li 2 CO 3 , LiF and Li 3 PO 4 , are computed by density function theory (DFT) approaches. Based on electron tunneling theory, it is estimated that sufficient to block electron tunneling. It is also found that the band gap decreases under tension while the work function remains the same, and thus the tunneling barrier decreases under tension and increases under compression. A new parameter, η, characterizing the average distances between anions, is proposed to unify the variation of band gap with strain under different loading conditions into a single linear function of η. An analytical model based on the tunneling results is developed to connect the irreversible capacity loss, due to the Li ions consumed in forming these SEI component layers on the surface of negative electrodes. As a result, the agreement between the model predictions and experimental results suggests that only the initial irreversible capacity loss is due to the self-limiting electron tunneling property of the SEI.

25 ENERGY STORAGE↗

A phase-field model of quasi-brittle fracture for pressurized cracks: Application to UO 2 high-burnup microstructure fragmentation

Here, in this paper, we present a phase-field model of quasi-brittle fracture with pressurized cracks, with dedicated applications for polycrystalline materials. The model is formulated as a minimization problem within the variational framework. The external work done by pressure on the crack surfaces is included in the objective function. Several careful modeling choices lead to a regularization-length-independent critical strength. The model is constructed to give a softening response with an underlying linear traction-separation law. The pressure-dependent softening response and the regularization of the prescribed pressure are demonstrated with a (quasi) one-dimensional numerical analysis. In a two-dimensional numerical analysis under plane strain assumptions, the critical stress corresponding to crack propagation (as predicted by our quasi-brittle fracture model) is compared with linear elastic fracture mechanics (LEFM) analytical solutions. Our model is further utilized to simulate fission-gas-induced fragmentation of the UO 2 high-burnup structure (HBS). Simulation results show that pressurized bubbles can cause crack nucleation and propagation, and that the bubble size and the surrounding external pressure affect the critical pressure corresponding to crack nucleation. Simulations of a partial HBS at different recrystallization stages show that different grain structures (due to recrystallization) also influence crack paths and fragmentation morphology.

36 MATERIALS SCIENCE↗

Ultrasensitive Sample Quantitation via Selected Reaction Monitoring Using CITP/CZE-ESI-Triple Quadrupole MS

We demonstrate the direct coupling of transient capillary isotachophoresis/ capillary zone electrophoresis (CITP/CZE) with a high sensitivity triple quadrupole mass spectrometer operating in selected reaction monitoring (SRM) mode for sample quantitation. The capability of CITP/CZE for in situ sample enrichment and separation has been shown to significantly improve the analytical figures of merit. A linear dynamic range spanning more than 4 orders of magnitude was observed. An average signal-to-noise ratio (S/N) of 49.6 was observed for 50 attomoles of targeted peptide in the presence of a complex and much more abundant bovine serum albumin (BSA) digest products. A correlation of variation (CV) less than 10 % for peak area was measured from triplicate sample analyses at 50 pM peptide concentration, showing good reproducibility of this online CITP/CZE-SRM mass spectrometry (MS) platform, and with limit of quantitation (LOQ) demonstrated to be well below 50 pM.

Environmental Molecular Sciences Laboratory↗

Active control of the resistive wall mode with power saturation

An analytic model of non-linear feedback stabilization of the resistive wall mode is presented. The non-linearity comes from either the current or the voltage saturation of the control coil power supply. For the so-called flux-to-current control, the current saturation of active coils always results in the loss of control. On the contrary, the flux-to-voltage control scheme tolerates certain degree of the voltage saturation. The minimal voltage limit is calculated, below which the control will be lost.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A current-driven resistive instability and its nonlinear effects in simulations of coaxial helicity injection in a tokamak

An instability observed in whole-device, resistive magnetohydrodynamic simulations of the driven phase of coaxial helicity injection in the National Spherical Torus eXperiment is identified as a current-driven resistive mode in an unusual geometry that transiently generates a current sheet. The mode consists of plasma flow velocity and magnetic field eddies in a tube aligned with the magnetic field at the surface of the injected magnetic flux. At low plasma temperatures (∼10–20 eV), the mode is benign, but at high temperatures (∼100 eV) its amplitude undergoes relaxation oscillations, broadening the layer of injected current and flow at the surface of the injected toroidal flux and background plasma. The poloidal-field structure is affected and the magnetic surface closure is generally prevented while the mode undergoes relaxation oscillations during injection. This study describes the mode and uses linearized numerical computations and an analytic slab model to identify the unstable mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A theoretical model of a thinning current sheet in the low-β plasmas

Magnetic reconnection is an important physical process in various explosive phenomena in the universe. In previous studies, it was found that fast reconnection takes place when the thickness of a current sheet becomes on the order of a microscopic length such as the ion Larmor radius or the ion inertial length. In this study, we investigated the pinching process of a current sheet by the Lorentz force in a low-β plasma using one-dimensional magnetohydrodynamics (MHD) simulations. It is known that there is an exact self-similar solution for this problem that neglects gas pressure. We compared the non-linear MHD dynamics with the analytic self-similar solution. From the MHD simulations, we found that with the gas pressure included the implosion process deviates from the analytic self-similar solution as t→t{sub 0}, where t{sub 0} is the explosion time when the thickness of a current sheet of the analytic solution becomes 0. We also found that a pair of MHD fast-mode shocks is generated and propagates after the formation of the pinched current sheet as t→t{sub 0}. On the basis of the Rankine–Hugoniot relations, we derived the scaling law of the physical quantities with respect to the initial plasma beta in the pinched current sheet. Our study could help us estimate the physical quantities in the pinched current sheet formed in a low-β plasma.

79 ASTRONOMY AND ASTROPHYSICS↗

Parametrization for the scale dependent growth in modified gravity

We propose a scale dependent analytic approximation to the exact linear growth of density perturbations in Scalar-Tensor (ST) cosmologies. In particular, we show that on large subhorizon scales, in the Newtonian gauge, the usual scale independent subhorizon growth equation does not describe the growth of perturbations accurately, as a result of scale-dependent relativistic corrections to the Poisson equation. A comparison with exact linear numerical analysis indicates that our approximation is a significant improvement over the standard subhorizon scale independent result on large subhorizon scales. A comparison with the corresponding results in the Synchronous gauge demonstrates the validity and consistency of our analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinematics and dynamics of disclination lines in three-dimensional nematics

An exact kinematic law for the motion of disclination lines in nematic liquid crystals as a function of the tensor order parameter Q is derived. Unlike other order parameter fields that become singular at their respective defect cores, the tensor order parameter remains regular. Following earlier experimental and theoretical work, the disclination core is defined to be the line where the uniaxial and biaxial order parameters are equal, or equivalently, where the two largest eigenvalues of Q cross. This allows an exact expression relating the velocity of the line to spatial and temporal derivatives of Q on the line, to be specified by a dynamical model for the evolution of the nematic. By introducing a linear core approximation for Q, analytical results are given for several prototypical configurations, including line interactions and motion, loop annihilation, and the response to external fields and shear flows. Behaviour that follows from topological constraints or defect geometry is highlighted. Finally, the analytic results are shown to be in agreement with three-dimensional numerical calculations based on a singular Maier–Saupe free energy that allows for anisotropic elasticity.

36 MATERIALS SCIENCE↗

Simulations of relativistic quantum plasmas using real-time lattice scalar QED

Real-time lattice quantum electrodynamics (QED) provides a unique tool for simulating plasmas in the strong-field regime, where collective plasma scales are not well separated from relativistic-quantum scales. As a toy model, we study scalar QED, which describes self-consistent interactions between charged bosons and electromagnetic fields. To solve this model on a computer, we first discretize the scalar-QED action on a lattice, in a way that respects geometric structures of exterior calculus and U(1)-gauge symmetry. The lattice scalar QED can then be solved, in the classical-statistics regime, by advancing an ensemble of statistically equivalent initial conditions in time, using classical field equations obtained by extremizing the discrete action. To demonstrate the capability of our numerical scheme, we apply it to two example problems. The first example is the propagation of linear waves, where we recover analytic wave dispersion relations using numerical spectrum. The second example is an intense laser interacting with a one-dimensional plasma slab, where we demonstrate natural transition from wakefield acceleration to pair production when the wave amplitude exceeds the Schwinger threshold. Our real-time lattice scheme is fully explicit and respects local conservation laws, making it reliable for long-time dynamics. The algorithm is readily parallelized using domain decomposition, and the ensemble may also be computed using quantum parallelism in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust Data-Driven Predictive Run-to-Run Control for Automated Serial Sectioning

This letter presents a one-step predictive run-to-run controller (R2R-MPC) for the automation of mechanical serial sectioning (MSS), a destructive material analysis process. To address the inherent uncertainty and disturbances in the MSS process, a robust closed-loop approach is presented. Here, the robust R2R-MPC models the uncertainty of the MSS process using a linear differential inclusion. As an analytical model of the MSS process is unavailable, the differential inclusion is identified from historical data. The R2R-MPC is posed as an optimization problem that computes incremental changes to the control input which minimize the worst-case material removal errors. This optimization-based controller is combined with a run-to-run controller to provide integral action that rejects constant disturbances and tracks constant reference removal rates. To demonstrate the efficacy of our robust R2R-MPC, we present simulation results which compare the presented controller with a conventional non-robust R2R.

42 ENGINEERING↗

Chance-Constrained AC Optimal Power Flow: Reformulations and Efficient Algorithms

Higher levels of renewable electricity generation increase uncertainty in power system operation. To ensure secure system operation, new tools that account for this uncertainty are required. Here, in this paper, we adopt a chance-constrained AC optimal power flow formulation, which guarantees that generation, power flows and voltages remain within their bounds with a pre-defined probability. We then discuss different chance-constraint reformulations and solution approaches for the problem. Additionally, we first discuss an analytical reformulation based on partial linearization, which enables us to obtain a tractable representation of the optimization problem. We then provide an efficient algorithm based on an iterative solution scheme which alternates between solving a deterministic AC OPF problem and assessing the impact of uncertainty. This more flexible computational framework enables not only scalable implementations, but also alternative chance-constraint reformulations. In particular, we suggest two sample based reformulations that do not require any approximation or relaxation of the AC power flow equations.

24 POWER TRANSMISSION AND DISTRIBUTION↗

SU-E-T-323: Dosimetric Evaluation of Small Fields for SBRT Treatment

Purpose: Stereotactic body radiation therapy (SBRT) is commonly employed to treat small targets for effective tumor control with radiation beams of small field sizes. The goal of this work was to evaluate dosimetrically a treatment planning system (TPS) by comparing the calculated dose for SBRT treatment with ion-chamber measurements. Methods: 3D images of a solid-water phantom with a pinpoint ion-chamber (0.015cm3) inside were acquired with a CT scanner. Active volume of the ion-chamber was delineated on CT images. Targets with a diameter of 1.5cm, 2cm, 3cm, 4cm and 5cm were drawn around the chamber. 3DCRT plans were generated for each target size with centrally opened 6MV beams and off-axis beams by changing the isocenter location, respectively, using a TPS with the Analytical Anisotropic Algorithm. A 21iX linear accelerator was employed for plan delivery. The measured and calculated doses were compared. To evaluate the dose calculations in heterogeneity for small fields SBRT treatment, similar plans were also generated and delivered on a heterogeneous thoracic phantom for 5 different size targets in the lung. Results: Dose comparisons between measurements and calculations showed 5.2%, 1.88%, 1.34%, 1.01% and 0.85% difference for SBRT plans with small central axis beams and 0.96%, 0.15%, 0.58%, 0.22% and 0.77% difference for plans with off-axis beams for five different size targets. For the thoracic phantom, the differences on dose between measurements and calculations are bigger, which are 8%, 5.9%, 4.5%, 3.9% and 4.5%, respectively. Conclusion: Dose verification for small fields used in the SBRT treatment has been performed based on ion-chamber measurements in both homogenous and heterogeneous phantoms. More than a 5% difference has been observed in the heterogeneous phantom, especially for very small fields. To meet the ICRU recommendation on a dose difference of no more than 5%, some corrections on the commissioning parameters of the TPS are needed.

60 APPLIED LIFE SCIENCES↗

Modifying the Asynchronous Jacobi Method for Data Corruption Resilience

Moving scientific computation from high-performance computing (HPC) and cloud computing (CC) environments to devices on the edge, i.e., physically near instruments of interest, has received tremendous interest in recent years. Such edge computing environments can operate on data in situ, offering enticing benefits over data aggregation to HPC and CC facilities that include avoiding costs of transmission, increased data privacy, and real-time data analysis. Because of the inherent unreliability of edge computing environments, new fault-tolerant approaches must be developed before the benefits of edge computing can be realized. Motivated by algorithm-based fault tolerance, a variant of the asynchronous Jacobi (ASJ) method is developed that achieves resilience to data corruption by rejecting solution approximations from neighbor devices according to a bound derived from convergence theory. Numerical results on a two-dimensional Poisson problem show that the new rejection criterion, along with a novel approximation to the shortest path length on which the criterion depends, restores convergence for the ASJ variant in the presence of certain types data corruption. Numerical results are obtained for when the singular values in the analytic bound are approximated. Additional linear systems are also explored, one with a more dense sparsity pattern and one that includes advection. All results indicate that successful resilience to data corruption depends on whether the bound tightens fast enough to reject corrupted data before the iteration evolution deviates significantly from that predicted by the convergence theory defining the bound. This observation generalizes to future work on algorithm-based fault tolerance for other asynchronous algorithms, including upcoming approaches that leverage Krylov subspaces.

97 MATHEMATICS AND COMPUTING↗

Efficient simulation and model reformulation of two-dimensional electrochemical thermal behavior of lithium-ion batteries

Lithium-ion batteries are an important technology to facilitate efficient energy storage and enable a shift from petroleum based energy to more environmentally benign sources. Such systems can be utilized most efficiently if good understanding of performance can be achieved for a range of operating conditions. Mathematical models can be useful to predict battery behavior to allow for optimization of design and control. An analytical solution is ideally preferred to solve the equations of a mathematical model, as it eliminates the error that arises when using numerical techniques and is usually computationally cheap. An analytical solution provides insight into the behavior of the system and also explicitly shows the effects of different parameters on the behavior. However, most engineering models, including the majority of battery models, cannot be solved analytically due to non-linearities in the equations and state dependent transport and kinetic parameters. The numerical method used to solve the system of equations describing a battery operation can have a significant impact on the computational cost of the simulation. In this paper, a model reformulation of the porous electrode pseudo three dimensional (P3D) which significantly reduces the computational cost of lithium ion battery simulation, while maintaining high accuracy, is discussed. This reformulation enables the use of the P3D model into applications that would otherwise be too computationally expensive to justify its use, such as online control, optimization, and parameter estimation. Furthermore, the P3D model has proven to be robust enough to allow for the inclusion of additional physical phenomena as understanding improves. In this study, the reformulated model is used to allow for more complicated physical phenomena to be considered for study, including thermal effects.

25 ENERGY STORAGE↗