Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bar codes”

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 181 records · Page 10

Application of the coded long-pulse technique to plasma line studies of the ionosphere

Recently, the coded long-pulse radar technique was tested at Arecibo Observatory, Puerto Rico using photoelectron-enhanced plasma lines in the daytime ionosphere. The technique immediately proved to be a powerful diagnostic tool for studying natural ionospheric phenomena. Our initial observations indicate that extremely accurate measurements of absolute electron density (0.01 to 0.03% error bars) can be achieved with an altitude resolution of 150 m and a temporal resolution of approx. 2 s. In addition, the technique provides information about electron density structure within a 150-m altitude cell and yields parameters from which the energy spectrum of suprathermal electrons (equal to or greater than 5 eV) can be deduced. Our earliest measurements are used to illustrate applications of the coded long-pulse technique to several aeronomic/ionsospheric areas of current interest. These include studies of neutral wave motions in the lower thermosphere, measurements of ion composition in the F(sub 1) region/upper ionosphere, and investigations of electron-gas thermal balance and photoelectron energy loss processes. The technique can be utilized to examine irregularity formation in the F region, probe electron acceleration processes in ionospheric modification experiments, verify the magnetic field dependence of Langmuir wave damping, and more generally test higher order corrections suggested for the Langmuir dispersion relation. It is anticipated that the latter tests will facilitate measurements of ionospheric currents.

Djuth, Frank T.↗

Nuclear data uncertainty and sensitivity analysis of the VHTRC benchmark using SCALE

The Very High Temperature Critical Assembly (VHTRC) experiment represents one of the few data sets available for the validation of HTGR lattice physics. Within the framework of the IAEA Coordinated Research Project on HTGR Uncertainty Analysis in Modeling, uncertainty and sensitivity analyses of this graphite-moderated facility are performed as the validation reference to the prismatic MHTGR-350 fuel block calculations. Nominal multi-group and continuous-energy criticality calculations with the KENO-VI Monte Carlo code of the SCALE code package are compared with the continuous-energy Monte Carlo Code Serpent. Good agreement with a maximum difference of 250 pcm is obtained. The experimental data set, however, differs by several hundred pcm with the results of both codes when using the ENDF-VII.0 library. When using ENDF/B-VII.1 data, this difference is reduced to a few hundred pcm and the calculations lie within the experimental error bars. Uncertainties of the VHTRC multiplication factors due to uncertainties in nuclear data are determined with SAMPLER/KENO-VI of SCALE 6.2 and with continuous-energy TSUNAMI of SCALE 6.2. For all experimental configurations, the obtained uncertainty is found to be 0.58% when using ENDF/B-VII.0 data. The top contributor to this uncertainty is the average number of neutrons per fission event of U-235. With ENDF/B-VII.1 data, uncertainties of about 0.66% are obtained in particular due to an increased uncertainty of U-235 nu-bar in the latest ENDF release. In conclusion, when considering these nuclear data uncertainties, the obtained 1s uncertainty intervals are overlapping with the experimental error bars.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TankSIM: A Cryogenic Tank Performance Prediction Program

Developed for predicting the behavior of cryogenic liquids inside propellant tanks under various environmental and operating conditions. Provides a multi-node analysis of pressurization, ullage venting and thermodynamic venting systems (TVS) pressure control using axial jet or spray bar TVS. Allows user to combine several different phases for predicting the liquid behavior for the entire flight mission timeline or part of it. Is a NASA in-house code, based on FORTRAN 90-95 and Intel Visual FORTRAN compiler, but can be used on any other platform (Unix-Linux, Compaq Visual FORTRAN, etc.). The last Version 7, released on December 2014, included detailed User's Manual. Includes the use of several RefPROP subroutines for calculating fluid properties.

Bolshinskiy, L. G.↗

Computations of Axisymmetric Flows in Hypersonic Shock Tubes

