Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moment equations”

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 343 records · Page 19

A similarity hypothesis for the two-point correlation tensor in a temporally evolving plane wake

The analysis demonstrated that the governing equations for the two-point velocity correlation tensor in the temporally evolving wake admit similarity solutions, which include the similarity solutions for the single-point moment as a special case. The resulting equations for the similarity solutions include two constants, beta and Re(sub sigma), that are ratios of three characteristic time scales of processes in the flow: a viscous time scale, a time scale characteristic of the spread rate of the flow, and a characteristic time scale of the mean strain rate. The values of these ratios depend on the initial conditions of the flow and are most likely measures of the coherent structures in the initial conditions. The occurrences of these constants in the governing equations for the similarity solutions indicates that these solutions, in general, will only be the same for two flows if these two constants are equal (and hence the coherent structures in the flows are related). The comparisons between the predictions of the similarity hypothesis and the data presented here and elsewhere indicate that the similarity solutions for the two-point correlation tensors provide a good approximation of the measures of those motions that are not significantly affected by the boundary conditions caused by the finite extent of real flows. Thus, the two-point similarity hypothesis provides a useful tool for both numerical and physical experimentalist that can be used to examine how the finite extent of real flows affect the evolution of the different scales of motion in the flow.

Ewing, D. W.↗

On the modelling of non-reactive and reactive turbulent combustor flows

A study of non-reactive and reactive axisymmetric combustor flows with and without swirl is presented. Closure of the Reynolds equations is achieved by three models: kappa-epsilon, algebraic stress and Reynolds stress closure. Performance of two locally nonequilibrium and one equilibrium algebraic stress models is analyzed assuming four pressure strain models. A comparison is also made of the performance of a high and a low Reynolds number model for combustor flow calculations using Reynolds stress closures. Effects of diffusion and pressure-strain models on these closures are also investigated. Two models for the scalar transport are presented. One employs the second-moment closure which solves the transport equations for the scalar fluxes, while the other solves the algebraic equations for the scalar fluxes. In addition, two cases of non-premixed and one case of premixed combustion are considered. Fast- and finite-rate chemistry models are applied to non-premixed combustion. Both show promise for application in gas turbine combustors. However, finite rate chemistry models need to be examined to establish a suitable coupling of the heat release effects on turbulence field and rate constants.

Nikjooy, Mohammad↗

Electromagnetic scattering from arbitrary shaped conducting bodies coated with lossy materials of arbitrary thickness

A simple and efficient numerical technique is presented to solve the electromagnetic scattering problem of coated conducting bodies of arbitrary shape. The surface equivalence principle is used to formulate the problem in terms of a set of coupled integral equations involving equivalent electric and magnetic surface currents which represent boundary fields. The conducting structures and the dielectric materials are modeled by planar triangular patches, and the method of moments is used to solve the integral equations. Numerical results for scattering cross sections are given for various structures and compared with other available data. These results are proved accurate by a number of representative examples.

Rao, Sadasiva M.↗

Gravity tectonics and seismic gaps in the mantle

The concept of gravity tectonics is applied to reveal the major clue as to the conditions which result in the correspondence of seismic and tectonic gaps in the mantle. An asymptotic theory is developed for the calculation of the thrust and moment when a descending lithospheric plate encounters resistance to its downward motion in the mesosphere. Dynamic analysis falls into two parts: (1) deriving equations for forces in the descending lithosphere, (2) deducing moment distribution which causes the detachment of lithosphere. For the analysis of forces a mathematical theory of shells is given. In order to determine the detachment mechanism, solutions of equations are obtained by asymptotic integration. It is found that a thrust N sub phi coupled with a moment M sub phi due to gravitational forces generated by density contrast may play a key role in the initial detachment of a piece of descending lithosphere. The results are in agreement with the observed seismic gaps beneath South America, Toga-Fiji, New Zealand and New Hebrides regions.

Liu, H. S.↗

An approximate method for solution to variable moment of inertia problems

An approximation method is presented for reducing a nonlinear differential equation (for the 'weather vaning' motion of a wind turbine) to an equivalent constant moment of inertia problem. The integrated average of the moment of inertia is determined. Cycle time was found to be the equivalent cycle time if the rotating speed is 4 times greater than the system's minimum natural frequency.

Beans, E. W.↗

Asymptotic Waveform Evaluation (AWE) Technique for Frequency Domain Electromagnetic Analysis

