Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model order reduction”

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 883 records · Page 49

Sensitivity Analysis of Genome-Scale Metabolic Flux Prediction

TRIMER, Transcription Regulation Integrated with MEtabolic Regulation, is a genome-scale modeling pipeline targeting at metabolic engineering applications. Using TRIMER, regulated metabolic reactions can be effectively predicted by integrative modeling of metabolic reactions with a Transcription Factor (TF)-gene regulatory network (TRN), which is modeled via a Bayesian network (BN). In this paper, we focus on sensitivity analysis of metabolic flux prediction for uncertainty quantification of BN structures for TRN modeling in TRIMER. We propose a computational strategy to construct the uncertainty class of TRN models based on the inferred regulatory order uncertainty given transcriptomic expression data. With that, we analyze the prediction sensitivity of the TRIMER pipeline for the metabolite yields of interest. The obtained sensitivity analyses can guide Optimal Experimental Design (OED) to help acquire new data that can enhance TRN modeling and achieve specific metabolic engineering objectives, including metabolite yield alterations. Here we have performed small- and large-scale simulated experiments, demonstrating the effectiveness of our developed sensitivity analysis strategy for BN structure learning to quantify the edge importance in terms of metabolic flux prediction uncertainty reduction and its potential to effectively guide OED.

59 BASIC BIOLOGICAL SCIENCES↗

Enabling off-highway diesel engine downsizing and performance improvement using electrically assisted turbocharging

Internal combustion engine (ICE) downsizing through various turbocharging configurations is generally known by the powertrain design community as an effective means to reduce frictional losses, increase waste heat recovery, and improve fuel efficiency while increasing engine power density. However, often is the case that turbocharging strategies, including variable geometry turbochargers, and regulated two-stage turbochargers, incur the performance tradeoff between transient response and fuel economy (pumping losses) at high engine speeds. For off-highway vehicles having particularly transient and high-powered duty cycles, efforts to improve this tradeoff and increase operational flexibility have turned to evaluating various electrified air system architectures. In this study, a 4.5 L diesel-ICE configured with an electrically driven compressor (eBooster ® ) is placed in series with a conventional turbocharger and integrated into a 48 V mild-hybrid powertrain architecture. The objective of this powertrain configuration is to enable engine downsizing by 34%, replacing the current 6.8 L ICE platform with the hybridized 4.5 L ICE concept. The commercial 1-D simulation software GT-SUITE is used for powertrain system development and system optimization. Development and validation of the GT-SUITE model and air system controls is concurrently supported through experimental data collection. The simulation model development includes using machine learning methods for optimizing exhaust gas dilution, injection timing, and eBooster ® power to improve steady-state and transient brake-specific fuel consumption while minimizing criteria pollutant emissions. It was found that total specific fluid consumption over the standardized non-road transient engine duty cycle could be reduced by 18% over the current 6.8 L engine by using an optimized eBoosted 4.5 L engine. The hybridized 4.5 L engine concept concurrently showed sufficient transient response capability and nearly an order of magnitude reduction in duty cycle total soot production.

25 ENERGY STORAGE↗

Plate measurement techniques and reduction methods used by the West German satellite observers, and resulting consequences for the observation

The West German tracking stations are equipped with ballistic cameras. Plate measurement and plate reduction must therefore follow photogrammetric methods. Approximately 100 star positions and 200 satellite positions are measured on each plate. The mathematical model for spatial rotation of the bundle of rays is extended by including terms for distortion and internal orientation of the camera as well as by providing terms for refraction which are computed for the measured coordinates of the star positions on the plate. From the measuring accuracy of the plate coordinates it follows that the timing accuracy for the exposures has to be about one millisecond, in order to obtain a homogeneous system.

Deker, H.↗

Beam dynamics corrections in the measurement of the anomalous precession frequency at the Muon $g-2$ experiment at Fermilab

