Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic functions”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Application of Contraction Mappings to the Control of Nonlinear Systems

The theoretical and applied aspects of successive approximation techniques are considered for the determination of controls for nonlinear dynamical systems. Particular emphasis is placed upon the methods of contraction mappings and modified contraction mappings. It is shown that application of the Pontryagin principle to the optimal nonlinear regulator problem results in necessary conditions for optimality in the form of a two point boundary value problem (TPBVP). The TPBVP is represented by an operator equation and functional analytic results on the iterative solution of operator equations are applied. The general convergence theorems are translated and applied to those operators arising from the optimal regulation of nonlinear systems. It is shown that simply structured matrices and similarity transformations may be used to facilitate the calculation of the matrix Green functions and the evaluation of the convergence criteria. A controllability theory based on the integral representation of TPBVP's, the implicit function theorem, and contraction mappings is developed for nonlinear dynamical systems. Contraction mappings are theoretically and practically applied to a nonlinear control problem with bounded input control and the Lipschitz norm is used to prove convergence for the nondifferentiable operator. A dynamic model representing community drug usage is developed and the contraction mappings method is used to study the optimal regulation of the nonlinear system.

Killingsworth, W. R., Jr.↗
DOE OSTI2019

Zeros of Riemann’s Zeta Functions in the Line z=1/2+it{sub 0}

It was found that, in addition to trivial zeros in points (z = − 2N, N = 1, 2…, natural numbers), the Riemann’s zeta function ζ(z) has zeros only on the line {z=(1/2)+it_0, t_0 is real}. All zeros are numerated, and for each number, N, the positions of the non-overlap intervals with one zero inside are found. The simple equation for the determination of centers of intervals is obtained. The analytical function η(z), leading to the possibility fix the zeros of the zeta function ζ(z), was estimated. To perform the analysis, the well-known phenomenon, phase-slip events, is used. This phenomenon is the key ingredient for the investigation of dynamical processes in solid-state physics, for example, if we are trying to solve the TDGLE (time-dependent Ginzburg-Landau equation).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2007

Representing Functions in n Dimensions to Arbitrary Accuracy

A method of approximating a scalar function of n independent variables (where n is a positive integer) to arbitrary accuracy has been developed. This method is expected to be attractive for use in engineering computations in which it is necessary to link global models with local ones or in which it is necessary to interpolate noiseless tabular data that have been computed from analytic functions or numerical models in n-dimensional spaces of design parameters.

Scotti, Stephen J.↗
NASA NTRS2014

Modeling Sound Propagation Through Non-Axisymmetric Jets

A method for computing the far-field adjoint Green's function of the generalized acoustic analogy equations under a locally parallel mean flow approximation is presented. The method is based on expanding the mean-flow-dependent coefficients in the governing equation and the scalar Green's function in truncated Fourier series in the azimuthal direction and a finite difference approximation in the radial direction in circular cylindrical coordinates. The combined spectral/finite difference method yields a highly banded system of algebraic equations that can be efficiently solved using a standard sparse system solver. The method is applied to test cases, with mean flow specified by analytical functions, corresponding to two noise reduction concepts of current interest: the offset jet and the fluid shield. Sample results for the Green's function are given for these two test cases and recommendations made as to the use of the method as part of a RANS-based jet noise prediction code.

Leib, Stewart J.↗
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 OSTI2023

Empirical Study of the Effect of Thermal Loading on the Heating Efficiency of Variable-Speed Air Source Heat Pumps

Heating buildings with air source heat pumps (ASHPs) has the potential to save energy compared to utilizing conventional heat sources. Accurate understanding of the efficiency of ASHPs is important to maximize the energy savings. While it is well understood that, in general, ASHP efficiency decreases with decreasing outdoor temperature, it is not well understood how the ASHP efficiency changes with different levels of thermal loading, even though it is an important consideration for sizing and controlling ASHPs. The goal of this study was to create an empirical model of the ASHP efficiency as a function of two independent variables–outside temperature and level of thermal loading. Four ductless mini-split ASHPs were evaluated in a cold chamber where the temperature (representing the outdoor temperature) was varied over a wide range. For each temperature, the ASHP performance data were collected at several levels of thermal loading. The data for all four ASHPs were combined and approximated with an analytical function that can be used as a general model for the ASHP steady-state efficiency as a function of the outside temperature and level of thermal loading. To the knowledge of the authors, no such empirical model that is solely based on third-party test data has been published before. While limitations exist, the model can be used to help guide future selection and operation of ASHPs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1974

Theoretical study of production of unique glasses in space