The Asymptotic Waveform Evaluation (AWE) technique is applied to a generalized frequency domain electromagnetic problem. Most of the frequency domain techniques in computational electromagnetics result in a matrix equation, which is solved at a single frequency. In the AWE technique, the Taylor series expansion around that frequency is applied to the matrix equation. The coefficients of the Taylor's series are obtained in terms of the frequency derivatives of the matrices evaluated at the expansion frequency. The coefficients hence obtained will be used to predict the frequency response of the system over a frequency range. The detailed derivation of the coefficients (called 'moments') is given along with an illustration for electric field integral equation (or Method of Moments) technique. The radar cross section (RCS) frequency response of a square plate is presented using the AWE technique and is compared with the exact solution at various frequencies.

Cockrell, C. R.↗

Uniaxial aerodynamic attitude control of artificial satellites

Within the context of a simple mechanical model the paper examines the movement of a satellite with respect to the center of masses under conditions of uniaxial aerodynamic attitude control. The equations of motion of the satellite take account of the gravitational and restorative aerodynamic moments. It is presumed that the aerodynamic moment is much larger than the gravitational, and the motion equations contain a large parameter. A two-parameter integrated surface of these equations is constructed in the form of formal series in terms of negative powers of the large parameter, describing the oscillations and rotations of the satellite about its lengthwise axis, approximately oriented along the orbital tangent. It is proposed to treat such movements as nominal undisturbed motions of the satellite under conditions of aerodynamic attitude control. A numerical investigation is made for the above integrated surface.

Sazonov, V. V.↗

Hinge Moments of Sealed-Internal-Balance Arrangements for Control Surfaces I : Theoretical Investigation

The results of a theoretical analysis of the hinge-moment characteristics of various sealed-internal-balance arrangements for control surfaces are presented. The analysis considered overhands sealed to various types of wing structure by flexible seals spanning gaps of various widths or sealed to the wing structure by a flexible system of linked plates. Leakage was not considered; the seal was assumed to extend the full spanwise length of the control surface. The effect of the developed width of the flexible seal and of the geometry of the structure to which the seal was anchored was investigated, as well as the effect of the gap width that is sealed. The results of the investigation indicated that the most nearly linear control-surface hinge-moment characteristics can probably be obtained from a flexible seal over a narrow gap (about 0.1 of the overhang chord), which is so installed that the motion of the seal is restricted to a region behind the point of attachment of the seal to the wing structure. Control-surface hinge moments that tend to be high at large deflections and low or overbalanced at small deflections will result if a very narrow seal is used.

CONTROL SURFACES, BALANCED - OVERHANG↗

A Polar Scaling Technique for the Regularization of Strongly Singular and Strongly Near-Singular Helmholtz Surface Integrals Evaluated Over 2-D Domains

The numerical integration of expressions containing strong singularities or strong near-singularities has long been a challenging problem in the electromagnetics community. Much attention has been paid to this problem, as strong $1/R^{{2}}$ singularities routinely appear when implementing electromagnetic simulation techniques like the method of moments (MoM). To date, several techniques, from singularity extraction to singularity cancellation (SC), have been employed to deal with problems that require the evaluation of 2-D strongly singular integrals. However, no single technique has been proposed that can deal with both strong singularities and strong near-singularities in a fully numerical manner for arbitrary 2-D domains. Moreover, it has been claimed that the Helmholtz-type strongly singular integral found in the MoM is convergent in a principal value sense, but this convergence value has yet to be proven mathematically. In this work, we will conduct the convergence proof and introduce a “polar scaling” change of variables method that may be used to evaluate Helmholtz integrals with both strong and weak singularities/near-singularities. The technique is fully numerical and can in principle be applied to any planar or curved polygon and any nonsingular basis function. We will also provide numerical results showing useful convergence behavior for integrals involving both exact and near-singularities.

47 OTHER INSTRUMENTATION↗

Prediction of heat release effects on a mixing layer

A fully second-order closure model for turbulent reacting flows is suggested based on Favre statistics. For diffusion flames the local thermodynamic state is related to single conserved scalar. The properties of pressure fluctuations are analyzed for turbulent flows with fluctuating density. Closure models for pressure correlations are discussed and modeled transport equations for Reynolds stresses, turbulent kinetic energy dissipation, density-velocity correlations, scalar moments and dissipation are presented and solved, together with the mean equations for momentum and mixture fraction. Solutions of these equations are compared with the experimental data for high heat release free mixing layers of fluorine and hydrogen in a nitrogen diluent.

