Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytic solution”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Solution of D dimensional Dirac equation for coulombic potential using NU method and its thermodynamics properties

The analytical solution of Ddimensional Dirac equation for Coulombic potential is investigated using Nikiforov-Uvarov method. The D dimensional relativistic energy spectra are obtained from relativistic energy eigenvalue equation by using Mat Lab software.The corresponding D dimensional radial wave functions are formulated in the form of generalized Jacobi and Laguerre Polynomials. In the non-relativistic limit, the relativistic energy equation reduces to the non-relativistic energy which will be applied to determine some thermodynamical properties of the system. The thermodynamical properties of the system are expressed in terms of error function and imaginary error function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Ultrasonic guided waves in eccentric annular pipes

This paper studies the feasibility of using ultrasonic guided waves to rapidly inspect tubes and pipes for possible eccentricity. While guided waves are well established in the long range inspection of structures such as pipes and plates, studies for more complex cross sections are limited and analytical solutions are often difficult to obtain. Recent developments have made the Semi Analytical Finite Element (SAFE) method widely accessible for researchers to study guided wave properties in complex structures. Here the SAFE method is used to study the effect of eccentricity on the modal structures and velocities of lower order guided wave modes in thin pipes of diameters typically of interest to the industry. Results are validated using experiments. The paper demonstrates that even a small eccentricity in the pipe can strongly affect guided wave mode structures and velocities and hence shows potential for pipe eccentricity inspection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2011

Experimental and modeling analysis of fast ionization wave discharge propagation in a rectangular geometry

Fast ionization wave (FIW), nanosecond pulse discharge propagation in nitrogen and helium in a rectangular geometry channel/waveguide is studied experimentally using calibrated capacitive probe measurements. The repetitive nanosecond pulse discharge in the channel was generated using a custom designed pulsed plasma generator (peak voltage 10-40 kV, pulse duration 30-100 ns, and voltage rise time {approx}1 kV/ns), generating a sequence of alternating polarity high-voltage pulses at a pulse repetition rate of 20 Hz. Both negative polarity and positive polarity ionization waves have been studied. Ionization wave speed, as well as time-resolved potential distributions and axial electric field distributions in the propagating discharge are inferred from the capacitive probe data. ICCD images show that at the present conditions the FIW discharge in helium is diffuse and volume-filling, while in nitrogen the discharge propagates along the walls of the channel. FIW discharge propagation has been analyzed numerically using quasi-one-dimensional and two-dimensional kinetic models in a hydrodynamic (drift-diffusion), local ionization approximation. The wave speed and the electric field distribution in the wave front predicted by the model are in good agreement with the experimental results. A self-similar analytic solution of the fast ionization wave propagation equations has also been obtained. The analytic model of the FIW discharge predicts key ionization wave parameters, such as wave speed, peak electric field in the front, potential difference across the wave, and electron density as functions of the waveform on the high voltage electrode, in good agreement with the numerical calculations and the experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Ion acoustic shock waves in plasmas with warm ions and kappa distributed electrons and positrons

The monotonic and oscillatory ion acoustic shock waves are investigated in electron-positron-ion plasmas (e-p-i) with warm ions (adiabatically heated) and nonthermal kappa distributed electrons and positrons. The dissipation effects are included in the model due to kinematic viscosity of the ions. Using reductive perturbation technique, the Kadomtsev-Petviashvili-Burgers (KPB) equation is derived containing dispersion, dissipation, and diffraction effects (due to perturbation in the transverse direction) in e-p-i plasmas. The analytical solution of KPB equation is obtained by employing tangent hyperbolic (Tanh) method. The analytical condition for the propagation of oscillatory and monotonic shock structures are also discussed in detail. The numerical results of two dimensional monotonic shock structures are obtained for graphical representation. The dependence of shock structures on positron equilibrium density, ion temperature, nonthermal spectral index kappa, and the kinematic viscosity of ions are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Finite-element analysis of the deformation of thin Mylar films due to measurement forces

