SEARCH · Engineering Papers
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.
The present paper treats less of the effect of the union (of spars and ribs, namely, the reduction of the bending moment at the fixed ends of the spars) than of its influence on the torsional rigidity of the wing. The calculations are carried out for a two-spar wing of constant cross section, in which the ribs are replaced by a continuous member of constant rigidity.
The feasibility of using forebody tangential blowing to control the roll-yaw motion of a wind tunnel model is experimentally demonstrated. An unsteady model of the aerodynamics is developed based on the fundamental physics of the flow. Data from dynamic experiments is used to validate the aerodynamic model. A unique apparatus is designed and built that allows the wind tunnel model two degrees of freedom, roll and yaw. Dynamic experiments conducted at 45 degrees angle of attack reveal the system to be unstable. The natural motion is divergent. The aerodynamic model is incorporated into the equations of motion of the system and used for the design of closed loop control laws that make the system stable. These laws are proven through dynamic experiments in the wind tunnel using blowing as the only actuator. It is shown that asymmetric blowing is a highly non-linear effector that can be linearized by superimposing symmetric blowing. The effects of forebody tangential blowing and roll and yaw angles on the flow structure are determined through flow visualization experiments. The transient response of roll and yaw moments to a step input blowing are determined. Differences on the roll and yaw moment dependence on blowing are explained based on the physics of the phenomena.
The feasibility of using forebody tangential blowing to control the roll-yaw motion of a wind tunnel model is experimentally demonstrated. An unsteady model of the aerodynamics is developed based on the fundamental physics of the flow. Data from dynamic experiments is used to validate the aerodynamic model. A unique apparatus is designed and built that allows the wind tunnel model two degrees of freedom, roll and yaw. Dynamic experiments conducted at 45 degrees angle of attack reveal the system to be unstable. The natural motion is divergent. The aerodynamic model is incorporated into the equations of motion of the system and used for the design of closed loop control laws that make the system stable. These laws are proven through dynamic experiments in the wind tunnel using blowing as the only actuator. It is shown that asymmetric blowing is a highly non-linear effector that can be linearized by superimposing symmetric blowing. The effects of forebody tangential blowing and roll and yaw angles on the flow structure are determined through flow visualization experiments. The transient response of roll and yaw moments to a step input blowing are determined. Differences on the roll and yaw moment dependence on blowing are explained based on the physics of the phenomena.
The research investigates the effect of wind tunnel model system dynamics on measured aerodynamic data. During wind tunnel tests designed to obtain lift and drag data, the required aerodynamic measurements are the steady-state balance forces and moments, pressures, and model attitude. However, the wind tunnel model system can be subjected to unsteady aerodynamic and inertial loads which result in oscillatory translations and angular rotations. The steady-state force balance and inertial model attitude measurements are obtained by filtering and averaging data taken during conditions of high model vibrations. The main goals of this research are to characterize the effects of model system dynamics on the measured steady-state aerodynamic data and develop a correction technique to compensate for dynamically induced errors. Equations of motion are formulated for the dynamic response of the model system subjected to arbitrary aerodynamic and inertial inputs. The resulting modal model is examined to study the effects of the model system dynamic response on the aerodynamic data. In particular, the equations of motion are used to describe the effect of dynamics on the inertial model attitude, or angle of attack, measurement system that is used routinely at the NASA Langley Research Center and other wind tunnel facilities throughout the world. This activity was prompted by the inertial model attitude sensor response observed during high levels of model vibration while testing in the National Transonic Facility at the NASA Langley Research Center. The inertial attitude sensor cannot distinguish between the gravitational acceleration and centrifugal accelerations associated with wind tunnel model system vibration, which results in a model attitude measurement bias error. Bias errors over an order of magnitude greater than the required device accuracy were found in the inertial model attitude measurements during dynamic testing of two model systems. Based on a theoretical modal approach, a method using measured vibration amplitudes and measured or calculated modal characteristics of the model system is developed to correct for dynamic bias errors in the model attitude measurements. The correction method is verified through dynamic response tests on two model systems and actual wind tunnel test data.
Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.
A paper presents further results of continuing research reported in several previous NASA Tech Briefs articles, the two most recent being Stochastic Representations of Chaos Using Terminal Attractors (NPO-41519), [Vol. 30, No. 5 (May 2006), page 57] and Physical Principle for Generation of Randomness (NPO-43822) [Vol. 33, No. 5 (May 2009), page 56]. This research focuses upon a mathematical formalism for describing post-instability motions of a dynamical system characterized by exponential divergences of trajectories leading to chaos (including turbulence as a form of chaos). The formalism involves fictitious control forces that couple the equations of motion of the system with a Liouville equation that describes the evolution of the probability density of errors in initial conditions. These stabilizing forces create a powerful terminal attractor in probability space that corresponds to occurrence of a target trajectory with probability one. The effect in configuration space (ordinary three-dimensional space as commonly perceived) is to suppress exponential divergences of neighboring trajectories without affecting the target trajectory. As a result, the post-instability 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.
Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.
In this work, we investigate the main challenges to prediction of turbulent external flows of practical interest with Reynolds-Averaged Navier–Stokes equations (RANS) and Scale-Resolving Simulation (SRS) models. This represents a crucial step toward further developing and establishing these formulations so they can be confidently utilized in engineering problems without reference data. The study initiates by identifying the major challenges to prediction. A literature review is performed to illustrate their effects in RANS and SRS computations. Afterward, we evaluate the impact of the challenges to prediction by analyzing representative statistically steady and unsteady flows with prominent RANS and SRS methods. These include multiple turbulent viscosity and second-moment RANS closures, and hybrid and bridging SRS models. The results demonstrate the potential of the selected SRS models to predict engineering flows. Yet, they also show the importance of considering the challenges to prediction during the setup and conduction of numerical experiments. These can suppress the advantages of using SRS formulations. The data also indicate that only SRS models can confidently predict statistically unsteady flows. In contrast, the results demonstrate that mean-flow quantities of statistically steady flows can be efficiently calculated with RANS closures, especially second-moment closures. Among the selected SRS methods, bridging models reveal better suited for prediction due to their ability to prevent commutation errors and enable the robust evaluation of numerical and modeling errors. This last property allows the use of a new validation technique that does not require reference data.
Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude S have N = 2S + 1 levels. Alternatively, the spin state is fully characterized by a set of N 2 - 1 local physical observables, which we interpret as generalized spin components. For example, a spin with S = 1 has three dipole components and five quadrupole components. Further, these components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU(N) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a nonequilibrium relaxation process that creates CP 2 skyrmions, which are topological defects with both dipole and quadrupole character.
A set of empirical methods was developed to predict low-speed lift, drag and pitching-moment variations with angle of attack for a class of low aspect ratio irregular planform wings suitable for application to advanced aerospace vehicles. The data base, an extensive series of wind-tunnel tests accomplished by the Langley Research Center of the National Aeronautics and Space Administration, is summarized. The approaches used to analyze the wind tunnel data, the evaluation of previously existing methods, data correlation efforts, and the development of the selected methods are presented and discussed. A summary of the methods is also presented to document the equations, computational charts and design guides which have been programmed for digital computer solution. Comparisons of predictions and test data are presented which show that the new methods provide a significant improvement in capability for evaluating the landing characteristics of advanced aerospace vehicles during the preliminary design phase of the configuration development cycle.
By using the Rayleigh-Ritz method of minimizing the total potential energy of a structural system, combined load (mechanical or thermal load) buckling equations are established for orthotropic rectangular sandwich panels supported under four different edge conditions. Two-dimensional buckling interaction curves and three dimensional buckling interaction surfaces are constructed for high-temperature honeycomb-core sandwich panels supported under four different edge conditions. The interaction surfaces provide easy comparison of the panel buckling strengths and the domains of symmetrical and antisymmetrical buckling associated with the different edge conditions. Thermal buckling curves of the sandwich panels also are presented. The thermal buckling conditions for the cases with and without thermal moments were found to be identical for the small deformation theory. In sandwich panels, the effect of transverse shear is quite large, and by neglecting the transverse shear effect, the buckling loads could be overpredicted considerably. Clamping of the edges could greatly increase buckling strength more in compression than in shear.
USM3D-ME solutions for the High-Fidelity CFD Verification Workshop 2024 are presented. Solutions are computed for the Reynolds-Averaged Navier-Stokes equations using the Spalart-Allmaras one equation turbulence model, SA-neg-QCR2000-R, that is enhanced with a nonlinear correction to the turbulence stresses and a correction for flow rotation. USM3D-ME solutions have been computed for the three verification cases, namely, a Joukowski airfoil, a subsonic three-dimensional flow around an extruded NACA 0012 wing in tunnel, and subsonic flow around a high-lift wing-body configuration. Iterative convergence of USM3D-ME solutions on several grids within various grid families is demonstrated. Grid convergence of integrated forces, pitching moment and sectional variations of surface pressure and skin friction are also presented.
Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.
Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.
Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.
The data obtained at the time of the Thunderstorm Research International Project storm at the Kennedy Space Center on July 11, 1978 are analyzed in a model-independent manner. The data base included data from three Doppler radars, a lightning detection and ranging system and a network of 25 electric field mills, and rain gages. Electric field measurements were used to analyze the charge moments transferred by lightning flashes, and the data were fitted to Weibull distributions; these were used to estimate statistical parameters of the lightning for both intracloud and cloud-to-ground flashes and to estimate the fraction of the flashes which were below the observation threshold. The displacement and the conduction current densities were calculated from electric field measurements between flashes. These values were used to derive the magnitudes and the locations of dipole and monopole generators by least squares fitting the measured Maxwell current densities to the displacement-dominated equations.
The precession and nutation of the Earth's equator arise from solar, lunar, and planetary torques on the oblate Earth. The mean lunar orbit plane is nearly coincident with the ecliptic plane. A small tilt out of the ecliptic is caused by planetary perturbations and the Earth's gravitational harmonic J(sub 2). These planetary perturbations on the lunar orbit result in torques on the oblate Earth which contribute to precession, obliquity rate, and nutation while the J(sub 2) perturbations contribute to precession and nutation. Small additional contributions to the secular rates arise from tidal effects and planetary torques on the Earth's bulge. The total correction to the obliquity rate is -0.024 sec/century, it is an observable motion in space (the much larger conventional obliquity rate is wholly from the motion of the ecliptic, not the equator), and it is not present in the IAU-adopted expressions for the orientation of the Earth's equator. The effects have generally been allowed for in past nutation theories and some precession theories. For the planetary effect, the contributions to the 18.6 yr nutation are -0.03 mas (milliarcseconds) for the in-phase Delta(psi) plus out-of-phase contributions of 0.14 mas in Delta(psi) and -0.03 mas in Delta(sub epsilon). The latter terms demonstrate that out-of-phase contributions can arise by means other than dissipation. The sum of the contributions to the precession rate is considered and the inferred value of the moment of inertia combination (C-A)/C, which is used to scale the coefficients in the nutation series, is evaluated. Using an updated value for the precession rate, the rigid body (C-A)/C =0.003 273 763 4 which, in combination with a satellite-derived J(sub 2), gives a normalized polar moment of inertia C/MR(exp 2) = 0.330 700 7. The planetary contributions to the precession and obliquity rates are not constant for long times causing accelerations in both quantities. Acceleration in precession also arises from tides and changing J(sub 2) Contributions from the improved theory, masses, ecliptic motion, and measured values of the precession rate and obliquity are combined to give expressions (polynomials in time) for precession, obliquity, and Greenwich Mean Sidereal Time.
Germanium selenide (GeSe) is a van der Waals–bonded layered material with promising optoelectronic properties, which has been experimentally synthesized for two-dimensional (2D) semiconductor applications. In the monolayer, due to reduced dimensionality and, thus, screening environment, perturbations such as the presence of defects have a significant impact on its properties. We apply density functional theory and many-body perturbation theory to understand the electronic and optical properties of GeSe containing a single selenium vacancy in the –2 charge state. We predict that the vacancy results in midgap “trap states” that strongly localize the electron and hole density and lead to sharp, low-energy optical absorption peaks below the predicted pristine optical gap. Analysis of the exciton wave function reveals that the 2D Wannier-Mott exciton of the pristine material monolayer is strongly modified by the presence of the defect. The lowest-energy exciton wave function is highly localized around the defect, with the Bohr radius reduced by a factor of 4 when compared to the pristine material and a dipole moment along the out-of-plane axis due to the defect-induced symmetry breaking. Altogether, these results suggest that the vacancy significantly modifies the excitonic properties of the system, demonstrating the importance of considering defects in the context of material design.