Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

INVESTIGATION OF HELICITY AND ENERGY FLUX TRANSPORT IN THREE EMERGING SOLAR ACTIVE REGIONS

We report the results of an investigation of helicity and energy flux transport from three emerging solar active regions (ARs). Using time sequence vector magnetic field observations obtained from the Helioseismic Magnetic Imager, the velocity field of plasma flows is derived by the differential affine velocity estimator for vector magnetograms. In three cases, the magnetic fluxes evolve to pump net positive, negative, and mixed-sign helicity flux into the corona. The coronal helicity flux is dominantly coming from the shear term that is related to horizontal flux motions, whereas energy flux is dominantly contributed by the emergence term. The shear helicity flux has a phase delay of 5–14 hr with respect to absolute magnetic flux. The nonlinear curve of coronal energy versus relative helicity identifies the configuration of coronal magnetic fields, which is approximated by a fit of linear force-free fields. The nature of coronal helicity related to the particular pattern of evolving magnetic fluxes at the photosphere has implications for the generation mechanism of two kinds of observed activity in the ARs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Non-Markovian dynamics of a superconducting qubit in an open multimode resonator

In this work, we study the dynamics of a transmon qubit that is capacitively coupled to an open multimode superconducting resonator. Our effective equations are derived by eliminating resonator degrees of freedom while encoding their effect in the Green's function of the electromagnetic background. We account for the dissipation of the resonator exactly by employing a spectral representation for the Green's function in terms of a set of non-Hermitian modes and show that it is possible to derive effective Heisenberg-Langevin equations without resorting to the rotating-wave, two-level, Born, or Markov approximations. A well-behaved time-domain perturbation theory is derived to systematically account for the nonlinearity of the transmon. We apply this method to the problem of spontaneous emission, capturing accurately the non-Markovian features of the qubit dynamics, valid for any qubit-resonator coupling strength.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2016

Modified sparse regularization for electrical impedance tomography

Electrical impedance tomography (EIT) aims to estimate the electrical properties at the interior of an object from current-voltage measurements on its boundary. It has been widely investigated due to its advantages of low cost, non-radiation, non-invasiveness, and high speed. Image reconstruction of EIT is a nonlinear and ill-posed inverse problem. Therefore, regularization techniques like Tikhonov regularization are used to solve the inverse problem. A sparse regularization based on L{sub 1} norm exhibits superiority in preserving boundary information at sharp changes or discontinuous areas in the image. However, the limitation of sparse regularization lies in the time consumption for solving the problem. In order to further improve the calculation speed of sparse regularization, a modified method based on separable approximation algorithm is proposed by using adaptive step-size and preconditioning technique. Both simulation and experimental results show the effectiveness of the proposed method in improving the image quality and real-time performance in the presence of different noise intensities and conductivity contrasts.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Theoretical study of various nonlinear phenomena in plasma systems and scaling of magneto-inertial-fusion targets

Plasma physics is an exciting field of study with a wide variety of nonlinear processes that come into play. Examples of such processes include the interaction of small-scale turbulence with large-scale plasma structures and the nonlinear saturation of plasma instabilities, for example those of magneto-hydrodynamical nature. During this Truman LDRD project, I studied a collection of nonlinear problems that are of interest to the field of plasma physics. This LDRD report summarizes four main research accomplishments. First, a new statistical model for describing inhomogeneous drift-wave turbulence inter- acting with zonal flows was developed. This new model includes the effects of nonlinear wave-wave collisions, which are expected to change the spectrum of the underlying DW turbulence and therefore the generation of zonal flows. Second, a new mathematical formalism was proposed to systematically apply the non- linear WKB approximation to general field theories, including those often used in fluid dynamics. This formalism represents an interesting tool for studying physical systems that show an explicit scale separation. Third, a weakly nonlinear model was developed to describe the magneto-Rayleigh-Taylor instability. This instability is of paramount importance to understand as it can reduce the performance of magnetic-inertial-fusion (MIF) platforms. The developed models captures the effects of harmonic generation and saturation of the linear growth of the instability. Finally, a framework was proposed for scaling magneto-inertial fusion (MIF) targets to larger pulsed-power drivers. From this framework, a set of scaling rules were derived that conserve the physical regimes of MIF systems when scaling up in peak current. By doing so, deleterious nonlinear processes that affect MIF performance may be kept at bay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Implicit and coupled fluid plasma solver with adaptive Cartesian mesh and its applications to non-equilibrium gas discharges

