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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Leveraging Existing Mission Tools in a Re-Usable, Component-Based Software Environment

Emerging methods in component-based software development offer significant advantages but may seem incompatible with existing mission operations applications. In this paper we relate our positive experiences integrating existing mission applications into component-based tools we are delivering to three missions. In most operations environments, a number of software applications have been integrated together to form the mission operations software. In contrast, with component-based software development chunks of related functionality and data structures, referred to as components, can be individually delivered, integrated and re-used. With the advent of powerful tools for managing component-based development, complex software systems can potentially see significant benefits in ease of integration, testability and reusability from these techniques. These benefits motivate us to ask how component-based development techniques can be relevant in a mission operations environment, where there is significant investment in software tools that are not component-based and may not be written in languages for which component-based tools even exist. Trusted and complex software tools for sequencing, validation, navigation, and other vital functions cannot simply be re-written or abandoned in order to gain the advantages offered by emerging component-based software techniques. Thus some middle ground must be found. We have faced exactly this issue, and have found several solutions. Ensemble is an open platform for development, integration, and deployment of mission operations software that we are developing. Ensemble itself is an extension of an open source, component-based software development platform called Eclipse. Due to the advantages of component-based development, we have been able to vary rapidly develop mission operations tools for three surface missions by mixing and matching from a common set of mission operation components. We have also had to determine how to integrate existing mission applications for sequence development, sequence validation, and high level activity planning, and other functions into a component-based environment. For each of these, we used a somewhat different technique based upon the structure and usage of the existing application.

Greene, Kevin↗
NASA NTRS2012

Linked-List-Based Multibody Dynamics (MBDyn) Engine

This new release of MBDyn is a software engine that calculates the dynamics states of kinematic, rigid, or flexible multibody systems. An MBDyn multibody system may consist of multiple groups of articulated chains, trees, or closed-loop topologies. Transient topologies are handled through conservation of energy and momentum. The solution for rigid-body systems is exact, and several configurable levels of nonlinear term fidelity are available for flexible dynamics systems. The algorithms have been optimized for efficiency and can be used for both non-real-time (NRT) and real-time (RT) simulations. Interfaces are currently compatible with NASA's Trick Simulation Environment. This new release represents a significant advance in capability and ease of use. The two most significant new additions are an application programming interface (API) that clarifies and simplifies use of MBDyn, and a link-list infrastructure that allows a single MBDyn instance to propagate an arbitrary number of interacting groups of multibody top ologies. MBDyn calculates state and state derivative vectors for integration using an external integration routine. A Trickcompatible interface is provided for initialization, data logging, integration, and input/output.

MacLean, John↗
NASA NTRS1997

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗
NASA NTRS1982

Preliminary study of the effects of a reversible chemical reaction on gas bubble dissolution

A preliminary investigation is carried out of the effects of a reversible chemical reaction on the dissolution of an isolated, stationary gas bubble in a glass melt. The exact governing equations for the model system are formulated and analyzed. The approximate quasi-steady-state version of these equations is solved analytically, and a calculation is made of bubble dissolution rates. The results are then compared with numerical solutions obtained from the finite difference form of the exact governing equations. It is pointed out that in the microgravity condition of space, the buoyant rise of a gas bubble in a glass melt will be negligible on the time scale of most experiments. For this reason, a determination of the behavior of a stationary gas bubble in a melt is relevant for an understanding of glass refining in space.

Weinberg, M. C.↗
NASA NTRS2003

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗
NASA NTRS1987

Discontinuous solutions to hyperbolic systems under operator splitting

Two-dimensional systems of linear hyperbolic equations are studied with regard to their behavior under a solution strategy that in alternate time-steps solves exactly the component one-dimensional operators. The initial data is a step function across an oblique discontinuity. The manner in which this discontinuity breaks up under repeated applications of the split operator is analyzed, and it is shown that the split solution will fail to match the true solution in any case where the two operators do not share all their eigenvectors. The special case of the fluid flow equations is analyzed in more detail, and it is shown that arbitrary initial data gives rise to pseudo acoustic waves and a non-physical stationary wave. The implications of these findings for the design of high-resolution computing schemes are discussed.

Roe, P. L.↗
NASA NTRS1959

On Fully Developed Channel Flows: Some Solutions and Limitations, and Effects of Compressibility, Variable Properties, and Body Forces

