Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling coefficient”

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 523 records · Page 29

Integrated airframe propulsion control

Perturbation equations which describe flight dynamics and engine operation about a given operating point are combined to form an integrated aircraft/propulsion system model. Included in the model are the dependence of aerodynamic coefficients upon atmospheric variables along with the dependence of engine variables upon flight condition and inlet performance. An off-design engine performance model is used to identify interaction parameters in the model. Inclusion of subsystem interaction effects introduces coupling between flight and propulsion variables. To analyze interaction effects on control, consideration is first given to control requirements for separate flight and engine models. For the separate airframe model, feedback control provides substantial improvement in short period damping. For the integrated system, feedback control compensates for the coupling present in the model and provides good overall system stability. However, this feedback control law involves many non-zero gains. Analysis of suboptimal control strategies indicates that performance of the closed loop integrated system can be maintained with a feedback matrix in which the number of non-zero gains is small relative to the number of components in the feedback matrix.

Fennell, R. E.↗

New free-wake analysis of rotorcraft hover performance using influence coefficients

Free-wake analyses of helicopter rotor wakes in hover using time stepping have been shown to encounter instabilities which preclude convergence to valid free-vortex solutions for rotor-wake geometries. Previous work has demonstrated that these convergence difficulties can be overcome by implementing a new free-wake analysis method based on the use of influence coefficients. The present paper reviews this approach and documents its incorporation into a hover performance analysis called Evaluation of Hover Performance using Influence Coefficients (EHPIC). The technical principles underlying the EHPIC code are described with emphasis on steps taken to develop the single-filament wake models used in previous work into a multifilament wake valid for realistic hover performance predictions. The coupling of the wake model to a lifting surface loads analysis is described, and sample problems are solved that illustrate the robustness of the method. Performance calculations are also undertaken for hover to illustrate the utility of EHPIC in the analysis of rotorcraft performance.

Quackenbush, Todd R.↗

Recent Advances in Statistical Energy Analysis

Statistical Energy Analysis (SEA) has traditionally been developed using modal summation and averaging approach, and has led to the need for many restrictive SEA assumptions. The assumption of 'weak coupling' is particularly unacceptable when attempts are made to apply SEA to structural coupling. It is now believed that this assumption is more a function of the modal formulation rather than a necessary formulation of SEA. The present analysis ignores this restriction and describes a wave approach to the calculation of plate-plate coupling loss factors. Predictions based on this method are compared with results obtained from experiments using point excitation on one side of an irregular six-sided box structure. Conclusions show that the use and calculation of infinite transmission coefficients is the way forward for the development of a purely predictive SEA code.

K H Heron↗

Jig Twist Optimization of Mach 0.745 Transonic Truss-Braced Wing Aircraft and High-Fidelity CFD Validation

This paper presents a jig twist optimization study of Mach 0.745 Transonic Truss-Braced Wing (TTBW) aircraft using an in-house developed aero-structural analysis solver VSPAERO coupled to BEAM3D. A vortex-lattice model of the TTBW model is developed, and a transonic and viscous flow correction method is implemented in the VSPAERO model to account for transonic and viscous flow effects. A correction method for the wing-strut interference aerodynamics is developed and applied to the VSPAERO solver. Also, a structural dynamic finite-element model of the TTBW aircraft is developed. This finite-element model includes the geometric nonlinear effect due to the tension in the struts which causes a deflection-dependent nonlinear stiffness. The VSPAERO model is coupled to the corresponding finite-element model to provide a rapid aero-structural analysis. A flight condition corresponding to Mach 0.745 at 42000 ft is selected for the TTBW aircraft jig twist optimization to reduce the drag coefficient. After the design is implemented, the drag coefficient of the twist optimized TTBW aircraft is reduced about 8 counts. At the end, a high-fidelity CFD solver FUN3D is used to validate the design.

Xiong, Juntao↗

Investigation of Cloud Interior Structure Using a Tomographic Technique Based on Radon Transform. Application to the Measurements Made by the Research Scanning Polarimeter During CAMP2Ex

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP's uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid-water cloud droplet sizes using the rainbow structure observed in the polarized reflectance measurements. The RSP observations also provide for derivation of 2D field of extinction coefficient inside the cloud using a tomographic technique. This technique utilizes the family of cloud shapes derived using "cutout" technique and corresponding to a number of thresholds in total reflectance (separating bright cloud from its darker background). These shapes can be interpreted as level curves of an abstract "reflectance density", which is then used for derivation of the directional cloud optical thickness (dCOT) tomogram. The latter is a collection of dCOTs parameterized by the angles and offsets of the corresponding view rays (chords) relative to the cloud center. After this, the inverse Radon Transform (the mathematical basis of the X-ray com-puter tomography) is applied to the dCOT tomogram yielding 2D spatial distribution of the extinction coefficient. The later can be converted into droplet number concentration using the droplet size profiles derived from the RSP's polarized reflectance measurements. After tests on synthetic data this technique was applied to real RSP measurements from NASA's Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) conducted in the vicinity of the Philippines during the Southwest Monsoon (August–September 2019). We have investigated the interiors of a number of clouds observed during CAMP2Ex focusing on Cu and CuCg (Tcu) cases. Our retrievals were routinely compared with the correlative measurements by lidar (HSRL-2) and cloud radar (APR-3) that were deployed on the same airborne platform (NASA P-3B) during this field experiment.