Analytical functional relationships describing homogeneous nucleation and crystallization in various supercooled liquids were developed. The time and temperature dependent relationships of nucleation and crystallization (intrinsic properties) are being used to relate glass forming tendency to extrinsic parameters such as cooling rate through computer simulation. Single oxide systems are being studied initially to aid in developing workable kinetic models and to indicate the primary materials parameters affecting glass formation. The theory and analytical expressions developed for simple systems is then extended to complex oxide systems. A thorough understanding of nucleation and crystallization kinetics of glass forming systems provides a priori knowledge of the ability of a given system to form a glass.

Larsen, D. C.↗
DOE OSTI2012

Benchmarking a new closed-form thermal analysis technique against a traditional lumped parameter, finite-difference method

A benchmarking effort was conducted to determine the accuracy of a new analytic generic geology thermal repository model developed at LLNL relative to a more traditional, numerical, lumped parameter technique. The fast-running analytical thermal transport model assumes uniform thermal properties throughout a homogenous storage medium. Arrays of time-dependent heat sources are included geometrically as arrays of line segments and points. The solver uses a source-based linear superposition of closed form analytical functions from each contributing point or line to arrive at an estimate of the thermal evolution of a generic geologic repository. Temperature rise throughout the storage medium is computed as a linear superposition of temperature rises. It is modeled using the MathCAD mathematical engine and is parameterized to allow myriad gridded repository geometries and geologic characteristics [4]. It was anticipated that the accuracy and utility of the temperature field calculated with the LLNL analytical model would provide an accurate 'birds-eye' view in regions that are many tunnel radii away from actual storage units; i.e., at distances where tunnels and individual storage units could realistically be approximated as physical lines or points. However, geometrically explicit storage units, waste packages, tunnel walls and close-in rock are not included in the MathCAD model.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

The dependence of divertor power sharing on magnetic flux balance in near double-null configurations on Alcator C-Mod

We report management of power exhaust will be a crucial task for tokamak fusion reactors. Reactor concepts are often proposed with double-null divertors, i.e. having two magnetic separatrices in an up-down symmetric configuration. This arrangement is potentially advantageous since the majority of the tokamak exhaust power tends to flow to the outer pair of divertor legs at large major radius, where the geometry is favorable for spreading the heat over a large surface area and there is more room for advanced divertor configurations. Despite the importance, there have been relatively few studies of divertor power sharing in near double null configurations and no studies at the poloidal magnetic fields and scrape-off layer power widths anticipated for a reactor. Motivated by this need we have undertaken a systematic study on Alcator C-Mod, examining the effect of magnetic flux balance on the power sharing among the four divertor legs in near double-null plasmas. Ohmic L-modes at three values of plasma current and ICRF-heated enhanced D-alpha (EDA) H-modes and I-modes at a single value of plasma current are explored, producing poloidal magnetic fields of 0.42, 0.62 and 0.85 Tesla. For Ohmic L-modes and ICRF-heated EDA H-modes, we find that the point of equal power sharing between upper and lower divertors occurs remarkably close to a balanced double null. Power sharing amongst the outer (upper versus lower) and inner (upper versus lower) pairs of divertors can be described in terms of a logistic function of magnetic flux balance, consistent with heat flux mapping along magnetic field lines to the outer midplane. Power sharing between inner and outer legs is found to follow a Gaussian-like function of magnetic flux balance with non-zero power to the inner divertors at double null. The overall behavior of H-modes operated near double null and for I-modes operating to within one heat flux e-folding of double null are found similar to Ohmic L-modes, with a significant reduction of power on the inner divertor legs. The results are encapsulated in terms of empirically-informed analytic functions of magnetic flux balance. When combined with magnetic equilibrium control system specifications, these relationships can be used to specify the power flux handling requirements for each of the four divertor target plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

The dependence of divertor power sharing on magnetic flux balance in near double-null configurations on Alcator C-Mod

Management of power exhaust will be a crucial task for tokamak fusion reactors. Reactor concepts are often proposed with double-null divertors, i.e. having two magnetic separatrices in an up-down symmetric configuration. This arrangement is potentially advantageous since the majority of the tokamak exhaust power tends to flow to the outer pair of divertor legs at large major radius, where the geometry is favorable for spreading the heat over a large surface area and there is more room for advanced divertor configurations. Despite the importance, there have been relatively few studies of divertor power sharing in near double null configurations and no studies at the poloidal magnetic fields and scrape-off layer power widths anticipated for a reactor. Motivated by this need we have undertaken a systematic study on Alcator C-Mod, examining the effect of magnetic flux balance on the power sharing among the four divertor legs in near double-null plasmas. Ohmic L-modes at three values of plasma current and ICRF-heated enhanced D-alpha (EDA) H-modes and I-modes at a single value of plasma current are explored, producing poloidal magnetic fields of 0.42, 0.62 and 0.85 Tesla. For Ohmic L-modes and ICRF-heated EDA H-modes, we find that the point of equal power sharing between upper and lower divertors occurs remarkably close to a balanced double null. Power sharing amongst the outer (upper versus lower) and inner (upper versus lower) pairs of divertors can be described in terms of a logistic function of magnetic flux balance, consistent with heat flux mapping along magnetic field lines to the outer midplane. Power sharing between inner and outer legs is found to follow a Gaussian-like function of magnetic flux balance with non-zero power to the inner divertors at double null. The overall behavior of H-modes operated near double null and for I-modes operating to within one heat flux e-folding of double null are found similar to Ohmic L-modes, with a significant reduction of power on the inner divertor legs. The results are encapsulated in terms of empirically-informed analytic functions of magnetic flux balance. When combined with magnetic equilibrium control system specifications, these relationships can be used to specify the power flux handling requirements for each of the four divertor target plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1958

