Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPONENTIAL FUNCTION”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Modeling of Aircraft Unsteady Aerodynamic Characteristics/Part 3 - Parameters Estimated from Flight Data: Parameters Estimated from Flight Data - Part 3

A nonlinear least squares algorithm for aircraft parameter estimation from flight data was developed. The postulated model for the analysis represented longitudinal, short period motion of an aircraft. The corresponding aerodynamic model equations included indicial functions (unsteady terms) and conventional stability and control derivatives. The indicial functions were modeled as simple exponential functions. The estimation procedure was applied in five examples. Four of the examples used simulated and flight data from small amplitude maneuvers to the F-18 HARV and X-31A aircraft. In the fifth example a rapid, large amplitude maneuver of the X-31 drop model was analyzed. From data analysis of small amplitude maneuvers ft was found that the model with conventional stability and control derivatives was adequate. Also, parameter estimation from a rapid, large amplitude maneuver did not reveal any noticeable presence of unsteady aerodynamics.

Klein, Vladislav↗
NASA NTRS1998

A Solution to the Fundamental Linear Fractional Order Differential Equation

This paper provides a solution to the fundamental linear fractional order differential equation, namely, (sub c)d(sup q, sub t) + ax(t) = bu(t). The impulse response solution is shown to be a series, named the F-function, which generalizes the normal exponential function. The F-function provides the basis for a qth order "fractional pole". Complex plane behavior is elucidated and a simple example, the inductor terminated semi- infinite lossy line, is used to demonstrate the theory.

Hartley, Tom T.↗
DOE OSTI2018

Fractional heat conduction with heat absorption in a sphere under Dirichlet boundary condition

The time-fractional heat conduction equation with the Caputo derivative and with heat absorption term proportional to temperature is considered in a sphere in the case of central symmetry. The fundamental solution to the Dirichlet boundary value problem is found, and the solution to the problem under constant boundary value of temperature is studied. The integral transform technique is used. The solutions are obtained in terms of series containing the Mittag-Leffler functions being the generalization of the exponential function. The numerical results are illustrated graphically.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1991

Determination of a natural basis function set on thick, structured slabs using Prony's method

The calculation of currents on a thick, structured slab such as a thick slab of honeycomb is discussed. Unfortunately, for certain applications, the slab can be on the order of several wavelengths, so that the straightforward application of the method of moments using subdomain basis functions is too expensive. The authors discuss how to apply Prony's method to the currents calculated for a thin, structured slab to obtain a natural set of basis functions to represent the currents in the interior of a thick slab. Prony's method represents the currents as a series of complex exponential functions. The thick slab problem is then solved by the method of moments using subdomain basis functions near the slab interfaces and one or two of the complex exponentials as basis functions within the slab.

Jorgenson, Roy E.↗
DOE OSTI2025

Multioutput Convolutional Neural Network for Improved Parameter Extraction in Time-Resolved Electrostatic Force Microscopy Data

Time-resolved scanning probe microscopy methods, like time-resolved electrostatic force microscopy (trEFM), enable imaging of dynamic processes ranging from ion motion in batteries to electronic dynamics in microstructured thin film semiconductors for solar cells. Reconstructing the underlying physical dynamics from these techniques can be challenging due to the interplay of cantilever physics with the actual transient kinetics of interest in the resulting signal. Previously, quantitative trEFM used empirical calibration of the cantilever or feed-forward neural networks trained on simulated data to extract the physical dynamics of interest. Both these approaches are limited by interpreting the underlying signal as a single exponential function, which serves as an approximation but does not adequately reflect many realistic systems. Here, we present a multi-branched, multi-output convolutional neural network (CNN) that uses the trEFM signal in addition to the physical cantilever parameters as input. The trained CNN accurately extracts parameters describing both single-exponential and bi-exponential underlying functions, and more accurately reconstructs real experimental data in the presence of noise. This article demonstrates an application of physics-informed machine learning to complex signal processing tasks, enabling more efficient and accurate analysis of trEFM.

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

Diffusion-Controlled growth of Ge nanocrystals in SiO{sub 2} films under conditions of ion synthesis at high pressure

