Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Static analysis”

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 1,639 records · Page 91

Low-speed Aerodynamic Investigations of a Hybrid Wing Body Configuration

Two low-speed static wind tunnel tests and a water tunnel static and dynamic forced-motion test have been conducted on a hybrid wing-body (HWB) twinjet configuration. These tests, in addition to computational fluid dynamics (CFD) analysis, have provided a comprehensive dataset of the low-speed aerodynamic characteristics of this nonproprietary configuration. In addition to force and moment measurements, the tests included surface pressures, flow visualization, and off-body particle image velocimetry measurements. This paper will summarize the results of these tests and highlight the data that is available for code comparison or additional analysis.

Vicroy, Dan D.↗

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗

Optimized Off-Design Performance of Flexible Wings with Continuous Trailing-Edge Flaps

This work assesses the potential aerodynamic performance benefits of a variable-camber, continuous-trailing-edge flap system on a generic transport aircraft at off-design conditions. A process to optimize transport wings while addressing static aeroelastic effects is presented. To establish a proper baseline, a transport wing is first aerodynamically optimized at a mid-cruise flight condition using an inviscid, aeroelastic analysis tool. The optimized wing is then analyzed at off-design cruise conditions. The optimization is repeated at these off-design conditions to determine how much performance is lost by the wing optimized solely for the mid-cruise condition. The full-span flap system is then adapted to improve performance of the mid-cruise-optimized wing at these off-design conditions. The measured improvement is quantified by a comparison with wings designed specifically for the off-design conditions. To evaluate the effects of aeroelasticity on the effectiveness of the flap system, this entire process is performed on both a conventionally stiff wing and a modern, more flexible wing. The results indicate that the flap system allows for recovery of near-optimal performance throughout cruise. The flap system is found to be advantageous even for modern wings with increased flexibility.

Rodriguez, David L.↗

Improvements in Block-Krylov Ritz Vectors and the Boundary Flexibility Method of Component Synthesis

A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence relationship, proposed by Wilson, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is based upon the boundary flexibility vectors of the component. Improvements have been made in the formulation of the initial seed to the Krylov sequence, through the use of block-filtering. A method to shift the Krylov sequence to create Ritz vectors that will represent the dynamic behavior of the component at target frequencies, the target frequency being determined by the applied forcing functions, has been developed. A method to terminate the Krylov sequence has also been developed. Various orthonormalization schemes have been developed and evaluated, including the Cholesky/QR method. Several auxiliary theorems and proofs which illustrate issues in component mode synthesis and loss of orthogonality in the Krylov sequence have also been presented. The resulting methodology is applicable to both fixed and free- interface boundary components, and results in a general component model appropriate for any type of dynamic analysis. The accuracy is found to be comparable to that of component synthesis based upon normal modes, using fewer generalized coordinates. In addition, the block-Krylov recurrence algorithm is a series of static solutions and so requires significantly less computation than solving the normal eigenspace problem. The requirement for less vectors to form the component, coupled with the lower computational expense of calculating these Ritz vectors, combine to create a method more efficient than traditional component mode synthesis.

Carney, Kelly Scott↗

Friction Experiment Data Analysis

The Spent Fuel Waste Disposition (SFWD) program under the U.S. Department of Energy (DOE) is planning a seismic shake table test of full-scale dry storage systems of spent nuclear fuel (SNF) to close the gap related to the seismic loads on the fuel assemblies in dry storage systems. This test will allow for quantifying the strains and accelerations on surrogate fuel assembly during representative earthquakes. A concrete layer will be installed on the shake table before the test to simulate conditions representative of an ISFSI pad. In the shake table tests with the vertical cask, the cask will be free-standing because this is representative of all, except two, ISFSIs in the U.S. with vertical dry storage casks. The static and dynamic friction coefficients between the steel bottom of the cask and the concrete layer on the shake table are important parameters that will affect cask behavior during the test. These parameters must be known for the pre- and post-test modelling, data analysis, and model validation. The friction experiment was performed at the Engineering Department of the University of New Mexico (UNM) to determine the friction coefficients between a steel plate with the same finish as the bottom of the vertical cask manufactured for the test and different concrete surfaces. In this experiment the steel plate was fixed and the concrete sample was pulled over the plate with a constant displacement rate using an MTS machine. This allowed for collecting continuous horizontal force data over the length of the steel plate. Four displacement rates and three vertical loads were used. The tests were performed with four concrete blocks with different degrees of the surface roughness - light sandblast, light to medium sandblast, medium bush hammer, and heavy sandblast. The total number of tests was 48. The data were used to calculate static and dynamic friction coefficients.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probabilistic Analysis and Design of a Raked Wing Tip for a Commercial Transport

