Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Geometry”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

LivermorE Al Projector for Computed Tomography Tasks

With recent computed tomography (CT) efforts using Artificial Intelligence (AI) and Deep Learning (DL)techniques, there is a strong need for differentiable forward projection models that can be integrated into existing DL frameworks. We developed a pytorch-based package library providing differentiable forward and back projection functions and classes to facilitate forward and back propagation of CT operations in the training procedure. This forward projectors support three CT projection geometries: cone, parallel and modular beams. This package can be used with both CPU and GPU with CUDA.

Kim, Hyojin↗
NASA NTRS1983

Computation of shock wave/target interaction

Computational results of shock waves impinging on targets and the ensuing diffraction flowfield are presented. A number of two-dimensional cases are computed with finite difference techniques. The classical case of a shock wave/cylinder interaction is compared with shock tube data and shows the quality of the computations on a pressure-time plot. Similar results are obtained for a shock wave/rectangular body interaction. Here resolution becomes important and the use of grid clustering techniques tend to show good agreement with experimental data. Computational results are also compared with pressure data resulting from shock impingement experiments for a complicated truck-like geometry. Here of significance are the grid generation and clustering techniques used. For these very complicated bodies, grids are generated by numerically solving a set of elliptic partial differential equations.

Mark, A.↗
NASA NTRS1995

Radial alignment simulation of solar wind streams observed by Pioneers 10 and 11 in 1974

A particularly favorable radial lineup between spacecraft in Earth orbit and Pioneers 10 and 11 (near the ecliptic plane at approximately 5.5 AU and approximately 4.5 AU, respectively) occurred in mid-1974, when the solar corona was in a stable and well-defined warped-dipole configuration. The radial alignment study reported here differs from previous applications of the technique in two respects: (1) It is the first time a two-dimensional (2-D) MHD model has been tested over such a lengthy propagation interval; the 2-D capability is crucial for treatment of the nonradial shearing motions occurring across the stream interface. (2) The three-dimensional (3-D) structure observed in the white light corona at that time is related to the systematic patterns of nonradial flow deflections appearing at the Pioneer corotating interaction region (CIR) fronts. Comparison of predicted and observed flows for pairs of streams in two successive rotations reveals that when the parent coronal hole projects far across the solar equator in a predominantly north-south orientation (i.e., nearest the 2-D idealization), the mapping is accurate down to details of the flow structures. But where the spacecraft tracks along a latitudinal boundary of a hole or the associated stream front is inclined at a shallow angle to the equator, the numerical projections deviate systematically from the observations. Among the sources of error are 3-D dynamical interactions neglected in the model, differential rotation effects, and slow temporal evolution of the coronal structures. A characteristic pattern of north-south and east-west deflections indicative of the 3-D geometry appears across the CIRs, but these patterns evidently reflect primarily the local, as opposed to global, orientation of the 3-D stream fronts. Such patterns appear common in CIRs observed by Pioneer during this period. These findings thus hold special relevance for the analysis of Ulysses observations, since the present coronal configuration is similar to that of 1974.

Pizzo, V. J.↗
NASA NTRS1997

Compensation for Spherical Geometric and Absorption Effects on Lower Thermospheric Emission Intensities Derived from High Earth Orbit Images

Remote sensing of the atmosphere from high earth orbit is very attractive due to the large field of view obtained and a true global perspective. This viewpoint is complicated by earth curvature effects so that slant path enhancement and absorption effects, small from low earth orbit, become dominant even at small nadir view angles. The effect is further complicated by the large range of local times and solar zenith angles in a single image leading to a modulation of the image intensity by a significant portion of the diurnal height variation of the absorbing layer. The latter effect is significant in particular for mesospheric, stratospheric and auroral emissions due to their depth in the atmosphere. As a particular case, the emissions from atomic oxygen (130.4 and 135.6 nm) and molecular nitrogen (two LBH bands, LBHS from 140 to 160 nm and LBHL from 160 to 180 nm) as viewed from the Ultraviolet Imager (UVI) are examined. The LBH emissions are of particular interest since LBHS has significant 02 absorption while LBHL does not, In the case of auroral emissions this differential absorption, well examined in the nadir, gives information about the height of the emission and therefore the energy of the precipitating particles. Using simulations of the viewing geometry and images from the UVI we examine these effects and obtain correction factors to adjust to the nadir case with a significant improvement of the derived characteristic energy. There is a surprisingly large effect on the images from the 02 diurnal layer height changes. An empirical compensation to the nadir case is explored based on the local nadir and local zenith angles for each portion of the image. These compensations are demonstrated as applied to the above emissions in both auroral and dayglow images and compared to models. The extension of these findings to other instruments, emissions and spectral regions is examined.

