Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 145 records · Page 8

A viscous flow analysis of the tip vortex generation process

A viscous primary-secondary flow analysis designed for numerical solution by forward marching integration is applied to the tip vortex generation problem. Equations governing streamwise momentum, streamwise vorticity, and a gas law are solved in conjunction with equations governing scalar and vector surface potentials for the secondary velocities. The numerical method used combines a consistently split linearized block implicit (LBI) scheme for parabolic equations and a scalar iterative ADI scheme for elliptic equations. Computed results are presented for a wing of rectangular planform immersed in a high Reynolds number stream at 6 degrees incidence. The physical mechanism for generation of the tip vortex, as revealed by the analysis and computations, is discussed. Detailed flow results for the case of 6 degree incidence are included.

Shamroth, S. J.↗

Radiation transport modeling using extended quadrature method of moments

The radiative transfer equation describes the propagation of radiation through a material medium. While it provides a highly accurate description of the radiation field, the large phase space on which the equation is defined makes it numerically challenging. As a consequence, significant effort has gone into the development of accurate approximation methods. Recently, an extended quadrature method of moments (EQMOM) has been developed to solve univariate population balance equations, which also have a large phase space and thus face similar computational challenges. The distinct advantage of the EQMOM approach over other moment methods is that it generates moment equations that are consistent with a positive phase space density and has a moment inversion algorithm that is fast and efficient. The goal of the current paper is to present the EQMOM method in the context of radiation transport, to discuss advantages and disadvantages, and to demonstrate its performance on a set of standard one-dimensional benchmark problems that encompass optically thin, thick, and transition regimes. Special attention is given in the implementation to the issue of realizability—that is, consistency with a positive phase space density. Numerical results in one dimension are promising and lay the foundation for extending the same framework to multiple dimensions.

97 MATHEMATICS AND COMPUTING↗

Approach guidance logic for a tilt-rotor aircraft

The distinctive feature of a tilt-rotor aircraft is that the pilot can change the rotor mast angles to go from a helicopter configuration for take-off and landing to an airplane configuration for high cruise speeds and vice-versa. An approach path for such an aircraft is proposed and the logic required to fly along this path in the presence of wind is determined. The main contribution of this work is an efficient and, to my knowledge, new method for generating the nominal state and control histories taking into account an estimate of the mean wind velocity and direction. The method requires the solution of algebraic (mostly linear) equations to generate a 'universal nominal', and feedforward and feedback gains. Then, in flight the additional state and control corrections due to deviation in descent rate, deceleration, and flight in a steady wind are obtained by multiplying simple precalculated functions of time by descent rate, deceleration or sine and cosine components of the mean wind vector. Simulations of approach flights for different wind conditions, assuming perfect state information in the feedback signal, indicated satisfactory performance.

Beser, J.↗

Multibody dynamics model building using graphical interfaces

In recent years, the extremely laborious task of manually deriving equations of motion for the simulation of multibody spacecraft dynamics has largely been eliminated. Instead, the dynamicist now works with commonly available general purpose dynamics simulation programs which generate the equations of motion either explicitly or implicitly via computer codes. The user interface to these programs has predominantly been via input data files, each with its own required format and peculiarities, causing errors and frustrations during program setup. Recent progress in a more natural method of data input for dynamics programs: the graphical interface, is described.

Macala, Glenn A.↗

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗

Theoretical and numerical modelling of chaotic electrostatic ion cyclotron (EIC) oscillations by Jerk equation