An examination of the effects of compressibility, variable properties, and body forces on fully developed laminar flow has indicated several limitations on such streams. In the absence of a pressure gradient, but presence of a body force (e.g., gravity), an exact fully developed gas flow results. For a liquid this follows also for the case of a constant streamwise pressure gradient. These motions are exact in the sense of a Couette flow. In the liquid case two solutions (not a new result) can occur for the same boundary conditions. An approximate analytic solution was found which agrees closely with machine calculations.In the case of approximately exact flows, it turns out that for large temperature variations across the channel the effects of convection (due to, say, a wall temperature gradient) and frictional heating must be negligible. In such a case the energy and momentum equations are separated, and the solutions are readily obtained. If the temperature variations are small, then both convection effects and frictional heating can consistently be considered. This case becomes the constant-property incompressible case (or quasi-incompressible case for free-convection flows) considered by many authors. Finally there is a brief discussion of cases wherein streamwise variations of all quantities are allowed but only a such form that independent variables are separable. For the case where the streamwise velocity varies inversely as the square root distance along the channel a solution is given.

Maslen, Stephen H.↗
NASA NTRS1992

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗
NASA NTRS1992

Conservative-variable average states for equilibrium gas multi-dimensional fluxes

Modern split component evaluations of the flux vector Jacobians are thoroughly analyzed for equilibrium-gas average-state determinations. It is shown that all such derivations satisfy a fundamental eigenvalue consistency theorem. A conservative-variable average state is then developed for arbitrary equilibrium-gas equations of state and curvilinear-coordinate fluxes. Original expressions for eigenvalues, sound speed, Mach number, and eigenvectors are then determined for a general average Jacobian, and it is shown that the average eigenvalues, Mach number, and eigenvectors may not coincide with their classical pointwise counterparts. A general equilibrium-gas equation of state is then discussed for conservative-variable computational fluid dynamics (CFD) Euler formulations. The associated derivations lead to unique compatibility relations that constrain the pressure Jacobian derivatives. Thereafter, alternative forms for the pressure variation and average sound speed are developed in terms of two average pressure Jacobian derivatives. Significantly, no additional degree of freedom exists in the determination of these two average partial derivatives of pressure. Therefore, they are simultaneously computed exactly without any auxiliary relation, hence without any geometric solution projection or arbitrary scale factors. Several alternative formulations are then compared and key differences highlighted with emphasis on the determination of the pressure variation and average sound speed. The relevant underlying assumptions are identified, including some subtle approximations that are inherently employed in published average-state procedures. Finally, a representative test case is discussed for which an intrinsically exact average state is determined. This exact state is then compared with the predictions of recent methods, and their inherent approximations are appropriately quantified.

Iannelli, G. S.↗
NASA NTRS1986

An inverse method for subcritical flows

An exact method is presented for two-dimensional subsonic flow within the limitations of the tangent gas approximation. While Woods (1961) studied these equations and proposed iterative methods for the solution of both the iterative and inverse problems, the inverse method here presented is noniterative and exact. It is shown that the direct Euler solution over the designed airfoil is very close to the input speed distribution, and that the constraints necessitated by upstream condition and closure requirements are easily incorporated.

Daripa, P. K.↗
NASA NTRS1974

Multiple scattering from finite inhomogeneous media

Utilizing the characteristic information concerning the apparent phase constant difference between the electric and magnetic fields propagating in an inhomogeneous medium, a theoretical analysis of the multiple scattering of electromagnetic waves in finite inhomogeneous media is presented. The solution is obtained by first approximating the coefficients of a pair of exact coupled first-order differential equations and then solving the equations by first-order iteration. The present first-order approximate solution with multiple scattering considerations is shown to be more accurate than the WKB solution. Methods to improve the accuracy of the first-order solution further are discussed. Application of the solution to slowly varying finite media with periodic properties demonstrates the validity of the solution. The same approach can be extended to frequencies in the optical region by retaining additional terms in the coefficients of the coupled differential equations.

Tang, C. C. H.↗
NASA NTRS1985

Numerical solution of transonic wing flows using an Euler/Navier-Stokes zonal approach

Transonic flow fields about wing geometries are computed using an Euler/Navier-Stokes approach in which the flow field is divided into several zones. The grid zones immediately adjacent to the wing surface are suitably clustered and solved with the Navier-Stokes equations. Grid zones removed from the wing are less finely clustered and are solved with the Euler equations. Wind tunnel wall effects are easily and accurately modeled with the new grid-zoning algorithm because the wind tunnel grid is constructed as an exact subset of the corresponding free-air grid. Solutions are obtained that are in good agreement with experiment, including cases with significant wind tunnel wall effects and shock-induced separation on the upper wing surface.

Holst, T. L.↗
NASA NTRS1975

Spinning mode sound propagation in ducts with acoustic treatment

