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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Current drive in recombining plasma

The Langevin equations describing the average collisional dynamics of suprathermal particles in nonstationary plasma remarkably admit an exact analytical solution in the case of recombining plasma. The current density produced by arbitrary particle fluxes is derived including the influence of charge recombination. Since recombination has the effect of lowering the charge density of the plasma, thus reducing the charged particle collisional frequencies, the evolution of the current density can be modified substantially compared to plasma with fixed charge density. The current drive efficiency is derived and optimized for discrete and continuous pulses of current, leading to the discovery of a nonzero ''residual'' current density that persists indefinitely under certain conditions, a feature not present in stationary plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

THE PLASMA ENVIRONMENT IN COMETS OVER A WIDE RANGE OF HELIOCENTRIC DISTANCES: APPLICATION TO COMET C/2006 P1 (MCNAUGHT)

On 2007 January 12, comet C/2006 P1 (McNaught) passed its perihelion at 0.17 AU. Abundant remote observations offer plenty of information on the neutral composition and neutral velocities within 1 million kilometers of the comet nucleus. In early February, the Ulysses spacecraft made an in situ measurement of the ion composition, plasma velocity, and magnetic field when passing through the distant ion tail and the ambient solar wind. The measurement by Ulysses was made when the comet was at around 0.8 AU. With the constraints provided by remote and in situ observations, we simulated the plasma environment of Comet C/2006 P1 (McNaught) using a multi-species comet MHD model over a wide range of heliocentric distances from 0.17 to 1.75 AU. The solar wind interaction of the comet at various locations is characterized and typical subsolar standoff distances of the bow shock and contact surface are presented and compared to analytic solutions. We find the variation in the bow shock standoff distances at different heliocentric distances is smaller than the contact surface. In addition, we modified the multi-species model for the case when the comet was at 0.7 AU and achieved comparable water group ion abundances, proton densities, plasma velocities, and plasma temperatures to the Ulysses/SWICS and SWOOPS observations. We discuss the dominating chemical reactions throughout the comet-solar wind interaction region and demonstrate the link between the ion composition near the comet and in the distant tail as measured by Ulysses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Passive and active plasma deceleration for the compact disposal of electron beams

Plasma-based decelerating schemes are investigated as compact alternatives for the disposal of high-energy beams (beam dumps). Analytical solutions for the energy loss of electron beams propagating in passive and active (laser-driven) schemes are derived. These solutions, along with numerical modeling, are used to investigate the evolution of the electron distribution, including energy chirp and total beam energy. In the active beam dump scheme, a laser-driver allows a more homogeneous beam energy extraction and drastically reduces the energy chirp observed in the passive scheme. These concepts could benefit applications requiring overall compactness, such as transportable light sources, or facilities operating at high beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Acoustic scattering of a cylindrical quasi-Gaussian beam with arbitrary incidence focused on a rigid elliptical cylinder

Using the partial-wave series expansion method in cylindrical coordinates, a formal analytical solution for the acoustical scattering of a 2D cylindrical quasi-Gaussian beam with an arbitrary angle of incidence θ{sub i}, focused on a rigid elliptical cylinder in a non-viscous fluid, is developed. The cylindrical focused beam expression is an exact solution of the Helmholtz equation. The scattering coefficients for the elliptical cylinder are determined by forcing the expression of the total (incident + scattered) field to satisfy the Neumann boundary condition for a rigid immovable surface, and performing the product of matrices involving an inversion procedure. Computations for the matrices elements require a single numerical integration procedure for each partial-wave mode. Numerical results are performed with particular emphasis on the focusing properties of the incident beam and its angle of incidence with respect to the major axis a of the ellipse as well as the aspect ratio a/b where b is the minor axis (assuming a > b). The method is validated and verified against previous results obtained via the T-matrix for plane waves. The present analysis is the first to consider an acoustical beam on an elliptic cylinder of variable cross-section as opposed to plane waves of infinite extent. Other 2D non-spherical and Chebyshev surfaces are mentioned that may be examined throughout this analytical formalism assuming a small deformation parameter ε.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Time evolution of decay of two identical quantum particles

An analytical solution for the time evolution of decay of two identical noninteracting quantum particles seated initially within a potential of finite range is derived using the formalism of resonant states. It is shown that the wave function, and hence also the survival and nonescape probabilities, for factorized symmetric and entangled symmetric or antisymmetric initial states evolve in a distinctive form along the exponentially decaying and nonexponential regimes. Our findings show the influence of the Pauli exclusion principle on decay. We exemplify our results by solving exactly the s-wave {delta} shell potential model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

MAGNETOHYDRODYNAMICS OF THE WEAKLY IONIZED SOLAR PHOTOSPHERE