Significant deformation of thin films occurs when measuring thickness by mechanical means. This source of measurement error can lead to underestimating film thickness if proper corrections are not made. Analytical solutions exist for Hertzian contact deformation, but these solutions assume relatively large geometries. If the film being measured is thin, the analytical Hertzian assumptions are not appropriate. ANSYS is used to model the contact deformation of a 48 gauge Mylar film under bearing load, supported by a stiffer material. Simulation results are presented and compared to other correction estimates. Ideal, semi-infinite, and constrained properties of the film and the measurement tools are considered.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Approximation solution of Schrodinger equation for Q-deformed Rosen-Morse using supersymmetry quantum mechanics (SUSY QM)

The approximate analytical solution of Schrodinger equation for Q-Deformed Rosen-Morse potential was investigated using Supersymmetry Quantum Mechanics (SUSY QM) method. The approximate bound state energy is given in the closed form and the corresponding approximate wave function for arbitrary l-state given for ground state wave function. The first excited state obtained using upper operator and ground state wave function. The special case is given for the ground state in various number of q. The existence of Rosen-Morse potential reduce energy spectra of system. The larger value of q, the smaller energy spectra of system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An explicitly solvable energy-conserving algorithm for pitch-angle scattering in magnetized plasmas

In this work, we develop an Explicitly Solvable Energy-Conserving (ESEC) algorithm for the Stochastic Differential Equation (SDE) describing the pitch-angle scattering process in magnetized plasmas. The Cayley transform is used to calculate both the deterministic gyromotion and stochastic scattering, affording the algorithm to be explicitly solvable and exactly energy conserving. An unusual property of the SDE for pitch-angle scattering is that its coefficients diverge at the zero velocity and do not satisfy the global Lipschitz condition. Consequently, when standard numerical methods, such as the Euler-Maruyama (EM), are applied, numerical convergence is difficult to establish. For the proposed ESEC algorithm, its energy-preserving property enables us to overcome this obstacle. We rigorously prove that the ESEC algorithm is order 1/2 strongly convergent. This result is confirmed by detailed numerical studies. For the case of pitch-angle scattering in a magnetized plasma with a constant magnetic field, the numerical solution is benchmarked against the analytical solution, and excellent agreements are found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Blobs in the framework of drift wave dynamics

Here, analytic solutions for blob-like structures obtained in S. I. Krasheninnikov [Phys. Lett. A 380, 3905 (2016)], are examined in detail. Numerical solutions for the cases where analytic consideration is not possible are provided

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Correlation of turbulent trailing vortex decay data

A correlation function, derived on the basis of self similar variable eddy viscosity decay, is introduced and utilized to correlate aircraft trailing vortex velocity data from ground and flight experiments. The correlation function collapses maximum tangential velocity data from scale model and flight tests to a single curve. The resulting curve clearly shows both the inviscid plateau and the downstream decay regions. A comparison between experimental data and numerical solution shows closer agreement with the variable eddy viscosity solution than the constant viscosity analytical solution.

Iversen, J. D.↗
NASA NTRS1976

Correlation of turbulent trailing vortex decay data

A correlation function, derived on the basis of self-similar variable eddy-viscosity decay, is introduced and utilized to correlate aircraft trailing vortex velocity data from ground and flight experiments. The correlation function collapses maximum tangential velocity data from scale-model and flight tests to a single curve. The resulting curve clearly shows both the inviscid plateau and the downstream decay regions. A comparison between experimental data and numerical solution shows closer agreement with the variable eddy viscosity solution than the constant viscosity analytical solution.

Iversen, J. D.↗
NASA NTRS1985

Opening of the magnetic field lines in a fast rotating magnetosphere, with an application to Jupiter

A simple model of a magnetosphere rotating around a rotating Jupiter-like object with a spin-aligned dipolar moment mu is considered. A low-energy plasma is released by inner sources located beyond the corotation radius and is diffused outward through closed lines. The lines form a closed thin equatorial disk, maintaining a quasi-static balance between the centrifugal force and the magnetic tension. Solutions are obtained for the analytical equations characterizing the plasma in the closed region; the structure of the magnetic field outside the disk; and the self-consistent value of the critical radius r sub zero. The analytical solutions are compared to the observed properties of the Jovian magnetosphere, and the results are discussed in detail.

Aly, J. J.↗
NASA NTRS1991

