Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAXED”

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 163 records · Page 9

Additive evaluation criteria for aircraft noise

The occurrence of unexpected aircraft noise events will frequently evoke intense complaints about annoyance over such events. It is recognized that the relationship between the volume of complaints and the corresponding maximum noise levels does in fact depend on the circumstances of the complainants and the time of year. The frequency of occurrence of the respective noise events is also a factor. The possible practical value of the addition of the maximal noise level, to the well-known cumulative noise descriptors L(sub eq), L(sub dn), etc. may be considered. One might start by considering the difference between the L(sub eq) and the average maximal noise level of the twenty loudest single noise events on an average day. A somewhat less sharply focused consideration of the maximum noise levels was adopted. In that standard the scope of noise-mitigation measures is defined generally with reference to L(sub eq). If the average maximal noise level L(sub max) of the entire aircraft fleet mix exceeds the L(sub eq) by more than 20 dB(A) and if more than 20 daily aircraft noise events exceed the L(sub eq) than the difference L(sub max)-20 becomes the key criterion for noise-mitigation measures. A statistical correlation of a large number of data from aircraft-noise-monitoring sensors located both directly underneath and laterally disposed to an aircraft flight-path has supplied a basis for the determination of the distribution of the maximal noise levels about the average value, L(sub max), of each type of aircraft. This distribution is given for the takeoff climb and for the landing approach.

Meyer, Thomas J.↗

Effects of simulated weightlessness and sympathectomy on maximum VO2 of male rats

The effects of simulated weightlessness (hind-limb suspension) and chemical sympathectomy (by repeated injections with guanethidine sulfate) on the maximum oxygen consumption (VO2 max) of female rats were investigated in rats assigned for 14 days to one of three groups: a head-down hind-limb suspension, a horizontal suspension with hind limbs weight bearing, or the caged control. The VO2 max values were assessed by having rats run on a treadmill enclosed in an airtight chamber. The hind-limb-suspended sympathectomized rats were found to exhibit shorter run times and lower mechanical efficiencies, compared to their presuspension values or the values from saline-injected suspended controls. On the other hand, the suspended sympathectomized rats did not demonstrate a decrease in the VO2 max values that was observed in saline-injected controls.

Woodman, C. R.↗

Range optimization for a supersonic aircraft

Range optimal trajectories for an aircraft flying in the vertical plane are obtained from Pontryagin's Minimum Principle. Control variables are load factor n which appears nonlinearly in the equations of motion and throttle setting eta, which appears only linearly. Both controls are subject to fixed bounds, namely eta between values of 0 and 1 and absolute value of n not greater than n(max). Additionally, a dynamic pressure limit is imposed, which represents a first-order state-inequality constraint. For fixed flight time, fixed initial coordinates, and partially fixed final coordinates, the effect of the load factor limit absolute value of n not greater than n(max) is studied. Upon varying n(max), six different switching structures are obtained. All trajectories involve singular control along arcs with active dynamic pressure limit.

Seywald, Hans↗

A statistical analysis of gamma-ray bursts detected by the Konus Experiment on Venera 11 and 12

We discuss the advantages of using the V/V sub max method to test well-defined samples of gamma-ray bursts for the spatial uniformity of their parent population. We have applied the V/V sub max test to gamma ray bursts of duration longer than 1 second recorded by the Konus experiment aboard Venera 11 and 12. Based on a sample of 123 bursts, we find V/V sub max equals 0.46 plus or minus 0.03, consistent with a uniform distribution in space. We urge that experimenters give careful attention to the detection limit for each recorded gamma ray burst, and that quantitative data for burst properties and detection limits be published.

Schmidt, Maarten↗

Instantaneous heat transfer coefficient based upon two-dimensional analyses of Stirling space engine components