An approach for conducting reliability-based design and optimization (RBDO) of a Boeing 767 raked wing tip (RWT) is presented. The goal is to evaluate the benefits of RBDO for design of an aircraft substructure. A finite-element (FE) model that includes eight critical static load cases is used to evaluate the response of the wing tip. Thirteen design variables that describe the thickness of the composite skins and stiffeners are selected to minimize the weight of the wing tip. A strain-based margin of safety is used to evaluate the performance of the structure. The randomness in the load scale factor and in the strain limits is considered. Of the 13 variables, the wing-tip design was controlled primarily by the thickness of the thickest plies in the upper skins. The report includes an analysis of the optimization results and recommendations for future reliability-based studies.

Mason Brian H.↗

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗

Analysis of shell-type structures subjected to time-dependent mechanical and thermal loading

A general mathematical model and solution methodologies are being developed for analyzing structural response of thin, metallic shell-type structures under large transient, cyclic, or static thermomechanical loads. Among the system responses, which were associated with these load conditions, were thermal buckling, creep buckling, and ratcheting. Thus, geometric as well as material-type nonlinearities (of high order) can be anticipated and must be considered in the development of the mathematical model. Furthermore, this must also be accommodated in the solution process.

Simitses, G. J.↗

Coupled flow, thermal and structural analysis of aerodynamically heated panels

A finite element approach to coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.

Thornton, Earl A.↗

Coupled flow, thermal and structural analysis of aerodynamically heated panels

A finite element approach for coupling flow, thermal and structural analyses of aerodynamically heated panels is presented. The Navier-Stokes equations for laminar compressible flow are solved together with the energy equation and quasi-static structural equations of the panel. Interactions between the flow, panel heat transfer and deformations are studied for thin stainless steel panels aerodynamically heated by Mach 6.6 flow.

Thornton, Earl A.↗

Estimates of tropical analysis differences in daily values produced by two operational centers

To assess the uncertainty of daily synoptic analyses for the atmospheric state, the intercomparison of three First GARP Global Experiment level IIIb datasets is performed. Daily values of divergence, vorticity, temperature, static stability, vertical motion, mixing ratio, and diagnosed diabatic heating rate are compared for the period of 26 January-11 February 1979. The spatial variance and mean, temporal mean and variance, 2D wavenumber power spectrum, anomaly correlation, and normalized square difference are employed for comparison.

Kasahara, Akira↗

Performance analysis of a large-grain dataflow scheduling paradigm

A paradigm for scheduling computations on a network of multiprocessors using large-grain data flow scheduling at run time is described and analyzed. The computations to be scheduled must follow a static flow graph, while the schedule itself will be dynamic (i.e., determined at run time). Many applications characterized by static flow exist, and they include real-time control and digital signal processing. With the advent of computer-aided software engineering (CASE) tools for capturing software designs in dataflow-like structures, macro-dataflow scheduling becomes increasingly attractive, if not necessary. For parallel implementations, using the macro-dataflow method allows the scheduling to be insulated from the application designer and enables the maximum utilization of available resources. Further, by allowing multitasking, processor utilizations can approach 100 percent while they maintain maximum speedup. Extensive simulation studies are performed on 4-, 8-, and 16-processor architectures that reflect the effects of communication delays, scheduling delays, algorithm class, and multitasking on performance and speedup gains.

Young, Steven D.↗

COMDES-MELT: A Turbofan Engine Icing Risk Analysis Tool, User’s Manual

The computational tool COMDES-MELT was developed to predict the susceptibility of turbofan engines to ice accretion when flying in an ice crystal environment. COMDES-MELT is a first-generation computational tool that can estimate the conditions at which the accretion of ice can occur, location within the turbofan engine compression system, and at which engine operating points within the flight trajectory. This is accomplished by tracking the parameters of static wet bulb temperature, ice particle melt ratio, and the ice-water to airflow rates ratio as ice crystals are ingested into the engine compression system. It has been used successfully to predict engine operating points for simulation in an altitude wind tunnel with ice crystal ingestion. These engine tests occur in a simulated ice crystal cloud environment. This predictive analysis has been used to generate the test matrix in preparation for tests in the Propulsion Systems Laboratory, an altitude test facility at the NASA Glenn Research Center. Furthermore, due to the speed of the code during an engine icing test, adjustments to the test matrix can be provided to the test engineer in real time to focus on a particular area of ice accretion susceptibility. The use of the COMDES-MELT code will be examined in this manual.