We investigate the importance of ambipolar diffusion and Hall currents for high-resolution comprehensive ({sup r}ealistic{sup )} photospheric simulations. To do so, we extended the radiative magnetohydrodynamics code MURaM to use the generalized Ohm's law under the assumption of local thermodynamic equilibrium. We present test cases comparing analytical solutions with numerical simulations for validation of the code. Furthermore, we carried out a number of numerical experiments to investigate the impact of these neutral-ion effects in the photosphere. We find that, at the spatial resolutions currently used (5-20 km per grid point), the Hall currents and ambipolar diffusion begin to become significant-with flows of 100 m s{sup -1} in sunspot light bridges, and changes of a few percent in the thermodynamic structure of quiet-Sun magnetic features. The magnitude of the effects is expected to increase rapidly as smaller-scale variations are resolved by the simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

ON HYDROMAGNETIC STRESSES IN ACCRETION DISK BOUNDARY LAYERS

Detailed calculations of the physical structure of accretion disk boundary layers, and thus their inferred observational properties, rely on the assumption that angular momentum transport is opposite to the radial angular frequency gradient of the disk. The standard model for turbulent shear viscosity satisfies this assumption by construction. However, this behavior is not supported by numerical simulations of turbulent magnetohydrodynamic (MHD) accretion disks, which show that angular momentum transport driven by the magnetorotational instability (MRI) is inefficient in disk regions where, as expected in boundary layers, the angular frequency increases with radius. In order to shed light on physically viable mechanisms for angular momentum transport in this inner disk region, we examine the generation of hydromagnetic stresses and energy density in differentially rotating backgrounds with angular frequencies that increase outward in the shearing-sheet framework. We isolate the modes that are unrelated to the standard MRI and provide analytic solutions for the long-term evolution of the resulting shearing MHD waves. We show that, although the energy density of these waves can be amplified significantly, their associated stresses oscillate around zero, rendering them an inefficient mechanism to transport significant angular momentum (inward). These findings are consistent with the results obtained in numerical simulations of MHD accretion disk boundary layers and challenge the standard assumption of efficient angular momentum transport in the inner disk regions. This suggests that the detailed structure of turbulent MHD accretion disk boundary layers could differ appreciably from those derived within the standard framework of turbulent shear viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Scattering of electromagnetic waves from a magnetized plasma column at oblique incidence

Scattering of electromagnetic waves from a magnetized plasma column is investigated using Maxwell's equations and applying boundary conditions. Backscattering cross section is evaluated by analytic solution of electric fields inside and outside of plasma column. Plots of backscattering cross section versus frequency, for the range up to J band, reveal two main peaks and two sidebands. Effects of plasma density and radius, as main parameters determining the characteristics of plasma column, on backscattering are discussed. Furthermore, the effect of electromagnetic wave incidence angle on backscattering of plasma column is included in the analysis. The influence of wave incidence angle and frequency, as well as, plasma density and radius on scattering pattern, which is an indicator of the distribution of scattered power in different azimuthal angles, is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Symmetric and asymmetric solitons trapped in H-shaped potentials

We report results of numerical and analytical studies of the spontaneous symmetry breaking in solitons, both two- and one-dimensional (2D and 1D), which are trapped in H-shaped potential profiles built of two parallel potential troughs linked by a narrow rung in the transverse direction. This system can be implemented in self-attractive Bose-Einstein condensates (BECs), as well as in a nonlinear bulk optical waveguide. We demonstrate that the introduction of the transverse link changes the character of the symmetry-breaking bifurcation (SBB) in the system from subcritical to supercritical. (In terms of the corresponding phase transition, it is a change between the first and second kinds.) A noteworthy feature of the SBB in this setting is a nonmonotonous dependence of the soliton's norm at the bifurcation point on the strength of the transverse link. In the full 2D system, the results are obtained in a numerical form. An exact analytical solution is found for the bifurcation in the 1D version of the model, with the transverse rung modeled by the local linear coupling between the parallel troughs with the {delta}-functional longitudinal profile. By replacing the {delta} function by its finite-width Gaussian counterpart, similar results are obtained by means of the variational approximation (VA). The VA is also applied to the 1D system with a mixed linear and nonlinear transverse localized coupling. Comparison of the results produced by the different varieties of the system clearly reveals basic features of its the symmetry-breaking transition.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Electrons in Quasi-1D Waveguides with the Spin–Orbit Interaction: Manifestations of Additional Spin Symmetry

