Search NASA⌕ Search

SEARCH · Search NASA

Results for “generation equation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Finite Element Simulation of Smart Structures

Finite element equations representing the behavior of piezoelectric materials when bounded to a typical structure and used as sensors and actuators were developed. Emphasis was placed on generating sensor output equations of piezoelectric sensors and responses of a typical structure bonded with piezoelectric sensors and actuators on the basis of finite element formulation. The model can predict not only structural responses due to both mechanical and electrical loading but also electrical potential due to mechanical or thermal effects. The resulted finite element equations were then used for simple control design and performance evaluation. In the control algorithm, voltages coming out from piezoelectric sensors, which are proportional to strains at sensing locations, are taken as input. The voltages applied to the piezoelectric actuators are used as output. The feasibility of integrating control algorithm with the element routine developed herein and FEAP was demonstrated. In particular, optimal independent modal space control was implemented in a software package on the basis of finite element formulation. A rudimentary finite element-control algorithm package was also developed to evaluate the performance of candidate control laws. A few numerical simulations using the software package developed herein were given. The integrated software package will provide a design tool to address issues such as how adaptive smart systems will scale to a full size aircraft, the amount of piezoelectric materials and the powers needed to actuate it for desired performance. It will also provide a viable new structural control design concept for practical applications in large flexible structures such as aerospace vehicles and aircraft.

Cui, Y. Lawrence↗
NASA NTRS2002

Numerical Simulations of Vortex Generator Vanes and Jets on a Flat Plate

Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Allan, Brian G.↗
NASA NTRS1989

The numerical solution of the Navier-Stokes equations for 3-dimensional, unsteady, incompressible flows by compact schemes

The present numerical method for the solution of unsteady, incompressible three-dimensional flow Navier-Stokes equations using velocity-vorticity variables and irregular Cartesian grids proceeds by solving: (1) equations of Cauchy-Riemann type for the velocity; and (2) transport-diffusion equations for the vorticity, whose solenoidal vorticity components are generated by a Poisson equation for an appropriate scalar potential. Iterations are used to solve the finite difference equations, facilitating the use of vector and parallel-computing methods; numerical experiments have verified the method's second-order spatial and temporal accuracy.

Gatski, T. B.↗
NASA NTRS1980

An investigation into grid patching techniques

In the past decade significant advances were made using flow field methods in the calculation of external transonic flows over aerodynamic configurations. It is now possible to calculate inviscid transonic flow over three dimensional configurations by solving the potential equation. However, with the exception of the transonic small disturbance methods which have the advantage of a simple cartesian grid, the configurations over which it is possible to calculate such flows are relatively simple. The major reason for this is the difficulty of producing compatibility between grid generation and flow equation solutions. The main programs in use, use essentially analytic transformations for prescribed configurations and, as such, are not easy to extend. While there is work in progress to extend this type of system to a limited extent, the long term effort is directed towards a more general approach. This approach should not be restricted to producing grid systems in isolation but rather a consideration of the overall problem of flow field solution.

Forsey, C. R.↗
NASA NTRS1997

Experimental Investigation of Jet Impingement Heat Transfer Using Thermochromic Liquid Crystals

Jet impingement cooling of a hypersonic airfoil leading edge is experimentally investigated using thermochromic liquid crystals (TLCS) to measure surface temperature. The experiment uses computer data acquisition with digital imaging of the TLCs to determine heat transfer coefficients during a transient experiment. The data reduction relies on analysis of a coupled transient conduction - convection heat transfer problem that characterizes the experiment. The recovery temperature of the jet is accounted for by running two experiments with different heating rates, thereby generating a second equation that is used to solve for the recovery temperature. The resulting solution requires a complicated numerical iteration that is handled by a computer. Because the computational data reduction method is complex, special attention is paid to error assessment. The error analysis considers random and systematic errors generated by the instrumentation along with errors generated by the approximate nature of the numerical methods. Results of the error analysis show that the experimentally determined heat transfer coefficients are accurate to within 15%. The error analysis also shows that the recovery temperature data may be in error by more than 50%. The results show that the recovery temperature data is only reliable when the recovery temperature of the jet is greater than 5 C, i.e. the jet velocity is in excess of 100 m/s. Parameters that were investigated include nozzle width, distance from the nozzle exit to the airfoil surface, and jet velocity. Heat transfer data is presented in graphical and tabular forms. An engineering analysis of hypersonic airfoil leading edge cooling is performed using the results from these experiments. Several suggestions for the improvement of the experimental technique are discussed.

