SEARCH · Search NASA
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.
Study of weak solutions of simple wave equation, inviscid Burgers equation, and Euler equations has resulted in technique for accurate prediction of shock waves occurring in inviscid supersonic flows.
The reported investigation is concerned with the solution of systems of differential equations which are derived from a Hamiltonian function in the extended phase space. The problem selected involves a one-dimensional perturbed harmonic oscillator. The van der Pol equation considered has an exact asymptotic value for its amplitude. Comparisons are made between a numerical solution and a known analytical solution. In addition to the van der Pol problem, known solutions regarding the restricted problem of three bodies are used as examples for perturbed Keplerian motion. The extended phase space Hamiltonian discussed by Stiefel and Scheifele (1971) is considered. A description is presented of two canonical formulations of the perturbed harmonic oscillator.
A two-phase proportional navigation scheme is developed for the case of two rigid bodies engaged in a rendezvous/docking maneuver. The target vehicle is nonmaneuvering, but does have constant nonzero angular and linear velocities. Under these conditions, it is shown that previously obtained solutions are not applicable. Analytical solutions are obtained leading to relationships between the transverse and LOS navigation constants. It is shown that the transverse navigation constant for the second phase of the maneuver must be 2. Also, initial conditions necessary for rendezvous are presented.
We describe a method for computing the history of the luminosity function in a young cluster still forming within a molecular cloud complex. Our method, which utilizes detailed results from stellar evolution theory, assumes that clusters arise from the continuous collapse of dense cloud cores over a protracted period of time. It is also assumed that stars reaching the main sequence are distributed in mass according to a prescribed initial mass function (IMF). We keep track separately of the contributions to the luminosity function from the populations of protostars, pre-main-sequence stars, and main-sequence stars. We derive expressions for the fractional contribution of these populations to both the total number of stars produced and the total cluster luminosity. In our model, the number of protostars rises quickly at first, but then levels off to a nearly constant value, which it maintains until the dispersal of the cloud complex. The number fraction of protostars always decreases with time. Averaged over the life of the parent cloud, this fraction is typically a few percent. The protostar mass distribution can be expressed as an integral over the IMF.
We investigate the expansion dynamics of a Bose–Einstein condensate that consists of two components and is initially confined in a quasi-one-dimensional trap. We classify the possible initial states of the two-component condensate by taking into account the nonuniformity of the distributions of its components and construct the corresponding phase diagram in the plane of nonlinear interaction constants. The differential equations that describe the condensate evolution are derived by assuming that the condensate density and velocity depend on the spatial coordinate quadratically and linearly, respectively, which reproduces the initial equilibrium distribution of the condensate in the trap in the Thomas–Fermi approximation. We have obtained self-similar solutions of these differential equations for several important special cases and write out asymptotic formulas describing the condensate motion on long time scales, when the condensate density becomes so low that the interaction between atoms may be neglected. The problem on the dynamics of immiscible components with the formation of dispersive shock waves is considered. We compare the numerical solutions of the Gross–Pitaevskii equations with their approximate analytical solutions and numerically study the situations where the analytical method being used admits no exact solutions.
An exact analytic solution of the Vlasov equation for the plasma distribution in a Paul trap is known to be a Maxwellian and thus, immune to collisions under the assumption of infinitely fast relaxation [K. Shah and H. S. Ramachandran, Phys. Plasmas 15, 062303 (2008)]. In this paper, it is shown that even for a more realistic situation of finite time relaxation, solutions of the Fokker-Planck equation lead to an equilibrium solution of the form of a Maxwellian with oscillatory temperature. This shows that the rf heating observed in Paul traps cannot be caused due to collisional effects alone.
An analytic solution of transient wall motion in a bulk ferromagnetic material is obtained from the dynamic torque equations with a Gilbert damping term. The steady state solution reduces to previously obtained, well known solutions.
Analytical solutions for inviscid supersonic corner flows are virtually nonexistent due to the complexity of the interference geometry. In view of this, numerical solutions for compressive-compressive and expansive-compressive corner flows are obtained. The governing equations are written in strong conservation-law form and are solved iteratively in nonorthogonal conical coordinates by use of a second-order, shock-capturing, finite-difference technique. The computed wave structure and surface pressure distributions are compared with high Reynolds number (Re greater than 2,500,000 ft) experimental data. The results clearly show that the wave structure in the corner is dominated by the inviscid field.
We consider a mathematical model describing the steady motion of a viscoelastic medium of Oldroyd type under the Navier slip condition at the boundary. In the rheological relation, we use the objective regularized Jaumann derivative. We prove the solubility of the corresponding boundary-value problem in the weak setting. We obtain an estimate for the norm of a solution in terms of the data of the problem. We show that the solution set is sequentially weakly closed. Furthermore, we give an analytic solution of the boundary-value problem describing the flow of a viscoelastic fluid in a flat channel under a slip condition at the walls. Bibliography: 13 titles. (paper)
Presence of SU(n) or other Lie group symmetry in a physical system is its powerful, usually underutilized property. In many cases it allows for finding analytical solutions to nonlinear differential equations describing this system. Power of the method is presented on diversified examples from mathematical physics: Lie-group symmetries in finding solutions of generalized, multidimensional theory of gravity; analytical Dirac–equation solutions for description of conducting polymers; stability of qubit states in quantum computers; spatial defects in condensed matter; reconstruction of 3D object from its 2D tomographic image; significant improvement of numerical solutions stability for Euler equations. The next question after obtaining such Lie group symmetric solution is: does a generalized solution with appropriate quantum group symmetry exists for the given physical system, and if yes what is the physical meaning of the deformation parameter q introduced by such solution. In many cases it can be identified. Any SU(n) solution is by its nature singular, assuming a perfect symmetry of the physical system discussed. Such solution gives a powerful insight to theoretical physics, yet the assumption may be too demanding for experimental applications. Deformation parameter q from a quantum group symmetry allows for a continuum of solutions, more applicable to experiments.
In stochastic systems, one is often interested in finding the optimal path that maximizes the probability of escape from a metastable state or of switching between metastable states. Even for simple systems, it may be impossible to find an analytic form of the optimal path, and in high-dimensional systems, this is almost always the case. In this article, we formulate a constructive methodology that is used to compute the optimal path numerically. The method utilizes finite-time Lyapunov exponents, statistical selection criteria, and a Newton-based iterative minimizing scheme. The method is applied to four examples. The first example is a two-dimensional system that describes a single population with internal noise. This model has an analytical solution for the optimal path. The numerical solution found using our computational method agrees well with the analytical result. The second example is a more complicated four-dimensional system where our numerical method must be used to find the optimal path. The third example, although a seemingly simple two-dimensional system, demonstrates the success of our method in finding the optimal path where other numerical methods are known to fail. In the fourth example, the optimal path lies in six-dimensional space and demonstrates the power of our method in computing paths in higher-dimensional spaces.
In this report, we investigate the effects of conductor losses in a multipole-based cable braid magnetic penetration model. Our multipole model uses a mesh of the actual cable geometry, which enables us to model more complicated structures. After summarizing the first principles model formulation, we consider a one-dimensional array of wires, for which an analytical solution is known in the lossless case. We extend this solution to the lossy case by using a complex-valued radius. We also model this structure analytically using a conformal-mapping solution. We then compare both the self-impedance and the transfer impedance results from our first principles cable braid electromagnetic penetration model to those obtained using the analytical solutions. An analysis for various frequencies (and skin depths) usually encountered in cable modeling is reported. These results are found in good agreement up to a radius to half spacing ratio of about 0.7, demonstrating a robustness needed for many commercial and non-commercial cables.
Recently a new class of biocompatible elastic polymers loaded with small ferrous particles, a magnetoelastic polymer, has been developed. This engineered material is formed into a thin film using spin casting. An applied magnetic field will deform the film. The magnetic deformation of this film has many possible applications, particularly in microfluidic pumps and pressure regulators. In this paper a finite element method suitable for the transient simulation of arbitrarily shaped three-dimensional magnetoelastic polymers subjected to time-varying magnetic fields is developed. The approach is similar to that employed in finite elment magnetohydrodynamic simulations, the key difference is a more complex hyperelastic material model. In order to confirm the validity of the approach, finite element solutions for an axially symmetric thin film are compared to an analytical solution based on the membrane (infinitely thin) approximation. For this particular problem the two approaches give qualitatively similar results and converge as the film thickness approaches zero.
We analytically solve the position-dependent mass (PDM) 1D Schrödinger equation for a new class of hyperbolic potentials V{sub q}{sup p}(x)=−V{sub 0}(sinh{sup p}x/cosh{sup q}x), p=−2,0,⋯q [see C. A. Downing, J. Math. Phys. 54, 072101 (2013)] among several hyperbolic single- and double-wells. For a solitonic mass distribution, m(x)=m{sub 0} sech{sup 2}(x), we obtain exact analytic solutions to the resulting differential equations. For several members of the class, the quantum mechanical problems map into confluent Heun differential equations. The PDM Poschl-Teller potential is considered and exactly solved as a particular case.
This paper reports the space-charge limited current (SLC) and virtual cathode behaviors in one-dimensional grounded drift space. A simple general analytical solution and an approximate solution for the planar diode are given. Through a semi-analytical method, a general solution for SLC in one-dimensional drift space is obtained. The behaviors of virtual cathode in the drift space, including dominant frequency, electron transit time, position, and transmitted current, are yielded analytically. The relationship between the frequency of the virtual cathode oscillation and the injected current presented may explain previously reported numerical works. Results are significant in facilitating estimations and further analytical studies.
Recent observations near the Galactic Center (GC) have found several molecular filaments displaying striking helically wound morphology that are collectively known as molecular tornadoes. We investigate the equilibrium structure of these molecular tornadoes by formulating a magnetohydrodynamic model of a rotating, helically magnetized filament. A special analytical solution is derived where centrifugal forces balance exactly with toroidal magnetic stress. From the physics of torsional Alfvén waves we derive a constraint that links the toroidal flux-to-mass ratio and the pitch angle of the helical field to the rotation laws, which we find to be an important component in describing the molecular tornado structure. The models are compared to the Ostriker solution for isothermal, nonmagnetic, nonrotating filaments. We find that neither the analytic model nor the Alfvén wave model suffer from the unphysical density inversions noted by other authors. A Monte Carlo exploration of our parameter space is constrained by observational measurements of the Pigtail Molecular Cloud, the Double Helix Nebula, and the GC Molecular Tornado. Observable properties such as the velocity dispersion, filament radius, linear mass, and surface pressure can be used to derive three dimensionless constraints for our dimensionless models of these three objects. A virial analysis of these constrained models is studied for these three molecular tornadoes. We find that self-gravity is relatively unimportant, whereas magnetic fields and external pressure play a dominant role in the confinement and equilibrium radial structure of these objects.
On the basis of an approximate analytic solution of a Cauchy problem for a nonlinear Schroedinger (NLS) equation describing steady-state light beams in a medium with saturating nonlinearity by the method of renormgroup (RG) symmetries, a classification of self-focusing solutions is given depending on two control parameters: the relative contributions of diffraction and nonlinearity and the saturation strength of the nonlinearity. The existence of tube-type self-focusing solutions is proved for an entering beam with Gaussian radial distribution of intensity. Numerical simulation is carried out that allows one to verify the theory developed and to determine its applicability limits.
The classical Neumann boundary value problem of an isotropic, homogeneous elastic half-plane under plane strain conditions is readdressed as the limiting case of the fully three-dimensional problem. Analytical solutions of the stress and strain tensors are obtained by taking the limit from known three-dimensional solutions. It is shown that the displacement fields for the plane strain problem are not well defined. A small number of simple expressions are developed, which provide a general solution for linearly-varying traction over arbitrary regions on the boundary. A simple, efficient, and rapidly convergent algorithm is developed which uses these solutions as analytic elements and provides a solution approach to the general boundary value problem. The method is verified against known solutions for Hertzian contact between parallel cylinders. Two numerical examples are presented for the analysis of shallow foundation systems. In the first, the boundary conditions are informed by analytical elastoplastic calculations and a strain influence analysis is performed and compared with the Schmertmann method. Subsequently, empirical laboratory contact traction distributions measured by Bauer et al., in both the normal and tangential directions are employed as boundary conditions for an analysis of the underlying stress field.