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 793 records · Page 44

Nonlinear Dynamic Analysis of Deploying Flexible Space Booms

A fundamental investigation of the planar deployment and lock-up of two flexible boom type appendages which have attached tip masses and are connected to a central rigid body through a rotational spring is presented. Nondimensional parameters are identified and it is shown that, in general, the solution depends only on two mass ratios and one nondimensional stiffness parameter. Results are presented for boom tip deflections, deployment time and root moments at lock-up. A threshold value of the nondimensional stiffness parameter is identified beyond which boom deflections become large. Also, a thorough examination of the effect of nonlinear terms in the equations governing the deployment phase is performed. Nonlinear terms in the deployment equations due to kinematics and structural deformation are required to predict more accurately boom deflections, but retention of an inconsistent set of nonlinear terms leads to erroneous results. In particular, retaining nonlinear kinematic terms while neglecting nonlinear structural terms can produce inaccurate results even below the threshold stiffness value.

Mcgowan, P. E.↗

Aeroelastic analysis of prop fan blades with a semiempirical dynamic stall model

The time-history response of a propfan wind-tunnel model with dynamic stall was studied analytically. The response obtained from the analysis was compared with available experimental data. The governing equations of motion were formulated in terms of blade normal modes calculated using the COSMIC-NASTRAN computer code. The response analysis considered the blade plunging and pitching motions. The lift, drag, and moment coefficients for angles of attack below the static stall angle were obtained from a quasi-steady theory. For angles above static stall angles, a semiempirical dynamic stall model based on a correction to the angle of attack was used to obtain lift, drag, and moment coefficients. Using these coefficients, the aerodynamic forces were calculated at a selected number of strips, and integrated to obtain the total generalized forces. The combined momentum-blade element theory was used to calculate the induced velocity. The semiempirical stall model predicted a limit cycle oscillation near the setting angle at which large vibratory stresses were observed in an experiment. The predicted mode and frequency of oscillation also agreed with those measured in the experiment near this setting angle. The results also correlated well with the other published data that used a semiempirical dynamic stall model based on a synthesized procedure.

Reddy, T. S. R.↗

Vibration and wave propagation characteristics of multisegmented elastic beams

Closed form analytical solutions are derived for the vibration and wave propagation of multisegmented elastic beams. Each segment is modeled as a Timoshenko beam with possible inclusion of material viscosity, elastic foundation and axial forces. Solutions are obtained by using transfer matrix methods. According to these methods formal solutions are first constructed which relate the deflection, slope, moment and shear force of one end of the individual segment to those of the other. By satisfying appropriate continuity conditions at segment junctions, a global 4x4 matrix results which relates the deflection, slope, moment and shear force of one end of the beam to those of the other. If any boundary conditions are subsequently invoked on the ends of the beam one gets the appropriate characteristic equation for the natural frequencies. Furthermore, by invoking appropriate periodicity conditions the dispersion relation for a periodic system is obtained. A variety of numerical examples are included.

Nayfeh, Adnan H.↗

The Lateral Stability of Equal-flanged Aluminum-alloy I-beams Subjected to Pure Bending

Equal-flange beams of a special extruded I-section of 27ST aluminum alloy were tested in pure bending. Complete end fixity was not attained. Loading was continued until a definite maximum value had been reached. Tensile tests were made on specimens cut from the flanges and the web of each beam. Compressive stress-strain characteristics were determined by pack compression tests on specimens cut from the flanges. Values computed from an equation previously suggested by one of the authors for the critical stress at which such beams become unstable were found to be in good agreement with values computed from experimentally determined critically bending moments.

Dumont, C↗

The effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings

An analysis is presented of the effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings with partial-span constant-percent-chord ailerons extending inboard from the wing tip. The geometric variables considered are aspect ratio, taper ratio, aileron span, and aileron chord. The shape of the wing-torsional-stiffness curve is assumed and the twisting moment is considered to result solely from the pressure distribution caused by aileron deflection, so that the necessity of using a successive-approximation method is avoided. Because of the complexity of the equations resulting from the analysis, numerical calculations from the equations are presented in a series of figures. A computational form is provided to be used in conjunction with these figures so that calculations can be made without reference to the analysis.

Tucker, Warren A↗

Application of Model Based Parameter Estimation for Fast Frequency Response Calculations of Input Characteristics of Cavity-Backed Aperture Antennas Using Hybrid FEM/MoM Technique

Model Based Parameter Estimation (MBPE) is presented in conjunction with the hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique for fast computation of the input characteristics of cavity-backed aperture antennas over a frequency range. The hybrid FENI/MoM technique is used to form an integro-partial- differential equation to compute the electric field distribution of a cavity-backed aperture antenna. In MBPE, the electric field is expanded in a rational function of two polynomials. The coefficients of the rational function are obtained using the frequency derivatives of the integro-partial-differential equation formed by the hybrid FEM/ MoM technique. Using the rational function approximation, the electric field is obtained over a frequency range. Using the electric field at different frequencies, the input characteristics of the antenna are obtained over a wide frequency range. Numerical results for an open coaxial line, probe-fed coaxial cavity and cavity-backed microstrip patch antennas are presented. Good agreement between MBPE and the solutions over individual frequencies is observed.

Reddy C. J.↗

Ethnic Differences in Bending Stiffness of the Ulna and Tibia