In this work, we present a new fluid plasma solver with adaptive Cartesian mesh (ACM) based on a full-Newton (nonlinear, implicit) scheme for non-equilibrium gas discharge plasma. The electrons and ions are described using drift-diffusion approximation coupled to Poisson equation for the electric field. The electron-energy transport equation is solved to account for electron thermal conductivity, Joule heating, and energy loss of electrons in collisions with neutral species. The rate of electron-induced ionization is a function of electron temperature and could also depend on electron density (important for plasma stratification). The ion and gas temperature are kept constant. The transport equations are discretized using a non-isothermal Scharfetter-Gummel scheme to resolve possible large temperature gradients in the sheaths. We demonstrate the new solver for simulations of direct current (DC) and radiofrequency (RF) discharges. The implicit treatment of the coupled equations allows using large time steps. The full-Newton method (FNM) enables fast nonlinear convergence at each time step, offering significantly improved simulation efficiency. We discuss the selection of time steps for solving different plasma problems. The new solver enables solving several problems we could not solve before with existing software: two- and three-dimensional structures of the entire DC discharges including cathode and anode regions, electric field reversals and double-layer formation, the normal cathode spot and an anode ring, moving striations in diffuse and constricted DC discharges, and standing striations in RF discharges. The developed FNM-ACM technique offers many benefits for tackling the disparity of gas discharge plasma systems' time scales and nonlinearity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Landau ghost pole problem in quantum field theory: From 50th of last century to the present day

In this paper we present our results of the investigation of asymptotical behavior of amplitude at short distances in four-dimensional scalar field theory with ϕ{sup 4} interaction. To formulate of our calculating model – two-particle approximation of the mean-field expansion we have used an Rochev’s iteration scheme of solution of the Schwinger-Dyson equations with the fermion bilocal source. We have considered the nonlinear integral equations in deep-inelastic region of momenta. As result we have a non-trivial behavior of amplitude at large momenta.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Surrogate-Based Asynchronous Decomposition Technique for Realistic Security-Constrained Optimal Power Flow Problems

Here we present a decomposition approach for obtaining good feasible solutions for the security-constrained, alternating-current, optimal power flow (SC-AC-OPF) problem at an industrial scale and under real-world time and computational limits. The approach was designed while preparing and participating in ARPA-E’s Grid Optimization Competition (GOC) Challenge 1. The challenge focused on a near-real-time version of the SC-AC-OPF problem, where a base operating point is optimized, taking into account possible single-element contingencies, after which the system adapts its operating point following the response of automatic frequency droop controllers and voltage regulators. Our solution approach for this problem relies on state-of-the-art nonlinear programming algorithms, and it employs nonconvex relaxations for complementarity constraints, a specialized two-stage decomposition technique with sparse approximations of recourse terms and contingency ranking and prescreening. The paper describes and justifies our approach and outlines the features of its implementation, including functions and derivatives evaluation, warm-starting strategies, and asynchronous parallelism. We discuss the results of the independent benchmark of our approach by ARPA-E’s GOC team in Challenge 1, where it was found to consistently produce high-quality solutions across a wide range of network sizes and difficulty, and conclude by outlining future extensions of the approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Formation of plasmon pulses in the cooperative decay of excitons of quantum dots near a metal surface

The formation of pulses of surface electromagnetic waves at a metal–dielectric boundary is considered in the process of cooperative decay of excitons of quantum dots distributed near a metal surface in a dielectric layer. It is shown that the efficiency of exciton energy transfer to excited plasmons can, in principle, be increased by selecting the dielectric material with specified values of the complex permittivity. It is found that in the mean field approximation, the semiclassical model of formation of plasmon pulses in the system under study is reduced to the pendulum equation with the additional term of nonlinear losses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Parametric decay of an extraordinary electromagnetic wave in relativistic plasma