Farshchi, M.↗

Prediction of heat release effects on a mixing layer

A fully second-order closure model for turbulent reacting flows is suggested based on Favre statistics. For diffusion flames the local thermodynamic state is related to a single conserved scalar. The properties of pressure fluctuations are analyzed for turbulent flows with fluctuating density. Closure models for pressure correlations are discussed and modelled transport equations for Reynolds stresses, turbulent kinetic energy dissipation, density-velocity correlations, scalar moments and dissipation are presented and solved, together with the mean equations for momentum and mixture fraction. Solutions of these equations are compared with the experimental data for high heat release free mixing layers of fluorine and hydrogen in a nitrogen diluent.

Farshchi, M.↗

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗

Transition characteristics and electron kinetics in microhollow cathode discharges

We demonstrate the transition characteristics and electron kinetics of microdischarges in a microgap with a cathode having microstructures using two-dimensional particle-in-cell/Monte Carlo collision (2d3v) simulations. It is found that the microdischarge is inside the hollow cavity at lower pressures, forming microhollow cathode discharges (MHCDs), while the dominant discharge moves outside the cavity at higher pressures, locating above the cathode rectangular protrusion tip. The spatial distributions of the microdischarge parameters (e.g., electric potential and electron density) are presented by capturing the transition characteristics. The electron kinetics of the microdischarges are examined based on the moment analysis of the electron Boltzmann equation. The collisional and non-collisional components of the electron power absorption are decomposed, which were seldomly investigated for microdischarges previously. The results provide fundamental insights into MHCD formation with a structured electrode at varying pressure conditions, which could be beneficial for the design and fabrication of microplasma devices in practical applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

The aim of the study is to analyze the stability of the energetic particle modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. Additionally, the analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n = 1 EPM and n = 2 global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis ($\rlap{-} \iota_{0}$). In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given $\rlap{-} \iota_{0}$, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n = 1 EPM, n = 2 toroidal AE (TAE) and n = 3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP β threshold with respect to the co- and no-ECCD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sinusoidal reaction formulation for radiation and scattering from conducting surfaces

A piecewise-sinusoidal reaction technique is proposed for the study of scattering and radiation from perfectly conducting bodies of arbitrary shape. The reaction technique described is a moment method based on the integral equation of an electric field. A theory of sinusoidal reaction formulation for radiation and scattering problems is outlined, in which the scattered field is obtained by integrating the surface currents. Numerical results are presented for the scattering patterns of rectangular plates and corner reflectors and the radiation patterns of corner-reflector antennas.

Wang, N. N.↗

Relation of a contracting earth to the apparent accelerations of the sun and moon

It has been concluded that the theoretical ratio of the apparent secular accelerations of the moon and the sun implied by tidal theory is far higher than the ratio of observed values. This paper shows that the discrepancy can be entirely removed within the range of uncertainty of the various quantities involved if terms making allowance for a secularly decreasing terrestrial moment of inertia are included in the dynamical equations for the apparent secular accelerations. Such a decrease has already been predicted on the basis of the phase-change hypothesis for the nature of earth's core. The development of this hypothesis is reviewed, the modified dynamical equations are derived, and values of tidal couples implied by existing observational data are examined. The numerically determined rate of decrease of earth's moment of inertia is found to be in close agreement with that predicted by the phase-change hypothesis. The results indicate that a diminishing gravitational constant is unlikely, that the moon would have been in close proximity to earth only 1 billion years ago, and that earth's rotation period at that time would have been about 5 hr.

Lyttleton, R. A.↗

Predicting aerodynamic sound utilizing a two-point, two-time turbulence theory

The feasibility of determining two-point, two-time turbulent velocity correlations based on extending closure models for one-point, one-time turbulence correlations is demonstrated. The procedure used is based on a spatial moment integral formulation of the governing equations using approximate, parameterized trial functions for the two-point, two-time velocity correlations. Solution of the equations is shown to give a set of anisotropic length scales and the separation-time-dependent decorrelation of the ensemble averaged turbulent velocities. The analysis is simplified using the assumption of homogeneous stationary turbulence and a constant shear, unidirectional mean flow. It is shown that the anisotropic behavior of measured turbulence correlations can be characterized through this technique.

Hecht, A. M.↗