Search NASA⌕ Search

SEARCH · Search NASA

Results for “structural equation model”

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

Control/structure interaction - Effects of actuator dynamics

The present evaluation of the effects of including actuator dynamics in the formulation of simple control laws for vibration suppression leads, theoretically, to a prediction of substantial interaction between actuator dynamics, structural dynamics, and control-law choices. Attention is given to the problem of closed-loop instability, under the assumptions that (1) the structure is modeled by an asymptotically stable system, (2) vibration-suppression is the purpose of control, and (3) the structure is well modeled by a finite-dimensional linear set of equations with viscous damping, such as a FEM model.

Inman, Daniel J.↗

Multigrid calculations of 3-D turbulent viscous flows

Convergence properties of a multigrid algorithm, developed to calculate compressible viscous flows, are analyzed by a vector sequence eigenvalue estimate. The full 3-D Reynolds-averaged Navier-Stokes equations are integrated by an implicit multigrid scheme while a k-epsilon turbulence model is solved, uncoupled from the flow equations. Estimates of the eigenvalue structure for both single and multigrid calculations are compared in an attempt to analyze the process as well as the results of the multigrid technique. The flow through an annular turbine is used to illustrate the scheme's ability to calculate complex 3-D flows.

Yokota, Jeffrey W.↗

Control law synthesis and stability robustness improvement using constrained optimization techniques

The present generic optimization procedure for a continuous or discrete control law (of arbitrary order), which will be applicable to a multiinput-multioutput system, is upon constraining used to satisfy conflicting design requirements on the mean-square responses and stability robustness at the plant input and output. The synthesis procedure is especially suitable for flexible airframes and large space structures modeled by a high-order state-space system of equations. Analytical expressions are obtained for the gradients of the cost function, together with design constraints on the mean-square response and minimum singular value.

Mukhopadhyay, Vivekananda↗

Parallel execution of a three-dimensional, chemically reacting, Navier-Stokes code on distributed-memory machines

This paper describes the parallel version of the three-dimensional, chemically reacting, computational fluid dynamics (CFD) code, SPARK. This work was performed on the Intel iPSC/860-based parallel computers. The SPARK code utilizes relatively simple explicit numerical algorithms, but models complex chemical reactions. The code solves the equations over a regular structured mesh so a simple dam decomposition is used to assign work to the individual processors. The explicit nature of the algorithm, combined with the computational intensity of the chemistry calculations, results in a very low communication-to-computation ratio when compared to typical CFD codes. The efficiency of the parallel code is examined and shown to be about 65 percent when the problem size is scaled with the number of processors. Two low-angle wall-jet injection cases are solved to demonstrate the capability of the parallel code for solving large problems efficiently.

Otto, John C.↗

Finite-difference simulation and visualization of elastodynamics in time-evolving generalized curvilinear coordinates

Modeling and simulation of free and forced structural vibrations is essential to an overall structural health monitoring capability. In the various embodiments, a first principles finite-difference approach is adopted in modeling a structural subsystem such as a mechanical gear by solving elastodynamic equations in generalized curvilinear coordinates. Such a capability to generate a dynamic structural response is widely applicable in a variety of structural health monitoring systems. This capability (1) will lead to an understanding of the dynamic behavior of a structural system and hence its improved design, (2) will generate a sufficiently large space of normal and damage solutions that can be used by machine learning algorithms to detect anomalous system behavior and achieve a system design optimization and (3) will lead to an optimal sensor placement strategy, based on the identification of local stress maxima all over the domain.

Kaul, Upender K.↗

Development of computational methods for heavy lift launch vehicles

The research effort has been focused on the development of an advanced flow solver for complex viscous turbulent flows with shock waves. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. A new computer program named CENS3D has been developed for viscous turbulent flows with discontinuities. Details of the code are described in Appendix A and Appendix B. With the developments of the numerical algorithm and dissipation model, the simulation of three-dimensional viscous compressible flows has become more efficient and accurate. The results of the research are expected to yield a direct impact on the design process of future liquid fueled launch systems.

