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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Linear fractional differential equations and eigenfunctions of fractional differential operators

Eigenfunctions associated with Riemann–Liouville and Caputo fractional differential operators are obtained by imposing a restriction on the fractional derivative parameter. Those eigenfunctions can be used to express the analytical solution of some linear sequential fractional differential equations. As a first application, we discuss analytical solutions for the so-called fractional Helmholtz equation with one variable, obtained from the standard equation in one dimension by replacing the integer order derivative by the Riemann–Liouville fractional derivative. A second application consists of an initial value problem for a fractional wave equation in two dimensions in which the integer order partial derivative with respect to the time variable is replaced by the Caputo fractional derivative. The classical Mittag-Leffler functions are important in the theory of fractional calculus because they emerge as solutions of fractional differential equations. Starting with the solution of a specific fractional differential equation in terms of these functions, we find a way to express the exponential function in terms of classical Mittag-Leffler functions. A remarkable characteristic of this relation is that it is true for any value of the parameter n appearing in the definition of the functions, i.e., we have an infinite family of different expressions for e{sup x} in terms of classical Mittag-Leffler functions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Discovering charge density functionals and structure-property relationships with PROPhet: A general framework for coupling machine learning and first-principles methods

Modern ab initio methods have rapidly increased our understanding of solid state materials properties, chemical reactions, and the quantum interactions between atoms. However, poor scaling often renders direct ab initio calculations intractable for large or complex systems. There are two obvious avenues through which to remedy this problem: (i) develop new, less expensive methods to calculate system properties, or (ii) make existing methods faster. This paper describes an open source framework designed to pursue both of these avenues. PROPhet (short for PROPerty Prophet) utilizes machine learning techniques to find complex, non-linear mappings between sets of material or system properties. The result is a single code capable of learning analytical potentials, non-linear density functionals, and other structure-property or property-property relationships. These capabilities enable highly accurate mesoscopic simulations, facilitate computation of expensive properties, and enable the development of predictive models for systematic materials design and optimization. Here, this work explores the coupling of machine learning to ab initio methods through means both familiar (e.g., the creation of various potentials and energy functionals) and less familiar (e.g., the creation of density functionals for arbitrary properties), serving both to demonstrate PROPhet’s ability to create exciting post-processing analysis tools and to open the door to improving ab initio methods themselves with these powerful machine learning techniques.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2011

Over-reflection of slow magnetosonic waves by homogeneous shear flow: Analytical solution

We have analyzed the amplification of slow magnetosonic (or pseudo-Alfvenic) waves (SMW) in incompressible shear flow. As found here, the amplification depends on the component of the wave-vector perpendicular to the direction of the shear flow. Earlier numerical results are consistent with the general analytic solution for the linearized magnetohydrodynamic equations, derived here for the model case of pure homogeneous shear (without Coriolis force). An asymptotically exact analytical formula for the amplification coefficient is derived for the case when the amplification is sufficiently large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

The frequency-dependent response of single aerosol particles to vapour phase oscillations and its application in measuring diffusion coefficients

A new method for measuring diffusion in the condensed phase of single aerosol particles is proposed and demonstrated. The technique is based on the frequency-dependent response of a binary particle to oscillations in the vapour phase of one of its chemical components. Here, we discuss how this physical situation allows for what would typically be a non-linear boundary value problem to be approximately reduced to a linear boundary value problem. For the case of aqueous aerosol particles, we investigate the accuracy of the closed-form analytical solution to this linear problem through a comparison with the numerical solution of the full problem. Then, using experimentally measured whispering gallery modes to track the frequency-dependent response of aqueous particles to relative humidity oscillations, we determine diffusion coefficients as a function of water activity. The measured diffusion coefficients are compared to previously reported values found using the two common experiments: (i) the analysis of the sorption/desorption of water from a particle after a step-wise change to the surrounding relative humidity and (ii) the isotopic exchange of water between a particle and the vapour phase. The technique presented here has two main strengths: first, when compared to the sorption/desorption experiment, it does not require the numerical evaluation of a boundary value problem during the fitting process as a closed-form expression is available. Second, when compared to the isotope exchange experiment, it does not require the use of labeled molecules. Therefore, the frequency-dependent experiment retains the advantages of these two commonly used methods but does not suffer from their drawbacks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2012

Low-frequency linear-mode regimes in the tokamak scrape-off layer