Twelve different cases of multidimensional models of Stirling engine components for space applications have been numerically investigated for oscillating, incompressible laminar flow with heat transfer. The cases studied covered wide ranges of Valensi number (from 44 to 700), Re(max) number (from 8250 to 60,000), and relative amplitude of fluid motion of 0.686 and 1.32. The Nusselt numbers obtained from the present study indicate a very complex shape with respect to time and axial location in the channel. The results indicate that three parameters can be used to define the local Nusselt number variation, namely: time average, amplitude, and phase angle. These parameters could be correlated respectively using: Re(max), Va and Re(max), and the relative amplitude of fluid motion.

Ibrahim, Mounir↗

A study of subterahertz HEMT monolithic oscillators

A detailed study of monolithic InP-based HEMT oscillators for subterahertz operation is presented. InAlAs/InGaAs HEMT's have been optimized for high frequency operation and showed very high maximum oscillation frequencies (f(sub max)) of 310 GHz using offset self-aligned gamma-gate technology. Power characteristics of HEMT oscillators are reported. An oscillation power of more than 10 mW was evaluated by large-signal analysis at 320 GHz using HEMT's with f(sub max) = 450 GHz, V(sub br) = 10 V and a gate width (W(sub g)) of 8 x 22.5 microns. Oscillator topology studies showed that complex feedback schemes such as dual and active feedback enhance the negative resistance. Push-push oscillator designs based on harmonic signal generation can finally be used to overcome the frequency barrier imposed by f(sub max).

Kwon, Youngwoo↗

Investigation at supersonic speeds of the variation with Reynolds number and Mach number of the total, base, and skin-friction drag of seven boattail bodies of revolution designed for minimum wave drag

An investigation has been conducted in the Langley 9-inch supersonic tunnel to determine the effect of the variation with Reynolds number and Mach number of the total, base, and skin-friction drag at zero lift of seven boattail bodies of revolution designed for minimum wave drag according to NACA TN 2550. The tests covered a Reynolds number range from approximately 1.0 x 10(6) to 10.0 x 10(6) at Mach numbers of 1.62, 1.93, and 2.41, respectively. The results show that base drag and, in general, the total drag increase with increasing values of the ratio of base area to maximum area B/S(sub max), although the results reported in NACA TN 3054 showed that the wave drag decreased with increasing values of B/S(sub max). The laminar skin-friction drag is in agreement with the theoretical predictions used, and, within the Mach number range of these tests, the simple Blasius incompressible theory gives a satisfactory prediction. Except for values of B/S (sub max) near 1, the Reynolds number of transition increases with increasing Mach number and as this ratio approaches 1 this variation is seen to reverse. These variations in Reynolds number of transition with Mach number appear to be associated with the changes in pressure gradient over the rear of the bodies.

Bromm, August F , Jr↗

Surface Fatigue Life of High Temperature Gear Materials

Three high temperature gear materials were evaluated using spur gear surface fatigue tests. These materials were, VASCO max 350, VASCO matrix 2, and nitralloy N and were evaluated for possible use in high temperature gear applications. The fatigue life of the three high temperature gear materials were compared with the life of the standard AISI 9310 aircraft gear material. Surface fatigue tests were conducted at a lubricant inlet temperature of 321 K (120 F), a lubricant outlet temperature of 350 K (170 F), a maximum Hertz stress of 1.71 GPa (248 ksi), a speed of 10,000 rpm, and with a synthetic paraffinic lubricant. The life of the nitralloy N was approximately the same as the AISI 9310, the life of the VASCO max 350 was much less than the AISI 9310 while the life of the VASCO matrix 2 was several times the life of the AISI 9310. The VASCO max 350 also showed very low fracture toughness with approximately half of the gears failed by tooth fracture through the fatigue spall. The VASCO matrix 2 had approximately 10-percent fracture failure through the fatigue spalls indicating moderate to good fracture toughness.

Townsend, Dennis P.↗

The ground vehicle manager's associate

An overview of MAX, a software framework for manager's associate systems, is presented. MAX is used to develop and execute a problem-solving strategy for the task planning of semi-autonomous agents with the assistance of human performance. This paper describes the use of MAX in the supervisory management of robotic vehicles as they explore a planetary surface.

