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 865 records · Page 48

Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation

A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments relies upon the global momentum conservation of the fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. Numerical examples illustrate the method's application to predicting bulk fluid motion including lateral propellant slosh in low-g conditions.

Orr, Jeb S.↗

Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation

A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments is based upon a novel approach that relies on the global momentum conservation of the closed fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. A numerical example illustrates the method's application to prediction of bulk fluid behavior during a spacecraft ullage settling maneuver.

Orr, Jeb S.↗

Analysis of Potential for Titanium Liner Buckling after Proof in a Large Kevlar/Epoxy COPV

We analyze the potential for liner buckling in a 40-in Kevlar49/epoxy overwrapped spherical pressure vessel (COPV) due to long, local depressions or valleys in the titanium liner, which appeared after proof testing (autofrettage). We begin by presenting the geometric characteristics of approximately 20 mil (0.02 in.) deep depressions measured by laser profilometry in several vessels. While such depths were more typical, depths of more than 40 mils (0.02 in.) were seen near the equator in one particular vessel. Such depressions are largely the result of overlap of the edges of overwrap bands (with rectangular cross-section prepreg tows) from the first or second wrap patterns particularly where they start and end. We then discuss the physical mechanisms of formation of the depressions during the autofrettage process in terms of uneven void compaction in the overwrap around the tow overlap lines and the resulting 10-fold increase in through-thickness stiffness of the overwrap. We consider the effects of liner plastic yielding mechanisms in the liner on residual bending moments and interface pressures with the overwrap both at the peak proof pressure (approx.6500 psi) and when reducing the pressure to 0 psi. During depressurization the Bauschinger phenomenon becomes very important whereby extensive yielding in tension reduces the magnitude of the yield threshold in compression by 30 to 40%, compared to the virgin annealed state of the liner titanium. In the absence of a depression, the liner is elastically stable in compression even at liner overwrap interface pressures nominally 6 times the approx. 1000 psi interface pressure that exists at 0 psi. Using a model based on a plate-on-an-elastic-foundation, we develop an extensive analysis of the possible destabilizing effects of a frozen-in valley. The analysis treats the modifying effects of the residual bending moments and interface pressures remaining after the proof hold as well as the Bauschinger effect on the compressive yield threshold. The key result is that depression depths of up to 40 mils can be tolerated, but above 40 mils, the Bauschinger effect drives destabilization, and buckling becomes increasingly likely depending on the details of depression formation during autofrettage. It is almost certain that destabilization and buckling will occur for depression depths beyond 55 mils. The main equations and formulas for treating the various phases of depression development and potential buckling, are only briefly outlined in the paper, but are available from the authors.

Phoenix, S. Leigh↗

Prediction of elastic-airplane lateral dynamics from rigid-body aerodynamics

Control-configured vehicle technology has increased the demand for detailed analysis of dynamic stability and control, handling and ride qualities, and control system dynamics at the early stages of preliminary design and development. For these early analyses an approximate, but reasonably accurate, set of equations of motion for elastic airplanes is needed. Such a formulation is developed for the lateral dynamics of elastic airplanes. It makes use of rigid-body aerodynamic stability derivatives and the antisymmetric elastic mode shapes and frequencies in formulating the forces and moments due to elastic motion. Verification of accuracy was made by comparison with B-1 airplane dynamics obtained by other methods. Frequencies and damping ratios of the coupled modes agree acceptably well with four antisymmetric elastic modes included.

Swaim, R. L.↗

On mixed and displacement finite element models of a refined shear deformation theory for laminated anisotropic plates

An improved plate theory that accounts for the transverse shear deformation is presented, and mixed and displacement finite element models of the theory are developed. The theory is based on an assumed displacement field in which the inplane displacements are expanded in terms of the thickness coordinate up to the cubic term and the transverse deflection is assumed to be independent of the thickness coordinate. The governing equations of motion for the theory are derived from the Hamilton's principle. The theory eliminates the need for shear correction factors because the transverse shear stresses are represented parabolically. A mixed finite element model that uses independent approximations of the displacements and moments, and a displacement model that uses only displacements as degrees of freedom are developed. A comparison of the numerical results for bending with the exact solutions of the new theory and the three-dimensional elasticity theory shows that the present theory (and hence the finite element models) is more accurate than other plate-theories of the same order.

Reddy, J. N.↗

Kinematics of the six-degree-of-freedom force-reflecting Kraft Master