Motivated by the wide range of physical parameters characterizing the scrape-off layer (SOL) of existing tokamaks, the regimes of low-frequency linear instabilities in the SOL are identified by numerical and analytical calculations based on the linear, drift-reduced Braginskii equations, with cold ions. The focus is put on ballooning modes and drift wave instabilities, i.e., their resistive, inertial, and ideal branches. A systematic study of each instability is performed, and the parameter space region where they dominate is identified. It is found that the drift waves dominate at high R/L{sub n}, while the ballooning modes at low R/L{sub n}; the relative influence of resistive and inertial effects is discussed. Electromagnetic effects suppress the drift waves and, when the threshold for ideal stability is overcome, the ideal ballooning mode develops. Our analysis is a first stage tool for the understanding of turbulence in the tokamak SOL, necessary to interpret the results of non-linear simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Analytic model of electron transport through and over non-linear barriers

Tunneling barriers are an essential component of electron sources, sensors, detectors, and vacuum nanoelectronics, and a pivotal factor in their performance, but the barriers themselves routinely depart from the analytic models used to model their behavior. Here, a new formalism is developed to analytically and accurately model emission through and over barriers associated with depletion layers, nanotip barriers, and MIM structures. The transmission probability for depletion layer and metal-insulator-metal (MIM) and metaloxide-semiconductor (MOS) barriers is accurately modeled as the electron energy exceeds the barrier height using approaches designed for rapid implementation demanded by simulation codes, and extensible to general barriers. The models supersede conventional thermal and field models in depletion and MIM/MOS barrier studies. Thermal-field methods are used to treat the transmission probability, and shape factor methods to treat the tunneling factor. Analytic formulae for current density are obtained. The methods ease device simulation and characterization of current-voltage relations for emerging technologically interesting barriers with better accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2012

Stabilization of the Rayleigh-Taylor instability in quantum magnetized plasmas

In this research, stabilization of the Rayleigh-Taylor instability (RTI) due to density gradients, magnetic fields, and quantum effects, in an ideal incompressible plasma, is studied analytically and numerically. A second-order ordinary differential equation (ODE) for the RTI including quantum corrections, with a continuous density profile, in a uniform external magnetic field, is obtained. Analytic expressions of the linear growth rate of the RTI, considering modifications of density gradients, magnetic fields, and quantum effects, are presented. Numerical approaches are performed to solve the second-order ODE. The analytical model proposed here agrees with the numerical calculation. It is found that the density gradients, the magnetic fields, and the quantum effects, respectively, have a stabilizing effect on the RTI (reduce the linear growth of the RTI). The RTI can be completely quenched by the magnetic field stabilization and/or the quantum effect stabilization in proper circumstances leading to a cutoff wavelength. The quantum effect stabilization plays a central role in systems with large Atwood number and small normalized density gradient scale length. The presence of external transverse magnetic fields beside the quantum effects will bring about more stability on the RTI. The stabilization of the linear growth of the RTI, for parameters closely related to inertial confinement fusion and white dwarfs, is discussed. Results could potentially be valuable for the RTI treatment to analyze the mixing in supernovas and other RTI-driven objects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Simple Expressions for the Design of Linear Tapers in Overmoded Corrugated Waveguides

In this paper, simple analytical formulae are presented for the design of linear tapers with very low mode conversion loss in overmoded corrugated waveguides. For tapers from waveguide radius a2 to a1, with a1 1 a 2 /λ. Here, λ is the wavelength of radiation. The fractional loss of the HE 11 mode in an optimized taper is 0.0293(a 2 -a 1 ) 4 /a$2\atop{1}$1a$2\atop{2}$. These formulae are accurate when a2≲2a 1 . Slightly more complex formulae, accurate for a 2 ≤4a 1 , are also presented in this paper. The loss in an overmoded corrugated linear taper is less than 1 % when a 2 ≤2.12a 1 and less than 0.1 % when a 2 ≤1.53a 1 . The present analytic results have been benchmarked against a rigorous mode matching code and have been found to be very accurate. The results for linear tapers are compared with the analogous expressions for parabolic tapers. Finally, parabolic tapers may provide lower loss, but linear tapers with moderate values of a 2 /a 1 may be attractive because of their simplicity of fabrication.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Investigating the Performance of a New Type of Preloaded Linear Stroke Length Magnetic Spring

This paper presents the analytic analysis, and proof-of-principle prototyping of a new type of magnetic spring with preload and a linear stroke length. An analytic based magnetic charge modeling approach is utilized to investigate the magnetic spring's energy density, stiffness characteristics and linearity. Here, it is shown that whilst the proposed magnetic spring has a lower mass and energy density than a mechanical spring, the magnetic spring offers several unique characteristics, such as contact-free operation, inherent preload as well as over-force failure protection. In addition, the operating principle of the presented magnetic spring can be extended to realize both positive and negative variable stiffness adjustment characteristics.

16 TIDAL AND WAVE POWER↗

DOE OSTI2019

Calibration of linear contact stiffnesses in discrete element models using a hybrid analytical-computational framework

