Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling 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 523 records · Page 29

Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry

Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the number of design variables is increased.

Nelson, Andrea↗

Compact Augmented Spark Igniters for Liquid Rocket Engines

The objective of this dissertation is to characterize augmented spark impinging (torch) igniters. Torch igniters inject propellants into a small prechamber where a spark igniter imparts the activation energy to initiate combustion. The flame stabilizes within the prechamber and exits into the main combustion chamber where it ignites the main flow directly or by mixing with a coaxial flow to create a secondary flame. The underlying phenomena that promote reliable ignition are multifaceted and required a broad scope of investigations. Initial efforts focused on spark discharges across an annular geometry. Gaseous propellent projects the electrical arc toward the prechamber, where the arc and the thermal energy is imparted to the surrounding fluid then meets with a combustible mixture. Discharges are parametrically examined against pressure, spark-gap size, and exciter types to emulate transient engine startup. The energy imparted to the flow, approximate velocity, and distance traveled by the exhaust plume are determined. Spark quenching (no energy deposition) and its effects on perceived versus actual measurements are shown. Results outline the prechamber target zone where a combustible mixture is necessary for ignition. Next, three-dimensional, time-accurate, and nonreacting computational fluid dynamics simulations assess effects of geometric and mass flow differences on the size, location, and composition of a combustible mixture. Fluid dynamic phenomena are elucidated. Correlations between nondimensional variables are found. Results are coupled with a conservation-law analysis to produce novel objective functions. Lastly, a full-scale, modular torch igniter is tested. Injector sizes, mixture ratios, and injector momentum ratios are systematically varied to map ignition probability and ignition delay. Results show the importance of local mixture ratios, provide injector momentum ratios guidelines, and demonstrate reliable ignition of core mixture ratios typically outside flammability limits. The deliverable of this project is not a torch igniter or derivative thereof. Instead, this dissertation produces novel objective-functions for igniter design. These equations are functions of variables that may be constrained by engine-level requirements, physics, and experimental results from this work. They provide a foundation for future igniter design and will minimize iterative processes from the design cycle for next-generation ignition hardware.

Torch Igniter↗

The coupling between flame surface dynamics and species mass conservation in premixed turbulent combustion

Current flamelot models based on a description of the flame surface dynamics require the closure of two inter-related equations: a transport equation for the mean reaction progress variable, (tilde)c, and a transport equation for the flame surface density, Sigma. The coupling between these two equations is investigated using direct numerical simulations (DNS) with emphasis on the correlation between the turbulent fluxes of (tilde)c, bar(pu''c''), and Sigma, (u'')(sub S)Sigma. Two different DNS databases are used in the present work: a database developed at CTR by A. Trouve and a database developed by C. J. Rutland using a different code. Both databases correspond to statistically one-dimensional premixed flames in isotropic turbulent flow. The run parameters, however, are significantly different, and the two databases correspond to different combustion regimes. It is found that in all simulated flames, the correlation between bar(pu''c'') and (u'')(sub S)Sigma is always strong. The sign, however, of the turbulent flux of (tilde)c or Sigma with respect to the mean gradients, delta(tilde)c/delta(x) or delta(Sigma)/delta(x), is case-dependent. The CTR database is found to exhibit gradient turbulent transport of (tilde)c and Sigma, whereas the Rutland DNS features counter-gradient diffusion. The two databases are analyzed and compared using various tools (a local analysis of the flow field near the flame, a classical analysis of the conservation equation for (tilde)(u''c''), and a thin flame theoretical analysis). A mechanism is then proposed to explain the discrepancies between the two databases and a preliminary simple criterion is derived to predict the occurrence of gradient/counter-gradient turbulent diffusion.

Trouve, A.↗

Control of Combustion-Instabilities Through Various Passive Devices