Presented here are kinematic equations for a six degree of freedom force-reflecting hand controller. The forward kinematics solution is developed and shown in simplified form. The Jacobian matrix, which uses terms from the forward kinematics solution, is derived. Both of these kinematic solutions require joint angle inputs. A calibration method is presented to determine the hand controller joint angles given the respective potentiometer readings. The kinematic relationship describing the mechanical coupling between the hand and controller shoulder and elbow joints is given. These kinematic equations may be used in an algorithm to control the hand controller as a telerobotic system component. The purpose of the hand controller is two-fold: operator commands to the telerobotic system are entered using the hand controller, and contact forces and moments from the task are reflected to the operator via the hand controller.

Williams, Robert L., II↗

Fuel injector: Air swirl characterization aerothermal modeling, phase 2, volume 1

A well integrated experimental/analytical investigation was conducted to provide benchmark quality relevant to a prefilming type airblast fuel nozzle and its interaction with the combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM), and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

Fuel Injector: Air swirl characterization aerothermal modeling, phase 2, volume 2

A well integrated experimental/analytical investigation was conducted to provide benchmark quality data relevant to prefilming type airblast fuel nozzle and its interaction with combustor dome air swirler. The experimental investigation included a systematic study of both single-phase flows that involved single and twin co-axial jets with and without swirl. A two-component Phase Doppler Particle Analyzer (PDPA) equipment was used to document the interaction of single and co-axial air jets with glass beads that simulate nonevaporating spray and simultaneously avoid the complexities associated with fuel atomization processes and attendant issues about the specification of relevant boundary conditions. The interaction of jets with methanol spray produced by practical airblast nozzle was also documented in the spatial domain of practical interest. Model assessment activities included the use of three turbulence models (k-epsilon, algebraic second moment (ASM) and differential second moment (DSM)) for the carrier phase, deterministic or stochastic Lagrangian treatment of the dispersed phase, and advanced numerical schemes. Although qualitatively good comparison with data was obtained for most of the cases investigated, the model deficiencies in regard to modeled dissipation rate transport equation, single length scale, pressure-strain correlation, and other critical closure issues need to be resolved before one can achieve the degree of accuracy required to analytically design combustion systems.

Nikjooy, M.↗

Steady-State Computation of Constant Rotational Rate Dynamic Stability Derivatives

Dynamic stability derivatives are essential to predicting the open and closed loop performance, stability, and controllability of aircraft. Computational determination of constant-rate dynamic stability derivatives (derivatives of aircraft forces and moments with respect to constant rotational rates) is currently performed indirectly with finite differencing of multiple time-accurate computational fluid dynamics solutions. Typical time-accurate solutions require excessive amounts of computational time to complete. Formulating Navier-Stokes (N-S) equations in a rotating noninertial reference frame and applying an automatic differentiation tool to the modified code has the potential for directly computing these derivatives with a single, much faster steady-state calculation. The ability to rapidly determine static and dynamic stability derivatives by computational methods can benefit multidisciplinary design methodologies and reduce dependency on wind tunnel measurements. The CFL3D thin-layer N-S computational fluid dynamics code was modified for this study to allow calculations on complex three-dimensional configurations with constant rotation rate components in all three axes. These CFL3D modifications also have direct application to rotorcraft and turbomachinery analyses. The modified CFL3D steady-state calculation is a new capability that showed excellent agreement with results calculated by a similar formulation. The application of automatic differentiation to CFL3D allows the static stability and body-axis rate derivatives to be calculated quickly and exactly.

Park, Michael A.↗

Some applications of aerodynamic formulations to problems in aircraft dynamics

Two applications of mathematical modeling to aerodynamic problems are cussed. The first application is an investigation of the capacity of a nonlinear aerodynamic mathematical model to describe the aerodynamic reactions on an airfoil with a deflecting flap in transonic flow. Flow field computational methods are used to evaluate the nonlinear, unsteady aerodynamic data in terms of characteristic motions called for by the model. Histories of unconstrained motions of the flap are generated from the flap equations of motion, with the aerodynamic reactions specified by the mathematical model. In the second application wing rock is investigated. The most recent model accommodates experimental results wing rock by admitting the existence of aerodynamic hysteresis in the variation of the steady state rolling moment coefficient with roll angle is described.

Schiff, L. B.↗

Theoretical aspects of sounder antenna operation

The fundamental problem of calculating the radiation from an electric dipole immersed in a magnetoplasma is one of solving Maxwell's equations coupled to the self-consistent field-particle equations describing the plasma. For very small applied antenna voltages, the plasma dynamics may be linearized, and the plasma reaction currents computed self-consistently in terms of the voltage on the antenna. Thus, the short dipole can be idealized as a point dipole of moment P sub o = C sub a, V sub o el, where C sub a is the antenna element-to-plasma capacitance for a sheath of dimensions roughly several plasma Debye lengths and el is the dipole element's length. In actual practice, it is at this point that the first difficulty with the antenna theory arises. Even for small signals, the sheath around such an antenna will exhibit a complex impedance. For larger signals, nonlinear effects appear. Attemps are made to solve this problem.

