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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2010

Snakelike nonautonomous solitons in a graded-index grating waveguide

We present a series of analytical solutions which describe nonautonomous solitons in a planar waveguide with an additional periodical structure, that is, a long-period grating. The explicit functions which describe the evolution of the width, peak, and trajectory of the soliton's wave center are presented exactly. The gain parameter has no effects on the motion of the soliton's wave center or its width; it affects just the evolution of the soliton's peak. The grating term affects the motion of the soliton's wave center without changing its shape. The evolution of the soliton under the propagation-distance-dependent gain term is investigated too. It is reported that an arbitrary additional structure can be added on the graded-index waveguide to control the motion of the soliton without affecting its shape.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Peculiarities of elastic scattering of intermediate-energy electrons related to their capture at quasi-stationary levels of an atom in a continuous spectrum

A theoretical model of electron scattering on an atom is constructed to study elastic atomic scattering of intermediate-energy electrons. The proposed model is based upon the combined Mensing potential with two spheres of atomic electrons, which admits analytical solutions of the radial Schröbinger equation. A procedure for matching the parameters of this scatterer to an approximate electrostatic potential of an atom in the form of a screened Coulomb potential has been determined. The screening radius of the latter potential has been calculated proceeding from the properties corresponding to the Thomas-Fermi method. A model of a scatterer determined according to the aforementioned procedure can be used to calculate the energy dependence of the cross section of elastic electron scattering on some atoms with s, p, and d shells representing elements neighboring zirconium. The main result is the establishment of factors responsible for the appearance of maxima on the energy dependences of the cross section of elastic electron scattering. These maxima are related to the resonant trapping of impinging electrons by quasi-stationary levels in a continuous spectrum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

An Investigation of Machining Characteristics in Micro-scale Milling Process

In this paper, an analytical solution of micro-scale milling process is presented in order to suggest available machining conditions. The size effect should be considered to determine cutting characteristics in micro-scale cutting. The feed per tooth is the most dominant cutting parameter related to the size effect in micro-scale milling process. In order to determine the feed per tooth at which chips can be formed, the finite element method is used. The finite element method is employed by utilizing the Johnson-Cook (JC) model as a constitutive model of work material flow stress. Machining experiments are performed to validate the simulation results by using a micro-machining stage. The validation is conducted by observing cutting force signals from a cutting tool and the conditions of the machined surface of the workpiece.

36 MATERIALS SCIENCE↗
DOE OSTI2011

Kadomtsev-Petviashvili solitons propagation in a plasma system with superthermal and weakly relativistic effects