Dempsey, Brian Paul↗
NASA NTRS1986

A preliminary theory of the origin of the 40-50 day oscillation

A theory of the 40 to 50 day oscillation in the tropical atmosphere which encompasses predominantly eastward propagation along the equator, spatial scales of wavenumber 1 and 2, slow time scales (20 to 70 days), and strong amplitude modulation in associated large-scale convective activity as a function of geographical location, i.e., pronounced convection is confined to the Indian Ocean and western Pacific, is presented. An initial disturbance generated near the equator disperses away from the source region, eastward by Kelvin waves and westward by Rossby waves. When the eastward disturbance passes over warm water it is expected to amplify by CISK due to enhanced supply of moisture from the ocean surface. Over cold water, its amplitude is expected to reduce due to decreased moisture supply. The climatological sea surface temperature contrast between the Indian Ocean/western Pacific and the eastern Pacific sets the spatial scale of the wave and its modulation to low wavenumbers and confines pronounced convection to the warm sector. If the eastward propagating wave can maintain itself over dissipation one circuit around the globe, it has a good chance to continue the process indefinitely, provided it is rejunvenated over the warm water with enough energy to go around for the next cycle.

Lau, K. M.↗
NASA NTRS1980

Application of the multigrid method to grid generation

The multigrid method (MGM), used to numerically solve the pair of nonlinear elliptic equations commonly used to generate two dimensional boundary-fitted coordinate systems is discussed. Two different geometries are considered: one involving a coordinate system fitted about a circle and the other selected for an impinging jet flow problem. Two different relaxation schemes are tried: one is successive point overrelaxation and the other is a four-color scheme vectorizeable to take advantage of a parallel processor computer for greater computational speed. Results using MGM are compared with those using SOR (doing successive overrelaxations with the corresponding relaxation scheme on the fine grid only). It is found that MGM becomes significantly more effective than SOR as more accuracy is demanded and as more corrective grids, or more grid points, are used. For the accuracy required, it is found that MGM is two to three times faster than SOR in computing time. With the four-color relaxation scheme as applied to the impinging jet problem, the advantage of MGM over SOR is not as great. This may be due to the effect of a poor initial guess on MGM for this problem.

Ohring, S.↗
NASA NTRS1993

High-Pressure Equation of State for Partially Ionic Solids

Recently, we showed that the cohesive energy of partially ionic solids may be characterized by a two-term energy relationship consisting of a Coulomb term arising from the valence-charge transfer delta Z between the atoms, and a scaled universal energy function E(sup *)(a(sup *)), which accounts for the partially covalent character of the bond and for the repulsion between the atomic cores for small R; a(sup *) is a scaled length. Normalized cohesive-energy curves of alkali halide crystals and of Ti and Ag halide crystals were obtained, and the cohesive-energy-curve parameters were used to generate theoretical equation-of-state (EOS) curves for the Li, Na, K, Cs, and Ag halides. Good agreement was obtained with the experimental isothermal compression curves over a wide pressure range (0-90 kbar). In this paper we verify that the cohesive-energy relationship is valid for divalent partially ionic solids; physically reasonable charge-transfer values (1.80 less than delta Z less than 2.0) are obtained for MgO, CaO, and CaS. Next, EOS curves for LiF, NaF, Nal, CsCl, Csl, MgO, CaO, and CaS are generated in terms of the cohesive-energy parameters. These EOS's yield excellent fits to experimental isothermal-compression data and to shock-wave data to very high pressures (P(sub max)= 250-1350 kbar).

