Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 991 records · Page 55

CFD Study of Turbo-Ramjet Interactions in Hypersonic Airbreathing Propulsion System

Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.

Chang, Ing↗

Application of fuzzy logic-neural network based reinforcement learning to proximity and docking operations: Special approach/docking testcase results

As part of the RICIS project, the reinforcement learning techniques developed at Ames Research Center are being applied to proximity and docking operations using the Shuttle and Solar Maximum Mission (SMM) satellite simulation. In utilizing these fuzzy learning techniques, we use the Approximate Reasoning based Intelligent Control (ARIC) architecture, and so we use these two terms interchangeably to imply the same. This activity is carried out in the Software Technology Laboratory utilizing the Orbital Operations Simulator (OOS) and programming/testing support from other contractor personnel. This report is the final deliverable D4 in our milestones and project activity. It provides the test results for the special testcase of approach/docking scenario for the shuttle and SMM satellite. Based on our experience and analysis with the attitude and translational controllers, we have modified the basic configuration of the reinforcement learning algorithm in ARIC. The shuttle translational controller and its implementation in ARIC is described in our deliverable D3. In order to simulate the final approach and docking operations, we have set-up this special testcase as described in section 2. The ARIC performance results for these operations are discussed in section 3 and conclusions are provided in section 4 along with the summary for the project.

Jani, Yashvant↗

Preliminary Investigation of a Ceramic Lining for a Combustion Chamber for Gas-turbine Use

Combustion chamber liners for gas turbines and experimental set-ups were tested for failure. A ceramic-lined test chamber was operated at fuel/air ratios up to 0.050. Thermal-shock evaluation indicated that a ceramic lining, which expands after firing, would crack but would not fall apart during operation. Refractoriness of the lining and the resistance to mechanical shock were adequate. In general, shell temperature reductions of approximately 400 deg F were effected by the use of this lining at fuel/air ratios of 0.016 and 0.050. The mechanism of failure of the ceramic lining was induced by sudden heating and cooling during operation.

Woodward, William H↗

Stability and Phase Noise Tests of Two Cryo-Cooled Sapphire Oscillators

A cryocooled Compensated Sapphire Oscillator (CSO), developed for the Cassini Ka-band Radio Science experiment, and operating in the 8K - 10K temperature range was previously demonstrated to show ultra-high stability of sigma(sub y) = 2.5 x 10 (exp -15) for measuring times 200 seconds less than or equal to tau less than or equal to 600 seconds using a hydrogen maser as reference. We present here test results for a second unit which allows CSO short-term stability and phase noise to be measured for the first time. Also included are design details of a new RF receiver and an intercomparison with the first CSO unit. Cryogenic oscillators operating below about 10K offer the highest possible short term stability of any frequency sources. However, their use has so far been restricted to research environments due to the limited operating periods associated with liquid helium consumption. The cryocooled CSO is being built in support of the Cassini Ka-band Radio Science experiment and is designed to operate continuously for periods of a year or more. Performance targets are a stability of 3-4 x 10 (exp -15) (1 second less than or equal to tau less than or equal to 100 seconds) and phase noise of -73dB/Hz @ 1Hz measured at 34 GHz. Installation in 5 stations of NASA's deep space network (DSN) is planned in the years 2000 - 2002. In the previous tests, actual stability of the CSO for measuring times tau less than or equal to 200 seconds could not be directly measured, being masked by short-term fluctuations of the H-maser reference. Excellent short-term performance, however, could be inferred by the success of an application of the CSO as local oscillator (L.O.) to the JPL LITS passive atomic standard, where medium-term stability showed no degradation due to L.O. instabilities at a level of (sigma)y = 3 x 10 (exp -14)/square root of tau. A second CSO has now been constructed, and all cryogenic aspects have been verified, including a resonator turn-over temperature of 7.907 K, and Q of 7.4 x 10 (exp 8). These values compare to a turn-over of 8.821 K and Q of 1.0 x 10 (exp 9) for the first resonator. Operation of this second unit provides a capability to directly verify for the first time the short-term (1 second less than or equal to tau less than or equal to 200 seconds) stability and the phase noise of the CSO units. The RF receiver used in earlier tests was sufficient to meet Cassini requirements for tau greater than or equal to 10 seconds but had short-term stability limited to 2-4 x 10 (exp -14) at tau = 1 second, a value 10 times too high to meet our requirements. A new low-noise receiver has been designed to provide approximately equal to 10-15 performance at 1 second, and one receiver is now operational, demonstrating again short-term CSO performance with H maser-limited stability. Short-term performance was degraded in the old receiver due to insufficient tuning bandwidth in a 100MHZ quartz VCO that was frequency-locked to the cryogenic sapphire resonator. The new receivers are designed for sufficient bandwidth, loop gain and low noise to achieve the required performance.