A time-accurate two-dimensional fluid code is used to compute test times in shock tubes operated at supersonic speeds. Unlike previous studies, this investigation resolves the finer temporal details of the shock-tube flow by making use of modern supercomputers and state-of-the-art computational fluid dynamic solution techniques. The code, besides solving the time-dependent fluid equations, also accounts for the finite rate chemistry in the hypersonic environment. The flowfield solutions are used to estimate relevant shock-tube parameters for laminar flow, such as test times, and to predict density and velocity profiles. Boundary-layer parameters such as bar-delta(sub u), bar-delta(sup *), and bar-tau(sub w), and test time parameters such as bar-tau and particle time of flight t(sub f), are computed and compared with those evaluated by using Mirels' correlations. This article then discusses in detail the effects of flow nonuniformities on particle time-of-flight behind the normal shock and, consequently, on the interpretation of shock-tube data. This article concludes that for accurate interpretation of shock-tube data, a detailed analysis of flowfield parameters, using a computer code such as used in this study, must be performed.

Sharma, Surendra P.↗

A structural model of the long-term degradation of the concrete biological shield

The concrete biological shield (CBS) of light water reactors is exposed to high neutron radiation dose in the long term, which may lead to the degradation of the concrete’s mechanical properties. Given the important shielding role of the CBS, it is necessary to investigate the irradiation effects at the structural scale and provide estimates of the damage extent from the wall’s inner surface to study potential license renewals. For this purpose, we developed a mechanical model accounting for radiation-induced expansion, creep, and damage in concrete using the Grizzly finite element code, informed by ex-core neutron flux calculations using the VERA tool. The model was applied to a 3D CBS structure represented by the CBS wall, a steel liner, reinforcement bars, and a concrete base mat and evaluated damage at 40, 60, and 80 years of operation. The VERA model predicted a maximum fluence of approximately 2 x 10 19 ncm -2 at 80 years of operation. The results showed that damage is highest at the inner surface of the CBS wall and gradually decreases with depth. It extends beyond the rebar after 60 years and reaches a depth of approximately 12 cm at 80 years.

42 ENGINEERING↗

NASTRAN nonlinear vibration analysis of beam and frame structures

A capability for the nonlinear vibration analysis of beam and frame structures suitable for use with NASTRAN level 15.5 is described. The nonlinearity considered is due to the presence of axial loads induced by longitudinal end restraints and lateral displacements that are large compared to the beam height. A brief discussion is included of the mathematical analysis and the geometrical stiffness matrix for a prismatic beam (BAR) element. Also included are a brief discussion of the equivalent linearization iterative process used to determine the nonlinear frequency, the required modifications to subroutines DBAR and XMPLBD of the NASTRAN code, and the appropriate vibration capability, four example problems are presented. Comparisons with existing experimental and analytical results show that excellent accuracy is achieved with NASTRAN in all cases.

Mei, C.↗

Modification and evaluation of a Barnes-type objective analysis scheme for surface meteorological data

The Purdue Regional Objective Analysis of the Mesoscale (PROAM) is a Barness-type scheme for the analysis of surface meteorological data. Modifications are introduced to the original version in order to increase its flexibility and to permit greater ease of usage. The code was rewritten for an interactive computer environment. Furthermore, a multiple iteration technique suggested by Barnes was implemented for greater accuracy. PROAM was subjected to a series of experiments in order to evaluate its performance under a variety of analysis conditions. The tests include use of a known analytic temperature distribution in order to quantify error bounds for the scheme. Similar experiments were conducted using actual atmospheric data. Results indicate that the multiple iteration technique increases the accuracy of the analysis. Furthermore, the tests verify appropriate values for the analysis parameters in resolving meso-beta scale phenomena.

Smith, D. R.↗

Analysis of Photogrammetry Data from ISIM Mockup