Swift, W.↗
NASA NTRS2003

Low-Power SOI CMOS Transceiver

The work aims at developing a low-power Silicon on Insulator Complementary Metal Oxide Semiconductor (SOI CMOS) Transceiver for deep-space communications. RF Receiver must accomplish the following tasks: (a) Select the desired radio channel and reject other radio signals, (b) Amplify the desired radio signal and translate them back to baseband, and (c) Detect and decode the information with Low BER. In order to minimize cost and achieve high level of integration, receiver architecture should use least number of external filters and passive components. It should also consume least amount of power to minimize battery cost, size, and weight. One of the most stringent requirements for deep-space communication is the low-power operation. Our study identified that two candidate architectures listed in the following meet these requirements: (1) Low-IF receiver, (2) Sub-sampling receiver. The low-IF receiver uses minimum number of external components. Compared to Zero-IF (Direct conversion) architecture, it has less severe offset and flicker noise problems. The Sub-sampling receiver amplifies the RF signal and samples it using track-and-hold Subsampling mixer. These architectures provide low-power solution for the short- range communications missions on Mars. Accomplishments to date include: (1) System-level design and simulation of a Double-Differential PSK receiver, (2) Implementation of Honeywell SOI CMOS process design kit (PDK) in Cadence design tools, (3) Design of test circuits to investigate relationships between layout techniques, geometry, and low-frequency noise in SOI CMOS, (4) Model development and verification of on-chip spiral inductors in SOI CMOS process, (5) Design/implementation of low-power low-noise amplifier (LNA) and mixer for low-IF receiver, and (6) Design/implementation of high-gain LNA for sub-sampling receiver. Our initial results show that substantial improvement in power consumption is achieved using SOI CMOS as compared to standard CMOS process. Potential advantages of SOI CMOS for deep-space communication electronics include: (1) Radiation hardness, (2) Low-power operation, and (3) System-on-Chip (SOC) solutions.

Fujikawa, Gene↗
DOE OSTI2020

Monitoring the Possible Effects of Corroding Fuel Cladding in Nuclear Ponds: a Modified ROV for Storage Ponds - 20144

Nuclear ponds are used worldwide to store radioactive nuclear material, often for extended periods of time. The pond water conditions need to be closely monitored so that mitigation strategies can be put in place should any leakage of radioactive material occur. This paper describes the selection, development and on-site deployment of a bespoke system whose design was based on of-the-shelf components. The system is able to identify, locate and differentiate the leakage of radioactive material stored in nuclear ponds from the background activity of the water in the pond. The selection of an appropriate detector and monitoring geometry was supported by MCNP modelling which allowed a high degree of confidence in the proposed design. The development of the system allowed a real-time indication of the canister integrity, without incurring the increased operator doses from manual sampling operations that were previously used on site. The benefits of the developed system include the fact that the system is lightweight, easily deployable, retrievable, decontaminable and results in negligible waste products. Moreover, it is able to analyse the samples in-situ, in this way reducing the time needed to manually collect and send the samples to a laboratory for lengthy analysis. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

3D printing of large, complex metallic glass structures