Dick, G. John↗

Performance of an Optimized Eta Model Code on the Cray T3E and a Network of PCs

In the year 2001, NASA will launch the satellite TRIANA that will be the first Earth observing mission to provide a continuous, full disk view of the sunlit Earth. As a part of the HPCC Program at NASA GSFC, we have started a project whose objectives are to develop and implement a 3D cloud data assimilation system, by combining TRIANA measurements with model simulation, and to produce accurate statistics of global cloud coverage as an important element of the Earth's climate. For simulation of the atmosphere within this project we are using the NCEP/NOAA operational Eta model. In order to compare TRIANA and the Eta model data on approximately the same grid without significant downscaling, the Eta model will be integrated at a resolution of about 15 km. The integration domain (from -70 to +70 deg in latitude and 150 deg in longitude) will cover most of the sunlit Earth disc and will continuously rotate around the globe following TRIANA. The cloud data assimilation is supposed to run and produce 3D clouds on a near real-time basis. Such a numerical setup and integration design is very ambitious and computationally demanding. Thus, though the Eta model code has been very carefully developed and its computational efficiency has been systematically polished during the years of operational implementation at NCEP, the current MPI version may still have problems with memory and efficiency for the TRIANA simulations. Within this work, we optimize a parallel version of the Eta model code on a Cray T3E and a network of PCs (theHIVE) in order to improve its overall efficiency. Our optimization procedure consists of introducing dynamically allocated arrays to reduce the size of static memory, and optimizing on a single processor by splitting loops to limit the number of streams. All the presented results are derived using an integration domain centered at the equator, with a size of 60 x 60 deg, and with horizontal resolutions of 1/2 and 1/3 deg, respectively. In accompanying charts we report the elapsed time, the speedup and the Mflops as a function of the number of processors for the non-optimized version of the code on the T3E and theHIVE. The large amount of communication required for model integration explains its poor performance on theHIVE. Our initial implementation of the dynamic memory allocation has contributed to about 12% reduction of memory but has introduced a 3% overhead in computing time. This overhead was removed by performing loop splitting in some of the high demanding subroutines. When the Eta code is fully optimized in order to meet the memory requirement for TRIANA simulations, a non-negligeable overhead may appear that may seriously affect the efficiency of the code. To alleviate this problem, we are considering implementation of a new algorithm for the horizontal advection that is computationally less expensive, and also a new approach for marching in time.

Kouatchou, Jules↗

Observations of Small Debris from the Cosmos 1408 Anti-Satellite Test using the HUSIR and Goldstone Radars