During ground testing of the Integrated Science Instrument Module (ISIM) for the James Webb Space Telescope (JWST), the ISIM Optics group plans to use a Photogrammetry Measurement System for cryogenic calibration of specific target points on the ISIM composite structure and Science Instrument optical benches and other GSE equipment. This testing will occur in the Space Environmental Systems (SES) chamber at Goddard Space Flight Center. Close range photogrammetry is a 3 dimensional metrology system using triangulation to locate custom targets in 3 coordinates via a collection of digital photographs taken from various locations and orientations. These photos are connected using coded targets, special targets that are recognized by the software and can thus correlate the images to provide a 3 dimensional map of the targets, and scaled via well calibrated scale bars. Photogrammetry solves for the camera location and coordinates of the targets simultaneously through the bundling procedure contained in the V-STARS software, proprietary software owned by Geodetic Systems Inc. The primary objectives of the metrology performed on the ISIM mock-up were (1) to quantify the accuracy of the INCA3 photogrammetry camera on a representative full scale version of the ISIM structure at ambient temperature by comparing the measurements obtained with this camera to measurements using the Leica laser tracker system and (2), empirically determine the smallest increment of target position movement that can be resolved by the PG camera in the test setup, i.e., precision, or resolution. In addition, the geometrical details of the test setup defined during the mockup testing, such as target locations and camera positions, will contribute to the final design of the photogrammetry system to be used on the ISIM Flight Structure.

Nowak, Maria↗

Analysis of Photogrammetry Data from ISIM Mockup, June 1, 2007

During ground testing of the Integrated Science Instrument Module (ISIM) for the James Webb Space Telescope (JWST), the ISIM Optics group plans to use a Photogrammetry Measurement System for cryogenic calibration of specific target points on the ISIM composite structure and Science Instrument optical benches and other GSE equipment. This testing will occur in the Space Environmental Systems (SES) chamber at Goddard Space Flight Center. Close range photogrammetry is a 3 dimensional metrology system using triangulation to locate custom targets in 3 coordinates via a collection of digital photographs taken from various locations and orientations. These photos are connected using coded targets, special targets that are recognized by the software and can thus correlate the images to provide a 3 dimensional map of the targets, and scaled via well calibrated scale bars. Photogrammetry solves for the camera location and coordinates of the targets simultaneously through the bundling procedure contained in the V-STARS software.

Nowak, Maria↗

Influence of Real Gas Effects on Chemical Kinetics in Oxycombustion in Supercritical Carbon Dioxide

Oxycombustion in supercritical carbon dioxide is an integral part of the Allam Cycle, a technology that enables carbon-neutral use of fossil-fuels and carbon-negative use of biofuels. We simulate oxycombustion in a realistic combustor geometry for two sets of fuel conditions: pure methane, and a 40 percent methane/60 percent carbon dioxide blend, both at 343.15 K. The fuel jet mixes with a preheated swirler of 20 percent oxygen and 80 percent carbon dioxide at 1005.35 K, with a 100 percent carbon dioxide coflow at 783.15 K. The entire system operates at a pressure of 300 bar, putting the entire system above the critical temperature and pressure of carbon dioxide. Simulations are performed using PeleC, a compressible block-structured adaptive-mesh refinement (AMR) reacting flow code. Direct comparison of results using thermodynamically self-consistent implementations of the ideal gas equation of state and the Soave-Redlich-Kwong (SRK) equation of state allow quantification of the real gas impacts on combustion kinetics, flow temperatures, and pollutant formation.

carbon dioxide↗

Integral Reactor Containment Condensation Model and Experimental Validation