In the last few years, third order explicit autonomous differential equations, known as jerk equations, have generated great interest as they show features of regular and chaotic motion. In this paper, we have modelled chaotic electrostatic ion cyclotron oscillations using a third order nonlinear ordinary differential equation (ODE) and investigated its nonlinear dynamical properties. The nonlinear ODE has been derived for a plasma system from a two fluid model in the presence of a source term, under the influence of an external magnetic field, which is perpendicular to the direction of the wave vector. It is seen that the equation does not require an external forcing term to obtain chaotic behaviour. The stability of the solutions of the equation has been investigated analytically as well as numerically, and the bifurcation diagram obtained shows a number of interesting phenomena for various regimes of parameters. The coexisting attractors as well as their corresponding basins are shown and the phase space portraits at different conditions are obtained numerically and shown here. The results obtained here are in agreement with preliminary experiments conducted for a similar configuration of a plasma system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of turbulence by density waves

A nonlinear system describes the microdynamical state of turbulence that is excited by density waves. It consists of an equation of propagation and a master equation. A group-scaling generates the scaled equations of many interacting groups of distribution functions. The two leading groups govern the transport processes of evolution and eddy diffusivity. The remaining sub-groups represent the relaxation for the approach of diffusivity to equilibrium. In strong turbulence, the sub-groups disperse themselves and the ensemble acts like a medium that offers an effective damping to close the hierarchy. The kinetic equation of turbulence is derived. It calculates the eddy viscosity and identifies the effective damping of the assumed medium self-consistently. It formulates the coupling mechanism for the intensification of the turbulent energy at the expense of the wave energy, and the transfer mechanism for the cascade. The spectra of velocity and density fluctuations find the power law k sup-2 and k sup-4, respectively.

Tichen, C. M.↗

Three-dimensional elliptic grid generation for an F-16

A case history depicting the effort to generate a computational grid for the simulation of transonic flow about an F-16 aircraft at realistic flight conditions is presented. The flow solver for which this grid is designed is a zonal one, using the Reynolds averaged Navier-Stokes equations near the surface of the aircraft, and the Euler equations in regions removed from the aircraft. A body conforming global grid, suitable for the Euler equation, is first generated using 3-D Poisson equations having inhomogeneous terms modeled after the 2-D GRAPE code. Regions of the global grid are then designated for zonal refinement as appropriate to accurately model the flow physics. Grid spacing suitable for solution of the Navier-Stokes equations is generated in the refinement zones by simple subdivision of the given coarse grid intervals. That grid generation project is described, with particular emphasis on the global coarse grid.

Reese L. Sorenson↗

Evaluation of the Accuracy of the Load Prediction Equations of Low-Load Balance Calibration Data

Accuracy and reliability of the load prediction equations of a low-load calibration data set of a force balance were investigated. First, independent load prediction equations were generated from the data of a full-load and low-load machine calibration. Then, the low- load data set was processed as a set of precision check loads for the load prediction equations that were obtained from the full-load and low-load calibrations. Finally, the load prediction equations of the two calibrations were applied to four manual check load data sets that were recorded between 2009 and 2022. The residuals of the predicted calibration and check loads were compared in all cases. No systematic improvement of the load prediction accuracy was observed when the load prediction equations of the low-load calibration were applied to the check load data sets. Therefore, it is recommended to apply a full-load calibration to a six-component balance even if the balance is not used across its entire design envelope. This approach has the advantage that the resulting load prediction equations are less likely to be applied outside of the calibration load ranges. In addition, it is expected that the numerical estimates of the gage sensitivities are more reliable assuming that the uncertainties of the applied loads and measured outputs of the calibration data are more or less constant across the load range of the balance.

strain-gage balance↗

Rooftop Solar Photovoltaic Technical Potential in the United States