Turbofan engine↗

Estimating Spatial Distribution Impacts of Rooftops Solar PV on Dynamic Hosting Capacity Evaluation for a Real Distribution Feeder

This paper models the sensitivity of dynamic hosting capacity analysis to the spatial distribution of distributed PV scenarios, random, close and far from substation deployments, at various penetration levels on a real distribution feeder, in a quasi-static time series (QSTS) power flow simulation using high resolution time-aware metrics. These spatial distribution scenarios were chosen to capture a wide range of potential operation impacts of these deployments using a set of thermaland voltage-based metrics such as instantaneous violations and moving averages for both thermal and voltage constraints. This study uses actual load and PV data up-scaled from 1-hour time step to 1-minute resolution to fully characterize the interaction between the daily changes in load and PV output, and their impacts on distribution system operations.

distribution system↗

RAFT (Response Amplitudes of Floating Turbines) [SWR-22-40]

RAFT is a Python code for frequency-domain analysis of floating wind turbines. RAFT works in the frequency domain to provide efficient computations of a floating wind system's response spectra accounting for platform hydrodynamics, quasi-static mooring reactions, rotor aerodynamics, and linearized control. It can be used to calculate response amplitude operators (RAO's), power spectral densities, mean/static properties, and response metrics based on mean and root-mean-square values of each output. RAFT serves as the Level 1 model in the WEIS (Wind Energy with Integrated Servo-control) toolkit, and features an OpenMDAO wrapper.

Hall, Matthew↗

Wabash CarbonSAFE. Geologic Analysis of the Potosi Dolomite Reservoir Interval and Potential Confining Units, at the Wabash #1 Well, Vigo Co., Indiana, Subtask 7.2, Technical Report

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to evaluate the feasibility of commercial-scale CO 2 storage near the site. The Mt. Simon Sandstone was the initial target for storage evaluation and was found to have generally poor reservoir qualities in the Wabash #1 well. Focus was then placed on the secondary target, the Potosi Dolomite. For storage and containment evaluation an extensive suite of geophysical logs was collected throughout the borehole, and whole cores were obtained from the Mt. Simon Sandstone, Eau Claire Formation, and Maquoketa Group. The Maquoketa Group is identified as a regionally extensive secondary sealing interval above the Mt. Simon Sandstone and as a primary seal for the Potosi Dolomite storage complex. Due to concerns over potential lost circulation zones in the Potosi Dolomite during drilling, no core samples or Formation Micro Imager (FMI) logs were acquired in the Potosi Dolomite. However, cased hole well testing was performed over a 20 feet (6 m) interval within the Potosi Dolomite during which a fluid swab sample was obtained. Approximately 35 miles of 2D seismic information was acquired in the project area near Wabash #1 to aid in reservoir and caprock characterization. The highly porous and permeable zones in the Potosi Dolomite and lower Oneota Dolomite observed within the Wabash #1 well are, for the purpose of this report, considered the same reservoir package and referred to as the Potosi reservoir interval. This report presents a regional depositional context for the Potosi reservoir interval, the overlying confining zones in the Knox Group, and Maquoketa Group primary seal. This context is developed through detailed information obtained from the Wabash #1 well such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Potosi reservoir interval. Information discussed in this report was compiled as a basis for the development of static and dynamic geologic models of the Potosi Dolomite storage complex. The static and dynamic modeling of CO 2 injection in the Potosi reservoir interval are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Minimum weight design of structures via optimality criteria

The state of the art of automated structural design through the use of optimality criteria, with emphasis on aerospace applications is reviewed. Constraints on stresses, displacements, and buckling strengths under static loading, as well as lower bound limits on natural frequencies and flutter speeds are presented. It is presumed that the reader is experienced in finite element methods of analysis, but is not familiar with optimal design techniques.

Kiusalaas, J.↗

Aerodynamic characteristics of a 1/4 scale powered helicopter model with a V-type empennage

An investigation was made in the Langley V/STOL tunnel to determine rotor induced effects on a 1/4-scale helicopter model with a conventional empennage and also a V-type empennage with dihedral angles of 45 deg, 50 deg, 55 deg, and 60 deg. Static longitudinal and lateral directional stability data are presented for rotor advance ratios of 0.057, 0.102, and 0.192 in level flight and climb attitudes. The data are presented without analysis or discussion.

Freeman, C. E.↗