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,009 records · Page 56

Electromagnetic Scattering by a Morphologically Complex Object: Fundamental Concepts and Common Misconceptions

Following Keller(Proc Symp Appl Math 1962;13:227:46), we classify all theoretical treatments of electromagnetic scattering by a morphologically complex object into first- principle (or "honest" in Keller s terminology) and phenomenological (or "dishonest") categories. This helps us identify, analyze, and dispel several profound misconceptions widespread in the discipline of electromagnetic scattering by solitary particles and discrete random media. Our goal is not to call for a complete renunciation of phenomenological approaches but rather to encourage a critical and careful evaluation of their actual origin, virtues, and limitations. In other words, we do not intend to deter creative thinking in terms of phenomenological short-cuts, but we do want to raise awareness when we stray (often for practical reasons) from the fundamentals. The main results and conclusions are illustrated by numerically-exact data based on direct numerical solutions of the macroscopic Maxwell equations.

Mischenko, Michael I.↗

Exact image method for Gaussian beam problems involving a planar interface

Exact image method, recently introduced for the solution of electromagnetic field problems involving sources above a planar interface or two homogeneous media, is shown to be valid also for sources located in complex space, which makes its application possible for Gaussian beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular shift of a TE polarized beam are correctly given as asymptotic results by the exact reflection image theory. Also, the apparent image location giving the correct Gaussian beam transmitted through the interface is obtained as another asymptotic check. The present theory makes it possible to calculate the exact coupling from the Gaussian beam to the reflected and refracted beams, as well as to the surface wave.

Lindell, I. V.↗

Exact-image method for Gaussian-beam problems involving a planar interface

The exact-image method, recently introduced for the solution of electromagnetic field problems involving sources above a planar interface between two homogeneous media, is shown to be valid also for sources located in complex space, which makes its application possible for Gaussian-beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular shift of a TE-polarized beam are correctly given as asymptotic results by the exact-reflection-image theory. Also, the apparent-image location giving the correct Gaussian beam transmitted through the interface is obtained as another asymptotic check. The theory described here makes it possible to calculate the exact coupling from the Gaussian beam to the reflected and refracted beams as well as to the surface wave.

Lindell, Ismo V.↗

Adaptively Refined Euler and Navier-Stokes Solutions with a Cartesian-Cell Based Scheme

A Cartesian-cell based scheme with adaptive mesh refinement for solving the Euler and Navier-Stokes equations in two dimensions has been developed and tested. Grids about geometrically complicated bodies were generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells were created using polygon-clipping algorithms. The grid was stored in a binary-tree data structure which provided a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations were solved on the resulting grids using an upwind, finite-volume formulation. The inviscid fluxes were found in an upwinded manner using a linear reconstruction of the cell primitives, providing the input states to an approximate Riemann solver. The viscous fluxes were formed using a Green-Gauss type of reconstruction upon a co-volume surrounding the cell interface. Data at the vertices of this co-volume were found in a linearly K-exact manner, which ensured linear K-exactness of the gradients. Adaptively-refined solutions for the inviscid flow about a four-element airfoil (test case 3) were compared to theory. Laminar, adaptively-refined solutions were compared to accepted computational, experimental and theoretical results.

Coirier, William J.↗

Direct handling of equality constraints in multilevel optimization

In recent years there have been several hierarchic multilevel optimization algorithms proposed and implemented in design studies. Equality constraints are often imposed between levels in these multilevel optimizations to maintain system and subsystem variable continuity. Equality constraints of this nature will be referred to as coupling equality constraints. In many implementation studies these coupling equality constraints have been handled indirectly. This indirect handling has been accomplished using the coupling equality constraints' explicit functional relations to eliminate design variables (generally at the subsystem level), with the resulting optimization taking place in a reduced design space. In one multilevel optimization study where the coupling equality constraints were handled directly, the researchers encountered numerical difficulties which prevented their multilevel optimization from reaching the same minimum found in conventional single level solutions. The researchers did not explain the exact nature of the numerical difficulties other than to associate them with the direct handling of the coupling equality constraints. The coupling equality constraints are handled directly, by employing the Generalized Reduced Gradient (GRG) method as the optimizer within a multilevel linear decomposition scheme based on the Sobieski hierarchic algorithm. Two engineering design examples are solved using this approach. The results show that the direct handling of coupling equality constraints in a multilevel optimization does not introduce any problems when the GRG method is employed as the internal optimizer. The optimums achieved are comparable to those achieved in single level solutions and in multilevel studies where the equality constraints have been handled indirectly.

Renaud, John E.↗

Nonplanar Method for Predicting Incompressible Aerodynamic Coefficients of Rectangular Wings with Circular-Arc Camber

