SEARCH · Search NASA

Results for “Integral calculus”

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 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Formalization of the Integral Calculus in the PVS Theorem Prover

The PVS Theorem prover is a widely used formal verification tool used for the analysis of safety-critical systems. The PVS prover, though fully equipped to support deduction in a very general logic framework, namely higher-order logic, it must nevertheless, be augmented with the definitions and associated theorems for every branch of mathematics and Computer Science that is used in a verification. This is a formidable task, ultimately requiring the contributions of researchers and developers all over the world. This paper reports on the formalization of the integral calculus in the PVS theorem prover. All of the basic definitions and theorems covered in a first course on integral calculus have been completed.The theory and proofs were based on Rosenlicht's classic text on real analysis and follow the traditional epsilon-delta method. The goal of this work was to provide a practical set of PVS theories that could be used for verification of hybrid systems that arise in air traffic management systems and other aerospace applications. All of the basic linearity, integrability, boundedness, and continuity properties of the integral calculus were proved. The work culminated in the proof of the Fundamental Theorem Of Calculus. There is a brief discussion about why mechanically checked proofs are so much longer than standard mathematics textbook proofs.

Butler, Ricky W.↗

NASA NTRS1985

Boolean integral calculus for digital systems

The concept of Boolean integration is introduced and developed. When the changes in a desired function are specified in terms of changes in its arguments, then ways of 'integrating' (i.e., realizing) the function, if it exists, are presented. Boolean integral calculus has applications in design of logic circuits.

Tucker, J. H.↗

NASA NTRS1988

Boolean integral calculus

The concept of Boolean integration is developed, and different Boolean integral operators are introduced. Given the changes in a desired function in terms of the changes in its arguments, the ways of 'integrating' (i.e. realizing) such a function, if it exists, are presented. The necessary and sufficient conditions for integrating, in different senses, the expression specifying the changes are obtained. Boolean calculus has applications in the design of logic circuits and in fault analysis.

Tucker, Jerry H.↗

NASA NTRS1992

A new approach to the digital redesign of continuous-time controllers

This paper presents a new method for digital redesign of both a continuous-time state-feedback control law and a cascaded analog controller. The concept of the Law of Mean from the integral calculus is utilized for the development of the equivalent digital controllers from the available analog controllers. A tuning parameter has been introduced into the digitally redesigned control gains so that the digitally controlled states will closely match the original continuous-time states. The proposed digital controllers give excellent performance.

Shieh, Leang-San↗

DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Computation of the Biot–Savart line integral with higher-order convergence using straight segments

One common approach to computing the magnetic field produced by a filamentary current-carrying coil is to approximate the coil as a series of straight segments. The Biot–Savart field from each straight segment is analytically known. However, if the endpoints of the straight segments are chosen to lie on the coil, then the accuracy of the Biot–Savart computation is generally only the second order in the number of endpoints. In this work, we propose a simple modification: shift each end point of the coil in the outward normal direction by an amount proportional to the local curvature. With this modification, the Biot–Savart accuracy increases to the fourth order and the numerical error is dramatically reduced for a given number of discretization points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A mathematical framework for ejecta cloud dynamics with application to source models and piezoelectric mass measurements

We present a mathematical framework for describing the dynamical evolution of an ejecta cloud generated by a generic ejecta source model. We consider a piezoelectric sensor fielded in the path of an ejecta cloud, for experimental configurations in which the ejecta are created at a singly shocked planar surface and fly ballistically through vacuum to the stationary sensor. To do so, we introduce the concept of a time- and velocity-dependent ejecta “areal mass function.” We derive expressions for the analytic (“true”) accumulated ejecta areal mass at the sensor and the measured (“inferred”) value obtained via the standard method for analyzing piezoelectric voltages. In this way, we derive an exact expression and upper bound for the error imposed upon a piezoelectric ejecta mass measurement (in a perfect system) by the assumption of instantaneous creation, which is commonly required for momentum diagnostic analyses. This error term is zero for truly instantaneous source models; otherwise, the standard piezoelectric analysis is guaranteed to overestimate the true mass. When combined with a piezoelectric dataset, this framework provides a unique solution for the ejecta particle velocity distribution, subject to the assumptions inherent in the data analysis. The framework also leads to strong boundary conditions that any ejecta source model must satisfy in order to be consistent with apparently global properties of piezoelectric measurements from a wide range of experiments. We demonstrate this methodology by applying it to the Richtmyer–Meshkov instability+self-similar velocity distribution ejecta source model currently under development at Los Alamos National Laboratory.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Analytic solutions for Asay foil trajectories with implications for ejecta source models and mass measurements