There is considerable information about the variations in bone mass associated with different opportunity to compare a mechanical property of bone in young college women of Caucasian, Hispanic and Asian descent who gave informed consent to participate in an exercise study. The subjects were sedentary, in good health, eumenorrheic, non-smokers and had body mass indices (BMI) less than 30. Measurements acquired were body weight, kg, and height, cm, calcaneal and wrist bone density, g/square cm (PIXI, Lunar GE) and bending stiffness (EI, Nm(exp 2)) in the ulna and tibia. E1 was determined non-invasively with an instrument called the Mechanical Response Tissue Analyzer (MRTA) that delivers a vibratory stimulus to the center of the ulna or tibia and analyzes the response curve based on the equation E1 = k(sub b) L(exp 3)/48 where k, is lateral bending stiffness, L is the length of the bone, E is Young's modulus of elasticity and I, the bending moment of inertia. The error of the test (CV) based on measurements of an aluminum rod with a known E1 was 4.8%, of calcaneal BMD, 0.54%, and of wrist bone density, 3.45%.

Arnaud, S. B.↗

Paired Pulse Basis Functions for the Method of Moments EFIE Solution of Electromagnetic Problems Involving Arbitrarily-shaped, Three-dimensional Dielectric Scatterers

A pair of basis functions is presented for the surface integral, method of moment solution of scattering by arbitrarily-shaped, three-dimensional dielectric bodies. Equivalent surface currents are represented by orthogonal unit pulse vectors in conjunction with triangular patch modeling. The electric field integral equation is employed with closed geometries for dielectric bodies; the method may also be applied to conductors. Radar cross section results are shown for dielectric bodies having canonical spherical, cylindrical, and cubic shapes. Pulse basis function results are compared to results by other methods.

MacKenzie, Anne I.↗

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Three dimensional numerical scattering from ocean-like surfaces

A new method of calculating the electric field scattered from two-dimensional conducting rough surfaces is presented. This method uses the operator conjugate gradients technique and the extended boundary condition equation. This avoids the build-up of round-off errors encountered in the inversion of large matrices. It also avoids the singularities of the integral kernels usually encountered in the method of moments. Techniques for improving the rate of convergence by applying physical constraints are discussed. This method is presently being used to calculate the polarization signatures from simulated ocean-like surfaces.

Rodriguez, Ernesto↗

Sloshing dynamics on rotating helium dewar tank

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics were based on the non-inertia frame spacecraft bound coordinate, and solve time dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers were derived. Results were widely published in the open journals.

Hung, R. J.↗

Research activities at the Center for Modeling of Turbulence and Transition

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Shih, Tsan-Hsing↗

Numerical studies of the surface tension effect of cryogenic liquid helium

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose of performing scientific observation during the normal spacecraft operation is investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics has been based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers, have been derived.

Hung, R. J.↗

Estimation of Aircraft Nonlinear Unsteady Parameters From Wind Tunnel Data

Aerodynamic equations were formulated for an aircraft in one-degree-of-freedom large amplitude motion about each of its body axes. The model formulation based on indicial functions separated the resulting aerodynamic forces and moments into static terms, purely rotary terms and unsteady terms. Model identification from experimental data combined stepwise regression and maximum likelihood estimation in a two-stage optimization algorithm that can identify the unsteady term and rotary term if necessary. The identification scheme was applied to oscillatory data in two examples. The model identified from experimental data fit the data well, however, some parameters were estimated with limited accuracy. The resulting model was a good predictor for oscillatory and ramp input data.

Klein, Vladislav↗

Turbulence and the Stabilization Principle

Further results of research, reported in several previous NASA Tech Briefs articles, were obtained on a mathematical formalism for postinstability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence). To recapitulate: Fictitious control forces are introduced to couple the dynamical equations with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in ordinary perceived three-dimensional space is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. Con sequently, the postinstability motion is represented by a set of functions describing the evolution of such statistical quantities as expectations and higher moments, and this representation is stable. The previously reported findings are analyzed from the perspective of the authors Stabilization Principle, according to which (1) stability is recognized as an attribute of mathematical formalism rather than of underlying physics and (2) a dynamical system that appears unstable when modeled by differentiable functions only can be rendered stable by modifying the dynamical equations to incorporate intrinsic stochasticity.

Zak, Michail↗

A $C^1$-Conforming Arbitrary-Order Two-Dimensional Virtual Element Method for the Fourth-Order Phase-Field Equation

We present a two-dimensional conforming virtual element method for the fourth-order phase-field equation. Our proposed numerical approach to the solution of this high-order phase-field (HOPF) equation relies on the design of an arbitrary-order accurate, virtual element space with $C^1$ global regularity. Such regularity is guaranteed by taking the values of the virtual element functions and their full gradient at the mesh vertices as degrees of freedom. Attaining high-order accuracy requires also edge polynomial moments of the trace of the virtual element functions and their normal derivatives. In this work, we detail the scheme construction, and prove its convergence by deriving error estimates in different norms. A set of representative test cases allows us to assess the behavior of the method.

97 MATHEMATICS AND COMPUTING↗

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FORTRAN 4 computer program for calculating critical speeds of rotating shafts

A FORTRAN 4 computer program, written for the IBM DCS 7094/7044 computer, that calculates the critical speeds of rotating shafts is described. The shaft may include bearings, couplings, extra masses (nonshaft mass), and disks for the gyroscopic effect. Shear deflection is also taken into account, and provision is made in the program for sections of the shaft that are tapered. The boundary conditions at the ends of the shaft can be fixed (deflection and slope equal to zero) or free (shear and moment equal to zero). The fixed end condition enables the program to calculate the natural frequencies of cantilever beams. Instead of using the lumped-parameter method, the program uses continuous integration of the differential equations of beam flexure across different shaft sections. The advantages of this method over the usual lumped-parameter method are less data preparation and better approximation of the distribution of the mass of the shaft. A main feature of the program is the nature of the output. The Calcomp plotter is used to produce a drawing of the shaft with superimposed deflection curves at the critical speeds, together with all pertinent information related to the shaft.

Trivisonno, R. J.↗