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 685 records · Page 38

An isoparametric finite-element Stokes flow test problem

Two-dimensional Stokes flow of an incompressible Newtonian fluid about a horizontal cylinder situated next to a moving wall is presented and suggested as a nontrivial benchmark problem. Both the primitive variable Galerkin finite-element and the exact solutions are generated on a sequence of refined meshes. Error norms and convergence rates are presented for the C(o) bilinear-velocity/C(-1) piecewise constant-pressure and the C(o) biquadratic-velocity/C(-1) linear-pressure quadrilateral elements.

Maslanik, Miriam K.↗

Nonlinear static and dynamic finite element analysis of an eccentrically loaded graphite-epoxy beam

The Dynamic Crash Analysis of Structures (DYCAT) and NIKE3D nonlinear finite element codes were used to model the static and implulsive response of an eccentrically loaded graphite-epoxy beam. A 48-ply unidirectional composite beam was tested under an eccentric axial compressive load until failure. This loading configuration was chosen to highlight the capabilities of two finite element codes for modeling a highly nonlinear, large deflection structural problem which has an exact solution. These codes are currently used to perform dynamic analyses of aircraft structures under impact loads to study crashworthiness and energy absorbing capabilities. Both beam and plate element models were developed to compare with the experimental data using the DYCAST and NIKE3D codes.

Fasanella, Edwin L.↗

Multiple vs single scattering - Assessment of magnitudes

All orders of scattering are analyzed for two artifical canopies. The SHL canopy consists of Small Horizontal Leaves that are much smaller than the leaf-to-leaf spacing. The IHL canopy consists of Infinite Horizontal Layers, where each leaf is of infinite extent (a horizontal plane). Hemispheric leaf reflectances and transmittances independent of the direction of illumination lead to exact solutions for these models. Sunlight that penetrates to a given leaf area index level is much stronger in an SHL canopy than that in IHL; but the difference becomes muted when leaf transmittance is large. Multiple scattering enhances the hemispheric canopy reflectance more strongly in SHL than it does in IHL. The enhancement depends linearly on leaf transmittance in SHL and on the transmittance squared in IHL. Comparison with measured reflectances indicates that IHL model grossly underestimates multiple scattering in soybean canopies.

Otterman, Joseph↗

Verification of the Proteus two-dimensional Navier-Stokes code for flat plate and pipe flows

The Proteus Navier-Stokes Code is evaluated for two-dimensional/axisymmetric, viscous, incompressible, internal and external flows. The particular cases to be discussed are laminar and turbulent flows over a flat plate, laminar and turbulent dveloping pipe flows and turbulent pipe flow with swirl. Results are compared with exact solutions, empirical correlations and experimental data. A detailed description of the code set-up, including boundary conditions, intitial conditions, grid size and grid packing is given for each case.

Conley, Julianne M.↗

An inverse method with regularity condition for transonic airfoil design

It is known from Lighthill's exact solution of the incompressible inverse problem that in the inverse design problem, the surface pressure distribution and the free stream speed cannot both be prescribed independently. This implies the existence of a constraint on the prescribed pressure distribution. The same constraint exists at compressible speeds. Presented here is an inverse design method for transonic airfoils. In this method, the target pressure distribution contains a free parameter that is adjusted during the computation to satisfy the regularity condition. Some design results are presented in order to demonstrate the capabilities of the method.

Zhu, Ziqiang↗

Pressure wave propagation studies for oscillating cascades

The unsteady flow field around an oscillating cascade of flat plates is studied using a time marching Euler code. Exact solutions based on linear theory serve as model problems to study pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step is demonstrated. Results show that an approximate non-reflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer than cases using reflective boundary conditions. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates does a better job of capturing the unsteady flow field than cases using uniform grids as long as the CFL number is less than one for a sufficient portion of the grid. Results for various stagger angles and oscillation frequencies show good agreement with linear theory as long as the grid is properly resolved.

Huff, Dennis L.↗

Application of program LAURA to perfect gas shock tube flows: A parametric study

The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) was originally developed to solve steady-flow problems. The desire to validate the algorithm with shock tube experimental data motivated the development of a time-accurate version of the LAURA code. The current work presents a test of the Algorithm. Computational results are compared with the exact solution for a simple shock tube case. The parameters examined are Courant number, relaxation sweeps, grid spacing, and the inviscid relaxation factor. The results of the study indicate that LAURA is capable of producing accurate solutions when appropriate values are used for each parameter.

Mitterer, K. F.↗