Two dimensional (2D) solitons are studied in a plasma system comprising of relativistically streaming ions, kappa distributed electrons, and positrons. Kadomtsev-Petviashvili (KP) equation is derived through the reductive perturbation technique. Analytical solution of the KP equation has been studied numerically and graphically. It is noticed that kappa parameters of electrons and positrons as well as the ions relativistic streaming factor have an emphatic influence on the structural as well as propagation characteristics of two dimensional solitons in the considered plasma system. Our results may be helpful in the understanding of soliton propagation in astrophysical and laboratory plasmas, specifically the interaction of pulsar relativistic wind with supernova ejecta and the transfer of energy to plasma by intense electric field of laser beams producing highly energetic superthermal and relativistic particles [L. Arons, Astrophys. Space Sci. Lib. 357, 373 (2009); P. Blasi and E. Amato, Astrophys. Space Sci. Proc. 2011, 623; and A. Shah and R. Saeed, Plasma Phys. Controlled Fusion 53, 095006 (2011)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Role of the heat accumulation effect in the multipulse modes of the femtosecond laser microstructuring of silicon

The results of quantitative evaluation of the heat accumulation effect during the femtosecond laser microstructuring of the surface of silicon are presented for discussion. In the calculations, the numerical–analytical method is used, in which the dynamics of electronic processes and lattice heating are simulated by the numerical method, and the cooling stage is described on the basis of an analytical solution. The effect of multipulse irradiation on the surface temperature is studied: in the electronic subsystem, as the dependence of the absorbance on the excited carrier density and the dependence of the absorbance on the electron-gas temperature; in the lattice subsystem, as the variation in the absorbance from pulse to pulse. It was shown that, in the low-frequency pulse-repetition mode characteristic of the femtosecond microstructuring of silicon, the heat accumulation effect is controlled not by the residual surface temperature by the time of the next pulse arrival, which corresponds to conventional concepts, but by an increase in the maximum temperature from pulse to pulse, from which cooling begins. The accumulation of the residual temperature of the surface can affect the microstructuring process during irradiation near the evaporation threshold or with increasing pulse-repetition rate.

36 MATERIALS SCIENCE↗
DOE OSTI2015

JET LUMINOSITY OF GAMMA-RAY BURSTS: THE BLANDFORD–ZNAJEK MECHANISM VERSUS THE NEUTRINO ANNIHILATION PROCESS

A neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH) is one of the plausible candidates for the central engine of gamma-ray bursts (GRBs). Two mechanisms, i.e., the Blandford–Znajek (BZ) mechanism and the neutrino annihilation process, are generally considered to power GRBs. Using the analytic solutions from Xue et al. and ignoring the effects of the magnetic field configuration, we estimate the BZ and neutrino annihilation luminosities as functions of the disk masses and BH spin parameters to contrast the observational jet luminosities of GRBs. Our results show that although the neutrino annihilation processes could account for most GRBs, the BZ mechanism is more effective, especially for long-duration GRBs. Actually, if the energy of the afterglows and flares of GRBs is included, then the distinction between these two mechanisms is more significant. Furthermore, massive disk mass and high BH spin are beneficial for powering the high luminosities of GRBs. Finally, we discuss possible physical mechanisms that could enhance the disk mass or neutrino emission rate of NDAFs and the relevant difference between these two mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Spontaneous symmetry breaking in Schroedinger lattices with two nonlinear sites

We introduce discrete systems in the form of straight (infinite) and ring-shaped chains, with two symmetrically placed nonlinear sites. The systems can be implemented in nonlinear optics (as waveguiding arrays) and Bose-Einstein condensates (by means of an optical lattice). A full set of exact analytical solutions for symmetric, asymmetric, and antisymmetric localized modes is found, and their stability is investigated in a numerical form. The symmetry-breaking bifurcation, through which the asymmetric modes emerge from the symmetric ones, is found to be of the subcritical type. It is transformed into a supercritical bifurcation if the nonlinearity is localized in relatively broad domains around two central sites, and also in the ring of a small size, i.e., in effectively nonlocal settings. The family of antisymmetric modes does not undergo bifurcations and features both stable and unstable portions. The evolution of unstable localized modes is investigated by means of direct simulations. In particular, unstable asymmetric states, which exist in the case of the subcritical bifurcation, give rise to breathers oscillating between the nonlinear sites, thus restoring an effective dynamical symmetry between them.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Quantum information entropies for position-dependent mass Schrödinger problem

The Shannon entropy for the position-dependent Schrödinger equation for a particle with a nonuniform solitonic mass density is evaluated in the case of a trivial null potential. The position S{sub x} and momentum S{sub p} information entropies for the three lowest-lying states are calculated. In particular, for these states, we are able to derive analytical solutions for the S{sub x} entropy as well as for the Fourier transformed wave functions, while the S{sub p} quantity is calculated numerically. We notice the behavior of the S{sub x} entropy, namely, it decreases as the mass barrier width narrows and becomes negative beyond a particular width. The negative Shannon entropy exists for the probability densities that are highly localized. The mass barrier determines the stability of the system. The dependence of S{sub p} on the width is contrary to the one for S{sub x}. Some interesting features of the information entropy densities ρ{sub s}(x) and ρ{sub s}(p) are demonstrated. In addition, the Bialynicki-Birula–Mycielski (BBM) inequality is tested for a number of states and found to hold for all the cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Test particle motion in the space-time of a Kerr black hole pierced by a cosmic string

