Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution theory”

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 433 records · Page 24

Interactive real time flow simulations

An interactive real time flow simulation technique is developed for an unsteady channel flow. A finite-volume algorithm in conjunction with a Runge-Kutta time stepping scheme was developed for two-dimensional Euler equations. A global time step was used to accelerate convergence of steady-state calculations. A raster image generation routine was developed for high speed image transmission which allows the user to have direct interaction with the solution development. In addition to theory and results, the hardware and software requirements are discussed.

Sadrehaghighi, I.↗

Modulational stability of rotating-disk flow

The problem of modulational primary stability of rotating-disk boundary layer flow for stationary modes was studied by using the method of multiple scales. The system of equations for infinitesimal disturbances is solved by linear stability theory to determine the modulated solution which represents the time and azimuthal variations of the wave pattern. The solution for the amplitude of the most critical disturbance indicates that the disturbance spreads out as the flow time increases and grows exponentially at the center.

Riahi, Daniel N.↗

Variational formulation of hybrid problems for fully 3-D transonic flow with shocks in rotor

Based on previous research, the unified variable domain variational theory of hybrid problems for rotor flow is extended to fully 3-D transonic rotor flow with shocks, unifying and generalizing the direct and inverse problems. Three variational principles (VP) families were established. All unknown boundaries and flow discontinuities (such as shocks, free trailing vortex sheets) are successfully handled via functional variations with variable domain, converting almost all boundary and interface conditions, including the Rankine Hugoniot shock relations, into natural ones. This theory provides a series of novel ways for blade design or modification and a rigorous theoretical basis for finite element applications and also constitutes an important part of the optimal design theory of rotor bladings. Numerical solutions to subsonic flow by finite elements with self-adapting nodes given in Refs., show good agreement with experimental results.

Liu, Gao-Lian↗

H2-optimal control with generalized state-space models for use in control-structure optimization

Several advances are provided solving combined control-structure optimization problems. The author has extended solutions from H2 optimal control theory to the use of generalized state space models. The generalized state space models preserve the sparsity inherent in finite element models and hence provide some promise for handling very large problems. Also, expressions for the gradient of the optimal control cost are derived which use the generalized state space models.

Wette, Matt↗

Performance and state-space analyses of systems using Petri nets

