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 289 records · Page 16

Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver

A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-epsilon two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user-specified region. The latter approach is applied for the solutions obtained using other one-and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.

Pandya, Mohagna J.↗

Application of optical distributed sensing and computation to control of large space structures

A real time holographic sensing technique is introduced and its advantages are investigated from the filtering and control point of view. A feature of holographic sensing is its capability to make distributed measurements of the position and velocity of moving objects, such as a vibrating flexible space structure. This work is based upon the distributed parameter models of linear time invariant systems, particularly including the linear oscillator equations describing the vibration of large flexible space structures. The general conclusion is that application of optical distributed sensors bring gains in the situation where Kalman filtering is necessary for state estimation. In this case, both steady state and transient filtering error covariance become smaller. This in turn results in smaller cost in the LQG problem.

Balakrishnan, A. V.↗

Identification of time-varying structural dynamic systems - An artificial intelligence approach

An application of the artificial intelligence-derived methodologies of heuristic search and object-oriented programming to the problem of identifying the form of the model and the associated parameters of a time-varying structural dynamic system is presented in this paper. Possible model variations due to changes in boundary conditions or configurations of a structure are organized into a taxonomy of models, and a variant of best-first search is used to identify the model whose simulated response best matches that of the current physical structure. Simulated model responses are verified experimentally. An output-error approach is used in a discontinuous model space, and an equation-error approach is used in the parameter space. The advantages of the AI methods used, compared with conventional programming techniques for implementing knowledge structuring and inheritance, are discussed. Convergence conditions and example problems have been discussed. In the example problem, both the time-varying model and its new parameters have been identified when changes occur.

Glass, B. J.↗

User's Guide for ENSAERO_FE Parallel Finite Element Solver

A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaces ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multizonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRANKOSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft model. It is used to create the stiffness matrices for each substructure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the substructure boundary nodes.

Eldred, Lloyd B.↗

A complete model of turbulence

A set of constitutive equations suitable for a priori computation of turbulent shear flows has been developed. Since no properties of a given turbulent flow need be known in advance in order to obtain a solution, the equations comprise a complete model of turbulence. Perturbation analysis shows that the model predicts a composite five-layer structure for an incompressible turbulent boundary layer, viz, a defect layer, a law-of-the-wall layer, a viscous sublayer, a near-surface roughness layer, and a viscous superlayer at the boundary-layer edge. Analysis of the defect layer demonstrates the key improvement of the model over its predecessor, the Saffman-Wilcox two-equation model of turbulence. Examination of model-predicted sublayer structure yields model-parameter boundary conditions appropriate for surfaces with roughness and mass injection. Results of numerical computations of compressible and incompressible equilibrium boundary layers show that, for such flows, the model is as accurate as mixing-length theory. Applications to transitional boundary layers and to nonequilibrium relaxation of a boundary layer passing from a rough to a smooth surface indicate that the model's applicability extends far beyond that of mixing-length theory's.

Wilcox, D. C.↗

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗

GVIPS Models and Software

Two reports discuss, respectively, (1) the generalized viscoplasticity with potential structure (GVIPS) class of mathematical models and (2) the Constitutive Material Parameter Estimator (COMPARE) computer program. GVIPS models are constructed within a thermodynamics- and potential-based theoretical framework, wherein one uses internal state variables and derives constitutive equations for both the reversible (elastic) and the irreversible (viscoplastic) behaviors of materials. Because of the underlying potential structure, GVIPS models not only capture a variety of material behaviors but also are very computationally efficient. COMPARE comprises (1) an analysis core and (2) a C++-language subprogram that implements a Windows-based graphical user interface (GUI) for controlling the core. The GUI relieves the user of the sometimes tedious task of preparing data for the analysis core, freeing the user to concentrate on the task of fitting experimental data and ultimately obtaining a set of material parameters. The analysis core consists of three modules: one for GVIPS material models, an analysis module containing a specialized finite-element solution algorithm, and an optimization module. COMPARE solves the problem of finding GVIPS material parameters in the manner of a design-optimization problem in which the parameters are the design variables.