We consider the trajectory of an Asay foil ejecta diagnostic for scenarios where ejecta are produced at a singly shocked planar surface and fly ballistically through a perfect vacuum to the sensor. We do so by building upon a previously established mathematical framework derived for the analytic study of stationary sensors. First, we derive the momentum conservation equation for the problem, in a form amenable to accelerating sensors, in terms of a generic ejecta source model. The result is an integrodifferential equation of motion for the foil trajectory. This equation yields an easily calculable closed-form implicit solution for the foil trajectory in instant-production scenarios. From there, we derive a boundary condition that particle velocity distributions must satisfy if their associated foil trajectories are to exhibit a smooth initial acceleration, as occurs in some experiments. This condition is identical to one derived previously from a consideration of piezoelectric voltage data obtained in similar experiments. We also compare techniques for inferring accumulated ejecta masses from foil trajectories, first by deriving the exact solution, and then by quantifying the error imposed by a frequently used approximate solution (both subject to the assumption of instantaneous ejecta production). Finally, we examine the common practice of presenting inferred cumulative ejecta masses as a function of implied ejecta velocity, establishing the conditions under which this methodology is most meaningful.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Novel superposed kinklike and pulselike solutions for several nonlocal nonlinear equations

In this work, we show that a number of nonlocal nonlinear equations, including the Ablowitz–Musslimani and Yang variant of the nonlocal nonlinear Schrödinger (NLS) equation, the nonlocal modified Korteweg de Vries (mKdV) equation, and the nonlocal Hirota equation, admit novel kinklike and pulselike superposed periodic solutions. Furthermore, we show that the nonlocal mKdV equation also admits the superposed (hyperbolic) kink–antikink solution. In addition, we show that while the nonlocal Ablowitz–Musslimani variant of the NLS admits complex parity-time reversal-invariant kink and pulse solutions, neither the local NLS nor the Yang variant of the nonlocal NLS admits such solutions. Finally, except for the Yang variant of the nonlocal NLS, we show that the other three nonlocal equations admit both the kink and pulse solutions in the same model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Planar turbulent wakes under pressure gradient: Integral and self-similarity analyses

By using a combination of integral and self-similarity analyses, the generalized analytical solutions for the mean transverse velocity and Reynolds shear stress are rigorously derived for the first time for the far field of planar turbulent wakes under arbitrary pressure gradients. Specifically, by assuming self-similarity for the mean axial velocity, the analytical formulation for the mean transverse velocity is obtained from the integral of the mean continuity equation, and the analytical formulation for the Reynolds shear stress is obtained from the integral of the momentum equation. The generalized analytical formulations for the mean transverse velocity and Reynolds shear stress consist of multiple components, each with its unique scale and physical mechanism. In the zero pressure gradient limit, the generalized formulations recover the single-scale equations reported by Wei, Liu, and Livescu. Furthermore, simpler approximate formulations for the mean transverse velocity and Reynolds shear stress are also obtained, and show excellent agreement with the experimental measurements. Here, the findings provide new insights into the properties of planar turbulent wakes under pressure gradients, filling some long-standing gaps in the existing literature.

42 ENGINEERING↗

DOE OSTI2024

Nonlinear Poisson–Boltzmann solutions for charged parallel plates: When opposite charges repel

I present an exact solution of the Poisson–Boltzmann equation for two parallel plates and discuss the solution properties. I discuss in more detail plates with opposite charges: In this case, there are two critical separations, L c,1 < L c,2 . For separations less than L c,1 , the force between plates is repulsive. It switches to attractive at L c,1 , but with the electric potential having the same sign on both plates. For L > L c,2 , the force remains attractive, and the potential at the plates has the same sign as the charge on each plate. I also describe charge regulation, determined by pK a , and provide formulas for both the critical distance where oppositely charged plates repel and their charging process. Finally, the implications of these results for the nanoparticle assembly, as driven by electrostatic interactions, are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

X-point effects on the ideal MHD modes in tokamaks in the description of dual-poloidal-region safety factor

The flux coordinates with dual-region safety factor (q) in the poloidal direction are developed in this work. The X-point effects on the ideal MHD modes in tokamaks are then analyzed using this coordinate system. Since the X-point effects mainly affect the edge region, the modes localized at the tokamak edge are particularly examined. Two types of modes are studied. The first is related to the conventional peeling or peeling-ballooning modes. The mode existence aligned with the local magnetic field in the poloidally core region as observed experimentally is confirmed. The X points are shown to contribute to a stabilizing effect for the conventionally treated modes with the surface-averaged q and with the tokamak edge portion truncated. The other is the axisymmetric modes localized in the vicinity of X points, which can affect the cross-field-line transport near the X points. The existence of axisymmetric modes points to the possibility of applying a toroidally axisymmetric resonant magnetic perturbation (RMP) in the X-point area for mitigating the edge localized modes, which can be an alternative to the current RMP design. The dual q description also has important implications for the existing non-axisymmetric RMP concept. It helps to understand why the RMP suppression of edge localized modes is difficult to achieve in the double-null tokamak configurations and points to the possibility of further improving the current RMP concept by considering the alignment to the local q.

Fourier analysis↗

DOE OSTI2026