On 15 November 2021, Russia conducted a direct-ascent anti-satellite (ASAT) test against its Cosmos 1408 satellite, which had been in orbit since 1982. The test produced at least 1500 fragments trackable by the U.S. Space Surveillance Network (SSN). The test is a significant event because the resulting cloud has the potential to endanger the International Space Station and other satellites in low Earth orbit (LEO). For almost 30 years, the NASA Orbital Debris Program Office (ODPO) has used the Haystack Ultrawideband Satellite Imaging Radar (HUSIR), operated by the Massachusetts Institute of Technology’s Lincoln Laboratory to perform statistical measurements of debris in LEO too small to be tracked by the SSN, nominally down to 5.5 mm at 1000 km altitude. The ODPO also utilizes the Goldstone Orbital Debris Radar (Goldstone), operated by NASA’s Jet Propulsion Laboratory, to characterize the small debris environment in LEO down to approximately 3 mm at 1000 km altitude. To characterize the small debris component of the Cosmos 1408 ASAT test, a series of observation campaigns were conducted with the HUSIR and Goldstone radars. This paper discusses the observation planning, including beam overlap analysis needed for Goldstone bistatic operation, using a model of the debris cloud produced with the NASA Standard Satellite Breakup Model. A description of the radars and the data processing techniques used to analyze the data are also discussed. Finally, results of the measurement campaigns including cumulative count rate versus size and surface area flux versus altitude and inclination are presented

James Murray↗

Oxygen-Methane Thruster

An oxygen-methane thruster was conceived with integrated igniter/injector capable of nominal operation on either gaseous or liquid propellants. The thruster was designed to develop 100 lbf (approximately 445 N) thrust at vacuum conditions and use oxygen and methane as propellants. This continued development included refining the design of the thruster to minimize part count and manufacturing difficulties/cost, refining the modeling tools and capabilities that support system design and analysis, demonstrating the performance of the igniter and full thruster assembly with both gaseous and liquid propellants, and acquiring data from this testing in order to verify the design and operational parameters of the thruster. Thruster testing was conducted with gaseous propellants used for the igniter and thruster. The thruster was demonstrated to work with all types of propellant conditions, and provided the desired performance. Both the thruster and igniter were tested, as well as gaseous propellants, and found to provide the desired performance using the various propellant conditions. The engine also served as an injector testbed for MSFC-designed refractory combustion chambers made of rhenium.

Pickens, Tim↗

Multi-User Space Link Extension (SLE) System

The Multi-User Space (MUS) Link Extension system, a software and data system, provides Space Link Extension (SLE) users with three space data transfer services in timely, complete, and offline modes as applicable according to standards defined by the Consultative Committee for Space Data Systems (CCSDS). MUS radically reduces the schedule, cost, and risk of implementing a new SLE user system, minimizes operating costs with a lights-out approach to SLE, and is designed to require no sustaining engineering expense during its lifetime unless changes in the CCSDS SLE standards, combined with new provider implementations, force changes. No software modification to MUS needs to be made to support a new mission. Any systems engineer with Linux experience can begin testing SLE user service instances with MUS starting from a personal computer (PC) within five days. For flight operators, MUS provides a familiar-looking Web page for entering SLE configuration data received from SLE. Operators can also use the Web page to back up a space mission's entire set of up to approximately 500 SLE service instances in less than five seconds, or to restore or transfer from another system the same amount of data from a MUS backup file in about the same amount of time. Missions operate each MUS SLE service instance independently by sending it MUS directives, which are legible, plain ASCII strings. MUS directives are usually (but not necessarily) sent through a TCP-IP (Transmission Control Protocol Internet Protocol) socket from a MOC (Mission Operations Center) or POCC (Payload Operations Control Center) system, under scripted control, during "lights-out" spacecraft operation. MUS permits the flight operations team to configure independently each of its data interfaces; not only commands and telemetry, but also MUS status messages to the MOC. Interfaces can use single- or multiple-client TCP/IP server sockets, TCP/IP client sockets, temporary disk files, the system log, or standard in, standard out, or standard error as applicable. By defining MUS templates in ASCII, the flight operations team can include any MUS system variable in telemetry or command headers or footers, and/or in status messages. Data fields can be arranged within messages in different sequences, according to the mission s needs. The only constraints imposed are on the format of MUS directive strings, and some bare minimum logical requirements that must be met in order for MUS to read the mission control center's spacecraft command inputs. The MUS system imposes no limits or constraints on the numbers and combinations of missions and SLE service instances that it will support simultaneously. At any time, flight operators may add, change, delete, bind, connect, or disconnect.