The Muon $g-2$ experiment at Fermilab (E989) aims to measure the anomalous magnetic moment of the muon with an accuracy of 140 ppb (parts per billions). This accuracy, obtained by adding in quadrature a statistical and a systematic contribution of comparable value (100 ppb), will allow to reduce the experimental uncertainty (from the previous E821 experiment at BNL \cite{bnl}) of a factor of 4, and represents one of the most precise tests of the Standard Model (SM) theory of elementary particles. The first result on the Run-1 dataset \cite{prl} was released on April 7, 2021, showing a very good agreement with the previous result from BNL experiment, with a slightly better uncertainty. The corresponding experimental average increases the significance of the discrepancy between the measured and Standard Model prediction of 4.2$\sigma$ \cite{white_paper}.\\ \noindent The measured quantity is the muon magnetic anomaly a$_\mu$=$\frac{g_\mu-2}{2}$ where g$_\mu$ is the gyromagnetic factor of the muon. Dirac's equation predicts g$_\mu$ = 2, while radiative corrections, dominated by the QED contribution due to an exchange of a virtual photon, causes a per-mille correction on this quantity. By including all the SM contributions, a$_\mu$ is known at 370 ppb. The E989 experiment measures $a_\mu$ injecting positive muons with momentum of 3.1 GeV/$c$ polarized longitudinally in bunches (called $fills$) with an average rate of 12 Hz, in a storage ring of 14 meters diameter. Due to the parity violation in the weak muon decay, high energy positrons produced are emitted preferably in the muon's spin direction. By counting the number of positrons with energy greater then 1.7 GeV in function of the time, the frequency precession of the muon spin is measured, that together with the measurement of the magnetic field, allows to extract $a_\mu$. The positrons are detected with 24 electromagnetic calorimeters, that measure the energy and the arrival time of the positrons, each made of 54 crystals of lead fluoride (PbF$_2$) read by silicon photomultipliers (SiPM). Together with calorimeters, two tracking detectors are used to make non-destructive measurements of the muon beam characteristic by reconstructing the muon decay position extrapolating backward the decay positrons. The knowledge of the beam motion inside the ring plays a fundamental role in the analysis of $a_{\mu}$, where the measured anomalous precession frequency $\omega_a$ must be corrected for four main beam dynamics effects. A first correction is associated with the presence of an electric field responsible of the vertical focusing on the storage ring, where vertical direction is orthogonal to the orbit and horizontal direction is along the storage ring radius. Due to the oscillations in the vertical direction, the so called vertical betatron oscillations, a second correction is necessary to account for an average angle associated with the muon motion off the ideal orbit. A third correction is caused by lost muons in the ring which have a different spin phase at the injection respect to the decay ones. Finally, due to the correlation between the vertical and horizontal beam motion of the muons and the acceptance of the calorimeter, a correction (called ``phase-acceptance") arises. In Run1, due to the presence of two damaged resistors in one electrostatic quadrupole, this correction had a prominent role. The replacement of the damaged resistors before Run2 decreased this effect by one order of magnitude, and a further reduction in Run-3 was provided by an improved orbit.\\ \noindent The work of this Thesis focuses on the beam dynamics corrections on $\omega_a$. Due to the prominent role in Run-1 a special attention has been put to the phase-acceptance correction\footnote{The effect of this correction on the E821 BNL measurement of the $g-2$ has been evaluated to be within the quoted systematic error.}. Chapter \ref{ch:anomalous_magnetic_moment} introduces the anomalous magnetic moment of the muon. Chapter \ref{ch:early_experiment} describes the history of the Muon $g-2$ experiments. Chapter \ref{ch:theory} discusses the Standard Model prediction and possible new physics scenario. Chapter \ref{ch:muon_g2_experiment} describes the E989 experiment, whose experimental technique and the apparatus are discussed focusing on the improvements needed to reach the final goal on $a_\mu$ measurement. The original work of the Thesis is discussed in the last 5 chapters. Chapter \ref{ch:experiment_technique} presents the analysis technique to extract $\omega_a$, Chapter \ref{ch:beam_dynamics} describes the beam dynamics in the E989 experiment, Chapter \ref{ch:beam_dynamics_correction} discusses the beam dynamics correction to $\omega_a$, and Chapter \ref{sec:pa_corr} focuses on the phase acceptance correction. Finally Chapter \ref{ch:run23_analysis_improvements} contains the Run-2/3 improvements which are expected to allow for an increased precision ons the beam dynamics co...

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗

Processes and mechanisms of coastal woody-plant mortality

Observations of woody plant mortality in coastal ecosystems are globally widespread, but the overarching processes and underlying mechanisms are poorly understood. This knowledge deficiency, combined with rapidly changing water levels, storm surges, atmospheric CO 2 , and vapor pressure deficit, creates large predictive uncertainty regarding how coastal ecosystems will respond to global change. Here we synthesize the literature on the mechanisms that underlie coastal woody-plant mortality, with the goal of producing a testable hypothesis framework. The key emergent mechanisms underlying mortality include hypoxic, osmotic, and ionic-driven reductions in whole-plant hydraulic conductance and photosynthesis that ultimately drive the coupled processes of hydraulic failure and carbon starvation. The relative importance of these processes in driving mortality, their order of progression, and their degree of coupling depends on the characteristics of the anomalous water exposure, on topographic effects, and on taxa-specific variation in traits and trait acclimation. Greater inundation exposure could accelerate mortality globally; however, the interaction of changing inundation exposure with elevated CO 2 , drought, and rising vapor pressure deficit could influence mortality likelihood. Models of coastal forests that incorporate the frequency and duration of inundation, the role of climatic drivers, and the processes of hydraulic failure and carbon starvation can yield improved estimates of inundation-induced woody-plant mortality.

salinity↗

Image Registration: A Necessary Evil

Registration of test and reference images is a key component of nearly all PSP data reduction techniques. This is done to ensure that a test image pixel viewing a particular point on the model is ratioed by the reference image pixel which views the same point. Typically registration is needed to account for model motion due to differing airloads when the wind-off and wind-on images are taken. Registration is also necessary when two cameras are used for simultaneous acquisition of data from a dual-frequency paint. This presentation will discuss the advantages and disadvantages of several different image registration techniques. In order to do so, it is necessary to propose both an accuracy requirement for image registration and a means for measuring the accuracy of a particular technique. High contrast regions in the unregistered images are most sensitive to registration errors, and it is proposed that these regions be used to establish the error limits for registration. Once this is done, the actual registration error can be determined by locating corresponding points on the test and reference images, and determining how well a particular registration technique matches them. An example of this procedure is shown for three transforms used to register images of a semispan model. Thirty control points were located on the model. A subset of the points were used to determine the coefficients of each registration transform, and the error with which each transform aligned the remaining points was determined. The results indicate the general superiority of a third-order polynomial over other candidate transforms, as well as showing how registration accuracy varies with number of control points. Finally, it is proposed that image registration may eventually be done away with completely. As more accurate image resection techniques and more detailed model surface grids become available, it will be possible to map raw image data onto the model surface accurately. Intensity ratio data can then be obtained by a "model surface ratio," rather than an image ratio. The problems and advantages of this technique will be discussed.

Bell, James↗

Reduced Order Modeling conditioned on monitored features for response and error bounds estimation in engineered systems

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to provide numerical error bounds associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Conditional VAEs↗

Geothermal-integrated thermally anisotropic building envelope for energy and peak-demand reduction

Buildings consume large amounts of energy for heating and cooling, while peak electricity demand places significant stress on the power grid. This paper presents a reduced-order co-simulation framework and load-oriented supervisory control strategy for a geothermal-integrated thermally anisotropic building envelope with a ground loop (TABE+GL). In TABE+GL, a hydronic loop embedded in the building envelope is directly coupled with a geothermal ground loop, allowing for bidirectional heat exchange between the envelope, the ground, and the indoor environment. A hybrid co-simulation framework was established by coupling a reduced-order resistor–capacitor (RC) thermal network model with EnergyPlus augmented with GHEDesigner modules. The RC model generated feasible heat flux options under three operating modes, and EnergyPlus predicted sensible loads, energy use, and pump energy demand. At each simulation step, a supervisory control algorithm selected the optimal loop configuration and duty factor that maximizes useful TABE geothermal utilization without exceeding the predicted sensible load, thereby avoiding overheating or cooling. Case studies were conducted for Los Angeles, California, Charleston, South Carolina, and Denver, Colorado. Results showed that the proposed framework reduced HVAC electricity consumption by 43%–67%, natural gas use for space heating by 11%–38%, and peak electricity demand by 43%–88%. These results highlight the potential of combining reduced-order envelope modeling, direct geothermal coupling, and load-oriented supervisory control to improve whole building energy performance and reduce peak demand across diverse weather conditions.

Howard, Daniel [Southern Adventist University]↗

Impact of thermodynamic fluctuations and pattern size on the nucleation behavior during area selective deposition

In this work, the authors explore the impact that thermodynamic fluctuations have on the spontaneous appearance of defects during area selective deposition on patterned surfaces that are fully covered by adsorbates under local thermodynamic equilibrium, such as self-assembled monolayers. By using a simple lattice gas model for the adsorbed monolayer, the authors were able to track the spontaneous formation of defects in the monolayer as a function of the pattern width. The results indicate that, for pattern widths of the order of tens of nanometers, roughening effects at the pattern edge can be the leading source for the spontaneous appearance of nucleation defects. This leads to an enhancement of the density of defects that can be up to three orders of magnitude higher than those expected in uniform (not patterned) surfaces. The model also predicts a density of defects that is inversely proportional to the pattern width. Finally, if the dynamic nucleation of defects during area selective deposition is driven by thermodynamic fluctuations, the model predicts that the nucleation rate should be proportional to the total precursor fluency and independent of purge times. Moreover, a tight confinement of the monolayer through a high quality smooth interface in the patterned substrate and strong cohesive interactions between adsorbates should each contribute to a reduction of the overall defect density.

36 MATERIALS SCIENCE↗

Solar Receiver with Integrated Thermal Storage for a Supercritical Carbon Dioxide Power Cycle

This project endeavors to design and demonstrate an integrated CSP system that meets the Department of Energy’s 6¢/kW-hr LCOE goal [1]. This design is based on the 10 MWe Supercritical Transformational Electric Power (STEP) engine specifications, but adaptable to any sCO 2 power module up to 10 MW e . The enabling technology to be tested is a new class of metal hydride based thermal energy storage materials which can store thermal energy in the hydrogen bonds between the metal hydride and hydrogen. This thermochemical energy storage solution provides a much greater energy density than sensible and latent energy storage enabling an integrated system design that closely linked the metal hydride to the power block. This allows the entire system to be mounted up tower to provide cost reductions. A novel metal hydride, capable of operating above 760 °C, was developed during the course of this project and thoroughly tested at a lab scale. The material’s hydrogen capacity, reaction enthalpy, thermal conductivity, and cyclic stability were measured. A demonstration test article was built that can store 2kWh of thermal energy and was coupled to Brayton Energy’s CO 2 loop where the performance of the metal hydride under commercial conditions was measured. An integrated model was created that uses the inputs from all the different sub-models in order to return an overall system LCOE using SAM. This project greatly advanced metal hydride based thermochemical energy storage technology and demonstrated its use. Previously there were no metal hydrides which were capable of operating at the proposed temperature limits (760°C), there is now a high temperature metal hydride with documented performance/stability and a path for low-cost production. Additionally, the metal hydride was tested at a scale >650x larger than what had previously been done in the lab. While not at commercial scale, it was an excellent demonstration of the technology and provided knowledge that can be leveraged for other applications.

08 HYDROGEN↗

Into the UV: A Precise Transmission Spectrum of HAT-P-41b Using Hubble’s WFC3/UVIS G280 Grism

The ultraviolet–visible wavelength range holds critical spectral diagnostics for the chemistry and physics at work in planetary atmospheres. To date, time-series studies of exoplanets to characterize their atmospheres have relied on several combinations of modes on the Hubble Space Telescope’s STIS/COS instruments to access this wavelength regime. Here for the first time, we apply the Hubble WFC3/UVIS G280 grism mode to obtain exoplanet spectroscopy from 200 to 800 nm in a single observation. We test the G280 grism mode on the hot Jupiter HAT-P-41b over two consecutive transits to determine its viability for the characterization of exoplanet atmospheres. We obtain a broadband transit depth precision of 29–33 ppm and a precision of on average 200 ppm in 10 nm spectroscopic bins. Spectral information from the G280 grism can be extracted from both the positive and negative first-order spectra, resulting in a 60% increase in the measurable flux. Additionally, the first Hubble Space Telescope orbit can be fully utilized in the time-series analysis. We present detailed extraction and reduction methods for use by future investigations with this mode, testing multiple techniques. We find the results to be fully consistent with STIS measurements of HAT-P-41b from 310 to 800 nm, with the G280 results representing a more observationally efficient and precise spectrum. HAT-P-41b’s transmission spectrum is best fit with a model with T {sub eq} = 2091 K, high metallicity, and significant scattering and cloud opacity. With these first-of-their-kind observations, we demonstrate that WFC3/UVIS G280 is a powerful new tool to obtain UV–optical spectra of exoplanet atmospheres, adding to the UV legacy of Hubble and complementing future observations with the James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

Reaction engineering for materials processing in space: Reduction of ilmenite by hydrogen and carbon monoxide

Oxygen is a consumable material which needs to be produced continuously in most space missions. Its use for propulsion as well as life support makes oxygen one of the largest volume chemicals to be produced in space. Production of oxygen from lunar materials is of particular interest and is very attractive possibility. The kinetics and mechanism of reduction of ilmenite by carbon monoxide and hydrogen at 800 to 1100 C were investigated. The temporal profiles of conversion for carbon monoxide have a sigmoidal shape and indicate the presence of three different stages (induction, acceleration, and deceleration) during the reduction reaction. The apparent activation energy decreases from 18 kcal/mole at 10 percent conversion to 10 kcal/mole at 50 percent conversion. The reaction is first order with respect to carbon monoxide under the experimental conditions studied. Both SEM and EDX analysis show that the diffusion of Fe product away from the reaction front and through the TiO2 phase, followed by the nucleation and growth of a separate Fe phase are important steps affecting the process kinetics. The results from hydrogen reduction show that the mechanism of ilmenite reduction by hydrogen is similar to that by carbon monoxide. However, the titanium dioxide can be further reduced by hydrogen at 800 to 1000 C. The detailed comparison and theoretical modeling of both reduction processes is presented.

Zhao, Y.↗

The Operation and Evolution of the Swift X-ray Telescope

The Swift X-ray Telescope (XRT) is a CCD based X-ray telescope designed for localization, spectroscopy and long term light curve monitoring of Gamma-Ray Bursts and their X-ray afterglows. Since the launch of Swift in November 2004, the XRT has undergone significant evolution in the way it is operated. Shortly after launch there was a failure of the thermo-electric cooler on the XRT CCD, which led to the XRT team being required to devise a method of keeping the XRT CCD temperature below 50C utilizing only passive cooling by minimizing the exposure of the XRT radiator to the Earth. We present in this paper an update on how the modeling of this passive cooling method has improved in first -1000 days since the method was devised, and the success rate of this method in day-to-day planning. We also discuss the changes to the operational modes and onboard software of the XRT. These changes include improved rapid data product generation in order to improve speed of rapid Gamma-Ray Burst response and localization to the community; changes to the way XRT observation modes are chosen in order to better fine tune data aquisition to a particular science goal; reduction of "mode switching" caused by the contamination of the CCD by Earth light or high temperature effects.

Kennea, Jamie↗

All-multiplicity one-loop amplitudes in Born-Infeld electrodynamics from generalized unitarity

We initiate a study of non-supersymmetric Born-Infeld electrodynamics in 4d at the quantum level. Explicit all-multiplicity expressions are calculated for the purely rational one-loop amplitudes in the self-dual (+ + … +) and next-to-self-dual (- + … +) helicity sectors. Using a supersymmetric decomposition, d-dimensional unitarity cuts of the integrand factorize into tree-amplitudes in a 4d model of Born-Infeld photons coupled to a massive complex scalar. The two-scalar tree-amplitudes needed to construct the Born-Infeld integrand are computed using two complimentary approaches: (1) as a double-copy of Yang-Mills coupled to a massive adjoint scalar with a dimensionally reduced form of Chiral Perturbation Theory, and (2) by imposing consistency with low-energy theorems under a reduction from 4d to 3d and T-duality. The Born-Infeld integrand is integrated in d = 4 - 2ϵ dimensions at order Ο(ϵ 0 ) using the dimension-shifting formalism. We comment on the implications for electromagnetic duality in quantum Born-Infeld theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gas-phase Radiolysis in Plutonium Dioxide Powder

A sample of PuO 2 powder undergoes self radiolysis, causing water present in the solid to generate hydrogen, oxygen, and minor species (e.g. H 2 O 2 ). This is a matter of concern for stored PuO 2 , since the gases could react exothermically and rupture the container. In experiments at LANL, most PuO 2 /salt samples generate mostly H 2 , but a few samples generate a mixture of H 2 and O 2 together. In these special cases, the mixture of H 2 and O 2 in the pore volumes of the powder will then be subject to radiolysis, causing the reverse reaction H 2 + O 2 $\Rightarrow$ H 2 O. In order to understand the steady-state quantities of these gases, it is necessary to estimate the amount of energy from the alpha irradiation that drives this reaction. This report outlines one possible approach for this process, yielding a value for the fractional amount of radiation deposited in the gas phase of 294 ppm. This corresponds to the case of a powder with porosity of 78%. We then use a comprehensive chemical kinetics model based on elementary gas phase reactions to calculate the rate of conversion back to water. We estimate a reduction in the H 2 pressure of about 12 kPa/month, which applies only when O 2 is simultaneously present in the gas phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A low-order global gravity field of Venus and dynamical implications

An estimate of sixth-degree and sixth-order harmonic coefficients of the global gravity field of Venus is obtained by processing the long periodic variations of the mean orbital elements of the Pioneer Venus orbiter. Approximately 220 days of data are included in this reduction, which provides almost complete longitudinal coverage. Oblateness is estimated to be -5.97 + or - 3.2 x 10 to the -6. It is noted that the amplitudes of other coefficients are similar to the predicted coefficients using Kaula's rule under equal stress assumption. Atmospheric density values as a function of altitude are obtained to help model drag perturbation. A radial acceleration map at 100 km above the Venus surface is generated, and correlation between gravity anomalies and major topographic features is observed. Spectral analysis of the harmonic model suggests that the interior density anomalies are like those on earth. The orientation angles of the principal axes of the moments of inertia are computed, and deviation of the maximum moment of inertia axis from the spin axis is observed to be small (less than 5 deg). It is found that the minimum moment of inertia axis passes through the Aphrodite and Beta regions.

Ananda, M. P.↗

A generic multi-flex-body dynamics, controls simulation tool for space station

An order (n) multiflex body Space Station simulation tool is introduced. The flex multibody modeling is generic enough to model all phases of Space Station from build up through to Assembly Complete configuration and beyond. Multibody subsystems such as the Mobile Servicing System (MSS) undergoing a prescribed translation and rotation are also allowed. The software includes aerodynamic, gravity gradient, and magnetic field models. User defined controllers can be discrete or continuous. Extensive preprocessing of 'body by body' NASTRAN flex data is built in. A significant aspect, too, is the integrated controls design capability which includes model reduction and analytic linearization.

London, Ken W.↗