Hybrid accretion disks in active galactic nuclei. I - Structure and spectra

A unified treatment is presented of the two distinct states of vertically thin AGN accretion disks: a cool (about 10 to the 6th K) optically thick solution, and a hot (about 10 to the 9th K) optically thin solution. A generalized formalism and a new radiative cooling equation valid in both regimes are introduced. A new luminosity limit is found at which the hot and cool alpha solutions merge into a single solution of intermediate optical depth. Analytic solutions for the disk structure are given, and output spectra are computed numerically. This is used to demonstrate the prospect of fitting AGN broadband spectra containing both the UV bump as well as the hard X-ray and gamma-ray tail, using a single accretion disk model. Such models are found to make definite predictions about the observed spectrum, such as the relation between the hard X-ray spectral index, the UV-to-X-ray luminosity ratio, and a feature of about 1 MeV.

Wandel, Amri↗
NASA NTRS2000

Prediction of Sound Waves Propagating Through a Nozzle Without/With a Shock Wave Using the Space-Time CE/SE Method

The benchmark problems in Category 1 (Internal Propagation) of the third Computational Aeroacoustics (CAA) Work-shop sponsored by NASA Glenn Research Center are solved using the space-time conservation element and solution element (CE/SE) method. The first problem addresses the propagation of sound waves through a nearly choked transonic nozzle. The second one concerns shock-sound interaction in a supersonic nozzle. A quasi one-dimension CE/SE Euler solver for a nonuniform mesh is developed and employed to solve both problems. Numerical solutions are compared with the analytical solution for both problems. It is demonstrated that the CE/SE method is capable of solving aeroacoustic problems with/without shock waves in a simple way. Furthermore, the simple nonreflecting boundary condition used in the CE/SE method which is not based on the characteristic theory works very well.

Wang, Xiao-Yen↗
NASA NTRSDate not supplied

Copernicus-LinCov (COPCOV) Software Integration in Support of Robust Trajectory Optimization

Robust trajectory optimization is the process of optimizing a trajectory while accounting for system uncertainty due to a variety of potential error sources. This work highlights the development and features of a novel tool known as CopCov to support robust trajectory optimization efforts. CopCov acts as an interface between Copernicus, a generalized trajectory design and optimization tool, and LinCov, a linear covariance analysis tool. By having a direct interface between these two software packages, Copernicus can receive covariance information from LinCov through a direct feedback loop, thus enabling optimization of a trajectory that is robust to trajectory dispersions and navigation errors. This paper details the architecture of CopCov and its flexibility to operate under varying configurations, including with both tools running locally or alternatively with the tools communicating via a remote connection. Additionally, the CopCov tool is demonstrated on a simple Hohmann transfer reference trajectory with varying numbers of Trajectory Correction Maneuvers (TCMs) and varying problem formulations. This example scenario is used to highlight how the inclusion of the CopCov interface affects burn placement of both major burns and minor burns (i.e., TCMs) in the optimized solution. Results are compared against analytical solutions and against a Genetic Algorithm (GA) optimizer for independent verification and validation.

Copernicus↗
DOE OSTI2022

Investigations about iso-geometric analysis for self-shielding calculations with the subgroup method

The implementation of a self-shielding method for a neutron transport calculation code based on the iso-geometric analysis (IGA) method that can resolve the multi-group neutron transport equation for arbitrary spatial domain, is presented. The self-shielding model based on the subgroup theory is adopted because the subgroup method can be used to perform calculations for any arbitrary geometrical domain which is the main purpose of our IGA code. Some basic theory of the subgroup method is given. A self-shielding calculation based on the PWR fuel pin composed of MOX fuel is presented and compared with a Monte Carlo calculation. The result is that the combination of SN transport theory, IGA and the subgroup method gives correctly shielded cross sections. Validation of the combination of the IGA method and SN neutron transport for subgroup calculations in two steps are presented: the first step is to ascertain the correctness of the IGA solutions compared to a known analytical solution in diffusion theory; the second step is a validation of the IGA solutions compared to a reference calculation in SN transport theory. The relation between the settings for the IGA solver and the accuracy of the results are elucidated. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