Edwards, Gary R.↗

An algorithm for solving the system-level problem in multilevel optimization

A multilevel optimization approach which is applicable to nonhierarchic coupled systems is presented. The approach includes a general treatment of design (or behavior) constraints and coupling constraints at the discipline level through the use of norms. Three different types of norms are examined: the max norm, the Kreisselmeier-Steinhauser (KS) norm, and the 1(sub p) norm. The max norm is recommended. The approach is demonstrated on a class of hub frame structures which simulate multidisciplinary systems. The max norm is shown to produce system-level constraint functions which are non-smooth. A cutting-plane algorithm is presented which adequately deals with the resulting corners in the constraint functions. The algorithm is tested on hub frames with increasing number of members (which simulate disciplines), and the results are summarized.

Balling, R. J.↗

Some requirements of a colliding comet source of gamma ray bursts

Colliding comets in the Solar System may be an important source of gamma ray bursts. The spherical gamma ray comet cloud required by the results of the Venera Satellites (Mazets and Golenetskii, 1987) and the Burst and Transient Source Experiment (BATSE) detector on the Compton Satellite (Meegan et al., 1992a, b) is neither the Oort Cloud nor the Kuiper Belt. To satisfy observations of N(greater than P(sub max)) vs P(sub max) for the maximum gamma ray fluxes, P(sub max) greater than 10(exp -5) ergs/sq cm/ s (about 30 bursts/yr), the comet density, n, should increase as n approximately a(exp 1) from about 40 to 100 AU where a is the comet heliocentric distance. The turnover above 100 AU requires n approximately a(exp -1/2) to 200 AU to fit the Venera results and n approximately a(exp 1/4) to 400 AU to fit the BATSE data. Then the masses of comets in the 3 regions are from: 40-100 AU, about 9 earth masses, m(sub E); 100-200 AU about 25 m(sub E); and 100-400 AU, about 900 m(sub E). The flux of 10(exp -5) ergs/sq cm/s corresponds to a luminosity at 100 AU of 3 x 10(exp 26) ergs/s. Two colliding spherical comets at a distance of 100 AU, each with nucleus of radius R of 5 km, density of 0.5 g/cu cm and Keplerian velocity 3 km/s have a combined kinetic energy of 3 x 10(exp 28) erg, a factor of about 100 greater tha required by the burst maximum fluxes that last for one second.

White, R. Stephen↗

Investigation at Mach Numbers of 0.20 to 3.50 of Blended Wing-Body Combinations of Sonic Design with Diamond, Delta, and Arrow Plan Forms

The models had aspect-ratio-2 diamond, delta, and arrow wings with the leading edges swept 45.00 deg, 59.04 deg, and 70.82 deg, respectively. The wing sections were computed by varying the section shape along with the body radii (blending process) to match the prescribed area distribution and wing plan form. The wing sections had an average value of maximum thickness ratio of about 4 percent of the local chords in a streamwise direction. The models were tested with transition fixed at Reynolds numbers of about 4,000,000 to 9,000,0000, based on the mean aerodynamic chord of the wings. The effect of varying Reynolds number was checked at both subsonic and supersonic speeds. The diamond model was superior to the other plan forms at transonic speeds ((L/D)max = 11.00 to 9.52) because of its higher lift-curve slope and near optimum wave drag due to the blending process. For the wing thickness tested with the diamond model, the marked body and wing contouring required for transonic conditions resulted in a large wave-drag penalty at the higher supersonic Mach numbers where the leading and trailing edges of the wing were supersonic. Because of the low sweep of the trailing edge of the delta model, this configuration was less adaptable to the blending process. Removing a body bump prescribed by the Mach number 1.00 design resulted in a good supersonic design. This delta model with 10 percent less volume was superior to the other plan forms at Mach numbers of 1.55 to 2.35 ((L/D)max = 8.65 to 7.24), but it and the arrow model were equally good at Mach numbers of 2.50 to 3.50 ((L/D)max - 6.85 to O.39). At transonic speeds the arrow model was inferior because of the reduced lift-curve slope associated with its increased sweep and also because of the wing base drag. The wing base-drag coefficients of the arrow model based on the wing planform area decreased from a peak value of 0.0029 at Mach number 1.55 to 0.0003 at Mach number 3.50. Linear supersonic theory was satisfactory for predicting the aerodynamic trends at Mach numbers from 1.55 to 3.50 of lift-curve slope, wave drag, drag due to lift, aerodynamic-center location, and maximum lift-drag ratios for each of the models.

