Search NASA⌕ Search

SEARCH · Search NASA

Results for “temperature advection”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

HSW-V v1.0: localized injections of interactive volcanic aerosols and their climate impacts in a simple general circulation model

Abstract. A new set of standalone parameterizations is presented for simulating the injection, evolution, and radiative forcing by stratospheric volcanic aerosols against an idealized Held–Suarez–Williamson (HSW) atmospheric background in the Energy Exascale Earth System Model version 2 (E3SMv2). In this model configuration (HSW with enabled volcanism, HSW-V), sulfur dioxide (SO2) and ash are injected into the atmosphere with a specified profile in the vertical, and they proceed to follow a simple exponential decay. The SO2 decay is modeled as a perfect conversion to a long-living sulfate aerosol which persists in the stratosphere. All three species are implemented as tracers in the model framework and are transported by the dynamical core's advection algorithm. The aerosols contribute simultaneously to local heating of the stratosphere and cooling of the surface by a simple plane-parallel Beer–Lambert law applied on two zonally symmetric radiation broadbands in the longwave and shortwave ranges. It is shown that the implementation parameters can be tuned to produce realistic temperature anomaly signatures of large volcanic events. In particular, results are shown for an ensemble of runs that mimic the volcanic eruption of Mt. Pinatubo in 1991. The design requires no coupling to microphysical subgrid-scale parameterizations and thus approaches the computational affordability of prescribed aerosol forcing strategies. The idealized simulations contain a single isolated volcanic event against a statistically uniform climate, where no background aerosols or other sources of externally forced variability are present. HSW-V represents a simpler-to-understand tool for the development of climate source-to-impact attribution methods.

Hollowed, Joseph P. (ORCID:0000000286581672)↗
NASA NTRS1993

Study of auroral dynamics with combined spacecraft and incoherent-scatter radar data

We have examined Sondrestrom incoherent-scatter radar observations of ionospheric plasma density and temperature distributions, as well as measurements of F-region ion drifts that were made during a prenoon pass by the DMSP-F7 satellite through the radar field of view. The spacecraft traversed a region of intense electron precipitation with a characteristic energy below approximately 200 eV. Particles with such low characteristic energies are believed to originate, either directly or indirectly, in the magnetosheath. The precipitation region had a width of about 2 deg invariant latitude. The corotating radar observed a patch of enhanced electron density and elevated electron temperature in the F2 region between about 10.5 and 12 magnetic local time in the same invariant latitude range where DMSP-F7 detected the soft-electron flux. The ion drift pattern, also obtained by radar, shows that it is unlikely that the plasma patch was produced by solar radiation and advected into the radar field of view. We suggest that the radar observed modifications of the ionospheric plasma distribution, which resulted from direct entry of magnetosheath electrons into the magnetosphere and down to ionospheric altitudes. Model calculations of the ionospheric response to the observed electron flux support our interpretation.

Watermann, Juergen↗
NASA NTRS1994

Ionospheric footprint of magnetosheathlike particle precipitation observed by an incoherent scatter radar

We have examined Sondrestrom incoherent scatter radar observations of ionospheric plasma density and temperature distributions and measurements of F region ion drifts that were made during a prenoon pass of the Defense Meteorological Satellite Program (DMSP)-F7 satellite through the radar field of view. The spacecraft traversed a region of intense electron precipitation with a characteristic energy below approximately 200 eV. Particles with such low characteristic energies are believed to be directly or indirectly of magnetosheath origin. The precipitation region had a width about 2 deg invariant latitude and covered the low-latitude boundary layer (LLBL), the cusp, and the equatorward section of the plasma mantle (PM). The corotating radar observed a patch of enhanced electron density and elevated electron temperature in the F2 region between about 10.5 and 12 magnetic local time in the same invariant latitude range where DMSP-F7 detected the soft-electron flux. The ion drift pattern, also obtained by radar, shows that it is unlikely that the plasma patch was produced by solar radiation and advected into the radar field of view. We suggest that the radar observed modifications of the ionospheric plasma distribution, which resulted from direct entry of magnetosheath electrons into the magnetosphere and down to ionospheric altitudes. Model calculations of the ionospheric response to the observed electron precipitation support our interpretation. The spectral characteristics of the electron flux in the LLBL, cusp, and equatorward section of the PM were in this case too similar to allow to distinguish between them by using incoherent scatter radar measurements only.

Watermann, Jurgen↗
NASA NTRS1993

Coupled ice-ocean variability in the Greenland Sea

