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 199 records · Page 11

Optimization of transformer ratio and beam loading in a plasma wakefield accelerator with a structure-exploiting algorithm

Plasma-based acceleration has emerged as a promising candidate as an accelerator technology for a future linear collider or a next-generation light source. We consider the plasma wakefield accelerator (PWFA) concept where a plasma wave wake is excited by a particle beam and a trailing beam surfs on the wake. For a linear collider, the energy transfer efficiency from the drive beam to the wake and from the wake to the trailing beam must be large, while the emittance and energy spread of the trailing bunch must be preserved. One way to simultaneously achieve this when accelerating electrons is to use longitudinally shaped bunches and nonlinear wakes. In the linear regime, there is an analytical formalism to obtain the optimal shapes. In the nonlinear regime, however, the optimal shape of the driver to maximize the energy transfer efficiency cannot be precisely obtained because currently no theory describes the wake structure and excitation process for all degrees of nonlinearity. In addition, the ion channel radius is not well defined at the front of the wake where the plasma electrons are not fully blown out by the drive beam. We present results using a novel optimization method to effectively determine a current profile for the drive and trailing beam in PWFA that provides low energy spread, low emittance, and high acceleration efficiency. We parameterize the longitudinal beam current profile as a piecewise-linear function and define optimization objectives. For the trailing beam, the algorithm converges quickly to a nearly inverse trapezoidal trailing beam current profile similar to that predicted by the ultrarelativistic limit of the nonlinear wakefield theory. For the drive beam, the beam profile found by the optimization in the nonlinear regime that maximizes the transformer ratio also resembles that predicted by linear theory. Furthermore, the current profiles found from the optimization method provide higher transformer ratios compared with the linear ramp predicted by the relativistic limit of the nonlinear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Single-mode perturbation growth in an idealized spherical implosion

This paper considers a dense imploding spherical shell, where perturbations on the inner surface grow due to hydrodynamic instabilities, geometric convergence and compression. A low-convergence implosion with a single spherical harmonic mode perturbation with mode numbers in the range from $\ell$ = 5 to $\ell$ = 100 and at three different amplitudes is considered. The linear theory of Epstein [1] is extended to incorporate reshock, making it applicable to the full implosion process (while perturbations are linear). Linear theory is employed to choose modes such that quantified geometric convergence and compression effects contribute significantly to perturbation growth at the lower mode number, while at the higher mode number the contribution from Rayleigh–Taylor instability dominates. Simulation results from four independent simulation codes (FLASH, HYDRA, Miranda and Flamenco) are presented. The simulation predictions are validated against linear theory pre-reshock and employed to validate the extended theory across the reshock. The simulations continue to substantially non-linear perturbation amplitudes, beyond the limits of the analytical approach, and the presented perturbation amplitudes can inform future non-linear modeling. In this work, the simulated perturbation amplitudes agree to within approximately 10% for most cases, with isolated cases having differences of greater than 50% from the simulated ensemble-mean. This occurs in the case of a low mode number perturbation with a very small amplitude, where the growth of secondary small-scale instabilities leads to substantial differences between the codes after re-shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