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

Stationary self-consistent distributions for a charged particle beam in the longitudinal magnetic field

A review of analytical solutions of the Vlasov equation for a beam of charged particles is given. These results are analyzed on the basis of a unified approach developed by the authors. In the context of this method, a space of integrals of motion is introduced in which the integrals of motion of particles are considered as coordinates. In this case, specifying a self-consistent distribution is reduced to defining a distribution density in this space. This approach allows us to simplify the construction and analysis of different self-consistent distributions. In particular, it is possible, in some cases, to derive new solutions by considering linear combinations of well-known solutions. This approach also makes it possible in many cases to give a visual geometric representation of self-consistent distributions in the space of integrals of motion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Charged particle tracking through electrostatic wire meshes using the finite element method

Wire meshes are used across many disciplines to accelerate and focus charged particles, however, analytical solutions are non-exact and few codes exist which simulate the exact fields around a mesh with physical sizes. A tracking code based in Matlab-Simulink using field maps generated using finite element software has been developed which tracks electrons or ions through electrostatic wire meshes. The fields around such a geometry are presented as an analytical expression using several basic assumptions, however, it is apparent that computational calculations are required to obtain realistic values of electric potential and fields, particularly when multiple wire meshes are deployed. The tracking code is flexible in that any quantitatively describable particle distribution can be used for both electrons and ions as well as other benefits such as ease of export to other programs for analysis. The code is made freely available and physical examples are highlighted where this code could be beneficial for different applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Numerical studies of nonlinear ultrasonic guided waves in uniform waveguides with arbitrary cross sections

Nonlinear guided waves have been investigated widely in simple geometries, such as plates, pipe and shells, where analytical solutions have been developed. This paper extends the application of nonlinear guided waves to waveguides with arbitrary cross sections. The criteria for the existence of nonlinear guided waves were summarized based on the finite deformation theory and nonlinear material properties. Numerical models were developed for the analysis of nonlinear guided waves in complex geometries, including nonlinear Semi-Analytical Finite Element (SAFE) method to identify internal resonant modes in complex waveguides, and Finite Element (FE) models to simulate the nonlinear wave propagation at resonant frequencies. Two examples, an aluminum plate and a steel rectangular bar, were studied using the proposed numerical model, demonstrating the existence of nonlinear guided waves in such structures and the energy transfer from primary to secondary modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Time-evolution of quantum systems via a complex nonlinear Riccati equation. I. Conservative systems with time-independent Hamiltonian

The sensitivity of the evolution of quantum uncertainties to the choice of the initial conditions is shown via a complex nonlinear Riccati equation leading to a reformulation of quantum dynamics. This sensitivity is demonstrated for systems with exact analytic solutions with the form of Gaussian wave packets. In particular, one-dimensional conservative systems with at most quadratic Hamiltonians are studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Are stealth scalar fields stable?

Nongravitating (stealth) scalar fields associated with Minkowski space in scalar-tensor gravity are examined. Analytical solutions for both nonminimally coupled scalar field theory and for Brans-Dicke gravity are studied and their stability with respect to tensor perturbations is assessed using a covariant and gauge-invariant formalism developed for alternative gravity. For Brans-Dicke solutions, the stability with respect to homogeneous perturbations is also studied. There are regions of parameter space corresponding to stability and other regions corresponding to instability.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

Efficient Modeling of Laser-Plasma Accelerators with INF&RNO

The numerical modeling code INF&RNO (INtegrated Fluid& paRticle simulatioN cOde, pronounced"inferno") is presented. INF&RNO is an efficient 2D cylindrical code to model the interaction of a short laser pulse with an underdense plasma. The code is based on an envelope model for the laser while either a PIC or a fluid description can be used for the plasma. The effect of the laser pulse on the plasma is modeled with the time-averaged poderomotive force. These and other features allow for a speedup of 2-4 orders of magnitude compared to standard full PIC simulations while still retaining physical fidelity. The code has been benchmarked against analytical solutions and 3D PIC simulations and here a set of validation tests together with a discussion of the performances are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Phase-space finite elements in a least-squares solution of the transport equation