It is well known that under some operating conditions, rocket engines (using solid or liquid fuels) exhibit unstable modes of operation that can lead to engine malfunction and shutdown. The sources of these instabilities are diverse and are dependent on fuel, chamber geometry and various upstream sources such as pumps, valves and injection mechanism. It is believed that combustion-acoustic instabilities occur when the acoustic energy increase due to the unsteady heat release of the flame is greater than the losses of acoustic energy from the system [1, 2]. Giammar and Putnam [3] performed a comprehensive study of noise generated by gasfired industrial burners and made several key observations; flow noise was sometimes more intense than combustion roar, which tended to have a characteristic frequency spectrum. Turbulence was amplified by the flame. The noise power varied directly with combustion intensity and also with the product of pressure drop and heat release rate. Karchmer [4] correlated the noise emitted from a turbofan jet engine with that in the combustion chamber. This is important, since it quantified how much of the noise from an engine originates in the combustor. A physical interpretation of the interchange of energy between sound waves and unsteady heat release rates was given by Rayleigh [5] for inviscid, linear perturbations. Bloxidge et al [6] extended Rayleigh s criterion to describe the interaction of unsteady combustion with one-dimensional acoustic waves in a duct. Solutions to the mass, momentum and energy conservation equations in the pre- and post-flame zones were matched by making several assumptions about the combustion process. They concluded that changes in boundary conditions affect the energy balance of acoustic waves in the combustor. Abouseif et al [7] also solved the one-dimensional flow equations, but they used a onestep reaction to evaluate the unsteady heat release rate by relating it to temperature and velocity perturbations. Their analysis showed that oscillations arise from coupling between entropy waves produced at the flame and pressure waves originating from the nozzle. Yang and Culick [8] assumed a thin flame sheet, which is distorted by velocity and pressure oscillations. Conservation equations were expressed in integral form and solutions for the acoustic wave equations and complex frequencies were obtained. The imaginary part of the frequency indicated stability regions of the flame. Activation energy asymptotics together with a one-step reaction were used by McIntosh [9] to study the effects of acoustic forcing and feedback on unsteady, one-dimensional flames. He found that the flame stability was altered by the upstream acoustic feedback. Shyy et al [10] used a high-accuracy TVD scheme to simulate unsteady, one-dimensional longitudinal, combustion instabilities. However, numerical diffusion was not completely eliminated. Recently, Prasad [11] investigated numerically the interactions of pressure perturbations with premixed flames. He used complex chemistry to study responses of pressure perturbations in one-dimensional combustors. His results indicated that reflected and transmitted waves differed significantly from incident waves.

Frendi, Kader↗

CAPRI: A Geometric Foundation for Computational Analysis and Design

CAPRI is a software building tool-kit that refers to two ideas; (1) A simplified, object-oriented, hierarchical view of a solid part integrating both geometry and topology definitions, and (2) programming access to this part or assembly and any attached data. A complete definition of the geometry and application programming interface can be found in the document CAPRI: Computational Analysis PRogramming Interface appended to this report. In summary the interface is subdivided into the following functional components: 1. Utility routines -- These routines include the initialization of CAPRI, loading CAD parts and querying the operational status as well as closing the system down. 2. Geometry data-base queries -- This group of functions allow all top level applications to figure out and get detailed information on any geometric component in the Volume definition. 3. Point queries -- These calls allow grid generators, or solvers doing node adaptation, to snap points directly onto geometric entities. 4. Calculated or geometrically derived queries -- These entry points calculate data from the geometry to aid in grid generation. 5. Boundary data routines -- This part of CAPRI allows general data to be attached to Boundaries so that the boundary conditions can be specified and stored within CAPRI s data-base. 6. Tag based routines -- This part of the API allows the specification of properties associated with either the Volume (material properties) or Boundary (surface properties) entities. 7. Geometry based interpolation routines -- This part of the API facilitates Multi-disciplinary coupling and allows zooming through Boundary Attachments. 8. Geometric creation and manipulation -- These calls facilitate constructing simple solid entities and perform the Boolean solid operations. Geometry constructed in this manner has the advantage that if the data is kept consistent with the CAD package, therefore a new design can be incorporated directly and is manufacturable. 9. Master Model access This addition to the API allows for the querying of the parameters and dimensions of the model. The feature tree is also exposed so it is easy to see where the parameters are applied. Calls exist to allow for the modification of the parameters and the suppression/unsuppression of nodes in the tree. Part regeneration is performed by a single API call and a new part becomes available within CAPRI (if the regeneration was successful). This is described in a separate document. Components 1-7 are considered the CAPRI base level reader.

Haimes, Robert↗