Schlosser, Herbert↗
NASA NTRS2010

Enhanced Modeling of First-Order Plant Equations of Motion for Aeroelastic and Aeroservoelastic Applications

A methodology is described for generating first-order plant equations of motion for aeroelastic and aeroservoelastic applications. The description begins with the process of generating data files representing specialized mode-shapes, such as rigid-body and control surface modes, using both PATRAN and NASTRAN analysis. NASTRAN executes the 146 solution sequence using numerous Direct Matrix Abstraction Program (DMAP) calls to import the mode-shape files and to perform the aeroelastic response analysis. The aeroelastic response analysis calculates and extracts structural frequencies, generalized masses, frequency-dependent generalized aerodynamic force (GAF) coefficients, sensor deflections and load coefficients data as text-formatted data files. The data files are then re-sequenced and re-formatted using a custom written FORTRAN program. The text-formatted data files are stored and coefficients for s-plane equations are fitted to the frequency-dependent GAF coefficients using two Interactions of Structures, Aerodynamics and Controls (ISAC) programs. With tabular files from stored data created by ISAC, MATLAB generates the first-order aeroservoelastic plant equations of motion. These equations include control-surface actuator, turbulence, sensor and load modeling. Altitude varying root-locus plot and PSD plot results for a model of the F-18 aircraft are presented to demonstrate the capability.

Pototzky, Anthony S.↗
DOE OSTI2025

Prediction of the Wave Normal Angle of Proton‐Band EMIC Waves Near Geosynchronous Orbit

We investigate how the wave normal angle (WNA) and polarization of proton-band electromagnetic ion cyclotron (EMIC) waves change as they travel from their source to Earth. This paper marks a significant milestone as the first full-wave simulation of proton-band EMIC waves reflecting from the ionosphere. Our findings show that the WNA can change rapidly during propagation, primarily due to plasma inhomogeneities, such as variations in the Alfvén speed. The wave polarization is strongly related to the WNA, consistent with theory. Newly generated EMIC waves near the equator propagate with a WNA of 0° , then the WNA gradually shifts to 90° as they move toward Earth. In contrast, reflecting waves having 90° of WNA at Earth maintain a relatively larger WNA even near the magnetic equator. As a result, only the newly generated waves close to the source, where the magnetic latitude is less than approximately 20° , show left-handed polarization, while linear polarization remains dominant throughout the rest of the propagation.

EMIC wave↗
NASA NTRS1988

Incompressible and anelastic flow simulations on numerically generated grids

In the numerical simulation of incompressible and anelastic flows, it is necessary to solve an elliptic equation at each time step. When the boundaries of such flows are nonrectangular, it may be advantageous to solve the equations on a new, numerically generated coordinate grid, in which the property of orthogonality has been preserved. Flow equations in general curvilinear coordinates maintaining the conservative form are given for both anelastic models using the momentum equations, and for incompressible models, using the vorticity equation. The general problem of grid-generation in two dimensions is presented, and a quasi-conformal transformation technique is discussed in detail. Some examples of grids generated by this technique are exhibited. Three examples of the flow of a stratified fluid over obstacles are presented, in which the grid-generation permits some new results to be obtained.

Sharman, R. D.↗
DOE OSTI2026

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗
DOE OSTI2026

Multibody for Everybody (M4E) - A Linearization Approach to Enable Frequency Domain Analysis, Time Integration and Control Co-Design

