Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

Measuring thermal profiles in high explosives using neural networks

We present a new method for calculating the temperature profile of high explosive (HE) material using a Convolutional Neural Network (CNN). To train/test the CNN, we have developed a hybrid experiment/simulation method for collecting acoustic and temperature data. We experimentally heat cylindrical containers of HE material until detonation/deflagration, where we continuously measure the acoustic bursts through the HE using multiple acoustic transducers lined around the exterior container circumference. However, measuring the temperature profile in the HE in an experiment would require inserting a large number of thermal probes, which would disrupt the heating process. Thus, we use two thermal probes, one at the HE center and one at the wall. We then use numerical simulation of the heating process to calculate the temperature distribution and correct the simulated temperatures based on the experimental center and wall temperatures. We calculate temperature errors on the order of 15 °C, which is ~12% of the range of temperatures in the experiment. We also investigate how the algorithm’s accuracy is affected by the number of acoustic receivers used to collect each measurement and the resolution of the temperature prediction. This work provides a means of assessing the safety status of HE material, which cannot be achieved using existing temperature measurement methods. In addition, it has implications for a range of other applications where internal temperature profile measurements would provide critical information. These applications include detecting chemical reactions, observing thermodynamic processes such as combustion, monitoring metal or plastic casting, determining the energy density in thermal storage capsules, and identifying abnormal battery operations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Convective penetrations and observable consequences of helium shell flashes in 3 and 6 solar-mass red supergiants

The results of 61 detailed numerical flash calculations are examined, with particular reference to 3 and 6 solar-mass stars with core masses ranging from 0.56 to 0.60 the solar mass and 0.94 to 0.95 the solar mass, respectively. First penetration of the convective envelope into pockets of carbon in the intershell region was found for the 15th flash for the 6 solar-mass case, with continuing deeper penetration for the flashes thereafter. However, no such type of envelope penetration could be found for the 3 solar-mass case.

Sackmann, I.-J.↗

Analysis of forced convection heat flow effects in horizontal ribbon growth from the melt

A heat transport analysis which considers forced convective fluid flow induced by the motion of a continuous solid ribbon over a melt has been done for horizontal ribbon growth. A model has been developed which treats both 'active' and 'passive' cooling at the ribbon surface. The results show that heat flow from the melt requires active cooling in the region of the leading growth edge or growth tip. Steady-state liquid-solid interface shape is analyzed and numerical results are given for steady-state pulling of silicon ribbon.

Zoutendyk, J. A.↗

A fully vectorized numerical solution of the incompressible Navier-Stokes equations

A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.

Patel, N.↗

Energy Flow Continuity in Solar Active Regions

The models for sunspots are combined into an active region model with consideration for the energy flow beneath active regions. An apparent average energy balance exists between the sunspot deficit and the facular excess, i.e., no 11 year variations in solar luminosity associated with the activity centers. This is seen as a consequence of the upper convection zone's inability to store these significant amounts of energy for periods greatly in excess of weeks. This view is supported by observed active region behavior and detailed numerical modelling. Increases in facular and spot brightness are nearly commensurate, with the faculae outlasting the spots on time scales of the order of weeks to a couple of months. Foukal finds the radiation (deficit from a sunspot blocking model) recovers slowly on a timescale of approximately 83 days.

Schatten, K. H.↗

Estimating net rainfall, evaporation and water storage of a bare soil from sequential L-band emissivities

A general method to use a time series of L-band emissivities as an input to a hydrological model for continuously monitoring the net rainfall and evaporation as well as the water content over the entire soil profile is proposed. The model requires a sufficiently accurate and general relation between soil emissivity and surface moisture content. A model which requires the soil hydraulic properties as an additional input, but does not need any weather data was developed. The method is shown to be numerically consistent.

Stroosnijder, L.↗

Data communication requirements for the advanced NAS network