The development of a nonplanar lifting surface method having a continuous distribution of singularities and satisfying the tangent flow boundary condition on the mean camber surface is given. The method predicts some incompressible longitudinal aerodynamic coefficients of rectangular wings which have circular-arc camber. The solution method is of the integral-equation type and the resulting surface integrals are evaluated by either using numerical or analytical techniques, as are appropriate. Applications are made and the results compared with those from an exact two-dimensional circular-arc camber solution, a three-dimensional flat-wing solution which represents the camber by a projected slope onto the flat surface, and a flat-wing experiment. From these comparisons, the present method is found to predict well the flat-wing experiment and limiting values, in addition to the center of pressure variation at an angle of attack of zero for any camber. For wings having camber ratios larger than about 1.25% and moderate to high aspect ratios, the results deterioriate due to the inadequacy of lifting pressure modes employed.

Lamar, J. E.↗

Synthesis of offset dual reflector antennas transforming a given feed illumination pattern into a specified aperture distribution

The problem of transforming a given primary feed pattern into a desired aperture field distribution through two reflections by an offset dual reflector system is investigated using the concepts of geometrical optics. A numerically rigorous solution for the reflector surfaces is developed which realizes an exact aperture phase distribution and an aperture amplitude distribution that is accurate to within an arbitrarily small numerical tolerance. However, this procedure does not always yield a smooth solution, i.e., the reflector surfaces thus realized may not be continuous or their slopes may vary too rapidly. In the event of nonexistence of a numerically rigorous smooth solution, an approximate solution that enforces the smoothness of the reflector surfaces can be obtained. In the approximate solution, only the requirement for the aperture amplitude distribution is relaxed, and the condition on the aperture phase distribution is continued to be satisfied exactly.

Mittra, R.↗

At-Least Version of the Generalized Minimum Spanning Tree Problem: Optimization Through Ant Colony System and Genetic Algorithms

The At-Least version of the Generalized Minimum Spanning Tree Problem (L-GMST) is a problem in which the optimal solution connects all defined clusters of nodes in a given network at a minimum cost. The L-GMST is NPHard; therefore, metaheuristic algorithms have been used to find reasonable solutions to the problem as opposed to computationally feasible exact algorithms, which many believe do not exist for such a problem. One such metaheuristic uses a swarm-intelligent Ant Colony System (ACS) algorithm, in which agents converge on a solution through the weighing of local heuristics, such as the shortest available path and the number of agents that recently used a given path. However, in a network using a solution derived from the ACS algorithm, some nodes may move around to different clusters and cause small changes in the network makeup. Rerunning the algorithm from the start would be somewhat inefficient due to the significance of the changes, so a genetic algorithm based on the top few solutions found in the ACS algorithm is proposed to quickly and efficiently adapt the network to these small changes.

Janich, Karl W.↗

Nonlinear analysis of a shock-loaded membrane.

Results from a computer method for analyzing the unsteady interaction of a fluid stream and a flat circular elastic membrane are presented. The loading on the membrane is assumed to be caused by the firing of a shock tube. The fluid pressures and velocities are determined from a scheme based on the numerical method of characteristics, and the membrane is analyzed using exact relations for membrane strain. The interactive solution is found to give peak stresses 40% lower than a solution which assumes a pressure invariant in space and time.

Madden, R.↗

Relaxation solutions for inviscid axisymmetric transonic flow over blunt or pointed bodies.

A finite-difference relaxation method is presented for numerical solution of the full potential equation and exact boundary conditions for general axisymmetric bodies is inviscid, steady transonic flow. Body-normal coordinates are used in the nose region and sheared cylindrical coordinates are used on the afterbody to accommodate corners such as boattails and flares. An improved difference scheme is used which does not require that the flow be nearly alined with a coordinate direction in supersonic regions, and which treats either subsonic or supersonic free streams. Numerical results are illustrated for some simple classical shapes such as spheres and ellipsoids, and for more practical shapes like tangent-ogives with boattails. Special attention is given to bodies which have been studied for area-rule applications. Agreement with available experimental results is good in cases where viscous effects and wind-tunnel wall interference are not important.

South, J. C., Jr.↗

Magnetic attitude control system for dual-spin satellites

A closed-loop control law is developed for a dual-spin satellite control system which utilizes the interaction of the geomagnetic field with the satellite dipole parallel to the spin axis. The control law consists of the linear combination of the pitch axis component of the rate of change of the geomagnetic field and the product of the roll angle and roll axis component of the geomagnetic field. Application of the method of multiple time scales yields approximate solutions for the feedback gains in terms of the system parameters. Approximate solutions are also obtained for the response of the system to disturbance torques. A comparison of the approximate solutions and numerical solutions obtained by numerical integration of the exact equations of motion is then given.

Alfriend, K. T.↗

Numerical computations of Orbiter flow fields and heating rates