Holdaway, George H.↗

Analyses of Fatigue Crack Growth and Closure Near Threshold Conditions for Large-Crack Behavior

A plasticity-induced crack-closure model was used to study fatigue crack growth and closure in thin 2024-T3 aluminum alloy under constant-R and constant-K(sub max) threshold testing procedures. Two methods of calculating crack-opening stresses were compared. One method was based on a contact-K analyses and the other on crack-opening-displacement (COD) analyses. These methods gave nearly identical results under constant-amplitude loading, but under threshold simulations the contact-K analyses gave lower opening stresses than the contact COD method. Crack-growth predictions tend to support the use of contact-K analyses. Crack-growth simulations showed that remote closure can cause a rapid rise in opening stresses in the near threshold regime for low-constraint and high applied stress levels. Under low applied stress levels and high constraint, a rise in opening stresses was not observed near threshold conditions. But crack-tip-opening displacement (CTOD) were of the order of measured oxide thicknesses in the 2024 alloy under constant-R simulations. In contrast, under constant-K(sub max) testing the CTOD near threshold conditions were an order-of-magnitude larger than measured oxide thicknesses. Residual-plastic deformations under both constant-R and constant-K(sub max) threshold simulations were several times larger than the expected oxide thicknesses. Thus, residual-plastic deformations, in addition to oxide and roughness, play an integral part in threshold development.

Newman, J. C., Jr.↗

Global to Microscale Evolution of the Pinatubo Volcanic Aerosol Derived from Diverse Measurements and Analyses

We assemble data on the Pinatubo aerosol from space, air, and ground measurements, develop a composite picture, and assess the consistency and uncertainties of measurement and retrieval techniques. Satellite infrared spectroscopy, particle morphology, and evaporation temperature measurements agree with theoretical calculations in showing a dominant composition of H2SO4-H20 mixture, with H2SO4 weight fraction of 65-80% for most stratospheric temperatures and humidities. Important exceptions are (1) volcanic ash, present at all heights initially and just above the tropopause until at least March 1992, and (2) much smaller H2SO4 fractions at the low temperatures of high-latitude winters and the tropical tropopause. Laboratory spectroscopy and calculations yield wavelength- and temperature-dependent refractive indices for the H2SO4-H20 droplets. These permit derivation of particle size information from measured optical depth spectra, for comparison to impactor and optical-counter measurements. All three techniques paint a generally consistent picture of the evolution of R(sub eff), the effective radius. In the first month after the eruption, although particle numbers increased greatly, R(sub eff) outside the tropical core was similar to preeruption values of approx. 0.1 to 0.2 microns, because numbers of both small (r < 0.2 microns) and large (r > 0.6 microns) particles increased. In the next 3-6 months, extracore R(sub eff) increased to approx. 0.5 microns, reflecting particle growth through condensation and coagulation. Most data show that R(sub eff) continued to increase for about 1 year after the eruption. R(sub eff) values up to 0.6 - 0.8 microns or more are consistent with 0.38 - 1 micron optical depth spectra in middle to late 1992 and even later. However, in this period, values from in situ measurements are somewhat less. The difference might reflect in situ undersampling of the very few largest particles, insensitivity of optical depth spectra to the smallest particles, or the inability of flat spectra to place an upper limit on particle size. Optical depth spectra extending to wavelengths lambda > 1 micron are required to better constrain R(sub eff), especially for R(sub eff) > 0.4 microns. Extinction spectra computed from in situ size distributions are consistent with optical depth measurements; both show initial spectra with lambda(sub max) <= 0.42 microns, thereafter increasing to 0.78 <= lambda(sub max) <= 1 micron. Not until 1993 do spectra begin to show a clear return to the preeruption signature of lambda(sub max) <= 0.42 microns. The twin signatures of large R(sub eff) (> 0.3 microns) and relatively flat extinction spectra (0.4 - 1 microns) are among the longest-lived indicators of Pinatubo volcanic influence. They persist for years after the peaks in number, mass, surface area, and optical depth at all wavelengths <= 1 microns. This coupled evolution in particle size distribution and optical depth spectra helps explain the relationship between global maps of 0.5- and 1.0-micron optical depth derived from the Advanced Very High Resolution Radiometer (AVHRR) and Stratospheric Aerosol and Gas Experiment (SAGE) satellite sensors. However, there are important differences between the AVHRR and SAGE midvisible optical thickness products. We discuss possible reasons for these differences and how they might be resolved.

