Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXACT SOLUTIONS”

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 919 records · Page 51

Dynamics of coorbital satellite rings

The dynamical behavior of a coorbital satellite ring is studied for N = 2-9 satellites, in terms of a simplified dynamical description, in which the motion is reduced to the separation angles between satellites. The number and stability of different kinds of stationary configurations is explored, revealing that equally spaced rings are not stable against small perturbations for N not greater than 6, while for N = 2-8 there exists another, stable compact solution. Integrations of exact equations confirm these results. Moreover, the systems are found to display chaotic characteristics for a certain range of energy. The behavior can be interpreted in terms of maximum velocity curves, defining the allowed region of the motion in the phase space.

Salo, H.↗

Radiation-driven winds of hot stars. VI - Analytical solutions for wind models including the finite cone angle effect

Analytical solutions for radiation-driven winds of hot stars including the important finite cone angle effect (see Pauldrach et al., 1986; Friend and Abbott, 1986) are derived which approximate the detailed numerical solutions of the exact wind equation of motion very well. They allow a detailed discussion of the finite cone angle effect and provide for given line force parameters k, alpha, delta definite formulas for mass-loss rate M and terminal velocity v-alpha as function of stellar parameters.

Kudritzki, R. P.↗

Ground-state energies of the nonlinear sigma model and the Heisenberg spin chains

A theorem on the O(3) nonlinear sigma model with the topological theta term is proved, which states that the ground-state energy at theta = pi is always higher than the ground-state energy at theta = 0, for the same value of the coupling constant g. Provided that the nonlinear sigma model gives the correct description for the Heisenberg spin chains in the large-s limit, this theorem makes a definite prediction relating the ground-state energies of the half-integer and the integer spin chains. The ground-state energies obtained from the exact Bethe ansatz solution for the spin-1/2 chain and the numerical diagonalization on the spin-1, spin-3/2, and spin-2 chains support this prediction.

Zhang, Shoucheng↗

The nonconvex multi-dimensional Riemann problem for Hamilton-Jacobi equations

Simple inequalities for the Riemann problem for a Hamilton-Jacobi equation in N space dimension when neither the initial data nor the Hamiltonian need be convex (or concave) are presented. The initial data is globally continuous, affine in each orthant, with a possible jump in normal derivative across each coordinate plane, x sub i = 0. The inequalities become equalities wherever a maxmin equals a minmax and thus an exact closed form solution to this problem is then obtained.

Osher, Stanley↗

Spectral element simulations of laminar and turbulent flows in complex geometries

Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.

Karniadakis, George EM↗

Buckling and vibration analysis of a simply supported column with a piecewise constant cross section

An analysis and sample results for the lateral buckling and vibration of a compressively loaded column is presented whose cross section is piecewise constant along its length. The column is symmetric about its mid-span and consists of three sections, the center section having a stiffer cross section than the two identical end sections. Buckling and vibration characteristics of the column are determined from numerical solution of the exact eigenvalue problems. Parametric structural efficiency analyses are performed using a nondimensionalized set of governing equations to determine the optimum ratio between the lengths of the center section and the outer sections based on both buckling load and vibration frequency requirements. In these analyses, two relationships exist. One is between cross-sectional mass and the cross section, and the other is a high-efficiency scheme. The effect of axial load on vibration frequency is also examined and compared with that of a uniform column.

Lake, Mark S.↗

Users manual for updated computer code for axial-flow compressor conceptual design

An existing computer code that determines the flow path for an axial-flow compressor either for a given number of stages or for a given overall pressure ratio was modified for use in air-breathing engine conceptual design studies. This code uses a rapid approximate design methodology that is based on isentropic simple radial equilibrium. Calculations are performed at constant-span-fraction locations from tip to hub. Energy addition per stage is controlled by specifying the maximum allowable values for several aerodynamic design parameters. New modeling was introduced to the code to overcome perceived limitations. Specific changes included variable rather than constant tip radius, flow path inclination added to the continuity equation, input of mass flow rate directly rather than indirectly as inlet axial velocity, solution for the exact value of overall pressure ratio rather than for any value that met or exceeded it, and internal computation of efficiency rather than the use of input values. The modified code was shown to be capable of computing efficiencies that are compatible with those of five multistage compressors and one fan that were tested experimentally. This report serves as a users manual for the revised code, Compressor Spanline Analysis (CSPAN). The modeling modifications, including two internal loss correlations, are presented. Program input and output are described. A sample case for a multistage compressor is included.