Parametric instability of an extraordinary electromagnetic wave in plasma preheated to a relativistic temperature is considered. A set of self-similar nonlinear differential equations taking into account the electron “thermal” mass is derived and investigated. Small perturbations of the parameters of the heated plasma are analyzed in the linear approximation by using the dispersion relation determining the phase velocities of the fast and slow extraordinary waves. In contrast to cold plasma, the evanescence zone in the frequency range above the electron upper hybrid frequency vanishes and the asymptotes of both branches converge. Theoretical analysis of the set of nonlinear equations shows that the growth rate of decay instability increases with increasing initial temperature of plasma electrons. This result is qualitatively confirmed by numerical simulations of plasma heating by a laser pulse injected from vacuum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Alkali-metal/alkaline-earth-metal fluorine beryllium borate NaSr{sub 3}Be{sub 3}B{sub 3}O{sub 9}F{sub 4} with large nonlinear optical properties in the deep-ultraviolet region

The linear optical response and second harmonic generation (SHG) in alkali-metal/alkaline-earth-metal fluorine beryllium borate NaSr{sub 3}Be{sub 3}B{sub 3}O{sub 9}F{sub 4} are investigated by means of density functional theory. Calculations are performed using four types of exchange correlations: Ceperley-Alder local density approximation, Perdew Burke and Ernzerhof general gradient approximation, Engel-Vosko generalized gradient approximation, and the recently modified Becke-Johnson potential (mBJ). The mBJ approach brings the calculated band gap (7.20 eV) in excellent agreement with the experimental one (7.28 eV). The calculated values of the uniaxial anisotropy δε=−0.076 and the birefringence Δn(0)=0.052 indicate considerable anisotropy in the linear optical properties, which makes it favorable for the second harmonic generation. The dominant component of the second harmonic generation is χ{sub 111}{sup (2)}(ω). The value of |χ{sub 111}{sup (2)}(ω)| is about 1.2 pm/V at λ = 1064 nm in agreement with previous calculations. To analyze the origin of the high SHG of NaSr{sub 3}Be{sub 3}B{sub 3}O{sub 9}F{sub 4} single crystals, we have correlated the features of |χ{sub 111}{sup (2)}(ω)| spectra with the features of ε{sub 2}(ω) spectra as a function of ω/2 and ω. From the calculated dominant component |χ{sub 111}{sup (2)}(ω)|, we find that the microscopic first hyperpolarizability, β{sub 111}, the vector components along the dipole moment direction is 0.5 × 10{sup −30} esu at static limit and 0.6 × 10{sup −30} esu at λ = 1064 nm.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Optomechanical interaction in fibre lasers with micro-optomechanical resonance structures

We report the results of studies of fibre lasers with passive modulators based on light-excited micro-optomechanical resonance structures (micro-oscillators). It is shown that in fibre lasers based on active fibres doped with rare-earth elements (Er, Er/Yb, Yb, Nd), the optomechanical interaction of laser radiation with micro-oscillators of various types (fibre-optic, microvolume) leads to self-oscillations of the characteristics of laser radiation at frequencies of relaxation oscillations and intermode beats synchronised with the frequencies of elastic eigenoscillations of micro-oscillators. It is found that in an ultra-long erbium – ytterbium fibre laser with a nonlinear mirror based on a microcantilever, laser photothermal excitation of the second mode of elastic transverse oscillations of the microcantilever makes it possible to perform passive mode locking exclusively due to Q-switching of the laser cavity. Pulsed lasing with a controlled repetition rate (∼76 kHz), a pulse duration of 2 – 5 μs and an output energy of 0.1 μJ per pulse is implemented. Based on a simplified physical model of the indicated fibre lasers with micro-oscillators, we have developed an approximate mathematical model describing the regimes of passive mode locking of fibre lasers with micro-oscillators that play the role of mirrors with a nonlinear reflection coefficient in the laser cavity. The prospects for the development and application of the considered laser systems are discussed. (control of laser radiation parameters)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Approximate method for stochastic chemical kinetics with two-time scales by chemical Langevin equations

The frequently used reduction technique is based on the chemical master equation for stochastic chemical kinetics with two-time scales, which yields the modified stochastic simulation algorithm (SSA). For the chemical reaction processes involving a large number of molecular species and reactions, the collection of slow reactions may still include a large number of molecular species and reactions. Consequently, the SSA is still computationally expensive. Because the chemical Langevin equations (CLEs) can effectively work for a large number of molecular species and reactions, this paper develops a reduction method based on the CLE by the stochastic averaging principle developed in the work of Khasminskii and Yin [SIAM J. Appl. Math. 56, 1766–1793 (1996); ibid. 56, 1794–1819 (1996)] to average out the fast-reacting variables. This reduction method leads to a limit averaging system, which is an approximation of the slow reactions. Because in the stochastic chemical kinetics, the CLE is seen as the approximation of the SSA, the limit averaging system can be treated as the approximation of the slow reactions. As an application, we examine the reduction of computation complexity for the gene regulatory networks with two-time scales driven by intrinsic noise. For linear and nonlinear protein production functions, the simulations show that the sample average (expectation) of the limit averaging system is close to that of the slow-reaction process based on the SSA. It demonstrates that the limit averaging system is an efficient approximation of the slow-reaction process in the sense of the weak convergence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2016