Perkins, Toby↗

Cubesat Proximity Operations Demonstration (CPOD)

The CubeSat Proximity Operations Demonstration (CPOD) project will demonstrate rendezvous, proximity operations and docking (RPOD) using two 3-unit (3U) CubeSats. Each CubeSat is a satellite with the dimensions 4 inches x 4 inches x 13 inches (10 centimeters x 10 centimeters x 33 centimeters) and weighing approximately 11 pounds (5 kilograms). This flight demonstration will validate and characterize many new miniature low-power proximity operations technologies applicable to future missions. This mission will advance the state of the art in nanosatellite attitude determination,navigation and control systems, in addition to demonstrating relative navigation capabilities.The two CPOD satellites are scheduled to be launched together to low-Earth orbit no earlier than Dec. 1, 2015.

Relative navigation capabilities↗

Die Attachment for -120 C to +20 C Thermal Cycling of Microelectronics for Future Mars Rovers: An Overview

Active thermal control for electronics on Mars Rovers imposes a serious penalty in weight, volume, power consumption, and reliability. Thus, we propose that thermal control be eliminated for future Rovers. From a functional standpoint there is no reason that the electronics could not operate over the entire temperature range of the Martian environment, which can vary from a low of approximately equal -90 C to a high of approximately equal +20 C during the Martian night and day. The upper end of this range is well within that for conventional electronics. Although the lower end is considerably below that for which conventional--even high-reliability electronics is designed or tested, it is well established that electronic devices can operate to such low temperatures. The primary concern is reliability of the overall electronic system, especially in regard to the numerous daily temperature cycles that it would experience over the duration of a mission on Mars. Accordingly, key reliability issues have been identified for elimination of thermal control on future Mars Rovers. One of these is attachment of semiconductor die onto substrates and into packages. Die attachment is critical since it forms a mechanical, thermal and electrical interface between the electronic device and the substrate or package. This paper summarizes our initial investigation of existing information related to this issue, in order to form an opinion whether die attachment techniques exist, or could be developed with reasonable effort, to withstand the Mars thermal environment for a mission duration of approximately I year. Our conclusion, from a review of literature and personal contacts. is that die attachment can be made sufficiently reliable to satisfy the requirements of future Mars Rovers. Moreover, it appears that there are several possible techniques from which to choose and that the requirements could be met by judicious selection from existing methods using hard solders, soft solders, or organic adhesives. Thus from the standpoint of die attachment. it appears feasible to eliminate thermal control for Rover electronics. We recommend that this be further investigated and verified for the specific hardware and thermal conditions appropriate to Mars Rovers.

Kirschman, Randall K.↗

Cryogenic refractor design techniques

The Infrared Astronomical Satellite (IRAS) was designed to operate at 2K, and over the spectral range of 8 to 120 micrometers. The focal plane is approximately 2 by 3 inches in size, and contains 62 individual field stop apertures, each with its own field lens, one or more filters and a detector. The design of the lenses involved a number of difficulties and challenges that are not usually encountered in optical design. Operating temperature is assumed during the design phase, which requires reliable information on dN/dT (Index Coefficient) for the materials. The optics and all supporting structures are then expanded to room temperature, which requires expansion coefficient data on the various materials, and meticulous attention to detail. The small size and dense packaging, as well as the high precision required, further contributed to the magnitude of the task.

Darnell, R. J.↗

Example MODIS Global Cloud Optical and Microphysical Properties: Comparisons between Terra and Aqua