Glassman, Arthur J.↗

A direct algorithm for convective adjustment of the vertical temperature profile for an arbitrary critical lapse rate

An efficient direct algorithm of convective adjustment for an arbitrary critical value of the vertical temperature lapse rate gamma is proposed. The algorithm provides an exact and unique solution of a standard convective adjustment problem for models with temperature specified either on nonuniformly spaced levels or for layers of different thicknesses in pressure, sigma, or other vertical coordinate related to pressure. The algorithm may be recommended for use either directly in atmospheric models not explicitly including a hydrologic cycle with prescribed gamma, or as a part of more complicated parameterizations of moist convection, where gamma may be calculated depending on relative humidity.

Akmaev, Rashid A.↗

Optimal and robust controllers for periodic and multirate systems

The problem of optimal rejection of bounded persistent disturbances is solved in the case of linear discrete-time periodic systems. The solution consists of solving an equivalent time-invariant standard l1 optimization problem subject to an additional constraint. This constraint assures the causality of the resulting periodic controller. By the duality theory, the problem is shown to be equivalent to a linear programming problem, which is no harder than the standard l1 problem. Also, it is shown that the method of solution presented applies exactly to the problem of disturbance rejection in the case of multirate sampled data systems. Finally, the results are applied to the problem of robust stabilization of periodic and multirate systems.

Dahleh, Munther A.↗

Nonlinear modes of the tensor Dirac equation and CPT violation

Recently, it has been shown that Dirac's bispinor equation can be expressed, in an equivalent tensor form, as a constrained Yang-Mills equation in the limit of an infinitely large coupling constant. It was also shown that the free tensor Dirac equation is a completely integrable Hamiltonian system with Lie algebra type Poisson brackets, from which Fermi quantization can be derived directly without using bispinors. The Yang-Mills equation for a finite coupling constant is investigated. It is shown that the nonlinear Yang-Mills equation has exact plane wave solutions in one-to-one correspondence with the plane wave solutions of Dirac's bispinor equation. The theory of nonlinear dispersive waves is applied to establish the existence of wave packets. The CPT violation of these nonlinear wave packets, which could lead to new observable effects consistent with current experimental bounds, is investigated.

Reifler, Frank J.↗

Critical comparison of second-order closures with direct numerical simulations of homogeneous turbulence