The goal of any modeling methodology is to develop a mathematical description of a system that is accurate in its representation and also permits analysis of structural and/or performance properties. Inherently, trade-offs exist between the level detail in the model and the ease with which analysis can be performed. Petri nets (PN's), a highly graphical modeling methodology for Discrete Event Dynamic Systems, permit representation of shared resources, finite capacities, conflict, synchronization, concurrency, and timing between state changes. By restricting the state transition time delays to the family of exponential density functions, Markov chain analysis of performance problems is possible. One major drawback of PN's is the tendency for the state-space to grow rapidly (exponential complexity) compared to increases in the PN constructs. It is the state space, or the Markov chain obtained from it, that is needed in the solution of many problems. The theory of state-space size estimation for PN's is introduced. The problem of state-space size estimation is defined, its complexities are examined, and estimation algorithms are developed. Both top-down and bottom-up approaches are pursued, and the advantages and disadvantages of each are described. Additionally, the author's research in non-exponential transition modeling for PN's is discussed. An algorithm for approximating non-exponential transitions is developed. Since only basic PN constructs are used in the approximation, theory already developed for PN's remains applicable. Comparison to results from entropy theory show the transition performance is close to the theoretic optimum. Inclusion of non-exponential transition approximations improves performance results at the expense of increased state-space size. The state-space size estimation theory provides insight and algorithms for evaluating this trade-off.

Watson, James Francis, III↗

The decay of a simple eddy

The principal result obtained in this report is a generalization of Taylor's formula for a simple eddy. The discussion of the properties of the eddy indicates that there is a slight analogy between the theory of eddies in a viscous fluid and the quantum theory of radiation. Another exact solution of the equations of motion of viscous fluid yields a result which reminds one of the well-known condition for instability in the case of a horizontally stratified atmosphere.

Bateman, H↗

Analytical description of the breakup of liquid jets in air

A viscous or inviscid cylindrical jet with surface tension in a vacuum tends to pinch due to the mechanism of capillary instability. Similarity solutions are constructed which describe this phenomenon as a critical time is encountered, for two physically distinct cases: inviscid jets governed by the Euler equations and highly viscous jets governed by the Stokes equations. In both cases the only assumption imposed is that at the time of pinching the jet shape has a radial length scale which is smaller than the axial length scale. For the inviscid case, we show that our solution corresponds exactly to one member of the one-parameter family of solutions obtained from slender jet theories and the shape of the jet is locally concave at breakup. For highly viscous jets our theory predicts local shapes which are monotonic increasing or decreasing indicating the formation of a mother drop connected to the jet by a thin fluid tube. This qualitative behavior is in complete agreement with both direct numerical simulations and experimental observations.

Papageorgiou, Demetrios T.↗

Scaling, elasticity, and CLPT

The first few viewgraphs describe the general solution properties of linear elasticity theory which are given by the following two statements: (1) for stress B.C. on S(sub sigma) and zero displacement B.C. on S(sub u) the altered displacements u(sub i)(*) and the actual stresses tau(sub ij) are elastically dependent on Poisson's ratio nu alone: thus the actual displacements are given by u(sub i) = mu(exp -1)u(sub i)(*); and (2) for zero stress B.C. on S(sub sigma) and displacement B.C. on S(sub u) the actual displacements u(sub i) and the altered stresses tau(sub ij)(*) are elastically dependent on Poisson's ratio nu alone: thus the actual stresses are given by tau(sub ij) = E tau(sub ij)(*). The remaining viewgraphs describe the minimum parameter formulation of the general classical laminate theory plate problem as follows: The general CLT plate problem is expressed as a 3 x 3 system of differential equations in the displacements u, v, and w. The eighteen (six each) A(sub ij), B(sub ij), and D(sub ij) system coefficients are ply-weighted sums of the transformed reduced stiffnesses (bar-Q(sub ij))(sub k); the (bar-Q(sub ij))(sub k) in turn depend on six reduced stiffnesses (Q(sub ij))(sub k) and the material and geometry properties of the k(sup th) layer. This paper develops a method for redefining the system coefficients, the displacement components (u,v,w), and the position components (x,y) such that a minimum parameter formulation is possible. The pivotal steps in this method are (1) the reduction of (bar-Q(sub ij))(sub k) dependencies to just two constants Q(*) = (Q(12) + 2Q(66))/(Q(11)Q(22))(exp 1/2) and F(*) - (Q(22)/Q(11))(exp 1/2) in terms of ply-independent reference values Q(sub ij); (2) the reduction of the remaining portions of the A, B, and D coefficients to nondimensional ply-weighted sums (with 0 to 1 ranges) that are independent of Q(*) and F(*); and (3) the introduction of simple coordinate stretchings for u, v, w and x,y such that the process is neatly completed.

Brunelle, Eugene J.↗

On the Stability of Three-Dimensional Boundary Layers: Secondary Instability - Pt. 2

The secondary instability of three-dimensional incompressible boundary layers is studied using Floquet theory. Starting from the equilibrium solutions that we obtained from the PSE computations documented in Part 1, we investigate the region where a purely stationary crossflow disturbance saturates for its secondary instability characteristics utilizing developed global and local eigenvalue solvers that are based on the Implicitly Restarted Arnoldi Method, and a Newton-Raphson technique, respectively. The main focuses of this study are on the existence of multiple roots in the eigenvalue spectrum that could explain experimental observations of time-dependent occurrences of an explosive growth of traveling disturbances, on the routes by which high-frequency disturbances enter the boundary layer, as well as on gaining more information about threshold amplitudes for the growth of secondary disturbances.

Janke, Erik↗

An Implicit LU/AF FDTD Method

There has been some recent work to develop two and three-dimensional alternating direction implicit (ADI) FDTD schemes. These ADI schemes are based upon the original ADI concept developed by Peaceman and Rachford and Douglas and Gunn, which is a popular solution method in Computational Fluid Dynamics (CFD). These ADI schemes work well and they require solution of a tridiagonal system of equations. A new approach proposed in this paper applies a LU/AF approximate factorization technique from CFD to Maxwell s equations in flux conservative form for one space dimension. The result is a scheme that will retain its unconditional stability in three space dimensions, but does not require the solution of tridiagonal systems. The theory for this new algorithm is outlined in a one-dimensional context for clarity. An extension to two and threedimensional cases is discussed. Results of Fourier analysis are discussed for both stability and dispersion/damping properties of the algorithm. Results are presented for a one-dimensional model problem, and the explicit FDTD algorithm is chosen as a convenient reference for comparison.

Beggs, John H.↗

Charged-particle transport in one dimension

A numerical solution to high-energy charged-particle transport is found by evaluation of the integral equations obtained by inverting the Boltzmann differential operator. An algorithm has been written for either continuous or discrete spectra at the boundary, which allows efficient and accurate evaluation of this rather complicated problem. Present results are compared with analytic solutions based on the perturbation theory.

NASA Discipline Radiation Health↗

Multiobjective Optimization of Low-Energy Trajectories Using Optimal Control on Dynamical Channels

We introduce a computational method to design efficient low-energy trajectories by extracting initial solutions from dynamical channels formed by invariant manifolds, and improving these solutions through variational optimal control. We consider trajectories connecting two unstable periodic orbits in the circular restricted 3-body problem (CR3BP). Our method leverages dynamical channels to generate a range of solutions, and approximates the areto front for impulse and time of flight through a multiobjective optimization of these solutions based on primer vector theory. We demonstrate the application of our method to a libration orbit transfer in the Earth-Moon system.

trajectory design↗

Polarization and Compressibility of Oblique Kinetic Alfven Waves

Even though solar wind, as a collisionless plasma, is properly described by the kineticMaxwell-Vlasov description, it can be argued that much of our understanding of solar wind observational data comes from an interpretation and numerical modeling which is based on a fluid description of magnetohydrodynamics. In recent years, there has been a significant interest in better understanding the importance of kinetic effects, i.e. the differences between the kinetic and usual fluid descriptions. Here we concentrate on physical properties of oblique kinetic Alfvn waves (KAWs), which are often recognized as one of the key ingredients in the solar wind turbulence cascade. We use three different fluid models with various degrees of complexity and calculate polarization and magnetic compressibility of oblique KAWs (propagation angle q = 88), which we compare to solutions derived from linear kinetic theory. We explore a wide range of possible proton plasma b = [0.1,10.0] and a wide range of length scales krL = [0.001,10.0]. It is shown that the classical isotropic two-fluid model is very compressible in comparison with kinetic theory and that the largest discrepancy occurs at scales larger than the proton gyroscale. We also show that the two-fluid model contains a large error in the polarization of electric field, even at scales krL 1. Furthermore, to understand these discrepancies between the two-fluid model and the kinetic theory, we employ two versions of the Landau fluid model that incorporate linear low-frequency kinetic effects such as Landau damping and finite Larmor radius (FLR) corrections into the fluid description. It is shown that Landau damping significantly reduces the magnetic compressibility and that FLR corrections (i.e. nongyrotropic contributions) are required to correctly capture the polarization.We also show that, in addition to Landau damping, FLR corrections are necessary to accurately describe the damping rate of KAWs. We conclude that kinetic effects are important even at scales which are significantly larger than the proton gyroscale krL 1.

Turbulence↗

Bifurcations from stationary to periodic solutions in a low-order model of forced, dissipative barotropic flow

The considered investigation is concerned with periodic solutions in the context of a forced, dissipative, barotropic spectral model truncated to three complex coefficients with constant forcing on only the intermediate scale. It is found that determining a periodic solution of this three-coefficient model also reduces to finding the algebraic roots of a real polynomial. In the derivation of this polynomial, a class of hydrodynamic spectral systems is described for which a periodic solution might be similarly specified. The existence of periodic solutions of the three-coefficient model is controlled by the roots of the stability polynomial of the basic stationary solution, which represents the simplest response to the constant forcing. When the forcing exceeds a critical value, the basic solution becomes unstable. Owing to the nature of the roots of the stability polynomial at critical forcing, bifurcation theory guarantees the existence of a periodic solution.

Mitchell, K. E.↗

Aerosol measurements in the winter/spring Antarctic stratosphere. I - Correlative measurements with ozone. II - Impact on polar stratospheric cloud theories

Aerosol measurements collected from August 25-November 3, 1986 at McMurdo Station using balloon-borne optical particle counters are examined in order to study the relationship between aerosol and ozone distribution and the formation of polar stratospheric clouds (PSCs). Ozone, aerosol, and condensation nuclei profiles, and pressure, temperature, and humidity measurements are analyzed. It is observed that the height of the stratospheric sulfate layer decreases over the period of measurement suggesting that upwelling in the votex is not important in the zone depletion process. Three theories on PSC formation are described, and the effects of the aerosol measurements on the theories are considered. The three theories are: (1) the original theory of water vapor pressure over a solution of H2SO4 of Steele et al. (1983) and Hamill and Mc Master (1984); (2) the nitric acid theory of PSCs of Toon et al. (1986) and Hamill et al. (1986); and (3) the quasi-cirrus cloud theory of Heymsfield (1986).

Hofmann, D. J.↗

Applications of fuzzy theories to multi-objective system optimization

Most of the computer aided design techniques developed so far deal with the optimization of a single objective function over the feasible design space. However, there often exist several engineering design problems which require a simultaneous consideration of several objective functions. This work presents several techniques of multiobjective optimization. In addition, a new formulation, based on fuzzy theories, is also introduced for the solution of multiobjective system optimization problems. The fuzzy formulation is useful in dealing with systems which are described imprecisely using fuzzy terms such as, 'sufficiently large', 'very strong', or 'satisfactory'. The proposed theory translates the imprecise linguistic statements and multiple objectives into equivalent crisp mathematical statements using fuzzy logic. The effectiveness of all the methodologies and theories presented is illustrated by formulating and solving two different engineering design problems. The first one involves the flight trajectory optimization and the main rotor design of helicopters. The second one is concerned with the integrated kinematic-dynamic synthesis of planar mechanisms. The use and effectiveness of nonlinear membership functions in fuzzy formulation is also demonstrated. The numerical results indicate that the fuzzy formulation could yield results which are qualitatively different from those provided by the crisp formulation. It is felt that the fuzzy formulation will handle real life design problems on a more rational basis.

Rao, S. S.↗

Application of Hamilton's law of varying action

The law of varying action enunciated by Hamilton in 1834-1835 permits the direct analytical solution of the problems of mechanics, both stationary and nonstationary, without consideration of force equilibrium and the theory of differential equations associated therewith. It has not been possible to obtain direct analytical solutions to nonstationary systems through the use of energy theory, which has been limited for 140 years to the principle of least action and to Hamilton's principle. It is shown here that Hamilton's law permits the direct analytical solution to nonstationary, initial value systems in the mechanics of solids without any knowledge or use of the theory of differential equations. Solutions are demonstrated for nonconservative, nonstationary particle motion, both linear and nonlinear.

Bailey, C. D.↗