Numerical computations of flow fields around an analytical description of the Space Shuttle Orbiter windward surface, including the root of the wing leading edge, are presented to illustrate the sensitivity of these calculations to several flow field modeling assumptions. Results of parametric flow field and boundary layer computations using the axisymmetric analogue concept to obtain three-dimensional heating rates, in conjunction with exact three-dimensional inviscid floe field solutions and two-dimensional boundary layer analysis - show the sensitivity of boundary layer edge conditions and heating rates to considerations of the inviscid flow field entropy layer, equilibrium air versus chemically and vibrationally frozen flow, and nonsimilar terms in the boundary layer computations. A cursory comparison between flow field predictions obtained from these methods and current Orbiter design methods has established a benchmark for selecting and adjusting these and future design methodologies.

Goodrich, W. D.↗

Beam area determination for multiple-beam satellite communication applications

An exact analytical model and a simple numerical solution are presented for determining the surface area of the earth illuminated by a conical beam, with an elliptical cross section, emanating from a satellite. The model is used to generate sets of parametric curves for determination of the beam area as a function of the different satellite locations, beam locations, and beamwidths. Numerical results show that for a circular beam with a 0.25 deg beamwidth, the beam illuminating the northern Continental United States (CONUS) can cover an area about five times larger than the area covered by the same size beam illuminating the southern CONUS. Comparisons are also made with the available approximate methods.

Salmasi, A. B.↗

Small-amplitude viscous motion on arbitrary potential flows

This paper is concerned with small-amplitude, unsteady, vortical and entropic motion imposed on steady potential flows. It is restricted to the case where the spatial scale of the unsteady motion is small compared to that of the mean flow. Under such conditions, the unsteady motion may be influenced by viscosity even if the mean flow is not. An exact high-frequency (small-wavelength) solution is obtained for the small-amplitude viscous motion imposed on a steady potential flow. It generalizes the one obtained by Pearson (1959) for the homogeneous-strain case to the case of quasi-homogeneous strain. This result is used to study the effect of viscosity on rapidly distorted turbulent flows. Specific numerical results are given for a turbulent flow near a two-dimensional stagnation point.

Goldstein, M. E.↗

Factors influencing the ultrasonic stress wave factor evaluation of composite material structures

To demonstrate that the finite-element model can be used to investigate some of the factors influencing the ultrasonic stress wave evaluation of materials, a hypothetical case was studied in which classical vibration theory was used. Vibration analysis and experiments for the undamaged case were conducted on an isotropic aluminum plate and a unidirectional graphite-epoxy plate, using a point source to excite the plates. The finite-element solution correlated within eight percent with the exact method. The frequencies predicted by the finite-element model were observed in the experiments in both plates, although in the composite plate, additional frequencies were observed which could not be accounted for. Damaged isotropic plates were also considered. The effects of increasing damage severity with constant damage area, and increasing damage area with constant severity on the resonant frequencies were analyzed.

Rebello, C. J.↗

Developing disturbances in shear

The origin and growth to moderate amplitude of disturbances in shear flow has been traditionally ascribed to the linear modal instability of the flow. Recent work on initial value problems has suggested that nonmodal growth of perturbations may be of equal and perhaps greater importance, at least in cases of rapid development. Examples of robust growth in model problems which support no instabilities and in baroclinic flows with realistic Ekman damping for which the exponential modes have zero or negative growth rates have been shown. Such examples have focused attention on the perturbations which are configured to tap the energy of mean flows. Here a critical examination of these favorably configured perturbations is given, making use of the simple constant free shear barotropic model which allows construction of exact two-dimensional isolated wave packet solutions.

Farrell, Brian↗

Investigations of the dynamical evolution of protoplanetary nebulae

The spectral energy distributions (SED's) of all the known Class I (protostellar) sources in the Taurus molecular cloud were modelled. The Tereby, Shu, & Cassen (1984, TSC) density distribution for a rotating, infalling envelope was adopted. The radiative equilibrium temperature distribution from the spherical average of the TSC density distribution was calculated. The resulting spherically-symmetric temperature distribution then provides the source function to obtain the emergent spectrum at a given inclination angle i from the formal solution of the transfer equation, using the exact density (opacity) distribution. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flat Spectrum T Tauri Stars were studied. The mid- to far-infrared fluxes of 'flat spectrum' T Tauri stars can be explained by radiative equilibrium emission from infalling dusty envelopes. This explanation indicated that models employing 'active' disks, in which the temperature distribution is a parameterized power law, should be invoked with caution. Infall also naturally explains the scattered light nebulae detected around many flat-spectrum sources. Results showed that the SED's of the protostar candidates in Taurus can be reproduced with TSC models having infall rates close to the values predicted by the theory of isothermal cloud collapse. Flattened infalling envelope models are also being studied.

Hartmann, Lee W.↗