Russell, P. B.↗

Dependence of Photosynthetic Capacity, Photosynthetic Pigment Allocation, and Carbon Storage on Nitrogen Levels in Foliage of Aspen Stands

The role of foliar nitrogen (N) in the seasonal dynamics and vertical canopy distribution of photosynthetic pigments, photosynthetic capacity, and carbon (C) storage was investigated in boreal broadleaved species. The study was conducted at two different aged stands (60 y and 15 y) in 1994 and 1996 in Saskatchewan, Canada as part of the Boreal Ecosystem-Atmosphere Study (BOREAS). Foliage in upper and lower strata was examined for aspen (Populus tremuloides Michx.) and its associated hazelnut shrub (Corylus americana Walt.). We determined that C accumulation, expressed as dry mass per unit leaf area (mg C cm (exp -2)), was linearly dependent on N content (approximately 0.3- 3.5 mg N cm (exp -2))(r (exp 2) = 0.93, n=383, P less than 0.001) when eleven foliage groups were defined according to species, site, and developmental stage. C assembly was greatest in the upper aspen strata of both sites (seasonal average, 40.1 plus or minus 0.6 mg C cm (exp -2)), intermediate in the lower aspen strata (32.7 plus or minus 0.6), and considerably lower, and similar, in the hazelnut shrub layers (23.7 plus or minus 0.6) and in expanding aspen leaves (23.8 plus or minus 0.5); the lowest C assembly per unit N occurred in the two youngest, emerging leaf groups (17.1 plus or minus 0.6). Other relationships among physiological and biochemical variables were typically non-linear and were confounded by inclusion of the three groups of young (i.e., emerging or expanding) leaves, unless these were separately identified. Net C uptake, measured as photosynthetic capacity (A (sub max), micromole CO2 m (exp -2) s (exp -1)), was greater in aspen throughout the season, and optimal in mid-summer at a C:N ratio of approximately 18 (approximately 2.3 %N). When young leaves were excluded and logarithms of both variables were used, A (sub max) was approximately linearly dependent on N (mg N cm (exp-2) (r (exp 2) = 0.85, n= 193, P less than 0.001), attributed to incorporation of N into photosynthetic complexes and enzymes. In mature leaves, differences in pigment content vs. N among canopy strata were accentuated when N was expressed per unit leaf area (Mg cm (exp -2)) . However, the simplest log-linear relationship between a pigment variable and N was obtained for a ratio describing the relative allocation of photosynthetic pigment to Chl a (Chl a/[Chl b + carotenoids], microgram cm (exp -2)/ microgram cm-2) vs. %N (r (exp 2) = 0.90, n=343, P less than 0.001). Attainment of comparable A (sub max) Chl a content and relative Chl a allocation per unit N (mg cm (exp -2)) was achieved at different foliar N levels per canopy group: the lowest N requirement was for hazelnut leaves in the lowest, shaded stratum at the older, closed canopy site; the highest N requirement was in aspen leaves of the upper-most stratum at the younger, more open canopy site. These results highlight the differences in physiological responses between young and fully expanded leaves and show that sustaining those foliar constituents and processes important to C balance may require higher foliar N levels in leaves of establishing vs. mature aspen stands. There may be implications for remote-sensing assessments made for carbon balance in springtime, or over a landscape mosaic comprised of different aged stands.

Middleton, Elizabeth M.↗

Flight Control System Development for the BURRO Autonomous UAV

Developing autonomous flying vehicles has been a growing field in aeronautical research within the last decade and will continue into the next century. With concerns about safety, size, and cost of manned aircraft, several autonomous vehicle projects are currently being developed; uninhabited rotorcraft offer solutions to requirements for hover, vertical take-off and landing, as well as slung load transportation capabilities. The newness of the technology requires flight control engineers to question what design approaches, control law architectures, and performance criteria apply to control law development and handling quality evaluation. To help answer these questions, this paper documents the control law design process for Kaman Aerospace BURRO project. This paper will describe the approach taken to design control laws and develop math models which will be used to convert the manned K-MAX into the BURRO autonomous rotorcraft. With the ability of the K-MAX to lift its own weight (6000 lb) the load significantly affects the dynamics of the system; the paper addresses the additional design requirements for slung load autonomous flight. The approach taken in this design was to: 1) generate accurate math models of the K-MAX helicopter with and without slung loads, 2) select design specifications that would deliver good performance as well as satisfy mission criteria, and 3) develop and tune the control system architecture to meet the design specs and mission criteria. An accurate math model was desired for control system development. The Comprehensive Identification from Frequency Responses (CIFER(R)) software package was used to identify a linear math model for unloaded and loaded flight at hover, 50 kts, and 100 kts. The results of an eight degree-of-freedom CIFER(R)-identified linear model for the unloaded hover flight condition are presented herein, and the identification of the two-body slung-load configuration is in progress.