Unsteady-flow-field predictions for oscillating cascades

The unsteady flow field around an oscillating cascade of flat plates with zero stagger was studied by using a time marching Euler code. This case had an exact solution based on linear theory and served as a model problem for studying pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step size was shown for a moderate reduced frequency. Results show that an approximate nonreflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer to the airfoils than when reflective boundaries are used. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates captures the unsteady flow field better than when uniform grids are used as long as the 'Courant Friedrichs Levy' (CFL) number is less than 1 for a sufficient portion of the grid. Finally, a solution based on an optimization of grid, CFL number, and boundary conditions shows good agreement with linear theory.

Huff, Dennis L.↗

The Proteus Navier-Stokes code

An effort is currently underway at NASA Lewis to develop two and three dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. Proteus solves the Reynolds-averaged, unsteady, compressible Navier-Stokes equations in strong conservation law form. Turbulence is modeled using a Baldwin-Lomax based algebraic eddy viscosity model. In addition, options are available to solve thin layer or Euler equations, and to eliminate the energy equation by assuming constant stagnation enthalpy. An extensive series of validation cases have been run, primarily using the two dimensional planar/axisymmetric version of the code. Several flows were computed that have exact solution such as: fully developed channel and pipe flow; Couette flow with and without pressure gradients; unsteady Couette flow formation; flow near a suddenly accelerated flat plate; flow between concentric rotating cylinders; and flow near a rotating disk. The two dimensional version of the Proteus code has been released, and the three dimensional code is scheduled for release in late 1991.

Towne, Charles E.↗

Electromagnetic scattering by coated convex surfaces and wedges simulated by approximate boundary conditions

Asymptotic/high-frequency solutions are developed for analyzing the non-specular scattering mechanisms associated with coated convex surfaces and edges simulated by approximate boundary conditions. In particular, the standard impedance boundary conditions (SIBC's) and the second order generalized impedance boundary conditions (GIBC's) are employed for a characterization of the edge diffraction, creeping wave, and surface diffracted wave contributions. To study the creeping wave and surface diffracted wave mechanisms, rigorous UTD (uniform geometrical theory of diffraction) diffraction coefficients are developed for a convex coated cylinder simulated with SIBC's and GIBC's. The ray solutions obtained remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficient for a convex cylinder are obtained via a generalization of the corresponding ones of the circular cylinder. To validate the asymptotic/high-frequency solution, integral equations are derived for both E and H-polarization and solved numerically using the method of moments. Results are presented for a single and three layered coated convex cylinder. Some insights are also provided on the accuracy of the employed GIBC's versus SIBC's for application to curved surfaces. To characterize the scattering by impedance wedges illuminated at skew incidence, diffraction coefficients are derived from an approximate solution of the governing functional difference equations. This solution exactly recovers the known ones for an impedance half plane or an arbitrary wedge at normal incidence, and to validate it for other wedge angles, a moment method code was used. Finally, to test the usefulness of the approximate skew incidence impedance wedge diffraction coefficient for three dimensional structures, equivalent currents are derived in the context of the physical theory of diffraction (PTD) for a finite length impedance wedge of arbitrary internal angles. These are incorporated in a standard general purpose PTD code and results are presented for a number of different impedance structures.

Syed, H. H.↗

Effect of an arcjet plume on satellite reflector performance

The effect of an arcjet plume on the performance of satellite reflector antennas is studied. The arcjet plume is modeled as a weakly ionized plasma. The spatial permittivity distribution of the plume is approximated using the measured electron density profile and a cold plasma model. Geometrical optics is applied to determine the ray paths as well as the transmitted fields through the inhomogeneous plume. The ray optics results are compared against several exact solutions for scattering from inhomogeneous dielectrics, and good agreement is observed for sufficiently large scatterer size. The far-field antenna patterns of the reflector in the presence of the plume are calculated from the transmitted ray fields using a ray-tube integration scheme. For arcjet prototypes in the 1-kW class, the plume effect on the antenna performance is small. As the electron density increases, the main beam and sidelobe level gradually degrade. The main beam also tends to squint away from the plume region.

Ling, Hao↗

Boundary formulations for shape sensitivity of temperature dependent conductivity problems

Used in concert with the Kirchhoff transformation, implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solids with temperature dependent thermal conductivity is shown to generate an accurate and economical approach for computation of shape sensitivities. For problems with specified temperature and heat flux boundary conditions, a linear problem results for both the analysis and sensitivity analysis. In problems with either convection or radiation boundary conditions, a nonlinear problem is generated. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems is presented that have closed form solutions. Exact analytical expressions for the shape sensitivities associated with these problems are developed and these are compared with the sensitivities computed using the boundary element formulation.