The linear Boltzmann transport equation is solved using a least-squares finite element approximation in the space, angular and energy phase-space variables. The method is applied to both neutral particle transport and also to charged particle transport in the presence of an electric field, where the angular and energy derivative terms are handled with the energy/angular finite elements approximation, in a manner analogous to the way the spatial streaming term is handled. For multi-dimensional problems, a novel approach is used for the angular finite elements: mapping the surface of a unit sphere to a two-dimensional planar region and using a meshing tool to generate a mesh. In this manner, much of the spatial finite-elements machinery can be easily adapted to handle the angular variable. The energy variable and the angular variable for one-dimensional problems make use of edge/beam elements, also building upon the spatial finite elements capabilities. The methods described here can make use of either continuous or discontinuous finite elements in space, angle and/or energy, with the use of continuous finite elements resulting in a smaller problem size and the use of discontinuous finite elements resulting in more accurate solutions for certain types of problems. The work described in this paper makes use of continuous finite elements, so that the resulting linear system is symmetric positive definite and can be solved with a highly efficient parallel preconditioned conjugate gradients algorithm. The phase-space finite elements capability has been built into the Sceptre code and applied to several test problems, including a simple one-dimensional problem with an analytic solution available, a two-dimensional problem with an isolated source term, showing how the method essentially eliminates ray effects encountered with discrete ordinates, and a simple one-dimensional charged-particle transport problem in the presence of an electric field. (authors)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Existence and classification of pseudo-asymptotic solutions for Tolman–Oppenheimer–Volkoff systems

The Tolman–Oppenheimer–Volkoff (TOV) equations are a partially uncoupled system of nonlinear and non-autonomous ordinary differential equations which describe the structure of isotropic spherically symmetric static fluids. Nonlinearity makes finding explicit solutions of TOV systems very difficult and such solutions are very rare. In this paper we introduce the notion of pseudo-asymptotic TOV systems and we show that the space of such systems is at least fifteen-dimensional. We also show that if the system is defined in a suitable domain (meaning the extended real line), then well-behaved pseudo-asymptotic TOV systems are genuine TOV systems in that domain, ensuring the existence of new fourteen analytic solutions for extended TOV equations. The pseudo-asymptotic solutions are classified according to the nature of the matter (ordinary or exotic) and to the existence of cavities and singularities. It is shown that at least three of them are realistic, in the sense that they are formed only by ordinary matter and contain no cavities or singularities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Stationary point analysis of the one-dimensional lattice Landau gauge fixing functional, aka random phase XY Hamiltonian

Research Highlights: > Exact results for all stationary points of some high-dimensional function are given. > They are interpreted as Gribov copies of a lattice Landau gauge fixing functional. > The Gribov ambiguity and the Neuberger problem in compact U(1) are illustrated. > Stationary points are used to discuss a criterion on the absence of phase transitions. - Abstract: We study the stationary points of what is known as the lattice Landau gauge fixing functional in one-dimensional compact U(1) lattice gauge theory, or as the Hamiltonian of the one-dimensional random phase XY model in statistical physics. An analytic solution of all stationary points is derived for lattices with an odd number of lattice sites and periodic boundary conditions. In the context of lattice gauge theory, these stationary points and their indices are used to compute the gauge fixing partition function, making reference in particular to the Neuberger problem. Interpreted as stationary points of the one-dimensional XY Hamiltonian, the solutions and their Hessian determinants allow us to evaluate a criterion which makes predictions on the existence of phase transitions and the corresponding critical energies in the thermodynamic limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Propagation of ion acoustic shock waves in negative ion plasmas with nonextensive electrons

Nonlinear ion acoustic shocks (monotonic as well as oscillatory) waves in negative ion plasmas are investigated. The inertialess electron species are assumed to be nonthermal and follow Tsallis distribution. The dissipation in the plasma is considered via kinematic viscosities of both positive and negative ion species. The Korteweg-de Vries Burgers (KdVB) equation is derived using small amplitude reductive perturbation technique and its analytical solution is presented. The effects of variation of density and temperature of negative ions and nonthermal parameter q of electrons on the strength of the shock structures are plotted for illustration. The numerical solutions of KdVB equation using Runge Kutta method are obtained, and transition from oscillatory to monotonic shock structures is also discussed in detail for negative ions nonthermal plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Generalized Bohm’s criterion and negative anode voltage fall in electric discharges

The value of the voltage fall across the anode sheath is found as a function of the current density. Analytic solutions are obtained in a wide range of the ratio of the directed velocity of plasma electrons v{sub 0} to their thermal velocity v{sub T}. It is shown that the voltage fall in a one-dimensional collisionless anode sheath is always negative. At the small values of v{sub 0}/v{sub T}, the obtained expression asymptotically transforms into the Langmuir formula. Generalized Bohm’s criterion for an electric discharge with allowance for the space charge density ρ(0), electric field E(0), ion velocity v{sub i}(0), and ratio v{sub 0}/v{sub T} at the plasma-sheath interface is formulated. It is shown that the minimum value of the ion velocity v{sub i}{sup *}(0) corresponds to the vanishing of the electric field at one point inside the sheath. The dependence of v{sub i}{sup *} (0) on ρ(0), E(0), and v{sub 0}/v{sub T} determines the boundary of the existence domain of stationary solutions in the sheath. Using this criterion, the maximum possible degree of contraction of the electron current at the anode is determined for a short high-current vacuum arc discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Study of tethered satellite active attitude control