We study the geodesic equation in the space-time of a Kerr black hole pierced by an infinitely thin cosmic string and give the complete set of analytical solutions of this equation for massive and massless particles in terms of Mino time that allows one to decouple the r and {theta} components of the geodesic equation. The solutions of the geodesic equation can be classified according to the particle's energy and angular momentum, the mass and angular momentum per mass of the black hole. We give examples of orbits showing the influence of the cosmic string. We also discuss the perihelion shift and the Lense-Thirring effect for bound orbits and show that the presence of a cosmic string enhances both effects. Comparing our results with experimental data from the LAGEOS satellites we find an upper bound on the energy per unit length of a string piercing the earth which is approximately 10{sup 16} kg/m. Our work has also applications to the recently suggested explanation of the alignment of the polarization vector of quasars using remnants of cosmic string decay in the form of primordial magnetic field loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Semi-analytical prediction of hydraulic resistance and heat transfer for pipe and channel flows of water at supercritical pressure

Within the Generation IV International Forum the Supercritical Water Reactor is investigated. For its core design and safety analysis the efficient prediction of flow and heat transfer parameters such as the wall-shear stress and the heat-transfer coefficient for pipe and channel flows is needed. For circular pipe flows a numerical model based on the one-dimensional conservation equations of mass, momentum end energy in the radial direction is presented, referred to as a 'semi-analytical' method. An accurate, high-order numerical method is employed to evaluate previously derived analytical solutions of the governing equations. Flow turbulence is modeled using the algebraic approach of Prandtl/van-Karman, including a model for the buffer layer. The influence of wall roughness is taken into account by a new modified numerical damping function of the turbulence model. The thermo-hydraulic properties of water are implemented according to the international standard of 1997. This method has the potential to be used within a sub-channel analysis code and as wall-functions for CFD codes to predict the wall shear stress and the wall temperature. The present study presents a validation of the method with comparison of model results with experiments and multi-dimensional computational (CFD) studies in a wide range of flow parameters. The focus is laid on forced convection flows related to reactor design and near-design conditions. It is found, that the method can accurately predict the wall temperature even under deterioration conditions as they occur in the selected experiments (Yamagata el al. 1972 at 24.5 MPa, Ornatski et al. 1971 at 25.5 and Swenson et al. 1963 at 22.75 MPa). Comparison of the friction coefficient under high heat flux conditions including significant viscosity and density reductions near the wall with various correlations for the hydraulic resistance will be presented; the best agreement is achieve with the correlation of Pioro et al. 2004. It is concluded, that the semi-analytical prediction method is accurate and efficient for forced convection flows even at high heat fluxes. The wall-roughness height is identified as an important model parameter. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

rHARM: ACCRETION AND EJECTION IN RESISTIVE GR-MHD

Turbulent magnetic diffusivity plays an important role for accretion disks and the launching of disk winds. We have implemented magnetic diffusivity and respective resistivity in the general relativistic MHD code HARM. This paper describes the theoretical background of our implementation, its numerical realization, our numerical tests, and preliminary applications. The test simulations of the new code rHARM are compared to an analytic solution of the diffusion equation and a classical shock tube problem. We have further investigated the evolution of the magnetorotational instability (MRI) in tori around black holes (BHs) for a range of magnetic diffusivities. We find an indication for a critical magnetic diffusivity (for our setup) beyond which no MRI develops in the linear regime and for which accretion of torus material to the BH is delayed. Preliminary simulations of magnetically diffusive thin accretion disks around Schwarzschild BHs that are threaded by a large-scale poloidal magnetic field show the launching of disk winds with mass fluxes of about 50% of the accretion rate. The disk magnetic diffusivity allows for efficient disk accretion that replenishes the mass reservoir of the inner disk area and thus allows for long-term simulations of wind launching for more than 5000 time units.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Oblique water wave damping by two submerged thin vertical porous plates of different heights