A major surface feature of the Greenland Sea during winter is the frequent eastward extension of sea ice south of 75 deg N and an associated embayment to the north. These features are nominally connected with the East Greenland Current, and both the promontory) and the embayment are readily apparent on climatic ice charts. However, there are significant changes in these features on time-scales as short as a few days. Using a combination of satellite microwave images (SSM/I) of ice cover, meteorological data and in situ velocity, temperature and salinity records, we relate the ice distribution and its changes to the developing structure and circulation of the upper ocean during winter 1988-1989. Our measurements illustrate the preconditioning that leads to convective overturn, which in turn brings warmer water to the surface and results in the rapid disappearance of ice. In particular, the surface was cooled to the freezing point by early December and the salinity then increased through ice formation (about 0.016 m/d) and brine rejection. Once the vertical density gradient was sufficiently eroded, a period of high heat flux (greater than 300 W/sq m) in late January provided enough buoyancy loss to convectively mix the upper water column to at least 200 m. We estimate vertical velocities at about 3 cm/s downward during the initial sinking. The deepening of the thermocline raised surface temperatures by over 1 C resulting in nearly 1.5 x 105 km of ice-melt within two days. Average rates of ice retreat are about 2 km/d southwestward, generally consistent with a wind-driven flow. Comparison of hydrographic surveys from before and after the overturning indicate the fresh water was advected out of the area, possibly to the south and east of our moorings.

Roach, A. T.↗
NASA NTRS2001

Evaluation of Stratospheric Transport in New 3D Models Using the Global Modeling Initiative Grading Criteria

The Global Modeling Initiative (GMI) Team developed objective criteria for model evaluation in order to identify the best representation of the stratosphere. This work created a method to quantitatively and objectively discriminate between different models. In the original GMI study, 3 different meteorological data sets were used to run an offline chemistry and transport model (CTM). Observationally-based grading criteria were derived and applied to these simulations and various aspects of stratospheric transport were evaluated; grades were assigned. Here we report on the application of the GMI evaluation criteria to CTM simulations integrated with a new assimilated wind data set and a new general circulation model (GCM) wind data set. The Finite Volume Community Climate Model (FV-CCM) is a new GCM developed at Goddard which uses the NCAR CCM physics and the Lin and Rood advection scheme. The FV-Data Assimilation System (FV-DAS) is a new data assimilation system which uses the FV-CCM as its core model. One year CTM simulations of 2.5 degrees longitude by 2 degrees latitude resolution were run for each wind data set. We present the evaluation of temperature and annual transport cycles in the lower and middle stratosphere in the two new CTM simulations. We include an evaluation of high latitude transport which was not part of the original GMI criteria. Grades for the new simulations will be compared with those assigned during the original GMT evaluations and areas of improvement will be identified.

Strahan, Susan E.↗
DOE OSTI2026

A worldwide climatology of extreme air masses

Extreme temperature events are among the most damaging weather phenomena. In a warming world, more heat extremes and fewer cold extremes are expected in most regions in the future, a trade-off that warrants further understanding of such events. Here, we track and analyze large, persistent areas of hot and cold extreme temperatures in parallel, relative to the location and time of year, to quantify overall regional exposure to extreme temperatures. To accomplish this, we compare the frequencies, movements, trends, and sources/sinks in each type of extreme air mass, calling them extreme cold or extreme hot air masses (ECAMs/EHAMs). For most land regions, ECAMs occur more often than EHAMs, and ECAMs are more common in each hemisphere’s winter, when cold-air advection is strongest and most widespread. Average movement of ECAMs has a stronger equatorward component in winter than in summer, while movement of EHAMs is eastward all year, with less meridional movement than ECAMs. EHAMs have become more common almost everywhere on land, and the reverse is true for ECAMs, with the strongest trends in the northern hemisphere occurring in autumn, especially in the Arctic. The number of EHAMs around the world is increasing at a higher rate (+ 2.32 events/year) than ECAM numbers are decreasing (-1.27 events/year), causing a net increase in extremes, especially in some midlatitude regions. These results help to document the different processes driving hot and cold extremes, and thus the asymmetric and regionally varying trends in frequency for extreme air masses.

Ryan, James M↗
DOE OSTI2026

Modeling glass degradation and release of radionuclides from vitrified waste for performance assessment simulations