The growth of Ge nanocrystals in SiO{sub 2} films is studied in relation to the dose of implanted Ge{sup +} ions and the annealing temperature at a pressure of 12 kbar. It is established that the dependences of the nanocrystal dimensions on the content of Ge atoms and the annealing time are described by the corresponding root functions. The nanocrystal radius squared is an exponential function of the inverse temperature. The dependences correspond to the model of the diffusion-controlled mechanism of nanocrystal growth. From the temperature dependence of the nanocrystal dimensions, the diffusion coefficient of Ge in SiO{sub 2} at a pressure of 12 kbar is determined: D = 1.1 × 10{sup –10} exp(–1.43/kT). An increase in the diffusion coefficient of Ge under pressure is attributed to the change in the activation volume of the formation and migration of point defects. Evidence in favor of the interstitial mechanism of the diffusion of Ge atoms to nanocrystal nuclei in SiO{sub 2} is reported.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Ocean feature recognition using genetic algorithms with fuzzy fitness functions (GA/F3)

A model for genetic algorithms with semantic nets is derived for which the relationships between concepts is depicted as a semantic net. An organism represents the manner in which objects in a scene are attached to concepts in the net. Predicates between object pairs are continuous valued truth functions in the form of an inverse exponential function (e sub beta lxl). 1:n relationships are combined via the fuzzy OR (Max (...)). Finally, predicates between pairs of concepts are resolved by taking the average of the combined predicate values of the objects attached to the concept at the tail of the arc representing the predicate in the semantic net. The method is illustrated by applying it to the identification of oceanic features in the North Atlantic.

Ankenbrandt, C. A.↗
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

Regression-based reduced-order models to predict transient thermal output for enhanced geothermal systems

Reduced-order modeling is a promising approach, as many phenomena can be described by a few parameters/mechanisms. An advantage and attractive aspect of a reduced-order model is that it is computational inexpensive to evaluate when compared to running a high-fidelity numerical simulation. A reduced-order model takes couple of seconds to run on a laptop while a high-fidelity simulation may take couple of hours to run on a high-performance computing cluster. The goal of this paper is to assess the utility of regression-based reduced-order models (ROMs) developed from high-fidelity numerical simulations for predicting transient thermal power output for an enhanced geothermal reservoir while explicitly accounting for uncertainties in the subsurface system and site-specific details. Numerical simulations are performed based on equally spaced values in the specified range of model parameters. Key sensitive parameters are then identified from these simulations, which are fracture zone permeability, well/skin factor, bottom hole pressure, and injection flow rate. We found the fracture zone permeability to be the most sensitive parameter. The fracture zone permeability along with time, are used to build regression-based ROMs for the thermal power output. The ROMs are trained and validated using detailed physics-based numerical simulations. Finally, predictions from the ROMs are then compared with field data. We propose three different ROMs with different levels of model parsimony, each describing key and essential features of the power production curves. The coefficients in the proposed regression-based ROMs are developed by minimizing a non-linear least-squares misfit function using the Levenberg–Marquardt algorithm. The misfit function is based on the difference between numerical simulation data and reduced-order model. ROM-1 is constructed based on polynomials up to fourth order. ROM-1 is able to accurately reproduce the power output of numerical simulations for low values of permeabilities and certain features of the field-scale data. ROM-2 is a model with more analytical functions consisting of polynomials up to order eight, exponential functions and smooth approximations of Heaviside functions, and accurately describes the field-data. At higher permeabilities, ROM-2 reproduces numerical results better than ROM-1, however, there is a considerable deviation from numerical results at low fracture zone permeabilities. ROM-3 consists of polynomials up to order ten, and is developed by taking the best aspects of ROM-1 and ROM-2. ROM-1 is relatively parsimonious than ROM-2 and ROM-3, while ROM-2 overfits the data. ROM-3 on the other hand, provides a middle ground for model parsimony. Based on R 2 -values for training, validation, and prediction data sets we found that ROM-3 is better model than ROM-2 and ROM-1. For predicting thermal drawdown in EGS applications, where high fracture zone permeabilities (typically greater than 10 –15 m 2 ) are desired, ROM-2 and ROM-3 outperform ROM-1. As per computational time, all the ROMs are 10 4 times faster when compared to running a high-fidelity numerical simulation. In conclusion, this makes the proposed regression-based ROMs attractive for real-time EGS applications because they are fast and provide reasonably good predictions for thermal power output.