Existing software was adapted for the study of tethered subsatellite rotational dynamics, an analytic solution for a stable configuration of a tethered subsatellite was developed, the analytic and numerical integrator (computer) solutions for this "test case' was compared in a two mass tether model program (DUMBEL), the existing multiple mass tether model (SKYHOOK) was modified to include subsatellite rotational dynamics, the analytic "test case,' was verified, and the use of the SKYHOOK rotational dynamics capability with a computer run showing the effect of a single off axis thruster on the behavior of the subsatellite was demonstrated. Subroutines for specific attitude control systems are developed and applied to the study of the behavior of the tethered subsatellite under realistic on orbit conditions. The effect of all tether "inputs,' including pendular oscillations, air drag, and electrodynamic interactions, on the dynamic behavior of the tether are included.

Colombo, G.↗
DOE OSTI2014

AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A Reduced-Boundary-Function Method for Convective Heat Transfer With Axial Heat Conduction and Viscous Dissipation

We introduce a new method of solution for the convective heat transfer under forced laminar flow that is confined by two parallel plates with a distance of 2a or by a circular tube with a radius of a. The advection-conduction equation is first mapped onto the boundary. The original problem of solving the unknown field T(x,r,t) is reduced to seek the solutions of T at the boundary (r = a or r = 0, r is the distance from the centerline shown in Fig. 1), i.e., the boundary functions T{sub a}(x,t) {triple_bond} T(x,r=a,t) and/or T{sub 0}(x,t) {triple_bond} T(x,r=0,t). In this manner, the original problem is significantly simplified by reducing the problem dimensionality from 3 to 2. The unknown field T(x,r,t) can be eventually solved in terms of these boundary functions. The method is applied to the convective heat transfer with uniform wall temperature boundary condition and with heat exchange between flowing fluids and its surroundings that is relevant to the geothermal applications. Analytical solutions are presented and validated for the steady-state problem using the proposed method.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2017

Spectra for the reemission of attosecond and shorter electromagnetic pulses by multielectron atoms

Based on the analytical solution of the Schrödinger equation, we have considered the reemission of attosecond and shorter electromagnetic pulses by multielectron atoms in the sudden perturbation approximation. We have developed a technique of calculating the spectra for the reemission of attosecond and shorter electromagnetic pulses by neutral multielectron atoms with nuclear charges from 1 to 92. The results are presented in the form of analytical formulas dependent on several coefficients and screening parameters tabulated for all of the atoms whose electron densities are described by the well-known Dirac–Hartree–Fock–Slater model. As examples we have calculated the spectra for the reemission by lithium, carbon, calcium, and iron atoms for two types of incident pulse: Gaussian and “sombrero.”.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Quantum theory of space charge limited current in solids

We present a quantum model of space charge limited current transport inside trap-free solids with planar geometry in the mean field approximation. We use a simple transformation which allows us to find the exact analytical solution for the steady state current case. We use our approach to find a Mott-Gurney like behavior and the mobility for single charge carriers in the quantum regime in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Potentials of the Heun class: The triconfluent case

We study special classes of potentials for which the one-dimensional (or radial) Schrödinger equation can be reduced to a triconfluent Heun equation by a suitable coordinate transformation together with an additional transformation of the wave function. In particular, we analyze the behaviour of those subclasses of the potential arising when the ordinary differential equation governing the coordinate transformation admits explicit analytic solutions in terms of the radial variable. Furthermore, we obtain formulae for solutions of the eigenvalue problem of the associated radial Schrödinger operator. Last but not least, using methods of supersymmetric quantum mechanics we relate the considered potentials to a new class of exactly solvable ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Theory of anisotropic whispering-gallery-mode resonators

An analytic solution for a uniaxial spherical resonator is presented using the method of Debye potentials. This serves as a starting point for the calculation of whispering gallery modes (WGMs) in such a resonator. Suitable approximations for the radial functions are discussed in order to best characterize WGMs. The characteristic equation and its asymptotic expansion for the anisotropic case is also discussed, and an analytic formula with a precision of the order O[{nu}{sup -1}] is also given. Our careful treatment of both boundary conditions and asymptotic expansions makes the present work a particularly suitable platform for a quantum theory of whispering gallery resonators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