Research Scanning Polarimeter↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

The pressure and deformation profile between two colliding lubricated cylinders

The pressure and deformation profiles between two colliding lubricated cylinders are obtained by solving the coupled, time-dependent elastohydrodynamic equations with an iterative procedure. The analysis includes effects which were not considered in a previous solution, namely, the effect of the lubricant compressibility and the effect of a lubricant with composite pressure-viscosity coefficients. It is found that the local approach velocity plays an important role during final stages of normal approach. It causes the lubricant to be entrapped within the contact region; neither the pressure nor the deformation profile converges to the Hertzian profile for a dry contact. The use of a smaller pressure-viscosity coefficient at high pressures reduces the sharp pressure gradient at the center of the contact and produces a much milder variation of load with respect to the film thickness. The effect of compressibility of the lubricant is found to be relatively small.

Lee, K.↗

Ultrasonic attenuation in superconducting molybdenum-rhenium alloys.

Investigation of longitudinal sound attenuation in superconducting Mo-Re alloys as a function of temperature, magnetic field, and frequency. Evaporated thin film CdS transducers were used for the measurements at frequencies up to 3 GHz. The normal state attenuation coefficient was found to be proportional to the square of frequency over this frequency range. Measurements in zero magnetic field yielded a value of the energy gap parameter close to the threshold value of 3.56 kTc, appropriate to a weakly coupled dirty limit superconductor.

Ashkin, M.↗

On the use of a finite difference method for solving anisotropic scattering problems

A new method of solving the radiative transfer equation is developed in which the scattering and absorption coefficients may have arbitrary variations with depth, and in which both internal (thermal) emission and incident radiation are allowed. Specular and diffuse reflection at both boundaries also is taken into account. The method begins by forming a paired set of coupled first-order differential equations for the symmetric and antisymmetric parts of the radiation field after writing the scattering integral as a numerical quadrature. These differential equations are broken into finite difference form, in which the symmetric and antisymmetric parts of the radiation field are found on alternate grid points. Numerical results for a number of test problems are shown, demonstrating that the method is very fast, that it returns specific intensities and fluxes that are accurate to at least a percent, and that it can be applied to optically thick problems.

Barkstrom, B. R.↗

A numerical simulation of seasonal stratospheric climate. I - Zonal temperatures and winds

Two numerical experiments were designed to isolate the effects of seasonal and latitudinal variations in longwave radiative transfer processes. The first experiment aimed at computing the rate of temperature change in the stratosphere due to longwave radiative transfer by CO2, H2O and O3 from Ramanathan's (1976) radiative-convective model applied to the troposphere-stratosphere circulation. The second experiment employs the Newtonian cooling approximation in which the rate of temperature change in the stratosphere by longwave radiative transfer is set equal to the product of the Newtonian cooling coefficient 'h' and the departure of the local temperature from a reference temperature. It is shown that the latitudinal temperature distribution of the lower stratosphere is maintained by the combined effects of dynamics, O3 solar heating and the longwave radiative coupling between troposphere and stratosphere. The latitudinal gradient in the troposphere-stratosphere longwave radiative coupling is maximum during winter and spring. The radiative response time (1/h) of the middle and upper stratosphere undergoes significant latitudinal and seasonal variations, largely due to the temperature dependence of h.

Ramanathan, V.↗

Airfoil design by numerical optimization using a minicomputer

A computer program developed for the automated design of low speed airfoils utilizes a generalized Joukowski method for aerodynamic analysis coupled with a conjugate gradient, penalty function, numerical optimization algorithm to give an efficient calculation technique for use with minicomputers. The program designs airfoils with a prescribed pressure distribution as well as those which minimize or maximize some aerodynamic force coefficient. At present the method is restricted to inviscid, incompressible flow. A typical design problem will execute in 4.5 hr on an HP 9830 minicomputer.

Hicks, R. M.↗

Numerical optimization of circulation control airfoils

A numerical procedure for optimizing circulation control airfoils, which consists of the coupling of an optimization scheme with a viscous potential flow analysis for blowing jet, is presented. The desired airfoil is defined by a combination of three baseline shapes (cambered ellipse, and cambered ellipse with drooped and spiralled trailing edges). The coefficients of these shapes are used as design variables in the optimization process. Under the constraints of lift augmentation and lift-to-drag ratios, the optimal airfoils are found to lie between those of cambered ellipse and the drooped trailing edge, towards the latter as the angle of attack increases. Results agree qualitatively with available experimental data.

Tai, T. C.↗

Self-consistent inclusion of space-charge in the traveling wave tube