1.1 Background/Objectives: Marine energy represents a promising yet underexploited source of power. To increase the harvested power, significant efforts have been made to improve wave energy converter (WEC) modeling capabilities and optimize power take-off (PTO) performance; however, these efforts have often treated WEC dynamics, PTO design, and controller development sequentially. In contrast, control co-design (CCD) is emerging as a promising strategy to address these issues directly, creating a growing need for fast analysis tools suitable for repeated simulation and parametric studies [1]. To support this need, this work presents the Multibody for Everybody (M4E) [2] linearization module, which employs a symbolic toolbox to provide deeper insight of WEC design parameters. The objective is to demonstrate that a minimal-coordinate linearization of articulated WEC dynamics can provide accurate wave response predictions and substantial computational savings relative to nonlinear time-domain simulation, while preserving compatibility with broader wave-energy analysis workflows, enabling CCD. 1.2 Approach/Activities: The proposed approach linearizes the equations of motion, generated by M4E, in minimal coordinates about a selected operating point and combines the resulting system with frequencydomain hydrodynamic terms to incorporate the reduced mass, damping, stiffness, and forcing operators. The linearized model is used for both impedance-based response amplitude operator (RAO) prediction and rapid regular-wave time integration. The methodology is demonstrated on a single-flap device and a FOSWEC configuration, with linearized M4E responses compared against the corresponding nonlinear M4E simulations and WEC-Sim results. Regular-wave time histories, RAO trends, and runtime differences are assessed. The framework is also compatible with broader wave-energy workflows, including coupling to WecOptTool, although that capability is not the focus of this work [3]. 1.3 Results/Lessons: The linearized M4E model reproduces key regularwave response characteristics such as integration and Response Amplitude over multiple frequencies. This module matches nonlinear M4E and WEC-Sim results while substantially reducing integration cost. Thus, the proposed framework can serve as a rapid analysis layer for articulated WEC design, parameter studies, and controls-oriented workflows. The analysis is most appropriate in the near-equilibrium regime, about the linearization point.

16 TIDAL AND WAVE POWER↗
NASA NTRS1972

Analysis and testing of stability augmentation systems

Testing and evaluation of stability augmentation systems for aircraft flight control were conducted. The flutter suppression system analysis of a scale supersonic transport wing model is described. Mechanization of the flutter suppression system is reported. The ride control synthesis for the B-52 aeroelastic model is discussed. Model analyses were conducted using equations of motion generated from generalized mass and stiffness data.

Sevart, F. D.↗
NASA NTRS1974

Chain pooling to minimize prediction error in subset regression

Monte Carlo studies using population models intended to represent response surface applications are reported. Simulated experiments were generated by adding pseudo random normally distributed errors to population values to generate observations. Model equations were fitted to the observations and the decision procedure was used to delete terms. Comparison of values predicted by the reduced models with the true population values enabled the identification of deletion strategies that are approximately optimal for minimizing prediction errors.

Holms, A. G.↗
NASA NTRS1978

Steering laws for double-gimbal control-moment gyros

Derivation of CMG steering laws was considered. The type of CMG considered was a double gimbal one which has a constant magnitude momentum vector. The equations for the generation of torque of an inertially fixed CMG are presented. To realize a particular commanded torque, the gimbal rates of the CMG must be specified. These gimbal rates are shown to be functions of gimbal angle orientation and the commanded torque. The conventional cross-product steering law was derived by minimizing a function which was the square of the magnitude of the error between the commanded torque and the realizable torque of the gyro. This steering law was severely dependent upon the inner gimbal angle. When the inner gimbal reaches 90 degrees, the cross product law requires infinite outer gimbal rates to realize finite torque values.

Winder, S. W.↗
NASA NTRS1980

Boundary-fitted coordinates for regions with highly curved boundaries and reentrant boundaries

A procedure was developed using the differential equation approach for generating boundary fitted coordinates for regions with highly curved boundaries as well as reentrant boundaries, such as those encountered in breaking surface waves. The resulting coordinates are nearly orthogonal and can provide adequate resolution even in the reentrant region. Consistent treatment of end boundaries and the use of a systematic initialization scheme and advanced implicit numerical solution techniques make the procedure highly efficient. The method developed for implicit enforcement of the periodicity boundary condition should be beneficial in the analysis of turbomachinery flow applications.

Ghia, U.↗