Oblique water wave damping by two fully submerged vertical parallel porous plates of different heights is investigated within the framework of linear water wave theory. A channel of infinite horizontal extent but of finite depth is divided into three regions in each of which three different velocity potentials are considered and subsequently the corresponding boundary value problems are formulated. By means of eigenfunction expansion and least square method, complete analytical solution in each of the regions is obtained. With the help of the matching conditions along the vertical boundaries between any two successive regions, the reflection and transmission coefficients, and hence the energy loss (in %), are obtained numerically by applying a matrix method. Numerical study is carried out for various relevant parameters and it is found that with an appropriate choice of the parameters and positioning of the porous plates, both of the reflection and transmission coefficients become significantly low and consequently the porous plates can be utilized as an effecting wave absorber. Comparison of this work with earlier available work is made, and the good agreement in the results suggests the applicability of this method in further engineering purposes. A particular case of two free surface piercing plates is also investigated on similar line and the results are presented in order to observe the difference with the earlier case.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Effect of quantum correction on nonlinear thermal wave of electrons driven by laser heating

In thermal interaction of laser pulse with a deuterium-tritium (DT) plane, the thermal waves of electrons are generated instantly. Since the thermal conductivity of electron is a nonlinear function of temperature, a nonlinear heat conduction equation is used to investigate the propagation of waves in solid DT. This paper presents a self-similar analytic solution for the nonlinear heat conduction equation in a planar geometry. The thickness of the target material is finite in numerical computation, and it is assumed that the laser energy is deposited at a finite initial thickness at the initial time which results in a finite temperature for electrons at initial time. Since the required temperature range for solid DT ignition is higher than the critical temperature which equals 35.9 eV, the effects of quantum correction in thermal conductivity should be considered. This letter investigates the effects of quantum correction on characteristic features of nonlinear thermal wave, including temperature, penetration depth, velocity, heat flux, and heating and cooling domains. Although this effect increases electron temperature and thermal flux, penetration depth and propagation velocity are smaller. This effect is also applied to re-evaluate the side-on laser ignition of uncompressed DT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Tkachenko modes and their damping in the vortex lattice regime of rapidly rotating bosons

We have found an exact analytical solution of the Bogoliubov-de Gennes equations for the Tkachenko modes of the vortex lattice in the lowest Landau level (LLL) in the thermodynamic limit (geometry of an infinite plane) at any momenta and calculated their damping rates. At finite temperatures both Beliaev and Landau damping leads to momentum-independent damping rates in the low-energy limit, which shows that, at sufficiently low energies, Tkachenko modes become strongly damped. We then found that the mean-square fluctuations of the density grow logarithmically at large distances, which indicates that the state is ordered in the vortex lattice only on a finite (although exponentially large) distance scale and introduces a low-momentum cutoff. Using this circumstance we showed that at finite temperatures the one-body density matrix undergoes an exponential decay at large distances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Fast iterative refinement method for mixed systems of integral and fractional integro-differential equations

The authors’ new method of iterative refinement has been successfully applied to the linear fractional Fredholm integro-differential equation. In this work, we adapt our method to obtain approximate analytical solutions for linear Volterra equations of both the first and the second kind of fractional type. A detailed convergence analysis is presented for each kind. The authors also apply their method to the challenging first kind linear Fredholm integral equation. Three different cases are considered to test the efficacy of our method. These include mixed systems of various forms of linear integral and fractional integro-differential equations. We compare our method with six well-established methods: Picard’s successive approximations, an accelerated form of the latter, the Adomian decomposition, the variational iteration, the homotopy perturbation and the homotopy analysis. The results show the versatility of the new method in solving a wide variety of equations, its high convergence speed and accuracy and its capability of directly dealing with equations of the first kind. We also prove that Picard’s method for linear equations is a special case of our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Magnetohydrostatic equilibrium. II. Three-dimensional multiple open magnetic flux tubes in the stratified solar atmosphere

