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.
A functional equation is derived for a spatially homogeneous distribution function and examples of spatially homogeneous gas motions are considered. Numerical calculation was conducted for one such motion with the use of a model equation.
This paper presents the numerical simulations of confined three-dimensional coaxial water jets. The objectives are to validate the newly proposed nonlinear turbulence models of momentum and scalar transport, and to evaluate the newly introduced scalar APDF and DWFDF equation along with its Eulerian implementation in the National Combustion Code (NCC). Simulations conducted include the steady RANS, the unsteady RANS (URANS), and the time-filtered Navier-Stokes (TFNS); both without and with invoking the APDF or DWFDF equation. When the APDF (ensemble averaged probability density function) or DWFDF (density weighted filtered density function) equation is invoked, the simulations are of a hybrid nature, i.e., the transport equations of energy and species are replaced by the APDF or DWFDF equation. Results of simulations are compared with the available experimental data. Some positive impacts of the nonlinear turbulence models and the Eulerian scalar APDF and DWFDF approach are observed.
Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field.
A typical semiconductor-based optoelectronic device, such as a diode laser, consists of three subsystems: an optical field, an electron-hole plasma (EHP), and a host crystal lattice. The physics of such a device involves the interplay of optical, electrical and thermal processes. A proper description of such a device requires that all three processes are treated on equal footing and in a self-consistent fashion. Furthermore, since a semiconductor laser has intrinsic spatial inhomogeneity, such a self-consistency naturally leads to a set of partial differential equations in space and time. There is a significant lacking of research interest and results on the transport aspects of optical devices in the literature with only a few exceptions. Even the most important carrier diffusion coefficient has not been properly derived and studied so far for optically excited plasma, while most of the work adopted results from electronics community where heavily doped semiconductors with mainly one type of carriers are dealt with. The corresponding transport equation for plasma energy or temperature has received even less attention. In this talk we describe our recent results on such a self-consistent derivation of temperature and carrier-density diffusion equations coupled with the lasing process. Starting from the microscopic semiconductor Bloch equations (SBEs) including the Boltzmann transport terms in the distribution function equations for electrons and holes, we derived a closed set of diffusion equations for carrier densities and temperatures with self-consistent coupling to Maxwell's equation and to an effective optical polarization equation. The coherent many-body effects are included within the screened Hartree-Fock approximation, while scatterings are treated within the second Born approximation including both the in- and out-scatterings. Microscopic expressions for electron-hole (e-h) and carrier-LO (c-LO) phonon scatterings are directly used to derive the momentum and energy relaxation rates. These rates expressed as functions of temperatures and densities lead to microscopic expressions for self- and mutual-diffusion coefficients in the coupled density-temperature diffusion equations. Approximations for reducing the general two-component description of the electron-hole plasma (EHP) to a single-component one are discussed. In particular, we show that a special single-component reduction is possible when e-h scattering dominates over c-LO phonon scattering. The ambipolar diffusion approximation is also discussed and we show that the ambipolar diffusion coefficients are independent of e-h scattering, even though the diffusion coefficients of individual components depend sensitively on the e-h scattering rates. Our discussions lead to new perspectives into the roles played in the single-component reduction by the electron-hole correlation in momentum space induced by scatterings and the electron-hole correlation in real space via internal static electrical field. Finally, the theory is completed by coupling the diffusion equations to the lattice temperature equation and to the effective optical polarization which in turn couples to the laser field. The equations derived above are implemented in various limiting cases to a typical diode laser to study the consequences of nonlinear diffusion and the cross diffusion terms on laser behavior, especially the dynamic behavior of a diode laser under modulation. Detailed results will be presented by comparing with the standard rate equation results.
A numerical technique to predict the viscous behavior in the three-dimensional domain of rotating and stationary blade passages is presented. The analysis is based upon the numerical integration of the incompressible Navier-Stokes equations. An approach is taken to combine the three momentum equations and the continuity equation into two sets of vorticity transport-stream function equations. These equations are expressed along families of arbitrarily defined orthogonally intersecting surfaces placed within the blade passage. Unlike the stream surfaces, the arbitrary surfaces allow fluid to flow across. The numerical code employs a nonorthogonal body fitted coordinate system on each of the intersecting surfaces. The solution process involves the iterative combination of the solutions on the two families of intersecting surfaces for a complete solution. The results of flow in thin blade stationary and rotating polar ducts is presented.
Differential equations soluble in finite terms of elementary functions
Nonlinear parabolic differential equation stability, boundedness and uniqueness of solutions obtained by Liapunov direct method
A comparative study of the efficiency of some finite difference methods for the solution of the Navier-Stokes equations was conducted. The study was restricted to the two-dimensional steady, uniform property vorticity-stream function equations. The comparisons were drawn by recording the CPU time required to obtain a solution as well as the accuracy of this solution using five numerical methods: central differences, first order upwind differences, second order upwind differences, exponential differences, and an ADI solution of the central difference equations. Solutions were obtained for two test cases: a recirculating eddy inside a square cavity with a moving top, and an impinging jet flow. The results show that whenever the central difference method is stable it generates results with a given accuracy for less CPU time than any other method.
The capabilities of the Marshall system for aerospace system simulation (MARSYAS) and how to use it are described. MARSYAS is a software system that allows easy setup and control of the simulation of the dynamics of large physical systems on a digital computer. The physical systems are modeled in the form of block diagrams or equations. The blocks can have multiple inputs and multiple outputs, and they can be nested to form hierarchies. The block diagrams can contain transfer functions, nonlinear and logical functions, equations, analog computer elements and FORTRAN programs. The input format of the equations can be combinations of nonlinear, time-varying differential equations and algebraic equations in their original format. MARSYAS could also serve as a storage and retrieval system for models as a basis for a model configuration control system on a central time-shared computer. The outputs of the simulation system can be not only time-responses but also other analysis data such as frequency response, power spectrum and stability parameters. The MARSYAS translator is written in FORTRAN running on the Univac 1108 computer under the EXEC 8 operating system.
Rapid iterative (or semidirect) computation methods are developed for the finite-difference solution of the nonlinear equations of subsonic and transonic aerodynamics. At each iteration, a fast, direct elliptic algorithm solves the entire computation field. In an application to subsonic flow over a lifting airfoil, the full nonlinear stream-function equation is solved. Finally, a direct Cauchy-Riemann solver is used for the nonlinear transonic small-disturbance equations for a biconvex airfoil. At M = 0.7, t/c = 0.1 (subcritical), three iterations on a 39 x 32 mesh (totaling 2.45 sec on an IBM 360/67 computer) obtain convergence within 0.1%. A slightly supercritical case requires seven iterations (6.75 sec) for convergence within 1%.
Liapunov function from quadratic polynomial for n-order nonlinear differential equations
Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.
The ability to fabricate additively manufactured laser waveguides with sharp dopant concentration interfaces is limited by diffusion of the dopants at the temperatures required to fully densify the material. Compositional analysis of bilayer samples, where each layer was either undoped YAG or YAG doped with Yb, Lu, or Nd, were fabricated such that all combinations were available for testing. Samples were fabricated at both 1750°C and 1850°C to determine the diffusion behavior of each dopant alone and also in the presence of a second dopant. It was found that the experimental concentration profiles exhibited both intragranular (bulk) and grain boundary contributions, and thus fitting to a complementary error function equation required the use of two diffusion coefficients. Nd always diffused further along grain boundaries than the other dopants, due in part to its small segregation coefficient in YAG. It is shown that the presence of Nd as a counter dopant inhibits intragranular diffusion of other dopants while enhancing their grain boundary diffusion. All the observed trends were attributed to a combination of intragranular lattice strain due to: dopant ions replacing yttrium substitutionally, the relative driving forces for the segregation of dopant ions to grain boundaries, and the ability of one dopant to affect the diffusion of a different dopant in the other direction to maintain charge neutrality.
Classical equations of elasticity in terms of stress components reduced to three stress functions and three equations
The stability of spinning flexible satellites in a force-free environment was analyzed. The satellite was modeled as a rigid core having attached to it a flexible appendage idealized as a collection of particles (point masses) interconnected by springs. Both Liapunov and Routh-Hurwitz stability procedures are used. In the former, the Hamiltonian of the system, constrained through the angular momentum integral so as to admit complete damping, is used as a testing function. Equations of motion are written using the hybrid coordinate formulation, which readily accepts a modal coordinate transformation ultimately allowing truncation to a level amenable to literal stability analysis. Closed form stability criteria are generated for the first mode of a restricted appendage model lying in a plane containing the system center of mass and orthogonal to the spin axis. The effects of spin on flexible bodies are discussed by considering a very elementary particle model. Control of passively unstable spacecraft is briefly considered.
It is shown that if the Hamiltonian constraint of general relativity is imposed as a restriction on the Hamilton principal functional in the classical theory, or on the state functional in the quantum theory, then the momentum constraints are automatically satisfied. This result holds both for closed and open spaces and it means that the full content of the theory is summarized by a single functional equation of the Tomonaga-Schwinger type.
The stream function equation is solved for steady two-dimensional (and axisymmetric) rotational flows. Both finite differences and finite volumes discretization techniques are studied, using generalized body fitted coordinates and unstructured staggered grids, respectively. For inviscid transonic flows, a new artificial viscosity scheme which does not produce any artificial vorticity is introduced, for the stability of the mixed flow calculations and for capturing shocks. The solution of Euler equations, in primitive variables, are also considered. The effects of the artificial viscosity and numerical boundary conditions on the total enthalpy and the vorticity distributions are demonstrated.
The algorithm employed in the present incompressible two-dimensional calculations of an impulsively-started lid-driven cavity has its basis in the time-dependent stream-function equation. While a Crank-Nicholson differencing scheme is used for the diffusion terms, the Adams-Bashforth scheme is used for the convection terms. The periodic asymptotic solutions obtained for Reynolds numbers of 5000 and 10,000 are found to be precisely periodic; it is demonstrated that they have reached asymptotic states. The indicators of that achievement are discussed.