We consider quasi-one-dimensional electron waveguides with spin–orbit interaction, which are formed in quantum wells grown along arbitrary crystallographic directions. An analytic solution to the Schrödinger equation is obtained for systems with Hamiltonians possessing additional spin symmetry. It is shown that the dispersion curves for electrons, which correspond to different size-quantization modes, can intersect only when such symmetry exists. We analyze the structure of dips appearing on the dependences of the conductance of an inhomogeneous waveguide on the energy of carriers. It is shown that the width of the dips substantially depends on the waveguide orientation in the plane of the quantum well. In particular, it vanishes when the waveguide is formed along the direction of the “magic” vector of the initial 2D system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2010

Formation of electrostatic solitons, monotonic, and oscillatory shocks in pair-ion plasmas

The nonlinear electrostatic structures in homogeneous, unmagnetized pair-ion plasma are studied. The dissipation in the system is taken through kinematic viscosities of both pair-ion species. The one dimensional (Korteweg-de Vries-Burgers) KdVB equation is derived using reductive perturbation method. The analytical solution of KdVB is obtained using tanh method. It is found that solitons and monotonic shocks structures are formed in such type of plasmas depending on the value of dissipation in the system. Both compressive and rarefactive structures of solitons and monotonic shocks are obtained depending on the temperatures of negative and positive ions. The oscillatory shock structure in pair-ion plasmas is obtained and its necessary conditions for formation are discussed. The numerical illustrations of potential structures for different values of dissipation in the system are also shown, which may have some relevance in the future experiments of laboratory produced pair-ion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

RADIATION TRANSPORT FOR EXPLOSIVE OUTFLOWS: OPACITY REGROUPING

Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) are methods used to stochastically solve the radiative transport and diffusion equations, respectively. These methods combine into a hybrid transport-diffusion method we refer to as IMC-DDMC. We explore a multigroup IMC-DDMC scheme that in DDMC, combines frequency groups with sufficient optical thickness. We term this procedure ''opacity regrouping''. Opacity regrouping has previously been applied to IMC-DDMC calculations for problems in which the dependence of the opacity on frequency is monotonic. We generalize opacity regrouping to non-contiguous groups and implement this in SuperNu, a code designed to do radiation transport in high-velocity outflows with non-monotonic opacities. We find that regrouping of non-contiguous opacity groups generally improves the speed of IMC-DDMC radiation transport. We present an asymptotic analysis that informs the nature of the Doppler shift in DDMC groups and summarize the derivation of the Gentile-Fleck factor for modified IMC-DDMC. We test SuperNu using numerical experiments including a quasi-manufactured analytic solution, a simple 10 group problem, and the W7 problem for Type Ia supernovae. We find that opacity regrouping is necessary to make our IMC-DDMC implementation feasible for the W7 problem and possibly Type Ia supernova simulations in general. We compare the bolometric light curves and spectra produced by the SuperNu and PHOENIX radiation transport codes for the W7 problem. The overall shape of the bolometric light curves are in good agreement, as are the spectra and their evolution with time. However, for the numerical specifications we considered, we find that the peak luminosity of the light curve calculated using SuperNu is ∼10% less than that calculated using PHOENIX.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics

We compare three thermodynamically consistent numerical fluxes known in the literature, appearing in a Voronoï finite volume discretization of the van Roosbroeck system with general charge carrier statistics. Our discussion includes an extension of the Scharfetter–Gummel scheme to non-Boltzmann (e.g. Fermi–Dirac) statistics. It is based on the analytical solution of a two-point boundary value problem obtained by projecting the continuous differential equation onto the interval between neighboring collocation points. Hence, it serves as a reference flux. The exact solution of the boundary value problem can be approximated by computationally cheaper fluxes which modify certain physical quantities. One alternative scheme averages the nonlinear diffusion (caused by the non-Boltzmann nature of the problem), another one modifies the effective density of states. To study the differences between these three schemes, we analyze the Taylor expansions, derive an error estimate, visualize the flux error and show how the schemes perform for a carefully designed p-i-n benchmark simulation. We present strong evidence that the flux discretization based on averaging the nonlinear diffusion has an edge over the scheme based on modifying the effective density of states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Subspace-based Inverse Uncertainty Quantification for Nuclear Data Assessment

Safety analysis and design optimization depend on the accurate prediction of various reactor attributes. Predictions can be enhanced by reducing the uncertainty associated with the attributes of interest. An inverse problem can be defined and solved to assess the sources of uncertainty, and experimental effort can be subsequently directed to further improve the uncertainty associated with these sources. In this work a subspace-based algorithm for inverse sensitivity/uncertainty quantification (IS/UQ) has been developed to enable analysts account for all sources of nuclear data uncertainties in support of target accuracy assessment-type analysis. An approximate analytical solution of the optimization problem is used to guide the search for the dominant uncertainty subspace. By limiting the search to a subspace, the degrees of freedom available for the optimization search are significantly reduced. A quarter PWR fuel assembly is modeled and the accuracy of the multiplication factor and the fission reaction rate are used as reactor attributes whose uncertainties are to be reduced. Numerical experiments are used to demonstrate the computational efficiency of the proposed algorithm. Our ongoing work is focusing on extending the proposed algorithm to account for various forms of feedback, e.g., thermal-hydraulics and depletion effects.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2014