Colbourne, Jason D.↗

Atlantic Tropical Cyclone Monitoring with AMSU-A: Estimation of Maximum Sustained Wind Speeds

The first Advanced Microwave Sounding Unit temperature sounder (AMSU-A) was launched on the NOAA-15 satellite on 13 May 1998. The AMSU-A's higher spatial and radiometric resolutions provide more useful information on the strength of the middle and upper tropospheric warm cores associated with tropical cyclones than have previous microwave temperature sounders. The gradient wind relationship suggests that the temperature gradient near the core of tropical cyclones increases nonlinearly with wind speed. We recast the gradient wind equation to include AMSU-A derived variables. Stepwise regression is used to determine which of these variables is most closely related to maximum sustained winds (V(sub max)). The satellite variables investigated include the radially averaged gradients at two spatial resolutions of AMSU-A channels 1 through 10 T(sub b) data (delta(sub r)T(sub b)), the squares of these gradients, a channel 15 based scattering index (SI-89), and area averaged T(sub b). Calculations of Tb and delta(sub r)T(sub b) from mesoscale model simulations of Andrew reveal the effects of the AMSU spatial sampling on the cyclone warm core presentation. Stepwise regression of 66 AMSU-A terms against National Hurricane Center (NHC) V(sub max) estimates from the 1998 and 1999 Atlantic hurricane season confirms the existence of a nonlinear relationship between wind speed and radially averaged temperature gradients near the cyclone warm core. Of six regression terms, four are dominated by temperature information, and two are interpreted as correcting for hydrometeor contamination. Jackknifed regressions were performed to estimate the algorithm performance on independent data. For the 82 cases that had in situ measurements of V(sub max), the average error standard deviation was 4.7 m/s. For 108 cases without in situ wind data, the average error standard deviation was 7.5 m/s. Operational considerations, including the detection of weak cyclones and false alarm reduction are also discussed.

Spencer, Roy↗

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