The High Mosaicity Illusion: Revealing the True Physical Characteristics of Macromolecular Crystals

An experimental system and software have been developed for simultaneously measuring the diffraction resolution and mosaic spread of macromolecular crystals. Hundreds of reflection profiles over a wide resolution range were rapidly measured by using a charge coupled device (CCD) area detector in combination with superfine phi slicing data collection. The contributions of the X-ray beam to the reflection widths were minimized by using a highly-parallel, highly-monochromatic synchrotron source. These contributions and Lorentz effects were evaluated and deconvoluted from the recorded data. Data collection and processing is described. From one degree of superfine phi slice data collected on a crystal of manganese superoxide dismutase the mosaicity of 261 reflections were measured. The average mosaicity was 0.0101 degrees (0.0035) at the full-width-at-half-maximum (FWHM) and ranged from 0.0011 degrees to 0.0188 degrees. Each reflection profile was individually fit with two gaussian profiles with the first gaussian contributing 55% and the second contributing 35% of the reflection. On average, the mosaicity of the first gaussian was 0.0054 degrees (0.0015) and the second was 0.0061 degrees (0.0023). The mosaicity of the crystal was anisotropic with fh, f k, and fl values of 0.0068 degrees, 0.0140 degrees and 0.0046 degrees, respectively at the FWHM. The anisotropic mosaicity analysis indicates that the crystal is the most perfect in the I direction which corresponds to the favored growth direction of the crystal.

Bellamy, Henry↗

Comparison of Different Global Information Sources Used in Surface Radiative Flux Calculation: Radiative Properties of the Surface

Direct estimates of surface radiative fluxes that resolve regional and weather-scale variabilty over the whole globe with reasonable accuracy have only become possible with the advent of extensive global, mostly satellite, datasets within the past couple of decades. The accuracy of these fluxes, estimated to be about 10-15 W per square meter is largely limited by the accuracy of the input datasets. The leading uncertainties in the surface fluxes are no longer predominantly induced by clouds but are now as much associated with uncertainties in the surface and near-surface atmospheric properties. This study presents a fuller, more quantitative evaluation of the uncertainties for the surface albedo and emissivity and surface skin temperatures by comparing the main available global datasets from the Moderate-Resolution Imaging Spectroradiometer product, the NASA Global Energy and Water Cycle Experiment Surface Radiation Budget project, the European Centre for Medium-Range Weather Forecasts, the National Aeronautics and Space Administration, the National Centers for Environmental Prediction, the International Satellite Cloud Climatology Project (ISCCP), the Laboratoire de Meteorologie Dynamique, NOAA/NASA Pathfinder Advanced Very High Resolution Radiometer project, NOAA Optimum Interpolation Sea Surface Temperature Analysis and the Tropical Rainfall Measuring Mission (TRMM) Microwave Image project. The datasets are, in practice, treated as an ensemble of realizations of the actual climate such that their differences represent an estimate of the uncertainty in their measurements because we do not possess global truth datasets for these quantities. The results are globally representative and may be taken as a generalization of our previous ISCCP-based uncertainty estimates for the input datasets. Surface properties have the primary role in determining the surface upward shortwave (SW) and longwave (LW) flux. From this study, the following conclusions are obtained. Although land surface albedos in the near near-infrared remain poorly constrained (highly uncertain), they do not cause too much error in total surface SW fluxes; the more subtle regional and seasonal variations associated with vegetation and snow are still on doubt. The uncertainty of the broadband black-sky SW albedo for land surface from this study is about 7%, which can easily induce 5-10 W per square meter uncertainty in (upwelling) surface SW flux estimates. Even though available surface (broadband) LW emissivity datasets differ significantly (3%-5% uncertainty), this disagreement is confined to wavelengths greater than 20 micrometers so that there is little practical effect (1-3 W per square meters) on the surface upwelling LW fluxes. The surface skin temperature is one of two leading factors that cause problems with surface LW fluxes. Even though the differences among the various datasets are generally only 2-4 K, this can easily cause 10-15 W per square meter uncertainty in calculated surface (upwelling) LW fluxes. Significant improvements could be obtained for surface LW flux calculations by improving the retrievals of (in order of decreasing importance): (1) surface skin temperature, (2) surface air and near-surface-layer temperature, (3) column precipitable water amount and (4) broadband emissivity. And for surface SW fluxes, improvements could be obtained (excluding improved cloud treatment) by improving the retrievals of (1) aerosols (from our sensitivity studies but not discussed in this work), and (2) surface (black-sky) albedo, of which, NIR part of the spectrum has much larger uncertainty.