A system of multiple open magnetic flux tubes spanning the solar photosphere and lower corona is modeled analytically, within a realistic stratified atmosphere subject to solar gravity. This extends results for a single magnetic flux tube in magnetohydrostatic equilibrium, described in Gent et al. Self-similar magnetic flux tubes are combined to form magnetic structures, which are consistent with high-resolution observations. The observational evidence supports the existence of strands of open flux tubes and loops persisting in a relatively steady state. Self-similar magnetic flux tubes, for which an analytic solution to the plasma density and pressure distribution is possible, are combined. We calculate the appropriate balancing forces, applying to the equations of momentum and energy conservation to preserve equilibrium. Multiplex flux tube configurations are observed to remain relatively stable for up to a day or more, and it is our aim to apply our model as the background condition for numerical studies of energy transport mechanisms from the solar surface to the corona. We apply magnetic field strength, plasma density, pressure, and temperature distributions consistent with observational and theoretical estimates for the lower solar atmosphere. Although each flux tube is identical in construction apart from the location of the radial axis, combinations can be applied to generate a non-axisymmetric magnetic field with multiple non-uniform flux tubes. This is a considerable step forward in modeling the realistic magnetized three-dimensional equilibria of the solar atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Application of matched asymptotic expansions to lunar and interplanetary trajectories. Volume 1: Technical discussion

Previously published asymptotic solutions for lunar and interplanetary trajectories have been modified and combined to formulate a general analytical solution to the problem on N-bodies. The earlier first-order solutions, derived by the method of matched asymptotic expansions, have been extended to second order for the purpose of obtaining increased accuracy. The derivation of the second-order solution is summarized by showing the essential steps, some in functional form. The general asymptotic solution has been used as a basis for formulating a number of analytical two-point boundary value solutions. These include earth-to-moon, one- and two-impulse moon-to-earth, and interplanetary solutions. The results show that the accuracies of the asymptotic solutions range from an order of magnitude better than conic approximations to that of numerical integration itself. Also, since no iterations are required, the asymptotic boundary value solutions are obtained in a fraction of the time required for comparable numerically integrated solutions. The subject of minimizing the second-order error is discussed, and recommendations made for further work directed toward achieving a uniform accuracy in all applications.

Lancaster, J. E.↗
DOE OSTI2016

KELVIN–HELMHOLTZ INSTABILITY IN SOLAR CHROMOSPHERIC JETS: THEORY AND OBSERVATION

Using data obtained by the high-resolution CRisp Imaging SpectroPolarimeter instrument on the Swedish 1 m Solar Telescope, we investigate the dynamics and stability of quiet-Sun chromospheric jets observed at the disk center. Small-scale features, such as rapid redshifted and blueshifted excursions, appearing as high-speed jets in the wings of the H α line, are characterized by short lifetimes and rapid fading without any descending behavior. To study the theoretical aspects of their stability without considering their formation mechanism, we model chromospheric jets as twisted magnetic flux tubes moving along their axis, and use the ideal linear incompressible magnetohydrodynamic approximation to derive the governing dispersion equation. Analytical solutions of the dispersion equation indicate that this type of jet is unstable to Kelvin–Helmholtz instability (KHI), with a very short (few seconds) instability growth time at high upflow speeds. The generated vortices and unresolved turbulent flows associated with the KHI could be observed as a broadening of chromospheric spectral lines. Analysis of the H α line profiles shows that the detected structures have enhanced line widths with respect to the background. We also investigate the stability of a larger-scale H α jet that was ejected along the line of sight. Vortex-like features, rapidly developing around the jet’s boundary, are considered as evidence of the KHI. The analysis of the energy equation in the partially ionized plasma shows that ion–neutral collisions may lead to fast heating of the KH vortices over timescales comparable to the lifetime of chromospheric jets.

79 ASTRONOMY AND ASTROPHYSICS↗