MODIS observations from the NASA EOS Terra spacecraft (launched in December 1999, 1030 local time equatorial crossing) have provided a unique data set of Earth observations. With the launch of the NASA Aqua spacecraft in May 2002 (1330 local time), two MODIS daytime (sunlit) and nighttime observations are now available in a 24 hour period, allowing for some measure of diurnal variability. We report on an initial analysis of several operational global (Level-3) cloud products from the two platforms. The MODIS atmosphere Level-3 products, which include clear-sky and aerosol products in addition to cloud products, are available as three separate files providing daily, eight-day, and monthly aggregations; each temporal aggregation is spatially aggregated to a 1 degree grid. The files contain approximately 600 statisitical datasets (from simple means and standard deviations to 1 - and 2-dimensional histograms). Operational cloud products include detection (cloud fraction), cloud-top properties, and daytimeonly cloud optical thickness and particle effective radius for both water and ice clouds. We will compare example global Terra and Aqua cloud fraction, optical thickness, and effective radius aggregations.

Hubanks, P. A.↗

Application of fuzzy logic-neural network based reinforcement learning to proximity and docking operations: Translational controller results

The reinforcement learning techniques developed at Ames Research Center are being applied to proximity and docking operations using the Shuttle and Solar Maximum Mission (SMM) satellite simulation. In utilizing these fuzzy learning techniques, we also use the Approximate Reasoning based Intelligent Control (ARIC) architecture, and so we use two terms interchangeable to imply the same. This activity is carried out in the Software Technology Laboratory utilizing the Orbital Operations Simulator (OOS). This report is the deliverable D3 in our project activity and provides the test results of the fuzzy learning translational controller. This report is organized in six sections. Based on our experience and analysis with the attitude controller, we have modified the basic configuration of the reinforcement learning algorithm in ARIC as described in section 2. The shuttle translational controller and its implementation in fuzzy learning architecture is described in section 3. Two test cases that we have performed are described in section 4. Our results and conclusions are discussed in section 5, and section 6 provides future plans and summary for the project.

Jani, Yashvant↗

Observations of Small Debris from the Cosmos 1408 Anti-Satellite Test using the HUSIR and Goldstone Radars

On 15 November 2021, Russia conducted a direct-ascent, anti-satellite (ASAT) test against its Cosmos 1408 satellite, which had been in orbit since 1982. The test produced at least 1500 fragments trackable by the U.S. Space Surveillance Network (SSN). This is a significant event because the resulting cloud has the potential to endanger the International Space Station and other satellites in low Earth orbit (LEO).For almost 30 years, the NASA Orbital Debris Program Office (ODPO) has been using the Haystack Ultrawideband Imaging Radar (HUSIR), located in Tyngsborough, Massachusetts and operated by the Massachusetts Institute of Technology’s Lincoln Laboratory, to perform statistical measurements of debris in LEO too small to be tracked by the SSN, nominally down to 5.5 mm at 1000 km altitude. The ODPO also utilizes the Goldstone Orbital Debris Radar (Goldstone), located near Barstow, California and operated by the NASA Jet Propulsion Laboratory, to characterize the small debris environment in LEO down to approximately 3 mm at 1000 km altitude. To characterize the small debris component of the Cosmos 1408 ASAT test, a series of observation campaigns were conducted with the HUSIR and Goldstone radars. This paper discusses the observation planning, including beam overlap analysis required for the Goldstone’s bistatic operation, using a model of the debris cloud produced using the NASA Standard Satellite Breakup Model. A description of the radars and the data processing techniques used to analyze the data are also discussed. Finally, results of the measurement campaigns including cumulative flux versus size, altitude, and inclination are presented.

James Murray↗

Conventional and Bimodal Nuclear Thermal Rocket (NTR) Artificial Gravity Mars Transfer Vehicle Concepts