Zhang, Yuanchong↗

Polarization of Directly Imaged Young Giant Planets as a Probe of Mass, Rotation, and Clouds

Young, hot gas giant planets at large separations from their primaries have been directly imaged around several nearby stars. More such planets will likely be detected by ongoing and new imaging surveys with instruments such as the Gemini Planet Imager (GPI). Efforts continue to model the spectra of these planets in order to constrain their masses, effective temperatures, composition, and cloud structure. One potential tool for analyzing these objects, which has received relatively less attention, is polarization. Linear polarization of gas giant exoplanets can arise from the combined influences of light scattering by atmospheric dust and a rotationally distorted shape. The oblateness of gas giant planet increases of course with rotation rate and for fixed rotation also rises with decreasing gravity. Thus young, lower mass gas giant planets with youthful inflated radii could easily have oblateness greater than that of Saturn s 10%. We find that polarizations of over 1% may easily be produced in the near-infrared in such cases. This magnitude of polarization may be measurable by GPI and other instruments. Thus if detected, polarization of a young Jupiter places constraints on the combination of its gravity, rotation rate, and degree of cloudiness. We will present results of our multiple scattering analysis coupled with a self-consistent dusty atmospheric models to demonstrate the range of polarizations that might be expected from resolved exoplanets and the range of parameter space that such observations may inform.

Marley, Mark Scott↗

Direct Multipoint Observations Capturing the Reformation of a Supercritical Fast Magnetosonic Shock

Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth’s bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the “foot” region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock’s rest frame. That conversion allows for a unique study of the relative spatial scales of the shock’s various features, including the shock’s growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Evidence is also presented of nonlinear wave decay from electromagnetic whistler-mode “lion roars” to electrostatic solitary waves in the downstream plasma regime. Results highlight the critical cross-scale coupling between electron-kinetic-and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks.

D L Turner↗

Discreditation of Bobdownsite and the Establishment of Criteria for the Identification of Minerals with Essential Monofluorophosphate (PO3F2-)

Bobdownsite, IMA number 2008-037, was approved as a new mineral by the Commission on New Minerals, Nomenclature and Classification (CNMNC) as the fluorine endmember of the mineral whitlockite. The type locality of bobdownsite is in Big Fish River, Yukon Canada, and bobdownsite was reported to be the first naturally occurring mineral with essential monofluorophosphate (PO3F2-). The type specimen of bobdownsite has been reinvestigated by electron probe microanalysis (EPMA), and our data indicate that fluorine abundances are below detection in the mineral. In addition, we conducted detailed analysis of bobdownsite from the type locality by gas chromatography isotope ratio mass spectrometry, Raman spectroscopy, EPMA, and NMR spectroscopy. These data were compared with previously published data on synthetic monofluorophosphate salts. Collectively, these data indicate that bobdownsite is indistinguishable from whitlockite with a composition along the whitlockite-merrillite solid solution. Bobdownsite is therefore discredited as a valid mineral species. An additional mineral, krásnoite, has been purported to have monofluorophosphate components in its structure, but reexamination of those data indicate that F- in krásnoite forms bonds with Al, similar to OH- bonded to Al in perhamite. Consequently, krásnoite also lacks monofluorophosphate groups, and there are currently no valid mineral species with monofluorophosphate in their structure. We recommend that any future reports of new minerals that contain essential monofluorophosphate anions be vetted by abundance measurements of fluorine, vibrational spectroscopy (both Raman and FTIR), and where paramagnetic components are permissibly low, NMR spectroscopy. Furthermore, we emphasize the importance of using synthetic compounds containing monofluorophosphate anions as a point of comparison in the identification of minerals with essential monofluorophosphate. Structural data that yield satisfactory P-F bond lengths determined by X-ray crystallography, coupled with direct chemical analyses of fluorine in a material do not constitute sufficient evidence alone to identify a new mineral with essential monofluorophosphate anions.