The release of radionuclides initially encapsulated in a slowly degrading solid waste form and contained in an eventually corroding canister defines the source term for numerical simulations for the assessment of a geologic repository for high-level radioactive waste. While the details of waste degradation, canister corrosion, and dissolution and mobilization of the radionuclides in pore water include complex chemical reaction and transport processes that are coupled to the thermal, hydrological, microbiological, and mechanical conditions in the repository, the source-term model suitable for use in a numerical performance assessment model should be a defensible abstraction of these mechanisms. We developed a radiological source-term model and implemented it into a non-isothermal flow and transport simulator. While the proposed source-term model is applicable to various waste forms, canister systems, and disposal concepts, we specifically considered radionuclide releases from vitrified high-level waste placed in a cylindrical canister disposed in a deep vertical borehole repository. In this model, waste degradation is a function of temperature, and it can be adjusted to evaluate the influence of and propagate uncertainties in pH, passivation reactions, and chemical conditions as well as geometrical factors. The time-dependent, congruent release of safety-relevant radionuclides present in the decaying inventory is then calculated. Finally, the radionuclides are mobilized by diffusive and advective transport according to the thermo-hydraulic conditions prevailing in the near field of the repository, from where they migrate through the geosphere to the accessible environment. We examine the influence of the source-term model’s parameters on performance assessment calculations through sensitivity and uncertainty propagation analyses, identifying influential factors and confirming the upper bound of their impact. These considerations align with the overarching goal of repository design, which is to demonstrate that engineered and natural barriers can collectively delay radionuclide migration for timescales far exceeding human planning, thereby providing multiple, redundant barriers against environmental contamination.

iTOUGH2↗
NASA NTRS1979

Frontal systems during passage of the Martian north polar hood over the Viking Lander 2 site prior to the first 1977 dust storm

Analysis of a 12-sol period of wind speed, wind direction, temperature, pressure, and optical depth at the Viking Lander 2 site presents the first in situ evidence of high- and low-pressure systems, complete with fronts, on the surface of Mars. The discovery of these systems in the Lander data occurred while analyzing a period during which the north polar hood was advected over the site at midday, dramatically decreasing the surface illumination and surface-to-atmospheric heat flux. This obscuration immediately preceded a global dust storm in the southern hemisphere and low latitudes of the northern hemisphere. The direct effects of the dust storm, reached 48 N, the Lander 2 latitude, later more gradually than they reached the 22 deg N latitude of Lander 1. The system responsible for the polar hood passages is a major disturbance, and it appears that radiational damping is inadequate to stop strong frontal formation. The front analyzed is characteristic of a repetitive series of systems that pass roughly every 3.3 sols. These systems are similar to those predicted by theoretical analyses and by general circulation modeling of the Martian atmosphere and those observed in laboratory experiments.

Tillman, J. E.↗
NASA NTRS2000

The Evolution of Tropospheric Temperature Field and its Relationship With The Onset of Asian Summer Monsoon

The mean state and year-to-year variations of the tropospheric temperature fields and their relationship with the establishment of the summertime East Asian monsoon (EAM) and the Indian monsoon (INM) are studied using the NCEP reanalysis data of 15 years (1982-1996). The results show that the seasonal shift of the South Asian High in the upper troposphere and the establishment of the EAM and the INM are closely related to the seasonal warming which causes a reversal of the meridional gradient of upper tropospheric mean temperature over the monsoon regions. On the average of 15 years, the reversal time of the temperature gradient in the EAM region (INM region) is concurrent with (one pentad earlier than) the onset time of the summer monsoon. In most years of the 15-year period, the reversal of temperature gradient coincides or precedes the onset time of the summer monsoon in both the EAM region and the INM region. The results suggest an important role of thermal processes on the establishment of the Asian monsoon. The contributors to the upper tropospheric warming over the EAM region are the strong horizontal warm advection and the diabetic heating against the adiabatic cooling due to upward motion. In the INM region, strong adiabatic heating by subsidence and the diabetic heating are major warming processes against the strong horizontal cold advection related to the persistent northwestlies to the southwestern periphery of the Tibetan Plateau. It appears that the early or late establishment of the Asian summer monsoon is not directly related to the differential warming near the surface.

He, H.↗
NASA NTRS2014

Modeling the Influences of Aerosols on Pre-Monsoon Circulation and Rainfall over Southeast Asia