58 GEOSCIENCES↗
DOE OSTI2014

The rigorous stochastic matrix multiplication scheme for the calculations of reduced equilibrium density matrices of open multilevel quantum systems

Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, i.e., how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernon influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Anisotropic inflation in Gauss-Bonnet gravity

We study anisotropic inflation with Gauss-Bonnet correction in presence of a massless vector field. In this scenario, exact anisotropic power-law inflation is realized when the inflaton potential, gauge coupling function and the Gauss-Bonnet coupling are exponential functions. We show that anisotropy becomes proportional to two slow-roll parameters of the theory and hence gets enhanced in presence of quadratic curvature corrections. The stability analysis reveals that anisotropic power-law solutions remain stable over a substantially large parameter region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

THE ARIZONA RADIO OBSERVATORY CO MAPPING SURVEY OF GALACTIC MOLECULAR CLOUDS. V. THE SH2-235 CLOUD IN CO J = 2 − 1, {sup 13}CO J = 2 − 1, AND CO J = 3 − 2

We present the results of a program to map the Sh2-235 molecular cloud complex in the CO and {sup 13}CO J = 2 − 1 transitions using the Heinrich Hertz Submillimeter Telescope. The map resolution is 38″ (FWHM), with an rms noise of 0.12 K brightness temperature, for a velocity resolution of 0.34 km s{sup −1}. With the same telescope, we also mapped the CO J = 3 − 2 line at a frequency of 345 GHz, using a 64 beam focal plane array of heterodyne mixers, achieving a typical rms noise of 0.5 K brightness temperature with a velocity resolution of 0.23 km s{sup −1}. The three spectral line data cubes are available for download. Much of the cloud appears to be slightly sub-thermally excited in the J = 3 level, except for in the vicinity of the warmest and highest column density areas, which are currently forming stars. Using the CO and {sup 13}CO J = 2 − 1 lines, we employ an LTE model to derive the gas column density over the entire mapped region. Examining a 125 pc{sup 2} region centered on the most active star formation in the vicinity of Sh2-235, we find that the young stellar object surface density scales as approximately the 1.6-power of the gas column density. The area distribution function of the gas is a steeply declining exponential function of gas column density. Comparison of the morphology of ionized and molecular gas suggests that the cloud is being substantially disrupted by expansion of the H ii regions, which may be triggering current star formation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

An approximate solution for interlaminar stresses in laminated composites: Applied mechanics program

An approximate solution for interlaminar stresses in finite width, laminated composites subjected to uniform extensional, and bending loads is presented. The solution is based upon the principle of minimum complementary energy and an assumed, statically admissible stress state, derived by considering local material mismatch effects and global equilibrium requirements. The stresses in each layer are approximated by polynomial functions of the thickness coordinate, multiplied by combinations of exponential functions of the in-plane coordinate, expressed in terms of fourteen unknown decay parameters. Imposing the stationary condition of the laminate complementary energy with respect to the unknown variables yields a system of fourteen non-linear algebraic equations for the parameters. Newton's method is implemented to solve this system. Once the parameters are known, the stresses can be easily determined at any point in the laminate. Results are presented for through-thickness and interlaminar stress distributions for angle-ply, cross-ply (symmetric and unsymmetric laminates), and quasi-isotropic laminates subjected to uniform extension and bending. It is shown that the solution compares well with existing finite element solutions and represents an improved approximate solution for interlaminar stresses, primarily at interfaces where global equilibrium is satisfied by the in-plane stresses, but large local mismatch in properties requires the presence of interlaminar stresses.

Rose, Cheryl A.↗
DOE OSTI2013

SECULAR BEHAVIOR OF EXOPLANETS: SELF-CONSISTENCY AND COMPARISONS WITH THE PLANET-PLANET SCATTERING HYPOTHESIS