Efficient selections of particle-scale contact parameters in discrete element modelling remain an open question. The purpose of this study is to provide a hybrid calibration framework to estimate linear contact stiffnesses (normal and tangential) for both two-dimensional and three-dimensional simulations. Analytical formulas linking macroscopic parameters (Young's modulus, Poisson's ratio) to mesoscopic particle parameters for granular systems are derived based on statistically isotropic packings under small-strain isotropic stress conditions. By taking the derived analytical solutions as initial approximations, the gradient descent algorithm automatically obtains a reliable numerical estimation. The introduced framework is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. The results show that this hybrid method practically reduces the time for artificial trials and errors to obtain reasonable stiffness parameters. The proposed framework can also be extended to other parameter calibration problems in DEM.

Calibration method↗

DOE OSTI2012

TWO-STAGE FRAGMENTATION FOR CLUSTER FORMATION: ANALYTICAL MODEL AND OBSERVATIONAL CONSIDERATIONS

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds shows that molecular clouds exhibit a preferred length scale for collapse that depends on the mass-to-flux ratio and neutral-ion collision time within the cloud. We extend this linear analysis to the context of clustered star formation. By combining the results of the linear analysis with a realistic ionization profile for the cloud, we find that a molecular cloud may evolve through two fragmentation events in the evolution toward the formation of stars. Our model suggests that the initial fragmentation into clumps occurs for a transcritical cloud on parsec scales while the second fragmentation can occur for transcritical and supercritical cores on subparsec scales. Comparison of our results with several star-forming regions (Perseus, Taurus, Pipe Nebula) shows support for a two-stage fragmentation model.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

Critical node treatment in the analytic function expansion method for Pin Power Reconstruction

Pin Power Reconstruction (PPR) was implemented in PARCS using the eight term analytic function expansion method (AFEN). This method has been demonstrated to be both accurate and efficient. However, similar to all the methods involving analytic functions, such as the analytic node method (ANM) and AFEN for nodal solution, the use of AFEN for PPR also has potential numerical issue with critical nodes. The conventional analytic functions are trigonometric or hyperbolic sine or cosine functions with an angular frequency proportional to buckling. For a critic al node the buckling is zero and the sine functions becomes zero, and the cosine function become unity. In this case, the eight terms of the analytic functions are no longer distinguishable from ea ch other which makes their corresponding coefficients can no longer be determined uniquely. The mode flux distribution of critical node can be linear while the conventional analytic functions can only express a uniform distribution. If there is critical or near critical node in a plane, the reconstructed pin power distribution is often be shown negative or very large values using the conventional method. In this paper, we propose a new method to avoid the numerical problem wit h critical nodes which uses modified trigonometric or hyperbolic sine functions which are the ratio of trigonometric or hyperbolic sine and its angular frequency. If there are no critical or near critical nodes present, the new pin power reconstruction method with modified analytic functions are equivalent to the conventional analytic functions. The new method is demonstrated using the L336C5 benchmark problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

Irradiation instability at the inner edges of accretion disks

An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around active galactic nuclei. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also present the linear growth rate as a function of β, the ratio between radiation force and gravity, and c{sub s}, the sound speed of the disk, obtained using two methods: a semi-analytic analysis of the linearized equations and a numerical simulation using the GPU-accelerated hydrodynamical code PEnGUIn. In particular, we find that IRI occurs at β ∼ 0.1 if the transition region extends as wide as ∼0.05r, and at higher β values if it is wider. This threshold value applies to c{sub s} ranging from 3% of the Keplerian orbital speed to 5%, and becomes higher if c{sub s} is lower. Furthermore, in the nonlinear evolution of the instability, disks with a large β and small c{sub s} exhibit 'clumping', extreme local surface density enhancements that can reach over 10 times the initial disk surface density.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

An analytical method as a preconditioning modeling for systems of linear equations

In this paper, we propose a simple and practical preconditioner based on the homotopy analysis method. Furthermore, using this modeling, we design an efficient algorithm to solve linear systems. Within the computational framework, convection–diffusion equation candidates are investigated and evaluated. Finally, computational comparisons of the standard methods against preconditioned methods based on experiments are presented which illustrate the merits of the simplicity, power, and effectiveness of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2012

Electron cyclotron resonance near the axis of the gas-dynamic trap