A variety of countermeasures have been developed to address the debilitating physiological effects of "zero-gravity" (0-g) experienced by cosmonauts and astronauts during their approximately 0.5-1.2 year long stays in LEO (Low Earth Orbit). Longer interplanetary flights, combined with possible prolonged stays in Mars orbit, could subject crewmembers to up to approximately 2.5 years of weightlessness. In view of known and recently diagnosed problems associated with 0-g, an artificial gravity spacecraft offers many advantages and may indeed be an enabling technology for human flights to Mars. A number of important human factors must be taken into account in selecting the rotation radius, rotation rate, and orientation of the habitation module or modules. These factors include the gravity gradient effect, radial and tangential Coriolis forces, along with cross-coupled acceleration effects. Artificial gravity (AG) Mars transfer vehicle (MTV) concepts are presented that utilize both conventional NTR, as well as, enhanced "bimodal" nuclear thermal rocket (BNTR) propulsion. The NTR is a proven technology that generates high thrust and has a specific impulse (I (sub sp)) capability of approximately 900 s - twice that of today's best chemical rockets. The AG/MTV concepts using conventional NTP carry twin cylindrical "ISS-type" habitation modules with their long axes oriented either perpendicular or parallel to the longitudinal spin axis of the MTV and utilize photovoltaic arrays (PVAs) for spacecraft power. The twin habitat modules are connected to a central operations hub located at the front of the MTV via two pressurized tunnels that provide the rotation radius for the habitat modules. For the BNTR AG/MTV option, each engine has its own "closed" secondary helium-xenon gas loop and Brayton rotating unit that can generate tens of kilowatts (kW (sub e)) of spacecraft electrical power during the mission coast phase eliminating the need for large PVAs. A single inflatable "TransHab-type" habitation module is also used with multiple vertical floors oriented radial to the MTV spin axis. The BNTR MTV's geometry - long and linear - is naturally compatible with AG operation. By rotating the vehicle about its center-of-mass and perpendicular to its flight vector at approximately 3.0 - 5.2 rpm, a centrifugal force and AG environment corresponding to approximately 0.38 - 1.0 g can be established to help maintain crew fitness out to Mars and back. Vehicles using NTP/BNTP can more readily accommodate the heavier payload mass and increased RCS propellant loading associated with AG operation, and can travel faster to and from Mars thereby reducing the crew's exposure to galactic cosmic radiation and solar flares. Mission scenario descriptions, key vehicle features and operational characteristics for each propulsion options are presented using the lift capability and payload volumes estimated for the SLS-1A and HLV.

Nuclear Propulsion↗

Low voltage 30cm ion thruster

The construction of an ion thruster module (including thruster, power conditioning, and control system) capable of operating for 10,000 hours over a five to one range at an effective specific impulse of approximately 2800 seconds is discussed. The several interrelated tasks involved in the construction of the engine are described. Performance tests of the engine and the effects of various modifications are reported. It was demonstrated that thruster performance and stability were not materially affected by reasonable changes from the nominal operating point.

Source record↗

Adiabatic Shock Capturing in Perfect Gas Hypersonic Flows

This paper considers the streamline-upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier-Stokes equations in conservation-variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock capturing operator, which is required to regularize the scheme in the presence of strong, shock-induced gradients. A standard shock capturing operator which has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator which are designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator.

Kirk, Benjamin S.↗

The 4th order GISS model of the global atmosphere

The new GISS 4th order model of the global atmosphere is described. It is based on 4th order quadratically conservative differences with the periodic application of a 16th order filter on the sea level pressure and potential temperature equations, a combination which is approximately enstrophy conserving. Several short range forecasts indicate a significant improvement over 2nd order forecasts with the same resolution (approximately 400 km). However the 4th order forecasts are somewhat inferior to 2nd order forecasts with double resolution. This is probably due to the presence of short waves in the range between 1000 km and 2000 km, which are computed more accurately by the 2nd order high resolution model. An operation count of the schemes indicates that with similar code optimization, the 4th order model will require approximately the same amount of computer time as the 2nd order model with the same resolution. It is estimated that the 4th order model with a grid size of 200 km provides enough accuracy to make horizontal truncation errors negligible over a period of a week for all synoptic scales (waves longer than 1000 km).

Kalnay-Rivas, E.↗