We conduct several sets of simulations with a version of NASA's Goddard Earth Observing System, version 5, (GEOS-5) Atmospheric Global Climate Model (AGCM) equipped with a two-moment cloud microphysical scheme to understand the role of biomass burning aerosol (BBA) emissions in Southeast Asia (SEA) in the pre-monsoon period of February-May. Our experiments are designed so that both direct and indirect aerosol effects can be evaluated. For climatologically prescribed monthly sea surface temperatures, we conduct sets of model integrations with and without biomass burning emissions in the area of peak burning activity, and with direct aerosol radiative effects either active or inactive. Taking appropriate differences between AGCM experiment sets, we find that BBA affects liquid clouds in statistically significantly ways, increasing cloud droplet number concentrations, decreasing droplet effective radii (i.e., a classic aerosol indirect effect), and locally suppressing precipitation due to a deceleration of the autoconversion process, with the latter effect apparently also leading to cloud condensate increases. Geographical re-arrangements of precipitation patterns, with precipitation increases downwind of aerosol sources are also seen, most likely because of advection of weakly precipitating cloud fields. Somewhat unexpectedly, the change in cloud radiative effect (cloud forcing) at surface is in the direction of lesser cooling because of decreases in cloud fraction. Overall, however, because of direct radiative effect contributions, aerosols exert a net negative forcing at both the top of the atmosphere and, perhaps most importantly, the surface, where decreased evaporation triggers feedbacks that further reduce precipitation. Invoking the approximation that direct and indirect aerosol effects are additive, we estimate that the overall precipitation reduction is about 40% due to the direct effects of absorbing aerosols, which stabilize the atmosphere and reduce surface latent heat fluxes via cooler land surface temperatures. Further refinements of our two-moment cloud microphysics scheme are needed for a more complete examination of the role of aerosol-convection interactions in the seasonal development of the SEA monsoon.

SE Asia Monsoon↗
NASA NTRS1981

Numerical modeling of diffusive-convective physical vapor transport in cylindrical vertical ampoules

Diffusive-convective physical vapor transport (PVT) in cylindrical, vertical ampoules of aspect ratio (length/radius) between 0.5 and 10 was modeled numerically. The transport of a crystal forming component through an inert component that undergoes zero net transport was considered. Systems were treated in which: (1) with unequal molecular weight of the components and with temperature gradients typically employed in PVT, convective flow arises dominantly from solutal density gradients; and (2) with equal molecular weight of the components, convective flow can arise only from thermal expansion. It was found that, due to the diffusion-induced horizontal density gradients, buoyancy-driven convective flows are superimposed on the diffusive-advective fluxes without threshold. Net recirculation sets in adjacent to the growing interface, in contrast to the corresponding monocomponent situation where marginally stable convective modes fill the whole fluid space. Depending on the orientation of the main transport direction with respect to gravity, convection can either reduce or enhance the diffusion-induced radial concentration gradients. Significant enhancement of the net transport rate was found to occur only when the whole vapor space between source and growing crystal is filled by a convective recirculation roll. Solutal and thermal convection results are similar; yet for quantitative discussions, thermal and solutal Rayleigh numbers are not interchangeable in contrast to convective situations that lack net mass transport across the fluid space.

Markham, B. L.↗
NASA NTRS1986

Large-scale short-period sea ice atmosphere interaction

Changes in the microwave brightness temperature measured by the Electrically Scanning Microwave Radiometer (ESMR) flown on board the Nimbus V satellite reveal large-scale sea ice fluctuations in the Antarctic marginal ice zone. These ice margin fluctuations arre predominantly wave numbers 1-4, with phase speeds of about 3 m/s independent of wave number. The spatial pattern and eastward advection of the sea ice anomalies match those of the atmospheric sea level pressure, and are consistent with sea ice displacement due to surface wind stress. Examination of the outgoing longwave radiation indicates that suppression of high clouds in regions of increased sea ice inc reases the radiative cooling which contributes to maintaining the ice. Data from three winter seasons indicate about a one-third probability of occurence of this large scale high frequency sea ice atmosphere interaction during any given 2-week period in winter.

Cahalan, Robert F.↗
DOE OSTI2021

The Implementation of Framework for Improvement by Vertical Enhancement Into Energy Exascale Earth System Model

Abstract The low cloud bias in global climate models (GCMs) remains an unsolved problem. Coarse vertical resolution in GCMs has been suggested to be a significant cause of low cloud bias because planetary boundary layer parameterizations cannot resolve sharp temperature and moisture gradients often found at the top of subtropical stratocumulus layers. This work aims to ameliorate the low cloud problem by implementing a new computational method, the Framework for Improvement by Vertical Enhancement (FIVE), into the Energy Exascale Earth System Model (E3SM). Three physics schemes representing microphysics, radiation, and turbulence as well as vertical advection are interfaced to vertically enhanced physics (VEP), which allows for these processes to be computed on a higher vertical resolution grid compared to the rest of the E3SM model. We demonstrate the better representation of subtropical boundary layer clouds with FIVE while limiting additional computational cost from the increased number of levels. When the vertical resolution approaches the large eddy simulation‐like vertical resolution in VEP, the climatological low cloud amount shows a significant increase of more than 30% in the southeastern Pacific Ocean. Using FIVE to improve the representation of low‐level clouds does not come with any negative side effects associated with the simulation of mid‐ and high‐level cloud and precipitation, that can occur when running the full model at higher vertical resolution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗
NASA NTRS1987