McCubbin, Francis M.↗

NASA's Modern Era Retrospective-Analysis for Research and Applications (MERRA): Early Results and Future Directions

This talk will review the status and progress of the NASA/Global Modeling and Assimilation Office (GMAO) atmospheric global reanalysis project called the Modern Era Retrospective-Analysis for Research and Applications (MERRA). An overview of NASA's emerging capabilities for assimilating a variety of other Earth Science observations of the land, ocean, and atmospheric constituents will also be presented. MERRA supports NASA Earth science by synthesizing the current suite of research satellite observations in a climate data context (covering the period 1979-present), and by providing the science and applications communities with of a broad range of weather and climate data with an emphasis on improved estimates of the hydrological cycle. MERRA is based on a major new version of the Goddard Earth Observing System Data Assimilation System (GEOS-5), that includes the Earth System Modeling Framework (ESMF)-based GEOS-5 atmospheric general circulation model and the new NOAA National Centers for Environmental Prediction (NCEP) unified grid-point statistical interpolation (GST) analysis scheme developed as a collaborative effort between NCEP and the GMAO. In addition to MERRA, the GMAO is developing new capabilities in aerosol and constituent assimilation, ocean, ocean biology, and land surface assimilation. This includes the development of an assimilation capability for tropospheric air quality monitoring and prediction, the development of a carbon-cycle modeling and assimilation system, and an ocean data assimilation system for use in coupled short-term climate forecasting.

Schubert, Siegfried↗

Whole earth telescope observations of the DBV white dwarf GD 358

We report on the analysis of 154 hours of early continuous high-speed photometry on the pulsating DB white dwarf (DBV) GD 358, obtained during the Whole Earth Telescope (WET) run of 1990 May. The power spectrum of the light curve is dominated by power in the range from 1000 to 2400 microHz with more than 180 significant peaks in the total spectrum. We identify all of the triplet frequencies as degree l = 1, and from the details of their spacings we derive the total stellar mass as 0.61 + or - 0.03 solar mass, the mass of the outer helium envelope as 2.0 + or - 1.0 x 10(exp -6) M(sub *), the absolute luminosity as 0.050 + or - 0.012 solar luminosity and the distance as 42 + or - 3 pc. We find strong evidence for differential rotation in the radial direction -- the outer envelope is rotating at least 1.8 times faster than the core -- and we detect the presence of a weak magnetic field with a strength of 1300 + or - 300 G. We also find a significant power at the sums and differences of the dominant frequencies, indicating nonlinear processes are significant, but they have a richness and complexity that rules out resonant mode coupling as a major cause.

Winget, D. E.↗

Space infrared telescope pointing control system. Infrared telescope tracking in the presence of target motion

The use of charge-coupled-devices, or CCD's, has been documented by a number of sources as an effective means of providing a measurement of spacecraft attitude with respect to the stars. A method exists of defocussing and interpolation of the resulting shape of a star image over a small subsection of a large CCD array. This yields an increase in the accuracy of the device by better than an order of magnitude over the case when the star image is focussed upon a single CCD pixel. This research examines the effect that image motion has upon the overall precision of this star sensor when applied to an orbiting infrared observatory. While CCD's collect energy within the visible spectrum of light, the targets of scientific interest may well have no appreciable visible emissions. Image motion has the effect of smearing the image of the star in the direction of motion during a particular sampling interval. The presence of image motion is incorporated into a Kalman filter for the system, and it is shown that the addition of a gyro command term is adequate to compensate for the effect of image motion in the measurement. The updated gyro model is included in this analysis, but has natural frequencies faster than the projected star tracker sample rate for dim stars. The system state equations are reduced by modelling gyro drift as a white noise process. There exists a tradeoff in selected star tracker sample time between the CCD, which has improved noise characteristics as sample time increases, and the gyro, which will potentially drift further between long attitude updates. A sample time which minimizes pointing estimation error exists for the random drift gyro model as well as for a random walk gyro model.

Powell, J. D.↗

Azimuthal Size Scales of Solar Wind Periodic Density Structures