Elliptic Functions and Integrals with Real Modulus in Fluid Mechanics

Advantage of the elliptic functions and of the more general functions of Schwarz for fluid mechanics. Flows outside and inside polygons. Application to the calculation of an elbow diffuser for a wind tunnel. Properties of the elliptic integrals of the first kind and of the elliptic functions. Properties of the theta functions and decomposition of the elliptic functions into products of theta functions. Properties of the zeta functions. Decomposition of the elliptic functions into sums of zeta functions and calculations of the elliptic integrals. Applications to the calculation of wing profiles, of compressor profiles, and to the study of the vibrations of airplane wings and of compressor vanes. The manuscript of the present paper was checked by Mr. Eichelbrenner who corrected several imperfections and suggested numerous improvements to make reading of the paper easier. However, the limited subject does not permit filling in more than an incomplete knowledge of the properties of analytic functions.

Legendre, Robert↗
NASA NTRS1999

Ocean Turbulence: One-Point Closure Model Momentum, Heat and Salt Vertical Diffusivities in the Presence of Shear - Paper 2

We develop and test a 1-point closure turbulence model with the following features: 1) we include the salinity field and derive the expression for the vertical turbulent diffusivities of momentum K(sub m) , heat K(sub h) and salt K(sub s) as a function of two stability parameters: the Richardson number R(sub i) (stratification vs. shear) and the Turner number R(sub rho) (salinity gradient vs. temperature gradient). 2) to describe turbulent mixing below the mixed layer (ML), all previous models have adopted three adjustable "background diffusivities" for momentum, heat and salt. We propose a model that avoids such adjustable diffusivities. We assume that below the ML, the three diffusivities have the same functional dependence on R( sub i) and R(sub rho) as derived from the turbulence model. However, in order to compute R(sub i) below the ML, we use data of vertical shear due to wave-breaking.measured by Gargett et al. The procedure frees the model from adjustable background diffusivities and indeed we employ the same model throughout the entire vertical extent of the ocean. 3) in the local model, the turbulent diffusivities K(sub m,h,s) are given as analytical functions of R(sub i) and R(sub rho). 5) the model is used in an O-GCM and several results are presented to exhibit the effect of double diffusion processes. 6) the code is available upon request.

Canuto, V. M.↗
DOE OSTI2022

Module-Fluidics: Building Blocks for Spatio-Temporal Microenvironment Control

Generating the desired solute concentration signal in micro-environments is vital to many applications ranging from micromixing to analyzing cellular response to a dynamic microenvironment. We propose a new modular design to generate targeted temporally varying concentration signals in microfluidic systems while minimizing perturbations to the flow field. The modularized design, here referred to as module-fluidics, similar in principle to interlocking toy bricks, is constructed from a combination of two building blocks and allows one to achieve versatility and flexibility in dynamically controlling input concentration. The building blocks are an oscillator and an integrator, and their combination enables the creation of controlled and complex concentration signals, with different user-defined time-scales. We show two basic connection patterns, in-series and in-parallel, to test the generation, integration, sampling and superposition of temporally-varying signals. All such signals can be fully characterized by analytic functions, in analogy with electric circuits, and allow one to perform design and optimization before fabrication. Such modularization offers a versatile and promising platform that allows one to create highly customizable time-dependent concentration inputs which can be targeted to the specific application of interest.

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

Fourier-transform microwave spectroscopy and determination of the three dimensional potential energy surface for Ar–CS

Pure rotational transitions of the Ar–CS van der Waals complex have been observed by Fourier Transform Microwave (FTMW) and FTMW-millimeter wave double resonance spectroscopy. Rotational transitions of v{sub s} = 0, 1, and 2 were able to be observed for normal CS, together with those of C{sup 34}S in v{sub s} = 0, where v{sub s} stands for the quantum number of the CS stretching vibration. The observed transition frequencies were analyzed by a free rotor model Hamiltonian, where rovibrational energies were calculated as dynamical motions of the three nuclei on a three-dimensional potential energy surface, expressed by analytical functions with 57 parameters. Initial values for the potential parameters were obtained by high-level ab initio calculations. Fifteen parameters were adjusted among the 57 parameters to reproduce all the observed transition frequencies with the standard deviation of the fit to be 0.028 MHz.

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