Recent acoustic data show larger noise attenuations than predicted for acoustically treated aircraft engine inlets without splitter rings. A theoretical study of the acoustic propagation of spinning modes in acoustically treated open circular ducts is presented, and a suppressor with splitter rings was modeled by using the rectangular approximation to the annular duct. Theoretical models were used to determine optimum impedance and maximum attenuation for several spinning lobe numbers from 0 to 50. Results of the analysis indicate that for circular ducts the maximum possible attenuation and the optimum wall impedance are strong functions of the lobe number, and for annular ducts the attenuation and optimum wall impedance are insensitive to the spinning lobe number for well cut-on modes. These results explain why suppressors with splitter rings were quite effective in spite of the lack of detailed information on the noise source modal structure. Conversely, effective use of outer wall treatment alone will require expanded knowledge of the noise source structure. Approximate solutions are presented to help interpret the more exact theoretical results.

Rice, E. J.↗
NASA NTRS1974

Spinning mode sound propagation in ducts with acoustic treatment

Recent acoustic data have shown larger noise attenuations than predicted for acoustically treated aircraft engine inlets without splitter rings. These data have stimulated a more detailed theoretical study of the acoustic propagation of spinning modes in acoustically treated open circular ducts. In addition, the suppressor with splitter rings was modeled by using the rectangular approximation to the annular duct. The theoretical models were used to determine optimum impedance and maximum attenuation for several spinning lobe numbers from 0 to 50. It is found that for circular ducts the maximum possible attenuation and the optimum wall impedance are strong functions of the lobe number. For annular ducts the attenuation and optimum wall impedance are insensitive to the spinning lobe number for well cut-on modes. The results help explain why suppressors with splitter rings have been quite effective in spite of the lack of detailed information on the noise-source modal structure. Conversely, effective use of outer-wall treatment alone will require expanded knowledge of the noise-source structure. Approximate solutions are presented to help interpret the more exact theoretical results.

Rice, E. J.↗
NASA NTRS1977

A numerical method for the exact calculation of airloads associated with impulsively started wings

A numerical method is developed to calculate three-dimensional potential flows due to the steady and impulsive motion of isolated wing and wing-wing interaction problems. The velocity potential is represented by a discrete set of constant-doublet quadrilaterals on wing and wake surfaces. The exact surface boundary condition is enforced, and the solution is obtained in a step-by-step fashion, configurations being impulsively started from rest. Free-wake geometries are generated for each time step with Rankine or Lamb viscous vortex segments used in wake-velocity calculations. Sample results include calculated performance to steady state for a thick wing and indicial lift of a wing-wing interaction problem.

Summa, J. M.↗
NASA NTRS1979

A study of the application of singular perturbation theory

A hierarchical real time algorithm for optimal three dimensional control of aircraft is described. Systematic methods are developed for real time computation of nonlinear feedback controls by means of singular perturbation theory. The results are applied to a six state, three control variable, point mass model of an F-4 aircraft. Nonlinear feedback laws are presented for computing the optimal control of throttle, bank angle, and angle of attack. Real Time capability is assessed on a TI 9900 microcomputer. The breakdown of the singular perturbation approximation near the terminal point is examined Continuation methods are examined to obtain exact optimal trajectories starting from the singular perturbation solutions.

Mehra, R. K.↗
NASA NTRS1978

Effect of noisy phase reference on coherent detection of unbalanced QPSK signals

A phase error in the reference signal used in coherent detection of PSK signals causes a degradation in performance; this effect on bit error probability has been much discussed in the literature for BPSK and balanced QPSK systems. For unbalanced QPSK (UQPSK) systems, where the power and data rates in the two channels are unequal, and where the channel data symbols might also be different, the problem of computing the error probability is more complex, and only approximate solutions have been previously given. In this paper, exact expressions are derived for the conditional probability of error for the I and Q channels as a function of the channel power and data rate ratio. This is done for all combinations of Manchester and NRZ signalling schemes, and it is shown that in some cases the low data rate channel can be more seriously affected by a phase offset. Numerical results are also given for the average probability of error, assuming a biphase Costas loop is used for carrier synchronization.

Osborne, H. C.↗
NASA NTRS1980

Approximate solution in gasdynamics

One-dimensional unsteady gasdynamics is considered. An approximation based mainly on the interaction of simple and entropy waves is adopted. A discussion supporting this approximation, based in part on shock expansion theory, is given. By the use of certain transformations the approximation leads to solution in terms of quadratures. Excellent agreement with exact numerical results is obtained over a wide range of cases.

Sirovich, L.↗