Propagation of an extraordinary electromagnetic wave in the vicinity of electron cyclotron resonance surface in an open linear trap is studied analytically, taking into account inhomogeneity of the magnetic field in paraxial approximation. Ray trajectories are derived from a reduced dispersion equation that makes it possible to avoid the difficulty associated with a transition from large propagation angles to the case of strictly longitudinal propagation. Our approach is based on the theory, originally developed by Zvonkov and Timofeev [Sov. J. Plasma Phys. 14, 743 (1988)], who used the paraxial approximation for the magnetic field strength, but did not consider the slope of the magnetic field lines, which led to considerable error, as has been recently noted by Gospodchikov and Smolyakova [Plasma Phys. Rep. 37, 768-774 (2011)]. We have found ray trajectories in analytic form and demonstrated that the inhomogeneity of both the magnetic field strength and the field direction can qualitatively change the picture of wave propagation and significantly affect the efficiency of electron cyclotron heating of a plasma in a linear magnetic trap. Analysis of the ray trajectories has revealed a criterion for the resonance point on the axis of the trap to be an attractor for the ray trajectories. It is also shown that a family of ray trajectories can still reach the resonance point on the axis if the latter generally repels the ray trajectories. As an example, results of general theory are applied to the electron cyclotron resonance heating experiment which is under preparation on the gas dynamic trap in the Budker Institute of Nuclear Physics [Shalashov et al., Phys. Plasmas 19, 052503 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Magnetic islands and singular currents at rational surfaces in three-dimensional magnetohydrodynamic equilibria

Using the recently developed multiregion, relaxed MHD (MRxMHD) theory, which bridges the gap between Taylor's relaxation theory and ideal MHD, we provide a thorough analytical and numerical proof of the formation of singular currents at rational surfaces in non-axisymmetric ideal MHD equilibria. These include the force-free singular current density represented by a Dirac δ-function, which presumably prevents the formation of islands, and the Pfirsch-Schlüter 1/x singular current, which arises as a result of finite pressure gradient. An analytical model based on linearized MRxMHD is derived that can accurately (1) describe the formation of magnetic islands at resonant rational surfaces, (2) retrieve the ideal MHD limit where magnetic islands are shielded, and (3) compute the subsequent formation of singular currents. The analytical results are benchmarked against numerical simulations carried out with a fully nonlinear implementation of MRxMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Closed bipolar electrode-enabled dual-cell electrochromic detectors for chemical sensing

Bipolar electrodes (BPE) are electrically floating metallic elements placed in electrified fluids that enable the coupling of anodic and cathodic redox reactions at the opposite ends by electron transfer through the electrode. One particularly compelling application allows electron transfer reactions at one end of a closed BPE to be read out optically by inducing a redox-initiated change in the optical response function of a reporter system at the other end. Here, a BPE-enabled method for electrochemical sensing based on the electrochromic response of a methyl viologen (MV) reporter is developed, characterized, and rendered in a field-deployable format. BPE-enabled devices based on two thin-layer-cells of ITO and Pt were fabricated to couple an analytical reaction in one cell with an MV reporter reaction, producing a color change in the complementary cell. Using Fe(CN) 6 3/4- as a model analyte, the electrochemically induced color change of MV was determined initially by measuring its absorbance via a CCD camera coupled to a microscope. Then, smartphone-based detection and RGB analysis were employed to further simplify the sensing scheme. In both methods a linear relationship between the analyte concentration was produced, the quantity of MV generated, and the colorimetric response, yielding a limit of detection of 1.0 μM. Similar responses were also observed in the detection of dopamine and acetaminophen. Further evolution of the device replaced the potentiostat with batteries to control potential, demonstrating the simplicity and portability of the device. Finally, the physical separation of the reporter and analytical cells renders the device competent to detect analytes in different (e.g. non-aqueous) phases, as demonstrated by using the electrochromic behavior of aqueous MV to detect ferrocene in acetonitrile in the analytical cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

RELAXATION OF WARPED DISKS: THE CASE OF PURE HYDRODYNAMICS

Orbiting disks may exhibit bends due to a misalignment between the angular momentum of the inner and outer regions of the disk. We begin a systematic simulational inquiry into the physics of warped disks with the simplest case: the relaxation of an unforced warp under pure fluid dynamics, i.e., with no internal stresses other than Reynolds stress. We focus on the nonlinear regime in which the bend rate is large compared to the disk aspect ratio. When warps are nonlinear, strong radial pressure gradients drive transonic radial motions along the disk's top and bottom surfaces that efficiently mix angular momentum. The resulting nonlinear decay rate of the warp increases with the warp rate and the warp width, but, at least in the parameter regime studied here, is independent of the sound speed. The characteristic magnitude of the associated angular momentum fluxes likewise increases with both the local warp rate and the radial range over which the warp extends; it also increases with increasing sound speed, but more slowly than linearly. The angular momentum fluxes respond to the warp rate after a delay that scales with the square root of the time for sound waves to cross the radial extent of the warp. These behaviors are at variance with a number of the assumptions commonly used in analytic models to describe linear warp dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