Yoon, Seokkwan↗

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗

Future Probes of the Neutron Star Equation of State Using X-ray Bursts

Observations with NASA s Rossi X-ray Timing Explorer (RXTE) have resulted in the discovery of fast (200 - 600 Hz), coherent X-ray intensity oscillations (hereafter, %urstoscillations ) during thermonuclear X-ray bursts from 12 low mass X-ray binaries (LMXBs). Although many of their detailed properties remain to be fully understood, it is now beyond doubt that these oscillations result from spin modulation of the thermonuclear burst flux from the neutron star surface. Among the new timing phenomena revealed by RXTE the burst oscillations are perhaps the best understood, in the sense that many of their properties can be explained in the framework of this relatively simple model. Because of this, detailed modelling of burst oscillations can be an extremely powerful probe of neutron star structure, and thus the equation of state (EOS) of supra-nuclear density matter. Both the compactness parameter beta = GM/c(sup 2)R, and the surface velocity, nu(sub rot) = Omega(sub spin)R, are encoded in the energy-dependent amplitude and shape of the modulation pulses. The new discoveries have spurred much new theoretical work on thermonuclear burning and propagation on neutron stars, so that in the near future it is not unreasonable to think that detailed physical models of the time dependent flux from burning neutron stars will be available for comparison with the observed pulse profiles from a future, large collecting area X-ray timing observatory. In addition, recent high resolution burst spectroscopy with XMM/Newton suggests the presence of redshifted absorption lines from the neutron star surface during bursts. This leads to the possibility of using large area, high spectral resolution measurements of X-ray bursts as a precise probe of neutron star structure. In this work I will explore the precision with which constraints on neutron star structure, and hence the dense matter EOS, can be made with the implementation of such programs.

Strohmayer, Tod E.↗

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗

Modeling of Sagging for 3D Printed Layers During the Curing Process

Here, a printed layer of silica-enforced poly(dimethylsiloxane)-co-(diphenylsiloxane) is modeled as a two-phase system consisting of air and polymer with an interface set up during the printing process. The structural geometry changes mostly due to the action of surface tension, while all material properties are strongly temperature dependent. Polymer flow is described using equations of the extended Herschel–Bulkley model, with parameters strongly dependent on temperature and degree of curing. Parameters of the model are determined using flow sweep measurements and separate experiments with vertical structure sagging at different temperatures. The curing process is modeled using dependencies between the curing rate, degree of curing, and temperature obtained in studies by differential scanning calorimetry. The developed model is used for simulations of printed structure deformations with different initial and boundary conditions.

36 MATERIALS SCIENCE↗

Effect of turbulence models on criticality conditions in swirling flows

The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Gatski, Thomas B.↗

Shock-wave structure in a partially ionized gas

The structure of a steady plane shock in a partially ionized gas has been investigated using the Boltzmann equation with a kinetic model as the governing equation and the discrete ordinate method as a tool. The effects of the electric field induced by the charge separation on the shock structure have also been studied. Although the three species of an ionized gas travel with approximately the same macroscopic velocity, the individual distribution functions are found to be very different. In a strong shock the atom distribution function may have double peaks, while the ion distribution function has only one peak. Electrons are heated up much earlier than ions and atoms in a partially ionized gas. Because the interactions of electrons with atoms and with ions are different, the ion temperature can be different from the atom temperature.

Lu, C. S.↗

Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas (Final Technical Report)