How much energy could we generate if PV modules were installed on all of the suitable roof area in the nation? To answer this question, we first use GIS methods to process a lidar dataset and determine the amount of roof area that is suitable for PV deployment in 128 cities nationwide, containing 23% of U.S. buildings, and provide PV-generation results for a subset of those cities. We then extend the insights from that analysis to the entire continental United States. We develop two statistical models -- one for small buildings and one for medium and large buildings -- and populate them with geographic variables that correlate with rooftop's suitability for PV. We simulate the productivity of PV installed on the suitable roof area, and present the technical potential of PV on both small buildings and medium/large buildings for every state in the continental US. Within the 128 cities covered by lidar data, 83% of small buildings have a location suitable for a PV installation, but only 26% of the total rooftop area of small buildings is suitable for development. The sheer number of buildings in this class, however, gives small buildings the greatest technical potential. Small building rooftops could accommodate 731 GW of PV capacity and generate 926 TWh/year of PV energy, approximately 65% of rooftop PV's total technical potential. We conclude by summing the PV-generation results for all building sizes and therefore answering our original question, estimating that the total national technical potential of rooftop PV is 1,118 GW of installed capacity and 1,432 TWh of annual energy generation. This equates to 39% of total national electric-sector sales.

14 SOLAR ENERGY↗

Evaluation of the Accuracy of the Load Prediction Equations of Low-Load Balance Calibration Data

Accuracy and reliability of the load prediction equations of a low-load calibration data set of a force balance were investigated. First, independent load prediction equations were generated from the data of a full-load and low-load machine calibration. Then, the low-load data set was processed as a set of precision check loads for the load prediction equations that were obtained from the full-load and low-load calibrations. Finally, the load prediction equations of the two calibrations were applied to four manual check load data sets that were recorded between 2009 and 2022. The residuals of the predicted calibration and check loads were compared in all cases. No systematic improvement of the load prediction accuracy was observed when the load prediction equations of the low-load calibration were applied to the check load data sets. Therefore, it is recommended to apply a full-load calibration to a six-component balance even if the balance is not used across its entire design envelope. This approach has the advantage that the resulting load prediction equations are less likely to be applied outside of the calibration load ranges. In addition, it is expected that the numerical estimates of the bridge sensitivities are more reliable assuming that the uncertainties of the applied loads and measured outputs of the calibration data are more or less constant across the load range of the balance.

strain-gage balance↗

A computer program to generate two-dimensional grids about airfoils and other shapes by the use of Poisson's equation

A method for generating two dimensional finite difference grids about airfoils and other shapes by the use of the Poisson differential equation is developed. The inhomogeneous terms are automatically chosen such that two important effects are imposed on the grid at both the inner and outer boundaries. The first effect is control of the spacing between mesh points along mesh lines intersecting the boundaries. The second effect is control of the angles with which mesh lines intersect the boundaries. A FORTRAN computer program has been written to use this method. A description of the program, a discussion of the control parameters, and a set of sample cases are included.

Sorenson, R. L.↗

Generalized Korteweg-de Vries equation induced from position-dependent effective mass quantum models and mass-deformed soliton solution through inverse scattering transform

We consider one-dimensional stationary position-dependent effective mass quantum model and derive a generalized Korteweg-de Vries (KdV) equation in (1+1) dimension through Lax pair formulation, one being the effective mass Schrödinger operator and the other being the time-evolution of wave functions. We obtain an infinite number of conserved quantities for the generated nonlinear equation and explicitly show that the new generalized KdV equation is an integrable system. Inverse scattering transform method is applied to obtain general solution of the nonlinear equation, and then N-soliton solution is derived for reflectionless potentials. Finally, a special choice has been made for the variable mass function to get mass-deformed soliton solution. The influence of position and time-dependence of mass and also of the different representations of kinetic energy operator on the nature of such solitons is investigated in detail. The remarkable features of such solitons are demonstrated in several interesting figures and are contrasted with the conventional KdV-soliton associated with constant-mass quantum model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Autonomous Noether Boundary-Value Problems not Solved with Respect to the Derivative