Periodic Density Structures (PDSs) are quasi-periodic variations of solar wind density ranging from a few minutes to a few hours. PDSs advect with the solar wind and have radial length scales (Lx) of tens to several thousand Mm, thus belonging to the class of “mesoscale structures”. Current interplanetary multi-spacecraft observations are not at spatial separations capable of directly measuring the 3D size scale of PDSs or other mesoscale structures. Instead, previous investigations estimated characteristic spatial scales in solar wind parameters using cross-correlation and/or coherence analysis applied to multi-spacecraft observations. For solar wind density and interplanetary magnetic field (IMF) intensity, the reported size scales perpendicular to the Sun-Earth line (Ly) ranged between ≈30 and ≈200 Earth Radii (RE). Here, we implemented a similar approach on the same parameters but focused on high density slow solar wind intervals with PDSs observed by the Wind and ARTEMIS-P1 spacecraft. Additionally, this is the first statistical study of the IMF intensity periodicities in relation to PDSs. We identified intervals in which the two spacecraft observed the same periodicity obtaining two PDS groups based on their radial length scale: Lx1 ≈86RE and Lx2 ≈35RE. Then, we classified the events based on the periodic variations’ coherence level. Reproducing the results with simulations of the PDSs’ transit, we inferred the Ly order of magnitudes for the two PDS groups: Ly1 ≈340RE and Ly2 ≈187RE. Knowing the PDSs’ size scales is fundamental to constrain models aimed at reproducing these structures and is critical to better understand the PDS-magnetosphere coupling.

Interplanetary physics↗

Changes in Thermospheric O/N2 Derived from UVI Auroral Images

A rigorous test of our understanding of the coupled ionosphere-thermosphere and its response to geomagnetic storms is the ability to reproduce observed storm effects as seen in the ionosphere and neutral atmosphere. The concept of compositional change is central to studies of thermosphere response to storm conditions. In particular, information about compositional change within the highly dynamic auroral region is limited. The Ultraviolet Imager (UVI) is designed to view the full auroral region using five filters to isolate emissions from atomic oxygen (1304 and 1356) and N2 LBH. This spectral resolution allows auroral energy characteristics to be derived by two separate methods from examining ratios of observed intensities (OI 1356/LBHL or LBHS/LBHL). The LBHS:LBHL ratio is typically used as the mean energy diagnostic since the OI 1356 emission is dependent on changes in the atomic oxygen density, and these changes relative to N2 can be large. However, once the mean energy has been specified by the LBH ratio, this variability in OI 1356 emission can be exploited as a direct diagnostic of total atomic oxygen column density. This opens the potential of using UVI images to monitor the temporal and spatial response of thermospheric O to high latitude forcing within the auroral regions. Initial results of this type of analysis will be presented along with discussion of its limitations and capabilities.

Germany, G. A.↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗

A Modern Load Relief Guidance Scheme for Space Launch Vehicles

Launch vehicle load relief algorithms are concerned with realizing a reduction of transient bending moments near maximum dynamic pressure. Traditional approaches to load relief typically use inner-loop acceleration feedback to reduce the wind-induced angle of attack. When implemented in the inner loop, load relief bandwidth is necessarily limited by the achievable stability margins, and when acceleration feedback is employed, by the uncertainty associated with structural modes that couple with the body-mounted accelerometer. The structure of inner loop load relief increases the dimensionality of the flight control gain and filter optimization problem. Most importantly, classical load relief laws do not take advantage of high-rate and high-accuracy GPS-aided inertial velocity data that is readily available from modern strap down IMUs. In this paper, a novel load relief guidance scheme is described that uses direct angle-of-attack feedback in a clever mechanization. The steering commands are determined by examining the wind-perturbed dynamics of a launch vehicle with respect to a gravity turn ascent trajectory. An angle of attack estimate is derived from GPS-aided inertial data and pre-launch range wind measurements, and it is shown that a reduction worst-case rigid-body loads can be realized without requiring air data. The algorithm also includes a high-rate navigation data preprocessing scheme that operates directly on the IMU delta-theta and delta-velocity measurements in order to produce a filtered acceleration estimate at the vehicle center of mass. The outer-loop guidance scheme simplifies the design process for the classical inner-loop autopilot. Algorithm performance is demonstrated using Monte Carlo analysis of a representative liquid booster in a production high fidelity launch vehicle simulation.

NESC↗