Fredricks, R. W.↗

Circulation control STOL aircraft design aspects

Since Davidson patented Circulation Control Airfoils in 1960, there have been only 2 aircraft designed and flown with circulation control (CC). Designing with CC is complex for the following reasons: the relation between lift increase and blowing momentum is nonlinear; for good cruise performance one must change the wing geometry in flight from a round to a sharp trailing edge. The bleed air from the propulsion engines or an auxiliary compressor, must be used efficiently. In designing with CC, the propulsion and control aspects are just as important as aerodynamics. These design aspects were examined and linearized equations are presented in order to facilitate a preliminary analysis of the performance potential of CC. The thrust and lift requirements for takeoff make the calculated runway length very sensitive to the bleed air ratio. Thrust vectoring improves performance and can offset nose down pitching moments. The choice of blowing jet to free stream velocity ratio determines the efficiency of applying bleed air power.

Loth, John L.↗

Equipotential transformation of multipole systems to dipole systems

The vector sum of fields produced by a set of magnetic dipoles with different magnetic moments that are both tilted from the planetary spin axis and offset from the planetary center by different amounts can be used to completely represent a planetary magnetic field in the conventional form of spherical harmonic expansion. The scalar equipotential transformation analytically results in 24 equations that may subsequently be solved for the 24 adjustable parameters in dipole systems with the predetermined main dipole. Attention is given to this method's application to the Jovian magnetic field.

Wang, W. X.↗

Periodic-disturbance accommodating control of the Space Station for asymptotic momentum management

Periodic-disturbance accommodating control is investigated for asymptotic momentum management of control moment gyros used as primary actuating devices for the Space Station. The proposed controller utilizes the concepts of quaternion feedback control and periodic-disturbance accommodation to achieve oscillations about the constant torque equilibrium attitude, while minimizing the control effort required. Three-axis coupled equations of motion, written in terms of quaternions, are derived for roll/yaw controller design and stability analysis. The quaternion feedback controller designed using the linear-quadratic regulator synthesis technique is shown to be robust for a wide range of pitch angles. It is also shown that the proposed controller tunes the open-loop unstable vehicle to a stable oscillatory motion which minimizes the control effort needed for steady-state operations.

Warren, Wayne↗

Periodic-disturbance accommodating control of the space station for asymptotic momentum management

Periodic-disturbance accommodating control is investigated for asymptotic momentum management of control moment gyros used as primary actuating devices for the Space Station. The proposed controller utilizes the concepts of quaternion feedback control and periodic-disturbance accommodation to achieve oscillations about the constant torque equilibrium attitude, while minimizing the control effort required. Three-axis coupled equations of motion, written in terms of quaternions, are derived for roll/yaw controller design and stability analysis. The quaternion feedback controller designed using the linear-quadratic regulator synthesis technique is shown to be robust for a wide range of pitch angles. It is also shown that the proposed controller tunes the open-loop unstable vehicle to a stable oscillatory motion which minimizes the control effort needed for steady-state operations.

Warren, Wayne↗

Computational attributes of the integral form of the equation of transfer

Difficulties can arise in radiative and neutron transport calculations when a highly anisotropic scattering phase function is present. In the presence of anisotropy, currently used numerical solutions are based on the integro-differential form of the linearized Boltzmann transport equation. This paper, departs from classical thought and presents an alternative numerical approach based on application of the integral form of the transport equation. Use of the integral formalism facilitates the following steps: a reduction in dimensionality of the system prior to discretization, the use of symbolic manipulation to augment the computational procedure, and the direct determination of key physical quantities which are derivable through the various Legendre moments of the intensity. The approach is developed in the context of radiative heat transfer in a plane-parallel geometry, and results are presented and compared with existing benchmark solutions. Encouraging results are presented to illustrate the potential of the integral formalism for computation. The integral formalism appears to possess several computational attributes which are well-suited to radiative and neutron transport calculations.

Frankel, J. I.↗

Broadband vertical interconnects using slot-coupled shielded microstrip lines

A full-wave space-domain integral equation analysis of aperture coupled shielded microstrip lines is presented based on the Equivalence Principle. The formulation includes the capability to model multilayered substrates through the derivation of the associated dyadic Green's functions which represent the layers through impedance boundary conditions. The method of moments is used to solve for the line currents and slot voltage with even and odd mode excitations which are then interpreted through transmission line analysis to determine the two-port scattering parameters. A parametric study together with experimental data is presented which demonstrates the behavior of the coupler and the accuracy of the technique.

Vandenberg, Norman L.↗