In monographs of N. V. Azbelev, A. M. Samoilenko, and A. A. Boichuk, constructive methods of study of Noether boundary-value problems have been developed. These methods continue the investigation of periodic problems stated by H. Poincaré, A. M. Lyapunov, N. M. Krylov, N. N. Bogolyubov, I. G. Malkin, and O. Veivoda by the methods of small parameter. We propose an improved scheme of study of autonomous Noether boundary-value problems for nonlinear systems in critical cases. In the case of multiple roots of the equation for generating constants, we obtain sufficient conditions of existence of solutions to an autonomous boundary-value problem not solved with respect to the derivative. The effectiveness of the scheme proposed is illustrated by an example of the periodic problem for the Liénard equation.

97 MATHEMATICS AND COMPUTING↗

Zonal-Mean Temperature Variations Inferred from SABER Measurements on TIMED Compared with UARS Observations

In the Numerical Spectral Model (NSM, Mayr et al., 2003), small-scale gravity waves propagating in the north/south direction can generate zonal mean (m = 0) meridional wind oscillations with periods between 2 and 4 months. These oscillations tend to be confined to low latitudes and have been interpreted to be the meridional counterpart of the wave-driven Quasi Biennial Oscillation in the zonal circulation. Wave driven meridional winds across the equator should generate, due to dynamical heating and cooling, temperature oscillations with opposite phase in the two hemispheres. We have analyzed SABER temperature measurements in the altitude range between 55 and 95 km to investigate the existence such variations. Because there are also strong tidal signatures (up to approximately 20 K) in the data, our algorithm estimates both mean values and tides together from the data. Based on SABER temperature data, the intra-annual variations with periods between 2 and 4 months can have amplitudes up to 5 K or more, depending on the altitude. Their amplitudes are in qualitative agreement with those inferred Erom UARS data (from different years). The SABER temperature variations also reveal pronounced hemispherical asymmetries, which are qualitatively consistent with wave driven meridional wind oscillations across the equator. Oscillations with similar periods have been seen in the meridional winds based on UARS data (Huang and Reber, 2003).

Huang, Frank T.↗

Rooftop Solar Photovoltaic Technical Potential in the United States. A Detailed Assessment

How much energy could be generated if PV modules were installed on all of the suitable roof area in the nation? To answer this question, we first use GIS methods to process a lidar dataset and determine the amount of roof area that is suitable for PV deployment in 128 cities nationwide, containing 23% of U.S. buildings, and provide PV-generation results for a subset of those cities. We then extend the insights from that analysis to the entire continental United States. We develop two statistical models--one for small buildings and one for medium and large buildings--and populate them with geographic variables that correlate with rooftop's suitability for PV. We simulate the productivity of PV installed on the suitable roof area, and present the technical potential of PV on both small buildings and medium/large buildings for every state in the continental US. Within the 128 cities covered by lidar data, 83% of small buildings have a location suitable for a PV installation, but only 26% of the total rooftop area of small buildings is suitable for development. The sheer number of buildings in this class, however, gives small buildings the greatest technical potential. Small building rooftops could accommodate 731 GW of PV capacity and generate 926 TWh/year of PV energy, approximately 65% of rooftop PV's total technical potential. We conclude by summing the PV-generation results for all building sizes and therefore answering our original question, estimating that the total national technical potential of rooftop PV is 1,118 GW of installed capacity and 1,432 TWh of annual energy generation. This equates to 39% of total national electric-sector sales.

14 SOLAR ENERGY↗

Approximate-factorization schemes for solving the transonic full-potential equation

The present paper provides a general discussion of approximate-factorization techniques applied to the transonic full-potential equation. Giving particular attention to the AF2 approximate-factorization scheme. This scheme was first introduced by Ballhaus and Steger (1975) for solving the low-frequency (unsteady), transonic small-disturbance equation. The full-potential equation algorithm is examined, taking into account the governing equations, grid generation, the artificial density scheme (spatial differencing), the alternating direction implicit scheme, the AF2 iteration scheme, temporal damping, and boundary conditions. Computed results are also presented. It is shown that fast, fully-implicit algorithms of the approximate-factorization variety are both efficient and reliable for solving the conservative full-potential equation.

Holst, T. L.↗