Metallic glasses (MGs) or amorphous alloys although have superior mechanical properties their products are limited to simple geometries such as foils/plates or rods with thin section-thickness due to the requirement of high cooling rates. In this study, 3D, large dimensions of amorphous structures with complex geometry are manufactured by our newly developed Laser Foil Printing (LFP) technology. Zr-based (LM105, provided by Liquidmetal Technologies, Inc.) amorphous foils of 100 μm thickness are used as feedstock, and they are laser welded, layer-by-layer, to become 3D amorphous structures. Test results by X-Ray diffraction (XRD), differential scanning calorimetry (DSC), and micro-hardness confirm that the printed structures at selected process parameters achieve the same or better degree of amorphization as the raw foils. A mathematical model was developed to calculate the heating and cooling rates during structure manufacturing which helps the selection of process parameters. Here this study expands MG products to 3D arbitrary geometries with large dimensions due to the inherited advantages of the LFP technology which would open many potential applications.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A set of parallel, implicit methods for a reconstructed discontinuous Galerkin method for compressible flows on 3D hybrid grids

A set of implicit methods are proposed for a third-order hierarchical WENO reconstructed discontinuous Galerkin method for compressible flows on 3D hybrid grids. An attractive feature in these methods are the application of the Jacobian matrix based on the P1 element approximation, resulting in a huge reduction of memory requirement compared with DG (P2). Also, three approaches -- analytical derivation, divided differencing, and automatic differentiation (AD) are presented to construct the Jacobian matrix respectively, where the AD approach shows the best robustness. A variety of compressible flow problems are computed to demonstrate the fast convergence property of the implemented flow solver. Furthermore, an SPMD (single program, multiple data) programming paradigm based on MPI is proposed to achieve parallelism. The numerical results on complex geometries indicate that this low-storage implicit method can provide a viable and attractive DG solution for complicated flows of practical importance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Rod-like plasmonic nanoparticles as optical building blocks: how differences in particle shape and structural geometry influence optical signal

Gold nanoparticles, particularly those with an anisotropic shape, have become a popular optical probe for experiments involving work on the nanoscale. However, to carry out such delicate and intricate experiments, it is first necessary to understand the detailed behavior of individual nanoparticles. In this series of experiments, optical and electron microscopy were utilized for the characterization of individual nanoparticles and small assemblies of nanoparticles. In the first experiment, gold nanorods were investigated. Single, isolated nanorods exhibit two maxima of localized surface plasmon resonance (LSPR), which are associated with the two nanorod axes. Upon the physical rotation of a nanorod at one of its LSPR wavelengths under polarized illumination, the optical behavior varies in a sinusoidal fashion. A dimer of nanorods exhibits optical behavior quite similar to a nanorod, except the LSPR maxima are shifted and broader. Under differential interference contrast (DIC) microscopy, a pair of nanorods separated by a distance below the diffraction limit can be distinguished from a single nanorod due to its optical behavior upon rotation. Dark field microscopy is unable to distinguish the two geometries. For the second set of experiments, the optical behavior of single gold nanorods at non-plasmonic wavelengths was investigated. The same nanorod was rotated with respect to a polarized light source under DIC, dark field, and polarized light microscopy. DIC microscopy was found to produce diffraction pattern peaks at non-plasmonic wavelengths, which could be altered by adjusting the setting of the polarizer. In the third set of experiments, the optical behavior of a single gold dumbbell and several simple dumbbell geometries were investigated with microscopy and simulations. The single dumbbell displayed behavior quite similar to that of a nanorod, but dumbbells exhibit a shift in both LSPR wavebands. Moreover, the shape of dumbbell particles allows them to interlock with one another quite easily. The dimers that form as a result display optical behavior that differs from what has been previously reported about nanorod dimers. Simulated surface charge density patterns reveal that hybridization of LSPR modes occurs readily along the lobes of individual dumbbells in some situations. A pentamer of dumbbells also displays hybridization of modes, and “hot spots” are observed at junctions between pairings of dumbbells. In the final set of experiments, the assembly behavior of nanoparticles in solution was observed in real time. In general, large assemblies of nanoparticles display backbone-like rigidity, but an interesting variety of movements is permitted within the larger structures.

36 MATERIALS SCIENCE↗
DOE OSTI2014

3D mapping of water in oolithic limestone at atmospheric and vacuum saturation using X-ray micro-CT differential imaging