Recently, several models have been proposed for closing the second-moment equations, in which the velocity-pressure gradient tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so-called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behaviors in different flows have not been fully examined because they are unmeasurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the possibility to do this kind of study. With direct numerical simulation, we may use the solution to calculate exactly the values of these correlation terms and then directly compare them with the values from their modeled formulations. In this paper, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of 45 homogeneous flows, which were done by Rogers et al. (1986) and Lee and Reynolds (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones that give the best performance. The paper also describes the modeling procedure, model constraints, and the various evaluated models. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

Applications of computer algebra to distributed parameter systems

In the analysis of vibrations of continuous elastic systems, one often encounters complicated transcendental equations with roots directly related to the system's natural frequencies. Typically, these equations contain system parameters whose values must be specified before a numerical solution can be obtained. The present paper presents a method whereby the fundamental frequency can be obtained in analytical form to any desired degree of accuracy. The method is based upon truncation of rapidly converging series involving inverse powers of the system natural frequencies. A straightforward method to developing these series and summing them in closed form is presented. It is demonstrated how Computer Algebra can be exploited to perform the intricate analytical procedures which otherwise would render the technique difficult to apply in practice. We illustrate the method by developing two analytical approximations to the fundamental frequency of a vibrating cantilever carrying a rigid tip body. The results are compared to the numerical solution of the exact (transcendental) frequency equation over a range of system parameters.

Storch, Joel A.↗

Discontinuities formed by steepening the nonlinear Alfven waves

Montgomery (1959) first studied how Alfven waves steepen to form collisionless shocks by use of a two-fluid model. He solved the coupled Lorentz equations of motion for electrons and ions and Maxwell equations by a perturbation technique and showed that linearized polarized Alfven waves in the course of time steepen to form shocks. We have extended this work and have obtained exact simple wave solutions and show that in the limit of slow fluid time scales, we recover the earlier result of Montgomery. For an arbitrary fluid time scale, the solutions include elliptically polarized Alfven waves which also steepen and evolve into shocks. These nonlinear Alfven waves are important in the study of the dynamic solar wind discontinuities.

Lee, N. C.↗

Two-Flux and Green's Function Method for Transient Radiative Transfer in a Semi-Transparent Layer

A method using a Green's function is developed for computing transient temperatures in a semitransparent layer by using the two-flux method coupled with the transient energy equation. Each boundary of the layer is exposed to a hot or cold radiative environment, and is heated or cooled by convection. The layer refractive index is larger than one, and the effect of internal reflections is included with the boundaries assumed diffuse. The analysis accounts for internal emission, absorption, heat conduction, and isotropic scattering. Spectrally dependent radiative properties are included, and transient results are given to illustrate two-band spectral behavior with optically thin and thick bands. Transient results using the present Green's function method are verified for a gray layer by comparison with a finite difference solution of the exact radiative transfer equations; excellent agreement is obtained. The present method requires only moderate computing times and incorporates isotropic scattering without additional complexity. Typical temperature distributions are given to illustrate application of the method by examining the effect of strong radiative heating on one side of a layer with convective cooling on the other side, and the interaction of strong convective heating with radiative cooling from the layer interior.

Siegel, Robert↗

Internal Radiation Effects in Zirconia Thermal Barrier Coatings

Using thermal barrier coatings on combustor liners, turbine vanes, and rotating blades is important for reducing metal temperatures in current and advanced aircraft engines. Zirconia is a common coating material, and it is partially transparent to thermal radiation. Radiation becomes more significant as temperatures are raised for higher efficiency in advanced engines. Calculations are often made with radiation effects neglected inside the coating. The effect of radiation is illustrated here, where an analytical procedure is provided by using the two-flux method for the radiative contribution. A detailed study was made of ceramic thermal barrier coatings for diesel engines, and a two-flux analysis was developed for radiation in semitransparent multilayer composites. These efforts provide the basis for the present analysis where illustrative solutions are obtained for typical conditions in an aircraft engine. The formulation and solution of the exact spectral radiative transfer equations including large scattering, as is characteristic of zirconia, are rather complicated. The two-flux method is used here to provide a simplified method.

Siegel, Robert↗

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

Siegel, Robert↗

Filtering Image Records Using Wavelets and the Zakai Equation

Consider the problem of detecting and localizing a faint object moving In an "essentially stationary" background, using a sequence of two-dimensional low-SNR images of the scene. A natural approach consists of "digitizing" each snapshot into a discrete set of observations, sufficiently (perhaps not exactly) matched to the object In question, then tracking the object using an appropriate stochastic filter. The tracking would be expected to make up for the low signal-to-noise ratio, this allowing one to "coherently" process successive images in order to beat down the noise and localize the object. Thus, "tracking" here does not refer to the usual notion of detecting then tracking: rather, we track in order to detect The problem then becomes one of choosing the appropriate image representation as well as the optimal (and necessarily nonlinear filter. We propose exact and approximate solutions using wavelets and the Zakai equation. The smoothness of the wavelets used is required in the derivation of the evolution equation for the conditional density giving the filter, and their orthogonality makes it possible to carry out actual computations of the Ito- and change-of-gauge-terms in the algorithm effectively.

Haddad, Ziad S.↗

Transition-Independent Decentralized Markov Decision Processes

There has been substantial progress with formal models for sequential decision making by individual agents using the Markov decision process (MDP). However, similar treatment of multi-agent systems is lacking. A recent complexity result, showing that solving decentralized MDPs is NEXP-hard, provides a partial explanation. To overcome this complexity barrier, we identify a general class of transition-independent decentralized MDPs that is widely applicable. The class consists of independent collaborating agents that are tied up by a global reward function that depends on both of their histories. We present a novel algorithm for solving this class of problems and examine its properties. The result is the first effective technique to solve optimally a class of decentralized MDPs. This lays the foundation for further work in this area on both exact and approximate solutions.

Becker, Raphen↗