The goal of the Numerical Aerodynamic Simulation (NAS) Program is to provide a powerful computational environment for advanced research and development in aeronautics and related disciplines. The present NAS system consists of a Cray 2 supercomputer connected by a data network to a large mass storage system, to sophisticated local graphics workstations, and by remote communications to researchers throughout the United States. The program plan is to continue acquiring the most powerful supercomputers as they become available. In the 1987/1988 time period it is anticipated that a computer with 4 times the processing speed of a Cray 2 will be obtained and by 1990 an additional supercomputer with 16 times the speed of the Cray 2. The implications of this 20-fold increase in processing power on the data communications requirements are described. The analysis was based on models of the projected workload and system architecture. The results are presented together with the estimates of their sensitivity to assumptions inherent in the models.

Levin, Eugene↗

Scaling of data communications for an advanced supercomputer network

The goal of NASA's Numerical Aerodynamic Simulation (NAS) Program is to provide a powerful computational environment for advanced research and development in aeronautics and related disciplines. The present NAS system consists of a Cray 2 supercomputer connected by a data network to a large mass storage system, to sophisticated local graphics workstations and by remote communication to researchers throughout the United States. The program plan is to continue acquiring the most powerful supercomputers as they become available. The implications of a projected 20-fold increase in processing power on the data communications requirements are described.

Levin, E.↗

Overview of numerical codes developed for predicted electrothermal deicing of aircraft blades

An overview of the deicing computer codes that have been developed at the University of Toledo under sponsorship of the NASA-Lewis Research Center is presented. These codes simulate the transient heat conduction and phase change occurring in an electrothermal deicier pad that has an arbitrary accreted ice shape on its surface. The codes are one-dimensional rectangular, two-dimensional rectangular, and two-dimensional with a coordinate transformation to model the true blade geometry. All modifications relating to the thermal physics of the deicing problem that have been incorporated into the codes will be discussed. Recent results of reformulating the codes using different numerical methods to increase program efficiency are described. In particular, this reformulation has enabled a more comprehensive two-dimensional code to run in much less CPU time than the original version. The code predictions are compared with experimental data obtained in the NASA-Lewis Icing Research Tunnel with a UH1H blade fitted with a B. F. Goodrich electrothermal deicer pad. Both continuous and cyclic heater firing cases are considered. The major objective in this comparison is to illustrate which codes give acceptable results in different regions of the airfoil for different heater firing sequences.

Keith, Theo G.↗

STRUTEX: A prototype knowledge-based system for initially configuring a structure to support point loads in two dimensions

Only recently have engineers begun making use of Artificial Intelligence (AI) tools in the area of conceptual design. To continue filling this void in the design process, a prototype knowledge-based system, called STRUTEX has been developed to initially configure a structure to support point loads in two dimensions. This prototype was developed for testing the application of AI tools to conceptual design as opposed to being a testbed for new methods for improving structural analysis and optimization. This system combines numerical and symbolic processing by the computer with interactive problem solving aided by the vision of the user. How the system is constructed to interact with the user is described. Of special interest is the information flow between the knowledge base and the data base under control of the algorithmic main program. Examples of computed and refined structures are presented during the explanation of the system.

Robers, James L.↗

On the continuous dependence with respect to sampling of the linear quadratic regulator problem for distributed parameter system

The convergence of solutions to the discrete- or sampled-time linear quadratic regulator problem and associated Riccati equation for infinite-dimensional systems to the solutions to the corresponding continuous time problem and equation, as the length of the sampling interval (the sampling rate) tends toward zero(infinity) is established. Both the finite-and infinite-time horizon problems are studied. In the finite-time horizon case, strong continuity of the operators that define the control system and performance index, together with a stability and consistency condition on the sampling scheme are required. For the infinite-time horizon problem, in addition, the sampled systems must be stabilizable and detectable, uniformly with respect to the sampling rate. Classes of systems for which this condition can be verified are discussed. Results of numerical studies involving the control of a heat/diffusion equation, a hereditary or delay system, and a flexible beam are presented and discussed.