If mutual gravitational scattering among exoplanets occurs, then it may produce unique orbital properties. For example, two-planet systems that lie near the boundary between circulation and libration of their periapses could result if planet-planet scattering ejected a former third planet quickly, leaving one planet on an eccentric orbit and the other on a circular orbit. We first improve upon previous work that examined the apsidal behavior of known multiplanet systems by doubling the sample size and including observational uncertainties. This analysis recovers previous results that demonstrated that many systems lay on the apsidal boundary between libration and circulation. We then performed over 12,000 three-dimensional N-body simulations of hypothetical three-body systems that are unstable, but stabilize to two-body systems after an ejection. Using these synthetic two-planet systems, we test the planet-planet scattering hypothesis by comparing their apsidal behavior, over a range of viewing angles, to that of the observed systems and find that they are statistically consistent regardless of the multiplicity of the observed systems. Finally, we combine our results with previous studies to show that, from the sampled cases, the most likely planetary mass function prior to planet-planet scattering follows a power law with index -1.1. We find that this pre-scattering mass function predicts a mutual inclination frequency distribution that follows an exponential function with an index between -0.06 and -0.1.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Statistics-Based Correction to Measurements of Hardness of Rough Surfaces by Instrumented Indentation

Abstract—A procedure for correcting hardness data measured on rough surfaces using instrumented indentation is presented. The correction procedure relies on the knowledge of surface topography, which can be obtained by atomic force microscopy or other appropriate instrumental technique. The analyzed surface is characterized by the standard deviation of height distribution and an autocorrelation function. The latter is proposed to be modeled with an exponential function with either a quadratic or linear exponent.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1977

A new algorithm for the integration of exponential and logarithmic functions

An algorithm for symbolic integration of functions built up from the rational functions by repeatedly applying either the exponential or logarithm functions is discussed. This algorithm does not require polynomial factorization nor partial fraction decomposition and requires solutions of linear systems with only a small number of unknowns. It is proven that if this algorithm is applied to rational functions over the integers, a computing time bound for the algorithm can be obtained which is a polynomial in a bound on the integer length of the coefficients, and in the degrees of the numerator and denominator of the rational function involved.

Rothstein, M.↗
NASA NTRSDate not supplied

Modeling of Response of Thermoplastic Composites Under High Energy Impact

This work investigates the response of thermoplastic composites under impact loading. Panels fabricated with the material system under consideration – carbon fiber reinforced low melt semi-crystalline resin TC1225 LMPAEK T700G (T700/LMPAEK) – were subject to impact with an Aluminum projectile, at low and high velocities. At lower velocities, the projectile rebounded from the panel after impact and caused damage e.g. delaminations. Finite element models were developed to model this test using LS-DYNA®, with material model MAT299, a continuum damage model based on the LaRC04 criteria, to capture the dynamic response and damage modes in the T700/LMPAEK panels. Interlaminar delaminations were modeled using the cohesive contact formulation available within LS-DYNA. Previous work by the authors showed that delaminations were consistently underpredicted in size, relative to the test. It was surmised that this was due to effects of strain rate on intralaminar and interlaminar response of the T700/LMPAEK material system. Therefore, the current work is focused on incorporating strain rate effects in modeling the interlaminar response. Rate dependence was incorporated using ‘one-way tie-break contact’ in LS-DYNA, wherein, the Mode-II fracture toughness of the LMPAEK resin, G IIc , was assumed to be an exponential function of the strain rate. In particular, G IIc , was modeled to exponentially decrease at higher strain rates. Models with this form of rate dependence favorably correlated with experimental data regarding delamination extent and rebound velocity, relative to prior models without strain rate effects.

Composite Materials↗
DOE OSTI2017

Effect of the link lifetime in a dynamical lattice on the properties of the avalanche processes on it

The properties of the avalanche processes that develop on a dynamical lattice, the structure of links in which changes due to a specific characteristic of each lattice node, namely, its “activity,” which determines the probability of connection of a certain node with neighboring nodes in one step of lattice evolution. The statistics of the sizes of the avalanches appearing in the lattice system is studied as a function of the node activity and the link lifetime (the lifetime of the links formed in the system). It is analytically and numerically shows that the type of avalanche dynamics in the system changes as a function of these parameters. The following three regimes can take place in the system: (1) avalanches of any sizes, from small to catastrophic, can appear, which is reflected in the power-law behavior of the probability density function of the appearance of avalanches of certain sizes; (2) avalanches of a certain average size mainly appear in the system, and the probability density is close to that of a normal distribution; and (3) transient regime, where the probability density function of the appearance of avalanches of certain sizes is close to an exponential function. These results open up the possibilities of controlling the behavior of a complex system; in particular, they can be used to prevent catastrophic avalanches by changing the link lifetime and the average node activity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