This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). In the three years of investigation duration, following the original proposal, the planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). A film flow model was developed for the scaling analysis, and the results suggest that the 1/3 scaled test facility covers large portion of laminar film flow, leading to a lower average heat transfer coefficient comparing to the prototypic value. Although it is conservative in reactor safety analysis, the significant reduction of heat transfer coefficient (50%) could under estimate the prototypic condensation heat transfer rate, resulting in inaccurate prediction of the decay heat removal capability. Further investigation is thus needed to quantify the scaling distortion for safety analysis code validation. Experimental investigations were performed in the existing MASLWR test facility at OrST with minor modifications. A total of 13 containment condensation tests were conducted for pressure ranging from 4 to 21 bar with three different static inventories of non-condensable gas. Condensation and heat transfer rates were evaluated employing several methods, notably from measured temperature gradients in the HTP as well as measured condensate formation rates. A detailed mass and energy accounting was used to assess the various measurement methods and to support simplifying assumptions required for the analysis. Condensation heat fluxes and heat transfer coefficients are calculated and presented as a function of pressure to satisfy the objectives of this investigation. The major conclusions for those tests are summarized below: (1) In the steam blow-down tests, the initial condensation heat transfer process involves the heating-up of the containment heat transfer plate. An inverse heat conduction model was developed to capture the rapid transient transfer characteristics, and the analysis method is applicable to SMR safety analysis. (2) The average condensation heat transfer coefficients for different pressure conditions and non-condensable gas mass fractions were obtained from the integral test facility, through the measurements of the heat conduction rate across the containment heat transfer plate, and from the water condensation rates measurement based on the total energy balance equation. 15 (3) The test results using the measured HTP wall temperatures are considerably lower than popular condensation models would predict mainly due to the side wall conduction effects in the existing MASLWR integral test facility. The data revealed the detailed heat transfer characteristics of the model containment, important to the SMR safety analysis and the validation of associated evaluation model. However this approach, unlike separate effect tests, cannot isolate the condensation heat transfer coefficient over the containment wall, and therefore is not suitable for the assessment of the condensation heat transfer coefficient against system pressure and noncondensable gas mass fraction. (4) The average condensation heat transfer coefficients measured from the water condensation rates through energy balance analysis are appropriate, however, with considerable uncertainties due to the heat loss and temperature distribution on the containment wall. With the consideration of the side wall conduction effects, the results indicate that the measured heat transfer coefficients in the tests is about 20% lower than the prediction of Dehbi’s correlation, mainly due to the side wall conduction effects. The investigation also indicates an increase in the condensation heat transfer coefficient at high containment pressure conditions, but the uncertainties invoked with this method appear to be substantial. (5) Non-condensable gas in the tests has little effects on the condensation heat transfer at high elevation measurement ports. It does affect the bottom measurements near the water level position. The results suggest that the heavier non-condensable gas is accumulated in the lower portion of the containment due to stratification in the narrow containment space. The overall effects of the non-condensable gas on the heat transfer process should thus be negligible for tall containments of narrow condensation spaces in most SMR designs. Therefore, the previous correlations with noncondensable gas effects are not appropriate to those small SMR containments due to the very poor mixing of steam and non-condensable gas. The MELCOR simulation results agree with the experimental data reasonably well. However, it is observed that the MELCOR overpredicts the heat flux for all analyzed tests. The MELCOR predicts that the heat fluxes for CCT’s approximately range from 30 to 45 kW/m2 whereas the experimental data (averaged) ranges from about 25 to 40 kW/m2. This may be due to the limited availability of liquid film models included in MELCOR. Also, it is believed that due to complex test geometry, measured temperature gradients across the heat transfer plate may have been underestimated and thus the heat flux had been underestimated. The MELCOR model predicts a film thickness on the order of 100 microns, which agrees very well with film flow model developed in this study for scaling analysis. However, the expected differences in film thicknesses for near vacuum and near atmospheric test conditions are not significant. Further study on the behavior of condensate film is expected to refine the simulation results. Possible refinements include but are not limited to, the followings: CFD simulation focusing on the liquid film behavior and benchmarking with experimental analyses for simpler geometries. 16 1 INTRODUCTION This NEUP funded project, NEUP 12-3630, is for experimental, numerical and analytical studies on high-pressure steam condensation phenomena in a steel containment vessel connected to a water cooling tank, carried out at Oregon State University (OrSU) and the University of Wisconsin at Madison (UW-Madison). The experimental results are employed to validate the containment condensation model in reactor containment system safety analysis code for integral SMRs. Such a containment condensation model is important to demonstrate the adequate cooling. In the three years of investigation, following the original proposal, the following planned tasks have been completed: (1) Performed a scaling study for the full pressure test facility applicable to the reference design for the condensation heat transfer process during design basis accidents (DBAs), modified the existing test facility to route the steady-state secondary steam flow into the high pressure containment for controllable condensation tests, and extended the operations at negative gage pressure conditions (OrSU). (2) Conducted a series of DBA and quasi-steady experiments using the full pressure test facility to provide a reliable high pressure condensation database (OrSU). (3) Analyzed experimental data and evaluated condensation model for the experimental conditions, and predicted the prototypic containment performance under accidental conditions (UW-Madison). The results are applicable to integral Small Modular Reactor (SMR) designs, including NuScale, mPower, Westinghouse SMR, Holtec-160 and other integral reactors with small containments of relatively high pressures under accidental conditions. Testing has been conducted at the OrSU laboratory in the existing MASLWR (Multi-Application Small Light Water Reactor) integral test facility sponsored by the US Department of Energy. Its highpressure stainless steel containment model (~2 MPa) is scaled to the NuScale SMR currently under development at NuScale Power, Inc.. Minor modifications to the model containment have been made to control the non-condensable gas fraction and to utilize the secondary loop stable steam flow for condensation testing. UW-Madison has developed a containment condensation model, which leveraged previous validated containment heat transfer work carried out at UW-Madison, and extended the range of applicability of the model to integral SMR designs that utilize containment vessels of high heat transfer efficiencies. In this final report, the research background and literature survey are presented in Chapter 2 and 3, respectively. The test facility description and modifications are summarized in Chapter 4, and the scaling analysis is introduced in Chapter 5. The tests description, procedures, and data analysis are presented in Chapter 6, while the numerical modeling is presented in Chapter 7, followed by a conclusion section in Chapter 8.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Gas Gun Model and Comparison to Experimental Performance of Pipe Guns Operating with Light Propellant Gases and Large Cryogenic Pellets