Rosen, I. G.↗

Demonstration of Plasma Arc Environmental Technology Applications for the Demilitrization of DOD Stockpiles

Since 1989 the US Army Construction Engineering Research Laboratories (USACERL) have been active participants in the research and development towards establishing Plasma Arc Technology (PAT) as an efficient, economical, and safe hazardous waste immobilization tool. A plasma torch capable of generating high temperatures makes this technology a viable and powerful tool for the thermal destruction of various military industrial waste streams into an innocuous ceramic material no longer requiring hazardous waste landfill disposal. The emerging plasma environmental thermal treatment process has been used to safely and efficiently meet the waste disposal needs for various demilitarized components disposal needs, such as: (1) pyrotechnic smoke assemblies, (2) thermal batteries, (3) proximity fuses, (4) cartridge actuated devices (CADs), and (5) propellant actuated devices (PADs). MSE Technology Applications, Inc., (MSE) has proposed and fabricated a Mobile Plasma Treatment System to be a technology demonstrator for pilotscale mobile plasma waste processing. The system is capable of providing small-scale waste remediation services, and conducting waste stream applicability demonstrations. The Mobile Plasma Treatment System's innovative concept provides the flexibility to treat waste streams at numerous sites and sites with only a limited quantity of waste, yet too hazardous to transport to a regional fixed facility. The system was designed to be operated as skid mounted modules; consisting of a furnace module, controls module, offgas module, and ancillary systems module. All system components have been integrated to be operated from a single control station with both semi-continuous feeding and batch slag-pouring capability.

Smith, Ed↗

Demonstration of Plasma Arc Environmental Technology Applications for the Demilitarization of DOD Stockpiles

Since 1989 the U.S. Army Construction Engineering Research Laboratories (USACERL) have been active participants in the research and development toward establishing Plasma Arc Technology (PAT) as an efficient, economical, and safe hazardous waste immobilization tool. A plasma torch capable of generating high temperatures makes this technology a viable and powerful tool for the thermal destruction of various military industrial waste streams into an innocuous ceramic material no longer requiring hazardous waste landfill (Class 1) disposal. The emerging pl asma environmental thermal treatment process, has been used to safely and efficiently meet the waste disposal needs for various demilitarized components disposal needs, such as: pyrotechnic smoke assemblies, thermal batteries, proximity fuses, cartridge actuated devices (CAD's), and propellant actuated devices (PAD's). MSE Technology Applications, Inc., (MSE) has proposed and fabricated a Mobile Plasma Treatment System to be a technology demonstrator for pilot-scale mobile plasma waste processing. The system is capable of providing small-scale waste remediation services, and conducting waste stream applicability demonstrations. The Mobile Plasma Treatment System's innovative concept provides the flexibility to treat waste streams at numerous sites and sites with only a limited quantity of waste, yet too hazardous to transport to a regional fixed facility. The system was designed to be operated as skid mounted modules; consisting of a furnace module, controls module, offgas module, and ancillary systems module. All system components have been integrated to be operated from a single control station with both semi-continuous feeding and batch slag-pouring capability.

Smith, Ed↗

Assimilation of MLS and OMI Ozone Data