Kane, James H.↗

Pressure wave propagation studies for oscillating cascades

The unsteady flowfield around an oscillating cascade of flat plates is studied using a time marching Euler code. Exact solutions based on linear theory serve as model problems to study pressure wave propagation in the numerical solution. The importance of using proper unsteady boundary conditions, grid resolution, and time step is demonstrated. Results show that an approximate non-reflecting boundary condition based on linear theory does a good job of minimizing reflections from the inflow and outflow boundaries and allows the placement of the boundaries to be closer than cases using reflective boundary conditions. Stretching the boundary to dampen the unsteady waves is another way to minimize reflections. Grid clustering near the plates does a better job of capturing the unsteady flowfield than cases using uniform grids as long as the CFL number is less than one for a sufficient portion of the grid. Results for various stagger angles and oscillation frequencies show good agreement with linear theory as long as the grid is properly resolved.

Huff, Dennis L.↗

Theoretical model for a thin cylindrical film optical fiber fluorosensor

The analytical treatment of power efficiency (P(eff) is undertaken for the case of a positively guiding optical fiber with a thin-film source distributed in the core-cladding interface. The approach adopts the exact solution of the cylindrical optical fiber with an infinite cladding to account for differences between the indices of refraction of the core and the cladding. The excitation of low-loss leaky modes by the cladding is ignored, and only the injection by the evanescent field is considered. The formulas permit the analysis of the power-injection efficiency of fibers with arbitrary differences in indices of refraction. P(eff) does not always increase with V number, but rather varies slightly with wavelength and fiber-core radius and varies significantly with the difference in the indices of refraction. The theoretical results of the work are of interest for designing an atomic-O chemical sensor based on evanescent-wave coupling.

Egalon, Claudio O.↗

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗

Fuel consumption in optimal control

A method has been developed for comparing three optimal control strategies based on fuel consumption. A general cost function minimization procedure was developed by applying two theorems associated with convex sets. Three cost functions associated with control saturation, pseudofuel, and absolute fuel are introduced and minimized. The first two cost functions led to the bang-bang and continuous control strategies, and the minimization of absolute fuel led to an impulsive strategy. The three control strategies were implemented on two elementary systems and a comparison of fuel consumption was made. The impulse control strategy consumes significantly less fuel than the continuous and bang-bang control strategies. This comparison suggests a potential for fuel savings in higher-order systems using impulsive control strategies. However, since exact solutions to fuel-optimal control for large-order systems are difficult if not impossible to achieve, the alternative is to develop near-optimal control strategies.

Redmond, Jim↗

Ultrasonic characterization of multilayered thick composite plates

Guided wave propagation in a periodically laminated plate is considered when each lamina is transversely isotropic with its symmetry axis lying in the plane of the graphite-epoxy plate. The dispersion behavior of propagation modes influences the dynamic response and can be used to ultrasonically characterize the mechanical properties of the laminated plate. It is assumed that fibers are oriented at 90 deg to one another in adjacent layers and all the laminae are of equal thickness. Results obtained are in good agreement with exact solutions of previous studies.

Karunasena, W.↗

Effect of contact stresses in four-point bend testing of graphite/epoxy and graphite/PMR-15 composite beams

The results of in-plane four-point bend experiments on unidirectionally reinforced composite beams are presented for graphite/epoxy (T300/934) and graphite/polyimide (G30-500/PMR-15) composites. The maximum load and the location of cracks formed during failure were measured for testpieces with fibers oriented at various angles to the beam axis. Since most of the beams failed near one or more of the load points, the strength of the beams was evaluated in terms of a proposed model for the local stress distribution. In this model, an exact solution to the problem of a localized contact force acting on a unidirectionally reinforced half plane is used to describe the local stress field. The stress singularity at the load points is treated in a manner similar to the stress singularity at a crack tip in fracture mechanisms problems. Using this approach, the effect of fiber angle and elastic material properties on the strength of the beam is described in terms of a load intensity factor. For fiber angles less than 45 deg from the beam axis, a single crack is initiated near one of the load points at a critical value of the load intensity factor. The critical load intensity factor decreases with increasing fiber angle. For larger fiber angles, multiple cracks occur at locations both near and away from the load points, and the load intensity factor at failure increases sharply with increasing fiber angle.

Binienda, W. K.↗