Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large grids”

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,243 records · Page 69

Large-eddy simulation of a backward facing step flow using a least-squares spectral element method

We report preliminary results obtained from the large eddy simulation of a backward facing step at a Reynolds number of 5100. The numerical platform is based on a high order Legendre spectral element spatial discretization and a least squares time integration scheme. A non-reflective outflow boundary condition is in place to minimize the effect of downstream influence. Smagorinsky model with Van Driest near wall damping is used for sub-grid scale modeling. Comparisons of mean velocity profiles and wall pressure show good agreement with benchmark data. More studies are needed to evaluate the sensitivity of this method on numerical parameters before it is applied to complex engineering problems.

Chan, Daniel C.↗

A Timescale for Major Events in Early Mars Crustal Evolution

The population of visible and buried impact basins greater than 200 km diameter revealed by high resolution gridded MOLA data and the cumulative frequency curves derived for these provide a basis for a chronology of major events in early martian history. The relative chronology can be given in terms of N(200) crater retention ages; absolute ages can be assigned using the Hartmann-Neukum (H&N) model chronology. In terms of billions of H&N years, the crustal dichotomy formed by large impact basins at 4.12 +/- 0.08 BYA [N(200) = 3.0 - 3.2] and the global magnetic field died at about or slightly before the same time (4.15 +/- 0.08 BYA) [N(200) = 3.5]. In this chronology, the buried lowlands are ~120 my younger than the buried highlands), approx. 160 my younger than the highlands overall and approx. 340 my younger than the oldest crater retention surface we see, defined by the largest impact basins.

Frey, H. V.↗

A Time Scale for Major Events in Early Mars Crustal Evolution