Injection of multiple large (~10 to 30 mm diameter) shattered pellets into ITER plasmas is presently part of the scheme planned to mitigate the deleterious effects of disruptions on the vessel components. To help in the design and optimize performance of the pellet injectors for this application, a model referred to as “the gas gun simulator” has been developed and benchmarked against experimental data. The computer code simulator is a Java program that models the gas-dynamics characteristics of a single-stage gas gun. Following a stepwise approach, the code utilizes a variety of input parameters to incrementally simulate and analyze the dynamics of the gun as the projectile is launched down the barrel. Using input data, the model can calculate gun performance based on physical characteristics, such as propellant-gas and fast-valve properties, barrel geometry, and pellet mass. Although the model is fundamentally generic, the present version is configured to accommodate cryogenic pellets composed of H2, D2, Ne, Ar, and mixtures of them and light propellant gases (H2, D2, and He). The pellets are solidified in situ in pipe guns that consist of stainless steel tubes and fast-acting valves that provide the propellant gas for pellet acceleration (to speeds ~200 to 700 m/s). The pellet speed is the key parameter in determining the response time of a shattered pellet system to a plasma disruption event. The calculated speeds from the code simulations of experiments were typically in excellent agreement with the measured values. With the gas gun simulator validated for many test shots and over a wide range of physical and operating parameters, it is a valuable tool for optimization of the injector design, including the fast valve design (orifice size and volume) for any operating pressure (~40 bar expected for the ITER application) and barrel length for any pellet size (mass, diameter, and length). Key design parameters and proposed values for the pellet injectors for the ITER disruption mitigation systems are discussed.

Combs, Stephen Kirk↗

A new multi-dimensional general relativistic neutrino hydrodynamics code for core-collapse supernovae. IV. The neutrino signal

Considering six general relativistic, two-dimensional (2D) supernova (SN) explosion models of progenitor stars between 8.1 and 27 M {sub ☉}, we systematically analyze the properties of the neutrino emission from core collapse and bounce to the post-explosion phase. The models were computed with the VERTEX-COCONUT code, using three-flavor, energy-dependent neutrino transport in the ray-by-ray-plus approximation. Our results confirm the close similarity of the mean energies, (E), of ν-bar {sub e} and heavy-lepton neutrinos and even their crossing during the accretion phase for stars with M ≳ 10 M {sub ☉} as observed in previous 1D and 2D simulations with state-of-the-art neutrino transport. We establish a roughly linear scaling of 〈E{sub ν-bar{sub e}}〉 with the proto-neutron star (PNS) mass, which holds in time as well as for different progenitors. Convection inside the PNS affects the neutrino emission on the 10%-20% level, and accretion continuing beyond the onset of the explosion prevents the abrupt drop of the neutrino luminosities seen in artificially exploded 1D models. We demonstrate that a wavelet-based time-frequency analysis of SN neutrino signals in IceCube will offer sensitive diagnostics for the SN core dynamics up to at least ∼10 kpc distance. Strong, narrow-band signal modulations indicate quasi-periodic shock sloshing motions due to the standing accretion shock instability (SASI), and the frequency evolution of such 'SASI neutrino chirps' reveals shock expansion or contraction. The onset of the explosion is accompanied by a shift of the modulation frequency below 40-50 Hz, and post-explosion, episodic accretion downflows will be signaled by activity intervals stretching over an extended frequency range in the wavelet spectrogram.