Ozone data from Aura Microwave Limb Sounder (MLS) and Ozone Monitoring Instrument (OMI) were assimilated into the ozone model at NASA's Global Modeling and Assimilation Office (GMAO). This assimilation produces ozone fields that are superior to those from the operational GMAO assimilation of Solar Backscatter Ultraviolet (SBUV/2) instrument data. Assimilation of Aura data improves the representation of the "ozone hole" and the agreement with independent Stratospheric Aerosol and Gas Experiment (SAGE) III and ozone sonde data. Ozone in the lower stratosphere is captured better: mean state, vertical gradients, spatial and temporal variability are all improved. Inclusion of OMI and MLS data together, or separately, in the assimilation system provides a way of checking how consistent OMI and MLS data are with each other, and with the ozone model. We found that differences between OMI total ozone column data and model forecasts decrease after MLS data are assimilated. This indicates that MLS stratospheric ozone profiles are consistent with OMI total ozone columns. The evaluation of error characteristics of OMI and MLS ozone will continue as data from newer versions of retrievals becomes available. We report on the initial step in obtaining global assimilated ozone fields that combine measurements from different Aura instruments, the ozone model at the GMAO, and their respective error characteristics. We plan to use assimilated ozone fields in estimation of tropospheric ozone. We also plan to investigate impacts of assimilated ozone fields on numerical weather prediction through their use in radiative models and in the assimilation of infrared nadir radiance data from NASA's Advanced Infrared Sounder (AIRS).

Stajner, I.↗

Using In-Situ Deposition of Metallic Thin Films on Mars to Monitor Atmospheric H20

We describe the use of a novel sensor, designed to characterize the reactive nature of the Martian atmosphere, to also characterize the instantaneous abundance of atmospheric H,O. The sensor deposits, in situ, a thin silver film onto a sapphire substrate and monitors oxidation by measuring resistance of the Ag film during both deposition and subsequent oxidation [1,2]. . The evaporation source is placed in the center of a hollow, open-ended tube. On flash heating the source, evaporated metallic silver rapidly reacts with oxidizing gases in the tube, depositing on the tube walls as a resisitive Ag oxide film. Unoxidized silver deposits over the Ag oxides, once the oxidizing gases in the vicinity of the source have been consumed. Since the ends of the tube are open, atmospheric gases in the tube cause a time delay as the evaporated silver initially reacts with the available O2, and H2O; once oxidizing gases in the tube are significantly depleted, a silver film closes the chemiresistor circuit. In principle this operation mode provides a repeatable measure of the variable H2O abundance, since O2, levels are constant in the atmosphere. To explore the utility of the instrument as an H2O sensor, predict its behavior, and identify possible failure modes, we developed a numerical model of the instrument, and exercised it for Mars conditions. The model predicts the behavior of the electrical circuit, the temperature dependence of resistivity for each component, and the resultant Joule heating. The model balances Joule heating against radiation, sensible heat loss to ambient CO2, and, latent heat loss from sublimating Ag. The flux of Ag atoms from the source is tracked continuously. In the gas phase, the number density of Ag, O2, and H2O is calculated, along with their reaction rates on collision. The model suggests that Ag substantially depletes H2O only near the center of the tube. The model supports the hypothesis that the onset of high conductivity between the A u electrodes is a function of the H2O abundance in the ambient gas, but only to a limit; when H2O abundances are buffered at or above 240K (an H2O abundance that would not have occurred on Mars), metallic Ag never deposits in excess of Ago, and the circuit does not close. A significant difference in the Mars simulation is that the Ag source sublimates more quickly, primarily because of differences in sensible heat fluxes from the source. This allows subsequent heating of the W filament used to support the Ag source. Had the originally-designed experimental sequence been carried out on Mars, the filament might have melted, making subsequent investigations impossible.

Zent, A. P.↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗

Data-Driven Adaptive Damping Controller for Wind Power Plants with Doubly-Fed Induction Generators: Preprint

This paper presents an adaptive damping controller for wind power plants in which the turbines are equipped with doubly-fed induction generators. The controller is designed to respond to an input control signal that is triggered according to the system operating conditions. A processing unit continuously estimates the electromechanical modes of oscillation based on real-time streaming data acquired from a phasor measurement unit that is strategically positioned on the grid. The decision to trigger (or not trigger) the control signal is automatic, based on the relative damping of the dominant mode. The modes are estimated using the dynamic mode decomposition algorithm with time-delay embedding. Numerical simulations performed on the two-area system demonstrate that the proposed controller enhances the rotor angle stability for both small-signal and large disturbances, and is adaptive to changing grid conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