This document details the work accomplished during the period 04/15/2012-01/31/2020 for the collaborative project “Interactions Between Microturbulence and Magnetohydrodynamic Instabilities in Magnetized Plasmas,” and serves as the final report on the University of California, San Diego part of this collaborative project for Department of Energy Office of Science - Fusion Energy Sciences Grant No. DE- SC0007783. This project focused on the development of analytic and computational tools, seeking to better understand the mutual interactions between small-scale temperature-gradient driven microturbulence and macroscopic magnetohydrodynamic (MHD) instabilities, particularly tearing modes, in magnetically confined plasmas relevant for the generation of fusion energy. This project advanced understanding through numerical simulations of different sets of model equations, corresponding to different microturbulence instabilities. Questions addressed included how the magnitudes, spatial structures, and scaling's of turbulent forcing's and drags impact magnetic island dynamics, as well as quantifying the “region of influence” of an island on the microturbulence dynamics. A key finding of the project is that the microturbulence fluctuations and associated cross-field particle and energy fluxes exhibit spatial asymmetries for sufficiently large islands. The results have improved our fundamental understanding of plasma dynamics, motivated new experimental investigations, and will help inform future fusion predictive modeling efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of Two Statistical Models for Low Boom Dose-response Relationships with Correlated Responses

This study compares two statistical models to construct summary dose-response curves for low boom community noise surveys. Data from two NASA field surveys are used that consist of multiple responses per survey participant. These data require an approach that accounts for the correlation among repeated annoyance observations from the same participant. A multilevel model accounts for the correlation by allowing estimated parameters to vary with each survey participant. On the other hand, a population average model utilizes generalized estimating equations and accounts for the correlation via a userspecified within-subject correlation structure. A visual comparison of the dose-response curves for these two methods reveals similar results. When comparing specific points along the summary curves, the multilevel model yields more precise confidence bounds than the population average model. The similarity between the summary curves derived from each model lends validity to both approaches for approximating a population representative summary curve, though modeling assumptions may lend favor to the multilevel logistic modeling approach over the population average model.

dose-response↗

Evolution of the solar nebula. I - Nonaxisymmetric structure during nebula formation

Numerical solutions of the equations of hydrodynamics, gravitation, and radiative transfer in three spatial dimensions are used to model the formation and time evolution of the early solar nebula in order to learn whether or not gravitational torques between nonaxisymmetric structures in the solar nebula can transport angular momentum rapidly enough to produce nebula clearing on astronomically indicated (10 to the 5 to 10 to the 7 yr) time scales. The models involve solutions for the collapse of spherical clouds with assumed initial density and rotation profiles onto protosuns of variable mass. Most of the models assume uniform initial density and rotation, and have variations in the initial parameters of cloud mass, cloud rotation rate, and protosun mass which are chosen to simulate a range of possible phases of early solar nebula evolution. The models show little tendency for directly forming small numbers of giant gaseous protoplanets through gaseous gravitational instability.

Boss, Alan P.↗

Aeroelastic oscillations of a pitching flexible wing with structural geometric nonlinearities: Theory and numerical simulation

In this paper, we focus on the derivation of an analytical model of the aeroelastic dynamics of an elastically mounted flexible wing. This equations of motion obtained serve to help understand the behaviour of the aeroelastic wind tunnel setup in question, which consists of a rectangular wing with a uniform NACA 0012 airfoil profile, whose base is free to rotate rigidly about a longitudinal axis. Of particular interest are the structural geometric nonlinearities primarily introduced by the coupling between the rigid body pitch degree-of-freedom and the continuous system. A coupled system of partial differential equations (PDEs) coupled with an ordinary differential equation (ODE) describing axial-bending-bending-torsion-pitch motion is derived using Hamilton's principle. A finite dimensional approximation of the system of coupled differential equations is obtained using the Galerkin method, leading to a system of coupled nonlinear ODEs. Subsequently, these nonlinear ODEs are solved numerically using Houbolt's method. The results that are obtained are verified by comparison with the results obtained by direct integration of the equations of motion using a finite difference scheme. Adopting a linear unsteady aerodynamic model, it is observed that the system undergoes coalescence flutter due to coupling between the rigid body pitch rotation dominated mode and the first flapwise bending dominated mode. Finally, the behaviour of the limit cycle oscillations is primarily influenced by the structural geometric nonlinear terms in the coupled system of PDEs and ODE.

42 ENGINEERING↗