79 ASTRONOMY AND ASTROPHYSICS↗

FY17 Status Report on the Initial EPP Finite Element Analysis of Grade 91 Steel

This report describes a modification to the elastic-perfectly plastic (EPP) strain limits design method to account for cyclic softening in Gr. 91 steel. The report demonstrates that the unmodified EPP strain limits method described in current ASME code case is not conservative for materials with substantial cyclic softening behavior like Gr. 91 steel. However, the EPP strain limits method can be modified to be conservative for softening materials by using softened isochronous stress-strain curves in place of the standard curves developed from unsoftened creep experiments. The report provides softened curves derived from inelastic material simulations and factors describing the transformation of unsoftened curves to a softened state. Furthermore, the report outlines a method for deriving these factors directly from creep/fatigue tests. If the material softening saturates the proposed EPP strain limits method can be further simplified, providing a methodology based on temperature-dependent softening factors that could be implemented in an ASME code case allowing the use of the EPP strain limits method with Gr. 91. Finally, the report demonstrates the conservatism of the modified method when applied to inelastic simulation results and two bar experiments.

36 MATERIALS SCIENCE↗

Neptune's Latitudinal Variations as Viewed with ALMA

We present spatially resolved millimeter maps of Neptune between 95 and 242 GHz taken with the Atacama Large Millimeter/submillimeter Array (ALMA) in 2016–2017. The millimeter weighting functions peak between 1 and 10 bar on Neptune, lying in between the altitudes probed at visible/infrared and centimeter wavelengths. Thus, these observations provide important constraints on the atmospheric structure and dynamics of Neptune. We identify seven well-resolved latitudinal bands of discrete brightness temperature variations, on the order of 0.5–3 K in all three observed ALMA spectral bands. We model Neptune’s brightness temperature using the radiative-transfer code Radio-BEAR and compare how various H{sub 2}S, CH{sub 4}, and ortho-/para-H{sub 2} abundance profiles can fit the observed temperature variations across the disk. We find that observed variations in brightness temperature with latitude can be explained by variations in the H{sub 2}S profile that range from sub- to supersaturations at altitudes above the 10 bar pressure level, while variations in CH{sub 4} improve the quality of fit near the equator. At the south polar cap, our best-fit model has a depleted deep atmospheric abundance of H{sub 2}S from 30 to only 1.5 times the protosolar value, while simultaneously depleting the CH{sub 4} abundance. This pattern of enhancement and depletion of condensible species is consistent with a global circulation structure where enriched air rises at the midlatitudes (32°–12°S) and north of the equator (2°–20°N), and dry air descends at the poles (90°–66°S) and just south of the equator (12°S–2°N). Our analysis finds more complex structure near the equator than accounted for in previous circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗

A new look at the simultaneous analysis and design of structures

The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.

Striz, Alfred G.↗

SINDA-NASTRAN interfacing program theoretical description and user's manual

The task of converting SINDA finite difference thermal model temperature results into NASTRAN finite element model thermal loads can be very labor intensive if there is not one node-to-one element, or systematic node-to-element. correlation between models. This paper describes the SINDA-NASTRAN Interfacing Program (SNIP), a FORTRAN computer code that generates NASTRAN structural model thermal load cards given by SINDA (or similar thermal model) temperature results and thermal model geometric data. SNIP generates NASTRAN thermal load cards for NASTRAN plate, shell, bar, and beam elements. The paper describes the interfacing procedures used by SNIP, and discusses set-up and operation of the program. Sample cases are included to demonstrate use of the program and show its performance under a variety of conditions. SNIP can provide structural model thermal loads that accurately reflect thermal model results while reducing the time required to interface thermal and structural models when compared to other methods.

Winegar, Steven R.↗