Cometary MHD and chemistry

An MHD and chemical comet-coma model was developed, applying the computer program of Huebner (1985) for the detailed chemical evolution of a spherically expanding coma and the program of Schmidt and Wegman (1982) and Wegman (1987) for the MHD flow of plasma and magnetic field in a comet to the Giotto-mission data on the ion abundances measured by the HIS ion mass spectrometer. The physics and chemistry of the coma are modeled in great detail, including photoprocesses, gas-phase chemical kinetics, energy balance with a separate electron temperature, multifluid hydrodynamics with a transition to free molecular flow, fast-streaming atomic and molecular hydrogen, counter and cross streaming of the ionized species relative to the neutral species in the coma-solar wind interaction region with momentum exchange by elastic collisions, mass-loading through ion pick-up, and Lorentz forces of the advected magnetic field. The results, both inside and outside of the contact surface, are discussed and compared with the relevant HIS ion mass spectra.

Wegmann, R.↗
NASA NTRS2001

The Extremely Warm Early Winter 2000 in Europe: What is the Forcing

High variability characterizes the winter climate of central Europe: interannual fluctuations in the surface-air temperature as large as 18 C over large areas are fairly common. The extraordinary early-winter 2000 in Europe appears to be a departure to an unprecedented extreme of the existing climate patterns. Such anomalous events affect agriculture, forestry, fuel consumption, etc., and thus deserve in-depth analysis. Our analysis indicates that the high anomalies of the surface-air temperature are predominantly due to the southwesterly flow from the eastern North Atlantic, with a weak contribution by southerly flow from the western Mediterranean. Backward trajectories based on the SSM/I and NCEP Reanalysis datasets traced from west-central Europe indicate that the warm air masses flowing into Europe originate in the southern North Atlantic, where the surface-air temperatures exceed by 15c or more the climatic norms in Europe for late-November or early-December. Because such large ocean-to-continent temperature differences characterize the winter conditions, we refer to this episode which started in late November as occurring in the early winter. In this season, with the sun low over the horizon in Europe, absorption of insolation by the surface has little significance. The effect of cloudiness, a corollary to the low-level maritime-air advection, is a warming by a reduction of heat loss (greenhouse effect). In contrast, in the summer, clouds, by reducing absorption of insolation, produce a cooling, effect at the surface.

Otterman, J.↗
NASA NTRS1994

Implications of summertime marine stratocumulus on the North American climate

This study focuses on the effects of summertime stratocumulus over the eastern Pacific. This cloud is linked to the semi-permanent sub-tropical highs that dominate the low-level circulation over the Pacific and Atlantic. Subsidence on the eastern flank of these highs creates an inversion based about 800 m above sea level that caps moist air near the surface. This air overlies cool waters driven by upwelling along the coastal regions of North America. Strong surface north-westerlies mix the boundary layer enough to saturate the air just below the capping inversion. Widespread stratocumulus is thus formed. All calculations were carried out using the GENESIS general circulation model that was run at MSFC. Among the more important properties of the model is that it includes radiative forcing due to absorption of solar radiation and the emission of infrared radiation, interactive clouds (both stratocumulus and cumulus types), exchanges of heat and moisture with the lower boundary. Clouds are interactive in the sense that they impact the circulation by modifying the fields of radiative heating and turbulent fluxes of heat and moisture in the boundary layer. In turn, clouds are modified by the winds through the advection of moisture. In order to isolate the effects of mid- and high-latitude stratocumulus, two runs were made with the model: one with and the other without stratocumulus. The runs were made for a year, but with perpetual July conditions, i.e., solar forcing was fixed. The diurnal solar cycle, however, was allowed for. The sea surface temperature distribution was fixed in both runs to represent climatological July conditions. All dependent variables were represented at 12 surfaces of constant sigma = p/p(sub O), where p is pressure and p(sub O) is surface pressure. To facilitate analysis, model output was transformed to constant pressure surfaces. Structures no smaller in size than 7.5 degrees longitude and 4.5 degrees in latitude were resolved. Smaller features of the circulation were parameterized. The model thus captures synoptic- and planetary-scale circulation features.

Clark, John H. E.↗
DOE OSTI2026

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