Determining the distribution of fluids in porous sedimentary rocks is of great importance in many geological fields. However, this is not straightforward, especially in the case of complex sedimentary rocks like limestone, where a multidisciplinary approach is often needed to capture its broad, multimodal pore size distribution and complex pore geometries. This paper focuses on the porosity and fluid distribution in two varieties of Massangis limestone, a widely used natural building stone from the southeast part of the Paris basin (France). The Massangis limestone shows locally varying post-depositional alterations, resulting in different types of pore networks and very different water distributions within the limestone. Traditional techniques for characterizing the porosity and pore size distribution are compared with state-of-the-art neutron radiography and X-ray computed microtomography to visualize the distribution of water inside the limestone at different imbibition conditions. X-ray computed microtomography images have the great advantage to non-destructively visualize and analyze the pore space inside of a rock, but are often limited to the larger macropores in the rock due to resolution limitations. In this paper, differential imaging is successfully applied to the X-ray computed microtomography images to obtain sub-resolution information about fluid occupancy and to map the fluid distribution in three dimensions inside the scanned limestone samples. The detailed study of the pore space with differential imaging allows understanding the difference in the water uptake behavior of the limestone, a primary factor that affects the weathering of the rock. - Highlights: • The water distribution in a limestone was visualized in 3D with micro-CT. • Differential imaging allowed to map both macro and microporous zones in the rock. • The 3D study of the pore space clarified the difference in water uptake behavior. • Trapped air is visualized in the moldic macropores at atmospheric saturation.

36 MATERIALS SCIENCE↗
DOE OSTI2016

DOUBLE POWER LAWS IN THE EVENT-INTEGRATED SOLAR ENERGETIC PARTICLE SPECTRUM

A double power law or a power law with exponential rollover at a few to tens of MeV nucleon{sup −1} of the event-integrated differential spectra has been reported in many solar energetic particle (SEP) events. The rollover energies per nucleon of different elements correlate with a particle's charge-to-mass ratio (Q/A). The probable causes are suggested as residing in shock finite lifetimes, shock finite sizes, shock geometry, and an adiabatic cooling effect. In this work, we conduct a numerical simulation to investigate a particle's transport process in the inner heliosphere. We solve the focused transport equation using a time-backward Markov stochastic approach. The convection, magnetic focusing, adiabatic cooling effect, and pitch-angle scattering are included. The effects that the interplanetary turbulence imposes on the shape of the resulting SEP spectra are examined. By assuming a pure power-law differential spectrum at the Sun, a perfect double-power-law feature with a break energy ranging from 10 to 120 MeV nucleon{sup −1} is obtained at 1 au. We found that the double power law of the differential energy spectrum is a robust result of SEP interplanetary propagation. It works for many assumptions of interplanetary turbulence spectra that give various forms of momentum dependence of a particle's mean free path. The different spectral shapes in low-energy and high-energy ends are not just a transition from the convection-dominated propagation to diffusion-dominated propagation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Al0 0.3 Ga 0.7 N PN diode with breakdown voltage >1600 V

Demonstration of Al0 0.3 Ga 0.7 N PN diodes grown with breakdown voltages in excess of 1600 V is reported. The total epilayer thickness is 9.1 μm and was grown by metal-organic vapour-phase epitaxy on 1.3-mm-thick sapphire in order to achieve crack-free structures. A junction termination edge structure was employed to control the lateral electric fields. A current density of 3.5 kA/cm 2 was achieved under DC forward bias and a reverse leakage current <3 nA was measured for voltages <1200 V. The differential on-resistance of 16 mΩ cm 2 is limited by the lateral conductivity of the n-type contact layer required by the front-surface contact geometry of the device. An effective critical electric field of 5.9 MV/cm was determined from the epilayer properties and the reverse current–voltage characteristics. To our knowledge, this is the first aluminium gallium nitride (AlGaN)-based PN diode exhibiting a breakdown voltage in excess of 1 kV. Finally, we note that a Baliga figure of merit (V br 2 /R spec,on ) of 150 MW/cm 2 found is the highest reported for an AlGaN PN diode and illustrates the potential of larger-bandgap AlGaN alloys for high-voltage devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗
DOE OSTI2017

Dose- and Ion-Dependent Effects in the Oxidative Stress Response to Space-Like Radiation Exposure in the Skeletal System