Near-infrared strong coupling between metamaterials and epsilon-near-zero modes in degenerately doped semiconductor nanolayers

Epsilon-near-zero (ENZ) modes provide a new path for tailoring light–matter interactions at the nanoscale. In this paper, we analyze a strongly coupled system at near-infrared frequencies comprising plasmonic metamaterial resonators and ENZ modes supported by degenerately doped semiconductor nanolayers. In strongly coupled systems that combine optical cavities and intersubband transitions, the polariton splitting (i.e., the ratio of Rabi frequency to bare cavity frequency) scales with the square root of the wavelength, thus favoring the long-wavelength regime. In contrast, we observe that the polariton splitting in ENZ/metamaterial resonator systems increases linearly with the thickness of the nanolayer supporting the ENZ modes. In this work, we employ an indium-tin-oxide nanolayer and observe a large experimental polariton splitting of approximately 30% in the near-infrared. As a result, this approach opens up many promising applications, including nonlinear optical components and tunable optical filters based on controlling the polariton splitting by adjusting the frequency of the ENZ mode.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Soft mean spherical approximation for dusty plasma liquids: Level of accuracy and analytic expressions

The soft mean spherical approximation is employed for the study of the thermodynamics of dusty plasma liquids, the latter treated as Yukawa one-component plasmas. Within this integral theory method, the only input necessary for the calculation of the reduced excess energy stems from the solution of a single non-linear algebraic equation. Consequently, thermodynamic quantities can be routinely computed without the need to determine the pair correlation function or the structure factor. The level of accuracy of the approach is quantified after an extensive comparison with numerical simulation results. The approach is solved over a million times with input spanning the whole parameter space and reliable analytic expressions are obtained for the basic thermodynamic quantities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Customizable two-species kinetic equilibria for nonuniform low-beta plasmas

Here, two methods are developed for constructing self-consistent two-species kinetic equilibria for low-beta electrostatic plasmas, in which the magnetic field is uniform and fixed in time. The first method uses Taylor series approximations to construct distribution functions that can be specified analytically. The second method relies on numerically solving a nonlinear ordinary differential equation and produces exact—to numerical precision—equilibria. In both methods, the equilibrium distribution functions for ions and electrons are expressed in terms of constants of motion and satisfy the steady-state Vlasov-Poisson equation system. Provided that the ion drift speed does not exceed the ion thermal speed, the equilibria can be specified with customizable density and electrostatic potential profiles. The methods can thereby be tailored to different applications and are successfully applied to construct kinetic equilibria for cross-field plasmas with sheared flows, large density variations, and different levels of magnetization. The equilibria are used to initialize fourth-order finite-volume Vlasov-Poisson simulations in (x, v x , v y ) coordinates and the associated temporal evolution is used to assess the accuracy of each method. Lastly, the low-amplitude deviations observed in these simulations demonstrate that the kinetic equilibria are robust and that they provide a valuable means of studying the dynamics of nonuniform magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Efficient calculation of cosmological neutrino clustering in the non-linear regime

We study in detail how neutrino perturbations can be followed in linear theory by using only terms up to l =2 in the Boltzmann hierarchy. We provide a new approximation to the third moment and demonstrate that the neutrino power spectrum can be calculated to a precision of better than ∼ 5% for masses up to ∼ 1 eV and k ∼< 10 h /Mpc. The matter power spectrum can be calculated far more precisely and typically at least a factor of a few better than with existing approximations. We then proceed to study how the neutrino power spectrum can be reliably calculated even in the non-linear regime by using the non-linear gravitational potential, sourced by dark matter overdensities, as it is derived from semi-analytic methods based on N -body simulations in the Boltzmann evolution hierarchy. Our results agree extremely well with results derived from N -body simulations that include cold dark matter and neutrinos as independent particles with different properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