Sum-rule corrections: A route to error cancellations in correlation matrix renormalisation theory

Here, we recently proposed the correlation matrix renormalisation (CMR) theory to efficiently and accurately calculate ground state total energy of molecular systems, based on the Gutzwiller variational wavefunction (GWF) to treat the electronic correlation effects. To help reduce numerical complications and better adapt the CMR to infinite lattice systems, we need to further refine the way to minimise the error originated from the approximations in the theory. This conference proceeding reports our recent progress on this key issue, namely, we obtained a simple analytical functional form for the one-electron renormalisation factors, and introduced a novel sum-rule correction for a more accurate description of the intersite electron correlations. Benchmark calculations are performed on a set of molecules to show the reasonable accuracy of the method.

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

The Symmetrical Quasi-Classical Model for Electronically Non-Adiabatic Processes Applied to Energy Transfer Dynamics in Site-Exciton Models of Light-Harvesting Complexes

In a recent series of papers, it has been illustrated that a symmetrical quasi-classical (SQC) windowing model applied to the Meyer–Miller (MM) classical vibronic Hamiltonian provides an excellent description of a variety of electronically non-adiabatic benchmark model systems for which exact quantum results are available for comparison. In this study, the SQC/MM approach is used to treat energy transfer dynamics in site-exciton models of light-harvesting complexes, and in particular, the well-known 7-state Fenna–Mathews–Olson (FMO) complex. Again, numerically “exact” results are available for comparison, here via the hierarchical equation of motion (HEOM) approach of Ishizaki and Fleming, and it is seen that the simple SQC/MM approach provides very reasonable agreement with the previous HEOM results. Finally, it is noted, however, that unlike most (if not all) simple approaches for treating these systems, because the SQC/MM approach presents a fully atomistic simulation based on classical trajectory simulation, it places no restrictions on the characteristics of the thermal baths coupled to each two-level site, e.g., bath spectral densities (SD) of any analytic functional form may be employed as well as discrete SD determined experimentally or from MD simulation (nor is there any restriction that the baths be harmonic), opening up the possibility of simulating more realistic variations on the basic site-exciton framework for describing the non-adiabatic dynamics of photosynthetic pigment complexes.

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

Comparing correlation components and approximations in Hartree–Fock and Kohn–Sham theories via an analytical test case study

The asymmetric Hubbard dimer is a model that allows for explicit expressions of the Hartree–Fock (HF) and Kohn–Sham (KS) states as analytical functions of the external potential, Δv, and of the interaction strength, U. We use this unique circumstance to establish a rigorous comparison between the individual contributions to the correlation energies stemming from the two theories in the {U, Δv} parameter space. Within this analysis of the Hubbard dimer, we observe a change in the sign of the HF kinetic correlation energy, compare the indirect repulsion energies, and derive an expression for the “traditional” correlation energy, i.e., the one that corrects the HF estimate, in a pure site-occupation function theory spirit [Eq. (45)]. Next, we test the performances of the Liu–Burke and the Seidl–Perdew–Levy functionals, which model the correlation energy based on its weak- and strong-interaction limit expansions and can be used for both the traditional and the KS correlation energies. Further, our results show that, in the Hubbard dimer setting, they typically work better for the HF reference, despite having been originally devised for KS. These conclusions are somewhat in line with prior assessments of these functionals on various chemical datasets. However, the Hubbard dimer model allows us to show the extent of the error that may occur in using the strong-interaction ingredient for the KS reference in place of the one for the HF reference, as has been carried out in most of the prior assessments.

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

Accelerating Markov Chain Monte Carlo sampling with diffusion models

Global fits of physics models require efficient methods for exploring high-dimensional and/or multimodal posterior functions. We introduce a novel method for accelerating Markov Chain Monte Carlo (MCMC) sampling by pairing a Metropolis-Hastings algorithm with a diffusion model that can draw global samples with the aim of approximating the posterior. We briefly review diffusion models in the context of image synthesis before providing a streamlined diffusion model tailored towards low-dimensional data arrays. We then present our adapted Metropolis-Hastings algorithm which combines local proposals with global proposals taken from a diffusion model that is regularly trained on the samples produced during the MCMC run. Our approach leads to a significant reduction in the number of likelihood evaluations required to obtain an accurate representation of the Bayesian posterior across several analytic functions, as well as for a physical example based on a global fit of parton distribution functions. Our method is extensible to other MCMC techniques, and we briefly compare our method to similar approaches based on normalising flows. A code implementation can be found at https://github.com/NickHunt-Smith/MCMC-diffusion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