Arnold, Steven M.↗

The effect of an unsteady three-dimensional wake on elastic blade-flapping eigenvalues in hover

This paper describes the formulation of a finite-state inflow model based on an acceleration potential and a helical wake geometry. The states of the model are coefficients of an inflow expansion in terms of a Fourier series (azimuthally) and of special polynomials (radially). The integrals over the wake are done in closed-form to obtain a set of ordinary differential equations for the inflow coefficients. The forcing functions for these equations are generalized forces which are integrals of the blade loading exactly as in structural dynamics. This model implicitly includes (for the hover case) Prandtl-Goldstein tip losses, dynamic inflow, and Theodorsen/Loewy lift deficiency. Thus, it is a fully three-dimensional unsteady wake model. Here, this model is coupled with elastic-blade equations in hover and eigenvalues are found. The results show that the three-dimensional wake has a large effect on the flap damping of all modes.

Peters, David A.↗

Aeroelastic analysis of launch vehicles in transonic flight

A numerical study of the aeroelastic stability of launch vehicles in transonic flight is performed. The unsteady aerodynanic flowfield about the vehicle is simulated using the three dimensional Reynolds-Averaged Navier-Stokes equations with the account for turbulent mixing done by a two-layer algebraic eddy viscosity model. The structural-dynamic formulation considers free-free flexural vibration of an elongated beam with variable properties and the equations are cast in modal form. The two sets of equations are integrated simultaneously in time, after a initial perturbation, to ascertain the aeroelastic stability of a given configuration. The method propossed is described in detail, and results that illustrate its application are presented.

Azevedo, Joao Luiz F.↗

Modal simulation of gearbox vibration with experimental correlation

A newly developed global dynamic model was used to simulate the dynamics of a gear noise rig at NASA Lewis Research Center. Experimental results from the test rig were used to verify the analytical model. In this global dynamic model, the number of degrees of freedom of the system are reduced by transforming the system equations of motion into modal coordinates. The vibration of the individual gear-shaft system are coupled through the gear mesh forces. A three-dimensional, axial-lateral coupled, bearing model was used to couple the casing structural vibration to the gear-rotor dynamics. The coupled system of modal equations is solved to predict the resulting vibration at several locations on the test rig. Experimental vibration data was compared to the predictions of the global dynamic model. There is excellent agreement between the vibration results from analysis and experiment.

Choy, Fred K.↗

Modal simulation of gearbox vibration with experimental correlation

A newly developed global dynamic model was used to simulate the dynamics of a gear noise rig at NASA Lewis Research Center. Experimental results from the test rig were used to verify the analytical model. In this global dynamic model, the number of degrees of freedom of the system are reduced by transforming the system equations of motion into modal coordinates. The vibration of the individual gear-shaft system are coupled through the gear mesh forces. A three-dimensional, axial-lateral coupled, bearing model was used to couple the casing structural vibration to the gear-rotor dynamics. The coupled system of modal equations is solved to predict the resulting vibration at several locations on the test rig. Experimental vibration data was compared to the predicitions of the global dynamic model. There is excellent agreement between the vibration results from analysis and experiment.

Choy, Fred K.↗

Analytic fits to atom-in-jellium shear modulus predictions

Fits to AJ shear modulus calculations from ρ0 to the limit of the corresponding AJ EOS table (usually 1000 ρ0), adjusted to reproduce STP shear modulus where possible. Fitting minimized the fractional difference, so low-pressure points would be reproduced accurately despite the wide-ranging tabulation. AJ calculations were at T=0. Isochoric variation G(T) has been small for all cases spot-tested so far, and may have either sign. The fitted equation usually matches the AJ data to within a few percent. Between numerical noise in the AJ calculation and probably-physical structure not captured by the equation, the deviation could be up to 20% in some places in most models, and 30% in a few. AJ does not capture structural phase transitions or phases stabilized by directional bonds, where the shear modulus may vary by a greater amount. Where possible, STP ρ,G were used as parameters, and low pressure AJ points were de-weighted or removed if necessary. Otherwise, such as where AJ fails to capture solid phases with a significantly different shear modulus, the AJ data were fitted as far down in pressure as possible, G0 was also fitted if necessary, and ρ0 was also adjusted if needed to keep G0>0. The resulting models are not intended for use at low pressure, though some are probably adequate for practical purposes