An ultra-fast method for generating synthetic down-scattered neutron data for inertial confinement fusion implosions

In inertial confinement fusion experiments at the National Ignition Facility, asymmetries are probed by a variety of neutron diagnostics, including neutron imaging systems, real-time neutron activation diagnostics (RTNADs), and neutron spectrometers. It is often useful to generate synthetic data based on these diagnostics to validate and tune models. However, current methods of doing so using Monte Carlo particle tracing are time-consuming. In this paper, an ultra-fast method is presented for generating synthetic neutron images, RTNAD data, and spectrometry data using line integrals and 3D convolutions. While it does not contain as much physics as particle tracing codes, it is thousands of times faster and produces nearly identical data. This enables analysis techniques that depend on generating large amounts of synthetic data, which will prove very useful for the study of asymmetries going forward.

Deuterium↗

NASA NTRS2009

Tactical Synthesis Of Efficient Global Search Algorithms

Algorithm synthesis transforms a formal specification into an efficient algorithm to solve a problem. Algorithm synthesis in Specware combines the formal specification of a problem with a high-level algorithm strategy. To derive an efficient algorithm, a developer must define operators that refine the algorithm by combining the generic operators in the algorithm with the details of the problem specification. This derivation requires skill and a deep understanding of the problem and the algorithmic strategy. In this paper we introduce two tactics to ease this process. The tactics serve a similar purpose to tactics used for determining indefinite integrals in calculus, that is suggesting possible ways to attack the problem.

Nedunuri, Srinivas↗

NASA NTRS2002

Trajectories for High Specific Impulse High Specific Power Deep Space Exploration

Flight times and deliverable masses for electric and fusion propulsion systems are difficult to approximate. Numerical integration is required for these continuous thrust systems. Many scientists are not equipped with the tools and expertise to conduct interplanetary and interstellar trajectory analysis for their concepts. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. An analytical method derived in the companion paper was also evaluated. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

NASA NTRS2002

Simulation of Trajectories for High Specific Impulse Deep Space Exploration

Difficulties in approximating flight times and deliverable masses for continuous thrust propulsion systems have complicated comparison and evaluation of proposed propulsion concepts. These continuous thrust propulsion systems are of interest to many groups, not the least of which are the electric propulsion and fusion communities. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. The analytical method derived in the companion paper was also used to simulate the trajectory. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

DOE OSTI2022

Nonlocal Kernel Network (NKN): a Stable and Resolution-Independent Deep Neural Network.

Neural operators have recently become popular tools for designing solution maps between function spaces in the form of neural networks. Differently from classical scientific machine learning approaches that learn parameters of a known partial differential equation (PDE) for a single instance of the input parameters at a fixed resolution, neural operators approximate the solution map of a family of PDEs [6, 7]. Despite their success, the uses of neural operators are so far restricted to relatively shallow neural networks and confined to learning hidden governing laws. In this work, we propose a novel nonlocal neural operator, which we refer to as nonlocal kernel network (NKN), that is resolution independent, characterized by deep neural networks, and capable of handling a variety of tasks such as learning governing equations and classifying images. Our NKN stems from the interpretation of the neural network as a discrete nonlocal diffusion reaction equation that, in the limit of infinite layers, is equivalent to a parabolic nonlocal equation, whose stability is analyzed via nonlocal vector calculus. The resemblance with integral forms of neural operators allows NKNs to capture long-range dependencies in the feature space, while the continuous treatment of node-to-node interactions makes NKNs resolution independent. The resemblance with neural ODEs, reinterpreted in a nonlocal sense, and the stable network dynamics between layers allow for generalization of NKN’s optimal parameters from shallow to deep networks. This fact enables the use of shallow-to-deep initialization techniques [8]. Our tests show that NKNs outperform baseline methods in both learning governing equations and image classification tasks and generalize well to different resolutions and depths.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Theoretical Analysis of Fractional Viscoelastic Flow in Circular Pipes: General Solutions

Fractional calculus is a relatively old yet emerging field of mathematics with the widest range of engineering and biomedical applications. Despite being an incredibly powerful tool, it, however, requires promotion in the engineering community. Rheology is undoubtedly one of the fields where fractional calculus has become an integral part of cutting-edge research. There exists extensive literature on the theoretical, experimental, and numerical treatment of various fractional viscoelastic flows in constraint geometries. However, the general theoretical approach that unites several most commonly used models is missing. Here we present exact analytical solutions for fractional viscoelastic flow in a circular pipe. We find velocity profiles and shear stresses for fractional Maxwell, Kelvin–Voigt, Zener, Poynting–Thomson, and Burgers models. The dynamics of these quantities are studied with respect to normalized pipe radius, fractional orders, and elastic moduli ratio. Three different types of behavior are identified: monotonic increase, resonant, and aperiodic oscillations. The models developed are applicable in the widest material range and allow for the alteration of the balance between viscous and elastic properties of the materials.

Gritsenko, Dmitry (ORCID:0000000248181512)↗