Space radiation may pose a risk to skeletal health during subsequent aging. Irradiation acutely stimulates bone remodeling in mice, although the long-term influence of space radiation on bone-forming potential (osteoblastogenesis) and possible adaptive mechanisms are not well understood. We hypothesized that ionizing radiation impairs osteoblastogenesis in an ion-type specific manner, with low doses capable of modulating expression of redox-related genes. 16-weeks old, male, C57BL6/J mice were exposed to low linear-energy-transfer (LET) protons (150 MeV/n) or high-LET 56Fe ions (600 MeV/n) using either low (5 or 10 cGy) or high (50 or 200 cGy) doses at NASA’s Space Radiation Lab. Five weeks or one year after irradiation, tissues were harvested and analyzed by microcomputed tomography for cancellous microarchitecture and cortical geometry. Marrow-derived, adherent cells were grown under osteoblastogenic culture conditions. Cell lysates were analyzed by RT-PCR during the proliferative or mineralizing phase of growth, and differentiation was analyzed by imaging mineralized nodules. As expected, a high dose (200 cGy), but not lower doses, of either 56Fe or protons caused a loss of cancellous bone volume/total volume. Marrow cells produced mineralized nodules ex vivo regardless of radiation type or dose; 56Fe (200 cGy) inhibited osteoblastogenesis by more than 90% (5 weeks and 1 year post-IR). After 5 weeks, irradiation (protons or 56Fe) caused few changes in gene expression levels during osteoblastogenesis, although a high dose 56Fe (200 cGy) increased Catalase and Gadd45. The addition of exogenous superoxide dismutase (SOD) protected marrow-derived osteoprogenitors from the damaging effects of exposure to low-LET (137Cs γ) when irradiated in vitro, but had limited protective effects on high-LET 56Fe-exposed cells. In sum, either protons or 56Fe at a relatively high dose (200 cGy) caused persistent bone loss, whereas only high-LET 56Fe increased redox-related gene expression, albeit to a limited extent, and inhibited osteoblastogenesis. Doses below 50 cGy did not elicit widespread responses in any parameter measured. We conclude that high-LET irradiation at 200 cGy impaired osteoblastogenesis and regulated steady-state gene expression of select redox-related genes during osteoblastogenesis, which may contribute to persistent bone loss.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1988

Rotating hot-wire investigation of the vortex responsible for blade-vortex interaction noise

This distribution of the circumferential velocity of the vortex responsible for blade-vortex interaction noise was measured using a rotating hot-wire rake synchronously meshed with a model helicopter rotor at the blade passage frequency. Simultaneous far-field acoustic data and blade differential pressure measurements were obtained. Results show that the shape of the measured far-field acoustic blade-vortex interaction signature depends on the blade-vortex interaction geometry. The experimental results are compared with the Widnall-Wolf model for blade-vortex interaction noise.

Fontana, Richard Remo↗
NASA NTRS2007

Summary of the Third AIAA CFD Drag Prediction Workshop

The workshop focused on the prediction of both absolute and differential drag levels for wing-body and wing-al;one configurations of that are representative of transonic transport aircraft. The baseline DLR-F6 wing-body geometry, previously utilized in DPW-II, is also augmented with a side-body fairing to help reduce the complexity of the flow physics in the wing-body juncture region. In addition, two new wing-alone geometries have been developed for the DPW-II. Numerical calculations are performed using industry-relevant test cases that include lift-specific and fixed-alpha flight conditions, as well as full drag polars. Drag, lift, and pitching moment predictions from previous Reynolds-Averaged Navier-Stokes computational fluid Dynamics Methods are presented, focused on fully-turbulent flows. Solutions are performed on structured, unstructured, and hybrid grid systems. The structured grid sets include point-matched multi-block meshes and over-set grid systems. The unstructured and hybrid grid sets are comprised of tetrahedral, pyramid, and prismatic elements. Effort was made to provide a high-quality and parametrically consistent family of grids for each grid type about each configuration under study. The wing-body families are comprised of a coarse, medium, and fine grid, while the wing-alone families also include an extra-fine mesh. These mesh sequences are utilized to help determine how the provided flow solutions fair with respect to asymptotic grid convergence, and are used to estimate an absolute drag of each configuration.

Vassberg, John C.↗
NASA NTRS2015

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.↗