It is shown how the complete field of the electron beam may be incorporated into the transmission line model theory of the traveling wave tube (TWT). The fact that the longitudinal component of the field due to the bunched beam is not used when formulating the beam-to-circuit coupling equation is not well-known. The fundamental partial differential equation for the traveling wave field is developed and compared with the older (now standard) one. The equation can be solved numerically using the same algorithms, but now the coefficients can be updated continuously as the calculation proceeds down the tube. The coefficients in the older equations are primarily derived from preliminary measurements and some trial and error. The newer coefficients can be found by a recursive method, since each has a well defined physical interpretation and can be calculated once a reasonable first trial solution is postulated. The results of the new expression were compared with those of the older forms, as well as to a field theory model to show the ease in which a reasonable fit to the field prediction is obtained. A complete summary of the existing transmission line modeling of the TWT is given to explain the somewhat vague ideas and techniques in the general area of drifting carrier-traveling circuit wave interactions. The basic assumptions and inconsistencies of the existing theory and areas of confusion in the general literature are examined and hopefully cleared up.

Freeman, Jon C.↗

Robust control of systems with real parameter uncertainty and unmodelled dynamics

During this research period we have made significant progress in the four proposed areas: (1) design of robust controllers via H infinity optimization; (2) design of robust controllers via mixed H2/H infinity optimization; (3) M-delta structure and robust stability analysis for structured uncertainties; and (4) a study on controllability and observability of perturbed plant. It is well known now that the two-Riccati-equation solution to the H infinity control problem can be used to characterize all possible stabilizing optimal or suboptimal H infinity controllers if the optimal H infinity norm or gamma, an upper bound of a suboptimal H infinity norm, is given. In this research, we discovered some useful properties of these H infinity Riccati solutions. Among them, the most prominent one is that the spectral radius of the product of these two Riccati solutions is a continuous, nonincreasing, convex function of gamma in the domain of interest. Based on these properties, quadratically convergent algorithms are developed to compute the optimal H infinity norm. We also set up a detailed procedure for applying the H infinity theory to robust control systems design. The desire to design controllers with H infinity robustness but H(exp 2) performance has recently resulted in mixed H(exp 2) and H infinity control problem formulation. The mixed H(exp 2)/H infinity problem have drawn the attention of many investigators. However, solution is only available for special cases of this problem. We formulated a relatively realistic control problem with H(exp 2) performance index and H infinity robustness constraint into a more general mixed H(exp 2)/H infinity problem. No optimal solution yet is available for this more general mixed H(exp 2)/H infinity problem. Although the optimal solution for this mixed H(exp 2)/H infinity control has not yet been found, we proposed a design approach which can be used through proper choice of the available design parameters to influence both robustness and performance. For a large class of linear time-invariant systems with real parametric perturbations, the coefficient vector of the characteristic polynomial is a multilinear function of the real parameter vector. Based on this multilinear mapping relationship together with the recent developments for polytopic polynomials and parameter domain partition technique, we proposed an iterative algorithm for coupling the real structured singular value.

Chang, Bor-Chin↗

Analysis of cascades using a two dimensional Euler aeroelastic solver

A two-dimensional unsteady aerodynamic Euler solver based on a flux differencing scheme is being developed to analyze oscillating cascades. The cascades can have subsonic, transonic, or supersonic flow with either subsonic or supersonic axial velocity. The aerodynamic solver is coupled with a typical section structural model for each blade of the cascade. Flutter analysis methods both in time and frequency domains are then implemented into the resulting aeroelastic solver. Methods that reduce computational time for calculating the unsteady aerodynamic coefficients, namely the influence coefficient method and the pulse response method, are also implemented and validated. The present solver showed good correlation with published results for all the flow regimes. It is shown that grid coarsening improved the accuracy of the predictions. A representative flutter calculation showed that both the frequency domain and time domain methods are implemented correctly into the aeroelastic solver.

Reddy, T. S. R.↗

Photochemistry of Saturn's Atmosphere: Hydrocarbon Chemistry and Comparisons with ISO Observations - 1

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Moses, Julianne I.↗

Galileo Heat Shield Response Experiment Data Reduction and Analyses

The Galileo heat shield response experiment consists of 10 Analog Resistance Ablation Detectors (ARADs) embedded in the carbon phenolic heat shield. As the vehicle descends into the atmosphere of Jupiter, these ARADs will in effect measure the recession as a function of time at various locations on the heat shield. This recession data will be used to reconstruct the time-dependent shape of the vehicle. The shape data is of critical importance to the Atmospheric Reconstruction Experiment which must know mass and drag coefficient of the vehicle in order to determine the structure of the Jovian atmosphere. The data is also intrinsically useful as a database for evaluating the accuracy of our numerical codes and methodology for predicting highly radiating reentry flowfields and coupled heat shield material response. This paper will document the reduction of the data from the ARAD sensors, reconstruction of the heat shield shape history, and comparison of the actual heat shield response with previous calculations and experiments. If possible, we will also include results from new calculations which couple the shock layer flow with the transient material response.

Milos, Frank S.↗