The population of visible and buried impact basins > 200 km diameter revealed by high resolution gridded MOLA data and the cumulative frequency curves derived for these pvide a basis for a chronology of major events in early martian history. The relative chronology can be given in terms of N(200) crater retention ages; 'absolute ages' can be assigued using the Hartmann-Neukum (H&N) model chronology. In terms of billions of H&N years, the crustal dichotomy formed by large impact basins at 4.12 +/- 0.08 BYA (N(200) = 3.0-3.2) and the global magnetic field died at about or slightly before the same time (4.15 +/- 0.08 BYA (N(200) = 3.5). In this chronology, the buried lowlands are approx. 120 my younger than the buried highlands, approx. 160 my younger than the highlands overall and approx. 340 my younger than the oldest crater retention surface we see, defined by the largest impact basins.

Frey, Herbert V.↗

Challenges to Computational Aerothermodynamic Simulation and Validation for Planetary Entry Vehicle Analysis

Challenges to computational aerothermodynamic (CA) simulation and validation of hypersonic flow over planetary entry vehicles are discussed. Entry, descent, and landing (EDL) of high mass to Mars is a significant driver of new simulation requirements. These requirements include simulation of large deployable, flexible structures and interactions with reaction control system (RCS) and retro-thruster jets. Simulation of radiation and ablation coupled to the flow solver continues to be a high priority for planetary entry analyses, especially for return to Earth and outer planet missions. Three research areas addressing these challenges are emphasized. The first addresses the need to obtain accurate heating on unstructured tetrahedral grid systems to take advantage of flexibility in grid generation and grid adaptation. A multi-dimensional inviscid flux reconstruction algorithm is defined that is oriented with local flow topology as opposed to grid. The second addresses coupling of radiation and ablation to the hypersonic flow solver - flight- and ground-based data are used to provide limited validation of these multi-physics simulations. The third addresses the challenges of retro-propulsion simulation and the criticality of grid adaptation in this application. The evolution of CA to become a tool for innovation of EDL systems requires a successful resolution of these challenges.

Gnoffo, Peter A.↗

Challenges to Computational Aerothermodynamic Simulation and Validation for Planetary Entry Vehicle Analysis

Challenges to computational aerothermodynamic (CA) simulation and validation of hypersonic flow over planetary entry vehicles are discussed. Entry, descent, and landing (EDL) of high mass to Mars is a significant driver of new simulation requirements. These requirements include simulation of large deployable, flexible structures and interactions with reaction control system (RCS) and retro-thruster jets. Simulation of radiation and ablation coupled to the flow solver continues to be a high priority for planetary entry analyses, especially for return to Earth and outer planet missions. Three research areas addressing these challenges are emphasized. The first addresses the need to obtain accurate heating on unstructured tetrahedral grid systems to take advantage of flexibility in grid generation and grid adaptation. A multi-dimensional inviscid flux reconstruction algorithm is defined that is oriented with local flow topology as opposed to grid. The second addresses coupling of radiation and ablation to the hypersonic flow solver--flight- and ground-based data are used to provide limited validation of these multi-physics simulations. The third addresses the challenges of retro-propulsion simulation and the criticality of grid adaptation in this application. The evolution of CA to become a tool for innovation of EDL systems requires a successful resolution of these challenges.

Gnoffo, Peter A.↗

Importance of the Tectonic Motions on Ganymede

Voyager close up views of Ganymede reveal some provinces of old, cratered, dark terrains that are broken into polygons and cut through by light stripes of grooved terrains. The groove terrain unit and its relationships with the other units was studied at various scales because these features seem to relate at least to two main scales of processes. Minor displacements concerning only local features, i.e., large scale processes and major displacements concerning usually the major evolution of the surface, i.e., small scale processes. It appears that tectonic processes, at least in the studied areas, if playing a role in the formation of local or regional features, are very limited. The conclusions reached agree with those of Thomas et al. who show that post grooves basins reveal the presence of a primitive grid that is not disturbed by the grooved terrains.

Forni, O. P.↗

High Reynolds number calculations using the dynamic subgrid-scale stress model

The dynamic subgrid-scale eddy viscosity model was used in the large-eddy simulation of the turbulent flow in a plane channel for Reynolds numbers based on friction velocity and channel half-width ranging between 200 and 2000, a range including values significantly higher than in previous simulations. The computed wall stress, mean velocity, and Reynolds stress profiles compare very well with experimental and direct simulation data. Comparison of higher moments is also satisfactory. Although the grid in the near-wall region is fairly coarse, the results are quite accurate: the turbulent kinetic energy peaks at y(+) of about 12, and the near-wall behavior of the resolved stresses is captured accurately. The model coefficient is o(0.001) in the buffer layer and beyond, where the cutoff wave numbers are in the decaying region of the spectra; in the near-wall region, the cutoff wave numbers are nearer the energy-containing range, and the resolved turbulent stresses become a constant fraction of the resolved stresses. This feature is responsible for the correct near-wall behavior of the model coefficient.

Piomelli, Ugo↗

Parallelization of the Physical-Space Statistical Analysis System (PSAS)

Atmospheric data assimilation is a method of combining observations with model forecasts to produce a more accurate description of the atmosphere than the observations or forecast alone can provide. Data assimilation plays an increasingly important role in the study of climate and atmospheric chemistry. The NASA Data Assimilation Office (DAO) has developed the Goddard Earth Observing System Data Assimilation System (GEOS DAS) to create assimilated datasets. The core computational components of the GEOS DAS include the GEOS General Circulation Model (GCM) and the Physical-space Statistical Analysis System (PSAS). The need for timely validation of scientific enhancements to the data assimilation system poses computational demands that are best met by distributed parallel software. PSAS is implemented in Fortran 90 using object-based design principles. The analysis portions of the code solve two equations. The first of these is the "innovation" equation, which is solved on the unstructured observation grid using a preconditioned conjugate gradient (CG) method. The "analysis" equation is a transformation from the observation grid back to a structured grid, and is solved by a direct matrix-vector multiplication. Use of a factored-operator formulation reduces the computational complexity of both the CG solver and the matrix-vector multiplication, rendering the matrix-vector multiplications as a successive product of operators on a vector. Sparsity is introduced to these operators by partitioning the observations using an icosahedral decomposition scheme. PSAS builds a large (approx. 128MB) run-time database of parameters used in the calculation of these operators. Implementing a message passing parallel computing paradigm into an existing yet developing computational system as complex as PSAS is nontrivial. One of the technical challenges is balancing the requirements for computational reproducibility with the need for high performance. The problem of computational reproducibility is well known in the parallel computing community. It is a requirement that the parallel code perform calculations in a fashion that will yield identical results on different configurations of processing elements on the same platform. In some cases this problem can be solved by sacrificing performance. Meeting this requirement and still achieving high performance is very difficult. Topics to be discussed include: current PSAS design and parallelization strategy; reproducibility issues; load balance vs. database memory demands, possible solutions to these problems.

Larson, J. W.↗

Session on High Speed Civil Transport Design Capability Using MDO and High Performance Computing

Since the inception of CAS in 1992, NASA Langley has been conducting research into applying multidisciplinary optimization (MDO) and high performance computing toward reducing aircraft design cycle time. The focus of this research has been the development of a series of computational frameworks and associated applications that increased in capability, complexity, and performance over time. The culmination of this effort is an automated high-fidelity analysis capability for a high speed civil transport (HSCT) vehicle installed on a network of heterogeneous computers with a computational framework built using Common Object Request Broker Architecture (CORBA) and Java. The main focus of the research in the early years was the development of the Framework for Interdisciplinary Design Optimization (FIDO) and associated HSCT applications. While the FIDO effort was eventually halted, work continued on HSCT applications of ever increasing complexity. The current application, HSCT4.0, employs high fidelity CFD and FEM analysis codes. For each analysis cycle, the vehicle geometry and computational grids are updated using new values for design variables. Processes for aeroelastic trim, loads convergence, displacement transfer, stress and buckling, and performance have been developed. In all, a total of 70 processes are integrated in the analysis framework. Many of the key processes include automatic differentiation capabilities to provide sensitivity information that can be used in optimization. A software engineering process was developed to manage this large project. Defining the interactions among 70 processes turned out to be an enormous, but essential, task. A formal requirements document was prepared that defined data flow among processes and subprocesses. A design document was then developed that translated the requirements into actual software design. A validation program was defined and implemented to ensure that codes integrated into the framework produced the same results as their standalone counterparts. Finally, a Commercial Off the Shelf (COTS) configuration management system was used to organize the software development. A computational environment, CJOPT, based on the Common Object Request Broker Architecture, CORBA, and the Java programming language has been developed as a framework for multidisciplinary analysis and Optimization. The environment exploits the parallelisms inherent in the application and distributes the constituent disciplines on machines best suited to their needs. In CJOpt, a discipline code is "wrapped" as an object. An interface to the object identifies the functionality (services) provided by the discipline, defined in Interface Definition Language (IDL) and implemented using Java. The results of using the HSCT4.0 capability are described. A summary of lessons learned is also presented. The use of some of the processes, codes, and techniques by industry are highlighted. The application of the methodology developed in this research to other aircraft are described. Finally, we show how the experience gained is being applied to entirely new vehicles, such as the Reusable Space Transportation System. Additional information is contained in the original.

Rehder, Joe↗

The Impact of Subsampling on MODIS Level-3 Statistics of Cloud Optical Thickness and Effective Radius

The MODIS Level-3 optical thickness and effective radius cloud product is a gridded l deg. x 1 deg. dataset that is derived from aggregation and subsampling at 5 km of 1 km, resolution Level-2 orbital swath data (Level-2 granules). This study examines the impact of the 5 km subsampling on the mean, standard deviation and inhomogeneity parameter statistics of optical thickness and effective radius. The methodology is simple and consists of estimating mean errors for a large collection of Terra and Aqua Level-2 granules by taking the difference of the statistics at the original and subsampled resolutions. It is shown that the Level-3 sampling does not affect the various quantities investigated to the same degree, with second order moments suffering greater subsampling errors, as expected. Mean errors drop dramatically when averages over a sufficient number of regions (e.g., monthly and/or latitudinal averages) are taken, pointing to a dominance of errors that are of random nature. When histograms built from subsampled data with the same binning rules as in the Level-3 dataset are used to reconstruct the quantities of interest, the mean errors do not deteriorate significantly. The results in this paper provide guidance to users of MODIS Level-3 optical thickness and effective radius cloud products on the range of errors due to subsampling they should expect and perhaps account for, in scientific work with this dataset. In general, subsampling errors should not be a serious concern when moderate temporal and/or spatial averaging is performed.

Oreopoulos, Lazaros↗

Volumetric Acoustic Vector Intensity Probe

A new measurement tool capable of imaging the acoustic intensity vector throughout a large volume is discussed. This tool consists of an array of fifty microphones that form a spherical surface of radius 0.2m. A simultaneous measurement of the pressure field across all the microphones provides time-domain near-field holograms. Near-field acoustical holography is used to convert the measured pressure into a volumetric vector intensity field as a function of frequency on a grid of points ranging from the center of the spherical surface to a radius of 0.4m. The volumetric intensity is displayed on three-dimensional plots that are used to locate noise sources outside the volume. There is no restriction on the type of noise source that can be studied. The sphere is mobile and can be moved from location to location to hunt for unidentified noise sources. An experiment inside a Boeing 757 aircraft in flight successfully tested the ability of the array to locate low-noise-excited sources on the fuselage. Reference transducers located on suspected noise source locations can also be used to increase the ability of this device to separate and identify multiple noise sources at a given frequency by using the theory of partial field decomposition. The frequency range of operation is 0 to 1400Hz. This device is ideal for the study of noise sources in commercial and military transportation vehicles in air, on land and underwater.

Klos, Jacob↗

Exercise Equipment Usability Assessment for a Deep Space Concept Vehicle

With international aspirations to send astronauts to deep space, the world is now faced with the complex problem of keeping astronauts healthy in unexplored hostile environments for durations of time never before attempted by humans. The great physical demands imparted by space exploration compound the problem of astronaut health, as the astronauts must not only be healthy, but physically fit upon destination arrival in order to perform the scientific tasks required of them. Additionally, future deep space exploration necessitates the development of environments conducive to long-duration habitation that would supplement propulsive vehicles. Space Launch System (SLS) core stage barrel sections present large volumes of robust structure that can be recycled and used for long duration habitation. This assessment will focus on one such conceptual craft, referred to as the SLS Derived Habitat (SLS-DH). Marshall Space Flight Center's (MSFC) Advanced Concepts Office (ACO) has formulated a high-level layout of this SLS-DH with parameters such as floor number and orientation, floor designations, grid dimensions, wall placement, etc. Yet to be determined, however, is the layout of the exercise area. Currently the SLS-DH features three floors laid out longitudinally, leaving 2m of height between the floor and ceilings. This short distance between levels introduces challenges for proper placement of exercise equipment such as treadmills and stationary bicycles, as the dynamic envelope for the 95th percentile male astronauts is greater than 2m. This study aims to assess the optimal equipment layout and sizing for the exercise area of this habitat. Figure 1 illustrates the layout of the DSH concept demonstrator located at MSFC. The exercise area is located on the lower level, seen here as the front half of the level occupied by a crew member. This small volume does not allow for numerous or bulky exercise machines, so the conceptual equipment has been limited to a treadmill and stationary bicycle. With the most current treadmill aboard the International Space Station (ISS), the Combined Operational Load-Bearing External Resistance Treadmill (COLBERT), being located in an International Standard Payload Rack (ISPR), the bottom of the conceptual treadmill features a height of 38in. Making the treadmill flush with the floor would be impossible in this rack configuration, as the distance from the outer wall of the spacecraft to the bottom floor would be too shallow. From preliminary sizing, the 38in required for the bottom of the treadmill combined with a 78in operational envelope for a 95th percentile may not be accommodated in the exercise area in a vertical orientation. Figure 2 demonstrates the volume required (in maroon) for an ISPR-bound treadmill in the concept demonstrator. Early indications as seen in this figure indicate that the crew members would contact the ceiling in such an arrangement. An assessment will be conducted to evaluate various orientations of exercise equipment in the concept demonstrator. Orientations to be tested include putting the bottom of the treadmill on the wall, having the treadmill at an angle in the floor both horizontally and vertically, and having a shorter (non-rack bound) treadmill in a vertical orientation on the floor. This assessment will yield findings regarding sizing of the area and how well participants feel they could exercise in such an environment. Due to the restrictions of assessing a microgravity vehicle in a normal-gravity environment, simulations in MSFC's Virtual Environments Lab (VEL) may be necessary. Final deliverables will include recommendations regarding the location and size of possible exercise equipment aboard the SLS-Derived DSH.

Rhodes, Brooke M.↗

Development of an Output-based Adaptive Method for Multi-Dimensional Euler and Navier-Stokes Simulations

The use of computational simulations in the prediction of complex aerodynamic flows is becoming increasingly prevalent in the design process within the aerospace industry. Continuing advancements in both computing technology and algorithmic development are ultimately leading to attempts at simulating ever-larger, more complex problems. However, by increasing the reliance on computational simulations in the design cycle, we must also increase the accuracy of these simulations in order to maintain or improve the reliability arid safety of the resulting aircraft. At the same time, large-scale computational simulations must be made more affordable so that their potential benefits can be fully realized within the design cycle. Thus, a continuing need exists for increasing the accuracy and efficiency of computational algorithms such that computational fluid dynamics can become a viable tool in the design of more reliable, safer aircraft. The objective of this research was the development of an error estimation and grid adaptive strategy for reducing simulation errors in integral outputs (functionals) such as lift or drag from from multi-dimensional Euler and Navier-Stokes simulations. In this final report, we summarize our work during this grant.

Darmofal, David L.↗

Wavenumber Domain Focusing of Squinted SAR Data with a Curved Orbit Geometry

Synthetic Aperture Radar systems provide raw data that need focusing to achieve full-resolution imaging. Current SAR applications, including interferometry, require accurate, phase-preserving, and precisely co-registered coherent images over large ground swaths with the highest achievable resolution. In addition to these challenges, stripmap SAR data may be acquired with an off-broadside (squinted) geometry, either by design or through platform motion. The precise batch focusing of these large aperture and wide bandwidth data sets is known to require a 2D frequency processing approach. The standard wave domain focusing algorithm, however, is only exact for data acquired on a rectilinear trajectory. We investigate a generalization of the standard omega-k focusing formulation that allows curved data acquisition tracks. The new formulation can be used in conjunction with a known extension for conical, squinted imaging grids. The approximations necessary to allow the generalized geometry are analyzed to determine the range of applicability of the proposed algorithm. The theory is validated using data simulated with parameters similar to the UAVSAR L-band SAR system.

curved orbit geometry↗

Impact of Satellite Viewing-Swath Width on Global and Regional Aerosol Optical Thickness Statistics and Trends

We use the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite aerosol optical thickness (AOT) product to assess the impact of reduced swath width on global and regional AOT statistics and trends. Alongtrack and across-track sampling strategies are employed, in which the full MODIS data set is sub-sampled with various narrow-swath (approximately 400-800 km) and single pixel width (approximately 10 km) configurations. Although view-angle artifacts in the MODIS AOT retrieval confound direct comparisons between averages derived from different sub-samples, careful analysis shows that with many portions of the Earth essentially unobserved, spatial sampling introduces uncertainty in the derived seasonal-regional mean AOT. These AOT spatial sampling artifacts comprise up to 60%of the full-swath AOT value under moderate aerosol loading, and can be as large as 0.1 in some regions under high aerosol loading. Compared to full-swath observations, narrower swath and single pixel width sampling exhibits a reduced ability to detect AOT trends with statistical significance. On the other hand, estimates of the global, annual mean AOT do not vary significantly from the full-swath values as spatial sampling is reduced. Aggregation of the MODIS data at coarse grid scales (10 deg) shows consistency in the aerosol trends across sampling strategies, with increased statistical confidence, but quantitative errors in the derived trends are found even for the full-swath data when compared to high spatial resolution (0.5 deg) aggregations. Using results of a model-derived aerosol reanalysis, we find consistency in our conclusions about a seasonal-regional spatial sampling artifact in AOT Furthermore, the model shows that reduced spatial sampling can amount to uncertainty in computed shortwave top-ofatmosphere aerosol radiative forcing of 2-3 W m(sup−2). These artifacts are lower bounds, as possibly other unconsidered sampling strategies would perform less well. These results suggest that future aerosol satellite missions having significantly less than full-swath viewing are unlikely to sample the true AOT distribution well enough to obtain the statistics needed to reduce uncertainty in aerosol direct forcing of climate.

remote sensing↗

“It’s Raining Bits”: Patterns in Directional Precipitation Persistence Across the United States

The spatial and temporal ordering of precipitation occurrence impacts ecosystems, streamflow, and water availability. For example, both large-scale climate patterns and local landscapes drive weather events, and the typical speeds and directions of these events moving across a basin dictate the timing of flows at its outlet. We address the predictability of precipitation occurrence at a given location, based on the knowledge of past precipitation at surrounding locations. We identify ‘‘dominant directions of precipitation influence’’ across the continental United States based on a gridded daily dataset. Specifically, we apply information theory–based measures that characterize dominant directions and strengths of spatial and temporal precipitation dependencies. On a national average, this dominant direction agrees with the prevalent direction of weather movement from west to east across the country, but regional differences reflect topographic divides, precipitation gradients, and different climatic drivers of precipitation. Trends in these information relationships and their correlations with climate indices over the past 70 years also show seasonal and spatial divides. This study expands upon a framework of information-based predictability to answer questions about spatial connectivity in addition to temporal persistence. The methods presented here are generally useful to understand many aspects of weather and climate variability.

North America↗

Determining True Sensor Spatial Resolution of Very High Resolution Optical Imagery

Some satellite data is delivered in images with gridded pixels. This gridded pixel size is often assumed to be the spatial resolution of the satellite sensor; however, this is not always the case. An image can be grided to any arbitrary pixel size, but the sensor resolution will remain constant. For example, an image with a pixel grid size much smaller than the sensor resolution will appear blurry along what should be sharp transitions. This discrepancy between an image’s pixel size and true sensor spatial resolution can be the source of much confusion and even misinformation among data users, which may lead them to waste time and resources on using images that do not suit their spatial resolution needs. This presentation will highlight our evaluation of the true spatial resolution of various government and commercial images in the pixel size range of 0.3 m to 60 m. Images evaluated include ESA’s Sentinel-2 (60 m, 20 m, 10 m pixels), USGS’s Landsat 8/9 (30 m & 15 m pixels), Planet’s SuperDoves (3 m pixels), BlackSky’s Globals (~1 m pixels), and the optical bands of Maxar’s WorldView-2 (2.4 m – 0.41 m pixels) and WorldView-3 (1.38 m – 0.31 m pixels). Our evaluation of true sensor spatial resolution, or ‘footprint size’ is based on the sensor’s line spread function (LSF). We calculate the width at half the height of the LSF to find the full width at half maximum (FWHM). The FWHM is how we report sensor spatial resolution. Different objects are examined for constructing the LSF depending on the sensor spatial resolution. Coarser resolution sensors in this evaluation such as Sentinel-2 and Landsat 8/9 are examined at bridges over a dark water background. The bright bridge acts as a line impulse, giving a sensor’s line spread function (LSF) in one direction. Additionally, we simulate the impacts of bridge width on the apparent LSF to obtain a true LSF without the effects of bridge width for these sensors. Finer resolution sensors will image the irregularities in bridges such as trusses, sidewalks, and in some cases painted lines, interfering with the LSF construction. Instead, these sensors are evaluated at large (60 m – 140 m) black and white checkerboards known as Cal/Val sites. At these locations, the image’s transition from black to white is extracted as an edge spread function (ESF). We calculate the derivative of this ESF to obtain the sensor’s LSF. From there, we find the FWHM as we do for the coarser resolution images. With the FWHM and pixel size, we determine how over- or under-sampled the images are. When the ratio of a sensor’s spatial resolution and the gridded image’s pixel size is less than 1, the image is considered under-sampled. In this case, each pixel’s information is unique but only a portion of that pixel’s ground area has been measured. On the other side, if the ratio is greater than 1, the image is considered over-sampled. That is, each pixel’s information is sourced from within the ground extent of the pixel and some extent outside additionally. We will show the true spatial resolution and the extent of over-/under-sampling in the imagery from ESA’s Sentinel-2 (60 m, 20 m, 10 m pixels), USGS’s Landsat 8/9 (30 m & 15 m pixels), Planet’s SuperDoves (3 m pixels), BlackSky’s Globals (~1 m pixels), and the optical bands of Maxar’s WorldView-2 (2.4 m – 0.41 m pixels) and WorldView-3 (1.38 m – 0.31 m pixels).

Alana Semple↗

MHD modelling of the earth's magnetosphere

Attention is given to the shock capturing technique, nonuniform grid system, and multiple time scale problem aspects of a global MHD model of the earth's magnetosphere that uses numerical methods. Because of the fast Alfven waves near the earth, the model is more difficult to solve than problems of hydrodynamic flow past bodies. The Rusanov (1962) scheme, which is a variant of the Lax (1954) method that maintains viscosity at a minimum value except where a large one is required, is used. The magnetosphere shape obtained by these means is noted.

Wu, C. C.↗