36 MATERIALS SCIENCE↗

Electron Beam Induced Luminescence of SiO2 Optical Coatings

Optical coatings of disordered thin film SiO,/SiO, samples undergoing electron beam bombardment exhibited a cathodoluminescence effect. Measurements were made of the absolute illuminance and emission spectra as functions of incident electron energy (0.2 keV to 30 keV). flux (approx. 0.1 nA/sq cm to > 100 nA/sq cm) and power (<50 nW/sq cm to > 3mW/sq cm) for coatings (approx.60 nm to approx.200 nm thick) on reflective metal substrates over a range of sample temperatures (approx.40 K to approx.400 K) and emission wavelengths (approx.260 nm to approx.5000 nm). Illuminance reached a plateau saturation intensity I(sub sat) approx. 10(exp -4)W/sq cm. scaling as [I(sub mc)/(I + I(sub mc)/I(sub sat)] as incident power I(sub mc) approached I(sub sat). Well below I(sub sat), illuminance scaled linearly with incident power for lower-energy non-penetrating radiation that deposited all of its energy in the dielectric film. At higher incident energies (> 1-5 keV) electrons penetrate the dielectric, depositing only a fraction of their energy in the dielectric film, where deposited energy scales as the ratio of film thickness to penetration range. Since electron range exhibits a peak below approx. 1 keV, the illuminance decreased with increasing energy for penetrating radiation. The temperature dependence of the illuminance. in general, increased with decreasing temperature. Changes in apparent color with temperature were also observed. Four bands in the UV/VIS/NIR range from 300 nm to 1000 nm were observed in spectral measurements; no significant luminescence was detected in a range from approx. 1000 nm to approx. 5000 nm. These different bands either increased or decreased in intensity and shifted slightly in energy with decreasing temperature. Our measurements are explained in terms of simple disordered band theory models based on atomic-scale stochastic transport of electrons in highly disordered insulating materials. The theoretical models provide a fundamental basis for understanding the dependence of cathodoluminescence on irradiation time, incident flux and energy. and sample thickness and temperature. The observed luminescence occurs when an incident high energy electron excites a valence band electron into the conduction band. Following the excitation. a complex series of transitions take place between the extended conduction states and the long-lived localized trapped states (or chromophores) associated with structural or compositional defects in the dielectric materials. Thc transport equations and models of thc spatial and energy distribution of the trapped states below the mobility edge, their occupancy, and transitions provide microscopic scale insights into this macroscopic phenomenon. Our observations and theoretical models are also compared with similar previous room temperature studies of similar SiO2 films.

Dennison, J. R.↗

Buckling of imperfect periodic lattice structures

A simplified buckling analysis is presented for a family of periodic lattice structures such as those proposed for large space structures. A transcendental 6 x 6 matrix of eigenvalues is shown to be sufficient for modeling buckling behavior because member stiffnesses are based on an exact solution of the beam-column equation. Exact stiffnesses are derived for a curved member, thus allowing modeling of imperfect lattice structures. Comparisons of predictions of the lattice model with those available from shell and beam theory underscore the inaccuracies introduced by treating the lattice structure as a continuum. Sample calculations are provided for an isogrid cylinder and a three element double-laced truss.

Anderson, M. S.↗

Large-scale processes in the solar nebula

Theoretical models of the structure of a minimum mass solar nebula should be able to provide the physical context to help evaluate the efficacy of any mechanism proposed for the formation of chondrules or Ca, Al-rich inclusions (CAI's). These models generally attempt to use the equations of radiative hydrodynamics to calculate the large-scale structure of the solar nebula throughout the planet-forming region. In addition, it has been suggested that chondrules and CAI's (=Ch&CAI's) may have been formed as a direct result of large-scale nebula processing such as passage of material through high-temperature regions associated with the global structure of the nebula. In this report we assess the status of global models of solar nebula structure and of various related mechanisms that have been suggested for Ch and CAI formation.

Boss, A. P.↗

The Numerical Analysis of a Turbulent Compressible Jet

A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).

DeBonis, James R.↗

Structural similitude and design of scaled down laminated models

The excellent mechanical properties of laminated composite structures make them prime candidates for wide variety of applications in aerospace, mechanical and other branches of engineering. The enormous design flexibility of advanced composites is obtained at the cost of large number of design parameters. Due to complexity of the systems and lack of complete design based informations, designers tend to be conservative in their design. Furthermore, any new design is extensively evaluated experimentally until it achieves the necessary reliability, performance and safety. However, the experimental evaluation of composite structures are costly and time consuming. Consequently, it is extremely useful if a full-scale structure can be replaced by a similar scaled-down model which is much easier to work with. Furthermore, a dramatic reduction in cost and time can be achieved, if available experimental data of a specific structure can be used to predict the behavior of a group of similar systems. This study investigates problems associated with the design of scaled models. Such study is important since it provides the necessary scaling laws, and the factors which affect the accuracy of the scale models. Similitude theory is employed to develop the necessary similarity conditions (scaling laws). Scaling laws provide relationship between a full-scale structure and its scale model, and can be used to extrapolate the experimental data of a small, inexpensive, and testable model into design information for a large prototype. Due to large number of design parameters, the identification of the principal scaling laws by conventional method (dimensional analysis) is tedious. Similitude theory based on governing equations of the structural system is more direct and simpler in execution. The difficulty of making completely similar scale models often leads to accept certain type of distortion from exact duplication of the prototype (partial similarity). Both complete and partial similarity are discussed. The procedure consists of systematically observing the effect of each parameter and corresponding scaling laws. Then acceptable intervals and limitations for these parameters and scaling laws are discussed. In each case, a set of valid scaling factors and corresponding response scaling laws that accurately predict the response of prototypes from experimental models is introduced. The examples used include rectangular laminated plates under destabilizing loads, applied individually, vibrational characteristics of same plates, as well as cylindrical bending of beam-plates.

Simitses, G. J.↗

Advanced Turbulence Models for Large-Scale Atmospheric Boundary Layer Flows

We present high-fidelity large-eddy-simulation (LES) modeling approaches for the turbulent atmospheric boundary layer (ABL) flows. Wind energy is a prime example of an application driven by ABL. Generation of electrical energy from farms of wind turbines at night in the stable ABL is a particularly interesting situation. In this report, we consider the well-known GEWEX (Global Energy and Water Cycle Experiment) Atmospheric Boundary Layer Study (GABLS) stably stratified benchmark LES case. We use a high-order spectral element code Nek5000/RS, which is supported under the DOE's Exascale Computing Project (ECP) Center for Efficient Exascale Discretizations (CEED) project, targeting application simulations on various acceleration-device based exascale computing platforms. In our earlier ANL report, we demonstrated our newly developed subgrid-scale (SGS) models based on high-pass filter (HPF), mean-field eddy viscosity (MFEV), and Smagorinsky (SMG) with no-slip and traction boundary conditions, provided with low-order statistics, convergence and turbulent structure analysis. In this report, we extend the range of our SGS modeling approaches in the context of the mean-field eddy viscosity (MFEV), to include the solution of an SGS turbulent kinetic energy equation (TKE). We demonstrate the model fidelity of Nek5000/RS in comparison to that of AMR-Wind, a block-structured second-order finite-volume code with adaptive-mesh-refinement capabilities, with which we studied scaling performance for both codes in comparison on DOE's leadership computing platforms.

17 WIND ENERGY↗