SEARCH · Engineering Papers
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.
In this work, a new optimized MacCormack-type scheme, which is 4th order accurate in time and space, is applied to Problems 1 and 2 of Category 1. The performance of this new scheme is compared to that of the 2-4 MacCormack scheme, and results for Problems 1 and 2 of Category 1 are presented and compared to the exact solutions.
An implementation of the Model Based Parameter Estimation (MBPE) technique is presented for obtaining the frequency response of the Radar Cross Section (RCS) of arbitrarily shaped, three-dimensional perfect electric conductor (PEC) bodies. An Electric Field Integral Equation (EFTE) is solved using the Method of Moments (MoM) to compute the RCS. The electric current is expanded in a rational function and the coefficients of the rational function are obtained using the frequency derivatives of the EFIE. Using the rational function, the electric current on the PEC body is obtained over a frequency band. Using the electric current at different frequencies, RCS of the PEC body is obtained over a wide frequency band. Numerical results for a square plate, a cube, and a sphere are presented over a bandwidth. Good agreement between MBPE and the exact solution over the bandwidth is observed.
During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.
Smart composite structures with distributed sensors and actuators have the capability to actively respond to a changing environment while offering significant weight savings and additional passive controllability through ply tailoring. Piezoelectric sensing and actuation of composite laminates is the most promising concept due to the static and dynamic control capabilities. Essential to the implementation of these smart composites are the development of accurate and efficient modeling techniques and experimental validation. This research addresses each of these important topics. A refined higher order theory is developed to model composite structures with surface bonded or embedded piezoelectric transducers. These transducers are used as both sensors and actuators for closed loop control. The theory accurately captures the transverse shear deformation through the thickness of the smart composite laminate while satisfying stress free boundary conditions on the free surfaces. The theory is extended to include the effect of debonding at the actuator-laminate interface. The developed analytical model is implemented using the finite element method utilizing an induced strain approach for computational efficiency. This allows general laminate geometries and boundary conditions to be analyzed. The state space control equations are developed to allow flexibility in the design of the control system. Circuit concepts are also discussed. Static and dynamic results of smart composite structures, obtained using the higher order theory, are correlated with available analytical data. Comparisons, including debonded laminates, are also made with a general purpose finite element code and available experimental data. Overall, very good agreement is observed. Convergence of the finite element implementation of the higher order theory is shown with exact solutions. Additional results demonstrate the utility of the developed theory to study piezoelectric actuation of composite laminates with pre-existing debonding. Significant changes in the modes shapes and reductions in the control authority result due to partially debonded actuators. An experimental investigation addresses practical issues, such as circuit design and implementation, associated with piezoelectric sensing and actuation of composite laminates. Composite specimens with piezoelectric transducers were designed, constructed and tested to validate the higher order theory. These specimens were tested with various stacking sequences, debonding lengths and gains for both open and closed loop cases. Frequency changes of 15% and damping on the order of more than 20% of critical damping, via closed loop control, was achieved. Correlation with the higher order theory is very good. Debonding is shown to adversely affect the open and closed loop frequencies, damping ratios, settling time and control authority.
Excimer laser photolysis in combination with time-resolved IR laser absorption detection of OH radicals has been used to study O3/OH(v = 0)/HO2 chain reaction kinetics at 298 K, (i.e.,(k(sub 1) is OH + 03 yields H02 + 02 and (k(sub 2) is H02 + 03 yields OH + 202). From time-resolved detection of OH radicals with high-resolution near IR laser absorption methods, the chain induction kinetics have been measured at up to an order of magnitude higher ozone concentrations ([03] less than or equal to 10(exp 17) molecules/cu cm) than accessible in previous studies. This greater dynamic range permits the full evolution of the chain induction, propagation, and termination process to be temporally isolated and measured in real time. An exact solution for time-dependent OH evolution under pseudo- first-order chain reaction conditions is presented, which correctly predicts new kinetic signatures not included in previous OH + 03 kinetic analyses. Specifically, the solutions predict an initial exponential loss (chain "induction") of the OH radical to a steady-state level ([OH](sub ss)), with this fast initial decay determined by the sum of both chain rate constants, k(sub ind) = k(sub 1) + k(sub 2). By monitoring the chain induction feature, this sum of the rate constants is determined to be k(sub ind) = 8.4(8) x 10(exp -14) cu cm/molecule/s for room temperature reagents. This is significantly higher than the values currently recommended for use in atmospheric models, but in excellent agreement with previous results from Ravishankara et al.
This paper deals with the buckling of the stiffened plate under uni-axial compression. The direct integration of the governing differential equation is performed and the exact solution to the problem is obtained. As examples, a square plate with single stiffener, and a stiffened three-span, continuous plate are investigated, with special attention given to the influence of stiffener misplacement on the buckling load and mode shape of the plate. It is found that a small misplacement of the stiffeners from the nominal configuration may change the buckling mode from a global one to a highly localized one.
Several algorithms to calculate a rain-rate profile from a single-frequency air-or spaceborne radar backscatter profile and a given path-integrated attenuation have been proposed. The accuracy of any such algorithm is limited by the ambiguities between the (multiple) exact solutions, which depend on the variability of the parameters in the Z-R and k-R relations used. In this study, coupled Z-R and k-R relations are derived based on the drop size distribution. It is then shown that, because of the coupling, the relative difference between the multiple mutually ambiguous rain-rate profiles solving the problem must remain acceptably low, provided the available path-integrated attenuation value is known to within 0.5 dB.
In this paper a new class of formations that maintain a constant shape as viewed from the Earth is introduced. An algorithm is developed to place n spacecraft in a constant shape formation spaced equally in time using the classical orbital elements. To first order, the dimensions of the formation are shown to be simple functions of orbit eccentricity and inclination. The performance of the formation is investigated over a Keplerian orbit using a performance measure based on a weighted average of the angular separations between spacecraft in formation. Analytic approximations are developed that yield optimum configurations for different values of n. The analytic approximations are shown to be in excellent agreement with the exact solutions.
A new, thermal-flow simulation code, called SFLOW. has been developed to model the gas dynamics, heat transfer, as well as O-ring and flow path erosion inside the space shuttle solid rocket motor joints by combining SINDA/Glo, a commercial thermal analyzer. and SHARPO, a general-purpose CFD code developed at Thiokol Propulsion. SHARP was modified so that friction, heat transfer, mass addition, as well as minor losses in one-dimensional flow can be taken into account. The pressure, temperature and velocity of the combustion gas in the leak paths are calculated in SHARP by solving the time-dependent Navier-Stokes equations while the heat conduction in the solid is modeled by SINDA/G. The two codes are coupled by the heat flux at the solid-gas interface. A few test cases are presented and the results from SFLOW agree very well with the exact solutions or experimental data. These cases include Fanno flow where friction is important, Rayleigh flow where heat transfer between gas and solid is important, flow with mass addition due to the erosion of the solid wall, a transient volume venting process, as well as some transient one-dimensional flows with analytical solutions. In addition, SFLOW is applied to model the RSRM nozzle joint 4 subscale hot-flow tests and the predicted pressures, temperatures (both gas and solid), and O-ring erosions agree well with the experimental data. It was also found that the heat transfer between gas and solid has a major effect on the pressures and temperatures of the fill bottles in the RSRM nozzle joint 4 configuration No. 8 test.
This work presents results from the study of a large volume of spacecraft flight data pertaining to gyroscope performance. We have examined long and short term trends of gyroscope biase of first, an exact solution for the time-dependence of the attitude part of the state error covariance, averaged over the three spacecraft axes, is given. This solution is more complete than the usual cubic polynomial for the variance in that it includes coupling among the three gyroscope axes and includes an important term arising from the initial correlations. Second, several continuous 24-hour spans of gyroscope data are examined to verify the short-term statistical model. This analysis demonstrates that in-flight data can be used to determine the strength of the white noise driving the random walk of the gyroscope bias. This may be useful for postlaunch improvement to the noise model and for diagnosing the health of the gyroscope. Third, the long-term trends in gyroscope biases show a nearly linear systematic variation over time scales of years. This has been found on three different missions. While the random walk model is adequate as a basis for onboard Kalman filters or for state estimation using relatively short time spans, these trends indicate that some applications could benefit by accounting for the secular changes in the biases. One example is a new gyroscope calibration method that is under development that allows for multi-epoch bias solutions.
This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques. we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.
This article considers a posteriori error estimation of specified functionals for first-order systems of conservation laws discretized using the discontinuous Galerkin (DG) finite element method. Using duality techniques, we derive exact error representation formulas for both linear and nonlinear functionals given an associated bilinear or nonlinear variational form. Weighted residual approximations of the exact error representation formula are then proposed and numerically evaluated for Ringleb flow, an exact solution of the 2-D Euler equations.
A numerical-analytic investigation of the impacts of vertical transport schemes on the model simulated age-of-air and the so-called 'tape recorder' will be presented using an idealized 1-D column transport model as well as a more realistic 3-D dynamical model. By comparing to the 'exact' solutions of 'age-of-air' and the 'tape recorder' obtainable in the 1-D setting, useful insight is gained on the impacts of numerical diffusion and dispersion of numerical schemes used in global models. Advantages and disadvantages of Eulerian, semi-Lagrangian, and Lagrangian transport schemes will be discussed. Vertical resolution requirement for numerical schemes as well as observing systems for capturing the fine details of the 'tape recorder' or any upward propagating wave-like structures can potentially be derived from the 1-D analytic model.
An axisymmetric Meshless Local Petrov-Galerkin (MLPG) algorithm is presented for the potential and elasticity problems. In this algorithm the trial and test functions are chosen from different spaces. By a judicious choice of these functions, the integrals involved in the weak form can be restricted to a local neighborhood. This makes the method truly meshless. The MLPG algorithm is used to study various potential and elasticity problems for which exact solutions are available. The sensitivity and effectiveness of the MLPG algorithm to various parameters such as the weight functions, basis functions and support domain radius, etc. was studied. The MLPG algorithm yielded accurate solutions for all weight functions, basis functions and support domain radii considered for all of the problems studied.
The proceedings of the Third Computational Aeroacoustics (CAA) Workshop on Benchmark Problems cosponsored by the Ohio Aerospace Institute and the NASA Glenn Research Center are the subject of this report. Fan noise was the chosen theme for this workshop with representative problems encompassing four of the six benchmark problem categories. The other two categories were related to jet noise and cavity noise. For the first time in this series of workshops, the computational results for the cavity noise problem were compared to experimental data. All the other problems had exact solutions, which are included in this report. The Workshop included a panel discussion by representatives of industry. The participants gave their views on the status of applying computational aeroacoustics to solve practical industry related problems and what issues need to be addressed to make CAA a robust design tool.
We describe all extension of the Markov decision process model in which a continuous time dimension is included ill the state space. This allows for the representation and exact solution of a wide range of problems in which transitions or rewards vary over time. We examine problems based on route planning with public transportation and telescope observation scheduling.
The case of isotropic compressible turbulence subjected to rapid isotropic compression is studied using inviscid rapid distortion theory and direct numerical simulation. An exact solution to the rapid distortion problem is given, and results are compared to those of direct numerical simulation. Implications for modelling turbulent flows are discussed.
This report documents the results of the numerical simulations of unsteady aerodynamic models. The results focus on numerical accuracy and efficiency, and the robustness of the numerical methods. The aerodynamic models includes the classical Wagner and Kussner functions and the Leishman-Beddoes dynamic stall model. The simulations includes the numerical approximations of the Duhamel's integrals using both indicial (step) and impulse responses, the numerical integrations of the state-space models, and the exact solutions. The report also presents the conversion among different model representations.