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 235 records · Page 13
All sourcesNASA NTRSDOE 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↗
NASA NTRS1977

Application of perturbation methods to improve analytical model correlation with test data

A linear perturbation method developed by White (1972) to improve the correlation between analytical models and experimental data is applied to modeling of the Viking Propulsion Subsystem (VPS), flown successfully as part of the Viking Spacecraft launch of 1975. Analytical models and modal test data developed for the VPS are employed in the application; a set of mathematically generated data and a corresponding perturbed analytical model, representative of actual flight hardware, are also subjected to the White method. It is found that the linear perturbation method does not yield a single unique improved model.

Garba, J. A.↗
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↗
NASA NTRS1987

Mesoscale lake-effect snowstorms in the vicinity of Lake Michigan - Linear theory and numerical simulations

Mesoscale lake-effect snowstorms in the vicinity of Lake Michigan are studied by a linear steady-state analytic model and a nonlinear time-dependent numerical model with parameterized subgrid-scale physics. The solutions of the linear model show that the orientation of the mean wind field to the surface heating pattern is crucial to the shapes of the disturbances. The results indicate that the relative warmth of the lake surface can induce three updraft centers under a westerly wind, two updraft centers/bands under a northwesterly wind, and a convergence band under a northerly wind. Such convergences are caused by the interaction between the mean wind and the local circulations forced by the curved thermal contrasts. The numerical results from the nonlinear model not only produce convergence patterns very similar to those found in the linear theory and in other numerical studies, but also capture the transient property in some of the lake-induced disturbances. All of these results are qualitatively confirmed by satellite images.

Hsu, Hsiao-Ming↗
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↗
NASA NTRS2017

A Review of the Proposed KIsi Offset-Secant Method for Size-Insensitive Linear-Elastic Fracture Toughness Evaluation

Recently proposed modifications to ASTM E399 would provide a new size-insensitive approach to analyzing the force-displacement test record. The proposed size-insensitive linear-elastic fracture toughness, KIsi, targets a consistent 0.5mm crack extension for all specimen sizes by using an offset secant that is a function of the specimen ligament length. The KIsi evaluation also removes the Pmax/PQ criterion and increases the allowable specimen deformation. These latter two changes allow more plasticity at the crack tip, prompting the review undertaken in this work to ensure the validity of this new interpretation of the force-displacement curve. This paper provides a brief review of the proposed KIsi methodology and summarizes a finite element study into the effects of increased crack tip plasticity on the method given the allowance for additional specimen deformation. The study has two primary points of investigation: the effect of crack tip plasticity on compliance change in the force-displacement record and the continued validity of linear-elastic fracture mechanics to describe the crack front conditions. The analytical study illustrates that linear-elastic fracture mechanics assumptions remain valid at the increased deformation limit; however, the influence of plasticity on the compliance change in the test record is problematic. A proposed revision to the validity criteria for the KIsi test method is briefly discussed.

James, Mark↗
NASA NTRS2017

A Review of the Proposed KIsi Offset-Secant Method for Size-Insensitive Linear-Elastic Fracture Toughness Evaluation

Recently proposed modifications to ASTM E399 would provide a new size-insensitive approach to analyzing the force-displacement test record. The proposed size-insensitive linear-elastic fracture toughness, KIsi, targets a consistent 0.5mm crack extension for all specimen sizes by using an offset secant that is a function of the specimen ligament length. The KIsi evaluation also removes the Pmax/PQ criterion and increases the allowable specimen deformation. These latter two changes allow more plasticity at the crack tip, prompting the review undertaken in this work to ensure the validity of this new interpretation of the force-displacement curve. This paper provides a brief review of the proposed KIsi methodology and summarizes a finite element study into the effects of increased crack tip plasticity on the method given the allowance for additional specimen deformation. The study has two primary points of investigation: the effect of crack tip plasticity on compliance change in the force-displacement record and the continued validity of linear-elastic fracture mechanics to describe the crack front conditions. The analytical study illustrates that linear-elastic fracture mechanics assumptions remain valid at the increased deformation limit; however, the influence of plasticity on the compliance change in the test record is problematic. A proposed revision to the validity criteria for the KIsi test method is briefly discussed.

James, Mark↗
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↗