Computational evaluation of the flexoelectric effect in dielectric solids

Flexoelectricity is a size-dependent electromechanical mechanism coupling polarization and strain gradient. It exists in a wide variety of materials, and is most noticeable for nanoscale objects, where strain gradients are higher. Simulations are important to understand flexoelectricity because experiments at very small scales are difficult, and analytical solutions are scarce. Here, we computationally evaluate the role of flexoelectricity in the electromechanical response of linear dielectric solids in two-dimensions. We deal with the higher-order coupled partial differential equations using smooth meshfree basis functions in a Galerkin method, which allows us to consider general geometries and boundary conditions. We focus on the most common setups to quantify the flexoelectric response, namely, bending of cantilever beams and compression of truncated pyramids, which are generally interpreted through approximate solutions. While these approximations capture the size-dependent flexoelectric electromechanical coupling, we show that they only provide order-of-magnitude estimates as compared with a solution fully accounting for the multidimensional nature of the problem. We discuss the flexoelectric mechanism behind the enhanced size-dependent elasticity in beam configurations. We show that this mechanism is also responsible for the actuation of beams under purely electrical loading, supporting the idea that a mechanical flexoelectric sensor also behaves as an actuator. The predicted actuation-induced curvature is in a good agreement with experimental results. The truncated pyramid configuration highlights the critical role of geometry and boundary conditions on the effective electromechanical response. Our results suggest that computer simulations can help understanding and quantifying the physical properties of flexoelectric devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Heat and Mass Transfer in the Drying of a Cylindrical Body in an Oscillating Magnetic Field

A problem on the heating of a cylindrical body of infinite length in an oscillating electromagnetic field in the process of its drying has been formulated and solved analytically with account of the intermittence of irradiation of the body defined by the Heaviside unit function, the exponential-law absorption of electromagnetic energy by it, and the convective heat and mass exchange between the surface of the body and the environment having constant parameters. The intensity of evaporation of moisture from the surface of the body was determined on the basis of analytical solution of the problem on the mass transfer (moisture diffusion) in it on the assumption that the phase transformations of the body proceed near its surface. Solutions of the problem on the heating of the cylindrical body have been obtained for the cases of nonuniform and uniform distributions of its local temperature, the temperature of the body averaged over its volume, and the temperature gradient near the surface of the body. The 'serviceability' of these solutions was verified on the basis of numerical simulation, with them, of the drying of a seed shaped as a cylinder under the action of an oscillating infrared radiation. As a result of the numerical simulation performed, a technological regime of drying of seeds at minimum and maximum temperatures of their heating by on oscillating infrared radiation for a definite period of time in a cycle, providing not only the drying of the seeds but also substantial improvement of their sowing properties (the sprouting energy and the germination power), has been found. It is shown that the oscillating infrared heating of seeds can be used for their drying in pseudofluidized and vibrofluidized beds.

42 ENGINEERING↗
DOE OSTI2017

Dissipation-preserving spectral element method for damped seismic wave equations

This article describes the extension of the conformal symplectic method to solve the damped acoustic wave equation and the elastic wave equations in the framework of the spectral element method. The conformal symplectic method is a variation of conventional symplectic methods to treat non-conservative time evolution problems, which has superior behaviors in long-time stability and dissipation preservation. To reveal the intrinsic dissipative properties of the model equations, we first reformulate the original systems in their equivalent conformal multi-symplectic structures and derive the corresponding conformal symplectic conservation laws. We thereafter separate each system into a conservative Hamiltonian system and a purely dissipative ordinary differential equation system. Based on the splitting methodology, we solve the two subsystems respectively. The dissipative one is cheaply solved by its analytic solution. While for the conservative system, we combine a fourth-order symplectic Nyström method in time and the spectral element method in space to cover the circumstances in realistic geological structures involving complex free-surface topography. The Strang composition method is adopted thereby to concatenate the corresponding two parts of solutions and generate the completed conformal symplectic method. A relative larger Courant number than that of the traditional Newmark scheme is found in the numerical experiments in conjunction with a spatial sampling of approximately 5 points per wavelength. A benchmark test for the damped acoustic wave equation validates the effectiveness of our proposed method in precisely capturing dissipation rate. The classical Lamb problem is used to demonstrate the ability of modeling Rayleigh wave in elastic wave propagation. More comprehensive numerical experiments are presented to investigate the long-time simulation, low dispersion and energy conservation properties of the conformal symplectic methods in both the attenuating homogeneous and heterogeneous media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