SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

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

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

NASA NTRS2004

Use of GOES, SSM/I, TRMM Satellite Measurements Estimating Water Budget Variations in Gulf of Mexico - Caribbean Sea Basins

This study presents results from a multi-satellite/multi-sensor retrieval system designed to obtain the atmospheric water budget over the open ocean. A combination of 3ourly-sampled monthly datasets derived from the GOES-8 5-channel Imager, the TRMM TMI radiometer, and the DMSP 7-channel passive microwave radiometers (SSM/I) have been acquired for the combined Gulf of Mexico-Caribbean Sea basin. Whereas the methodology has been tested over this basin, the retrieval system is designed for portability to any open-ocean region. Algorithm modules using the different datasets to retrieve individual geophysical parameters needed in the water budget equation are designed in a manner that takes advantage of the high temporal resolution of the GOES-8 measurements, as well as the physical relationships inherent to the TRMM and SSM/I passive microwave measurements in conjunction with water vapor, cloud liquid water, and rainfall. The methodology consists of retrieving the precipitation, surface evaporation, and vapor-cloud water storage terms in the atmospheric water balance equation from satellite techniques, with the water vapor advection term being obtained as the residue needed for balance. Thus, the intent is to develop a purely satellite-based method for obtaining the full set of terms in the atmospheric water budget equation without requiring in situ sounding information on the wind profile. The algorithm is validated by cross-checking all the algorithm components through multiple- algorithm retrieval intercomparisons. A further check on the validation is obtained by directly comparing water vapor transports into the targeted basin diagnosed from the satellite algorithms to those obtained observationally from a network of land-based upper air stations that nearly uniformly surround the basin, although it is fair to say that these checks are more effective m identifying problems in estimating vapor transports from a leaky operational radiosonde network than in verifying the transport estimates determined from the satellite algorithm system Total columnar atmospheric water budget results are presented for an extended annual cycle consisting of the months of October-97, January-98, April-98, July-98,October-98, and January 1999. These results are used to emphasize the changing relationship in E-P, as well as in the varying roles of storage and advection in balancing E-P both on daily and monthly time scales and on localized and basin space scales. Results from the algorithm-to-algorithm intercomparisons are also presented in the context of sensitivity testing to help understand the intrinsic uncertainties in evaluating the water budget terms by an all-satellite algorithm approach.

Smith, Eric A.↗

NASA NTRS2006

Simulated convective systems using a cloud resolving model: Impact of large-scale temperature and moisture forcing using observations and GEOS-3 reanalysis

The GCE (Goddard Cumulus Ensemble) model, which has been developed and improved at NASA Goddard Space Flight Center over the past two decades, is considered as one of the finer and state-of-the-art CRMs (Cloud Resolving Models) in the research community. As the chosen CRM for a NASA Interdisciplinary Science (IDS) Project, GCE has recently been successfully upgraded into an MPI (Message Passing Interface) version with which great improvement has been achieved in computational efficiency, scalability, and portability. By basically using the large-scale temperature and moisture advective forcing, as well as the temperature, water vapor and wind fields obtained from TRMM (Tropical Rainfall Measuring Mission) field experiments such as SCSMEX (South China Sea Monsoon Experiment) and KWAJEX (Kwajalein Experiment), our recent 2-D and 3-D GCE simulations were able to capture detailed convective systems typical of the targeted (simulated) regions. The GEOS-3 [Goddard EOS (Earth Observing System) Version-3] reanalysis data have also been proposed and successfully implemented for usage in the proposed/performed GCE long-term simulations (i.e., aiming at producing massive simulated cloud data -- Cloud Library) in compensating the scarcity of real field experimental data in both time and space (location). Preliminary 2-D or 3-D pilot results using GEOS-3 data have generally showed good qualitative agreement (yet some quantitative difference) with the respective numerical results using the SCSMEX observations. The first objective of this paper is to ensure the GEOS-3 data quality by comparing the model results obtained from several pairs of simulations using the real observations and GEOS-3 reanalysis data. The different large-scale advective forcing obtained from these two kinds of resources (i.e., sounding observations and GEOS-3 reanalysis) has been considered as a major critical factor in producing various model results. The second objective of this paper is therefore to investigate and present such an impact of large-scale forcing on various modeled quantities (such as hydrometeors, rainfall, and etc.). A third objective is to validate the overall GCE 3-D model performance by comparing the numerical results with sounding observations, as well as available satellite retrievals.

Shie, C.-L.↗

NASA NTRS2006

Monthly-Diurnal Water Budget Variability Over Gulf of Mexico-Caribbean Sea Basin from Satellite Observations

This study presents results from a multi-satellite/multi-sensor retrieval system design d to obtain the atmospheric water budget over the open ocean. A combination of hourly-sampled monthly datasets derived from the GOES-8 5-channel Imager, the TRMM TMI radiometer, and the DMSP 7-channel passive microwave radiometers (SSM/I) have been acquired for the combined Gulf of Mexico-Caribbean Sea basin. Whereas the methodology has been tested over this basin, the retrieval system is designed for portability to any open-ocean region. Algorithm modules using the different datasets to retrieve individual geophysical parameters needed in the water budget equation are designed in a manner that takes advantage of the high temporal resolution of the GOES-8 measurements, as well as the physical relationships inherent to the TRMM and SSM/I passive microwave measurements in conjunction with water vapor, cloud liquid water, and rainfall. The methodology consists of retrieving the precipitation, surface evaporation, and vapor-cloud water storage terms in the atmospheric water balance equation from satellite techniques, with the water vapor advection term being obtained as the residue needed for balance. Thus, the intent is to develop a purely satellite-based method for obtaining the full set of terms in the atmospheric water budget equation without requiring in situ sounding information on the wind profile. The algorithm is validated by cross-checking all the algorithm components through multiple-algorithm retrieval intercomparisons. A further check on the validation is obtained by directly comparing water vapor transports into the targeted basin diagnosed from the satellite algorithms to those obtained observationally from a network of land-based upper air stations that nearly uniformly surround the basin, although it is fair to say that these checks are more effective in identifying problems in estimating vapor transports from a "leaky" operational radiosonde network than in verifying the transport estimates determined from the satellite algorithm system. Total columnar atmospheric water budget results are presented for an extended annual cycle consisting of the months of October-97, January-98, April-98, July-98,October-98, and January- 1999. These results are used to emphasize the changing relationship in E-P, as well as in the varying roles of storage and advection in balancing E-P both on daily and monthly time scales and on localized and basin space scales. Results from the algorithm-to-algorithm intercomparisons are also presented in the context of sensitivity testing to help understand the intrinsic uncertainties in evaluating the water budget terms by an all-satellite algorithm approach.

Smith, E. A.↗

NASA NTRS2007

Stable Low Cloud Phase II: Nocturnal Event Study

This report describes the work done by the Applied Meteorology Unit (AMU) in developing a database of nights that experienced rapid (< 90 minutes) low cloud formation in a stable atmosphere, resulting in ceilings at the Shuttle Landing Facility (TTS) that violated Space Shuttle Flight Rules (FR). This work is the second phase of a similar AMU task that examined the same phenomena during the day. In the first phase of this work, the meteorological conditions favoring the rapid formation of low ceilings include the presence of any inversion below 8000 ft, high relative humidity (RH) beneath the inversion and a clockwise turning of the winds from the surface to the middle troposphere (-15000 ft). The AMU compared and contrasted the atmospheric and thermodynamic conditions between nights with rapid low ceiling formation and nights with low ceilings resulting from other mechanisms. The AMU found that there was little to discern between the rapidly-forming ceiling nights and other low ceiling nights at TTS. When a rapid development occurred, the average RH below the inversions was 87% while non-events had an average RH of 79%. One key parameter appeared to be the vertical wind profile in the Cape Canaveral, FL radiosonde (XMR) sounding. Eighty-three percent of the rapid development events had veering winds with height from the surface to the middle troposphere (-15,000 ft) while 61% of the non-events had veering winds with height. Veering winds indicate a warm-advection regime, which supports large-scale rising motion and ultimately cloud formation in a moist environment. However, only six of the nights (out of 86 events examined) with low cloud ceilings had an occurrence of rapidly developing ceilings. Since only 7% rapid development events were observed in this dataset, it is likely that rapid low cloud development is not a common occurrence during the night, or at least not as common as during the day. In the AMU work on the daytime rapid low cloud development (Case and Wheeler 2005), nearly 30% of the low cloud ceiling cases investigated were identified as rapidly developing events. Forecasters at the Space Meteorology Group (SMG) issue 30 to 90 minute forecasts for low cloud ceilings at TTS to support Space Shuttle landings. Mission verification statistics have shown ceilings to be the number one forecast challenge. More specifically, forecasters at SMG are concerned with any rapidly developing clouds/ceilings below 8000 R in a stable, capped thermodynamic environment. Therefore, the AMU was tasked to examine archived events of rapid stable cloud formation resulting in ceilings below 8000 ft, and document the atmospheric regimes favoring this type of cloud development. The AMU examined the cool season months of November to March during the years of 1994-2005 for nights that had low-level inversions and rapid, stable low cloud formation that resulted in ceilings violating the Space Shuttle FR. The AMU wrote and modified existing code to identify inversions from the evening and morning XMR radiosonde during the cool season and output pertinent sounding information. They parsed all days with cloud ceilings below 8000 ft at TTS, forming a database of possible rapidly-developing low ceiling events. Nights with precipitation or noticeable fog burn-off situations were excluded from the database. Only the nighttime hours were examined for possible ceiling development events since the daytime events were examined in the first phase of this work. The report presents one sample case of rapidly-developing low cloud ceilings. The case depicts the representative meteorological and thermodynamic characteristics of such events. The case also illustrates how quickly the cloud decks can develop, sometimes forming in 30 minutes or less. The report also summarizes the composite meteorological conditions for 6 event nights with rapid low cloud ceiling formation and 80 non-events nights consisting of advection or widespread low cloud ceilings. The teorological conditions were quite similar for both the event and non-event nights, since both types of nights experienced low cloud ceilings. Both types of nights had a relatively moist environment beneath an inversion based below 8000 ft.

Bauman, William H., III↗

NASA NTRS2007

A Stellar-mass Black Hole in the Ultra-luminous X-ray Source M82 X-1

We have analyzed the archival XMM-Newton data of the archetypal Ultra-Luminous X-ray Source (ULX) M82 X-1 with an LO5 ksec exposure when the source was in the steady state. Thanks to the high photon statistics from the large effective area and long exposure, we were able to discriminate different X-ray continuum spectral models. Neither the standard accretion disk model (where the radial dependency of the disk effective temperature is T(r) proportional to r(sup -3/4)) nor a power-law model gives a satisfactory fit. In fact, observed curvature of the M82 X-1 spectrum was just between those of the two models. When the exponent of the radial dependence (p in T(r) proportional to r(sup -P)) of the disk temperature is allowed to be free, we obtained p = 0.61 (sup +0.03)(sub -0.02). Such a reduction of p from the standard value 3/4 under extremely high mass accretion rates is predicted from the accretion disk theory as a consequence of the radial energy advection. Thus, the accretion disk in M82 X-1 is considered to be in the Slim disk state, where an optically thick Advection Dominant Accretion Flow (ADAF) is taking place. We have applied a theoretical slim disk spectral model to M82 X-1, and estimated the black hole mass approximately equal to 19 - 32 solar mass. We conclude that M82 X-1 is a stellar black hole which has been produced through evolution of an extremely massive star, shining at a several times the super-Eddington luminosity.

Okajima, Takashi↗

NASA NTRS2006

EOS Microwave Limb Sounder Observations of 'Frozen-in' Anticyclonic Air in Arctic Summer

A previously unreported phenomenon, a 'frozen-in' anticyclone (FrIAC) after the 2005 Arctic spring vortex breakup, was discovered in Earth Observing System (EOS) Microwave Limb Sounder (MLS) long-lived trace gas data. A tongue of low-latitude (high-N2O, low-H2O) air was drawn into high latitudes and confined in a tight anticyclone, then advected intact in the summer easterlies through late August. A similar feature in O3 disappeared by early April as a result of chemical processes. The FrIAC was initially advected upright at nearly the same speed at all levels from approx.660 to 1300 K (approx.25-45 km); increasing vertical wind shear after early June tilted the FrIAC and weakened it at higher levels. The associated feature in PV disappeared by early June; transport calculations fail to reproduce the remarkable persistence of the FrIAC, suggesting deficiencies in summer high-latitude winds. The historical PV record suggests that this phenomenon may have occurred several times before. The lack of a persistent signature in O3 or PV, along with its small size and rapid motion, make it unlikely that a FrIAC could have been reliably identified without hemispheric daily longlived trace gas profiles such as those from EOS MLS.

Arctic vortex↗

NASA NTRS2006

Forecast, Measurement, and Modeling of an Unprecedented Polar Ozone Filament Event over Mauna Loa Observatory, Hawaii

In mid-March 2005 the northern lower stratospheric polar vortex experienced a severe stretching episode, bringing a large polar filament far south of Alaska toward Hawaii. This meridional intrusion of rare extent, coinciding with the polar vortex final warming and breakdown, was followed by a zonal stretching in the wake of the easterly propagating subtropical main flow. This caused polar air to remain over Hawaii for several days before diluting into the subtropics. After being successfully forecasted to pass over Hawaii by the high-resolution potential vorticity advection model Modele Isentrope du transport Meso-echelle de l'Ozone Stratospherique par Advection (MIMOSA), the filament was observed on isentropic surfaces between 415 K and 455 K (17-20 km) by the Jet Propulsion Laboratory stratospheric ozone lidar measurements at Mauna Loa Observatory, Hawaii, between 16 and 19 March 2005. It was materialized as a thin layer of enhanced ozone peaking at 1.6 ppmv in a region where the climatological values usually average 1.0 ppmv. These values were compared to those obtained by the three dimensional Chemistry-Transport Model MIMOSA-CHIM. Agreement between lidar and model was excellent, particularly in the similar appearance of the ozone peak near 435 K (18.5 km) on 16 March, and the persistence of this layer at higher isentropic levels for the following three days. Passive ozone, also modeled by MIMOSA-CHIM, was at about 3-4 ppmv inside the filament while above Hawaii. A detailed history of the modeled chemistry inside the filament suggests that the air mass was still polar ozone- depleted when passing over Hawaii. The filament quickly separated from the main vortex after its Hawaiian overpass. It never reconnected and, in less than 10 days, dispersed entirely in the subtropics.

atmospheric processes↗

NASA NTRS2007

The alpha(3) Scheme - A Fourth-Order Neutrally Stable CESE Solver

The conservation element and solution element (CESE) development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To provide a solid foundation for a systematic CESE development of high order schemes, in this paper we describe a new 4th-order neutrally stable CESE solver of the advection equation Theta u/Theta + alpha Theta u/Theta x = 0. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables u(sup n) (sub j), (u(sub x))(sup n) (sub j) , and (uxz)(sup n) (sub j) (the numerical analogues of u, Theta u/Theta x, and Theta(exp 2)u/Theta x(exp 2), respectively) and four equations per mesh point, the new scheme is referred to as the alpha(3) scheme. As in the case of other similar CESE neutrally stable solvers, the alpha(3) scheme enforces conservation laws in space-time locally and globally, and it has the basic, forward marching, and backward marching forms. These forms are equivalent and satisfy a space-time inversion (STI) invariant property which is shared by the advection equation. Based on the concept of STI invariance, a set of algebraic relations is developed and used to prove that the alpha(3) scheme must be neutrally stable when it is stable. Moreover it is proved rigorously that all three amplification factors of the alpha(3) scheme are of unit magnitude for all phase angles if |v| <= 1/2 (v = alpha delta t/delta x). This theoretical result is consistent with the numerical stability condition |v| <= 1/2. Through numerical experiments, it is established that the alpha(3) scheme generally is (i) 4th-order accurate for the mesh variables u(sup n) (sub j) and (ux)(sup n) (sub j); and 2nd-order accurate for (uxx)(sup n) (sub j). However, in some exceptional cases, the scheme can achieve perfect accuracy aside from round-off errors.

Chang, Sin-Chung↗

NASA NTRS2005

The Closure of the Ocean Mixed Layer Temperature Budget using Level-Coordinate Model Fields

Entrainment is an important element of the mixed layer mass, heat, and temperature budgets. Conventional procedures to estimate entrainment heat advection often do not permit the closure of heat and temperature budgets because of inaccuracies in its formulation. In this study a rigorous approach to evaluate the effect of entrainment using the output of a general circulation model (GCM) that does not have an explicit prognostic mixed layer model is described. The integral elements of the evaluation are 1) the rigorous estimates of the temperature difference between mixed layer water and entrained water at each horizontal grid point, 2) the formulation of the temperature difference such that the budget closes over a volume greater than one horizontal grid point, and 3) the apparent warming of the mixed layer during the mixed layer shoaling to account for the weak vertical temperature gradient within the mixed layer. This evaluation of entrainment heat advection is compared with the estimates by other commonly used ad hoc formulations by applying them in three regions: the north-central Pacific, the Kuroshio Extension, and the Nino-3 areas in the tropical Pacific. In all three areas the imbalance in the mixed layer temperature budget by the ad hoc estimates is significant, reaching a maximum of about 4 K yr(exp -1).

temperature gradient↗

NASA NTRS2007

Multi-Layer Arctic Mixed-Phase Clouds Simulated by a Cloud-Resolving Model: Comparison with ARM Observations and Sensitivity Experiments

A cloud-resolving model (CRM) is used to simulate the multiple-layer mixed-phase stratiform (MPS) clouds that occurred during a three-and-a-half day subperiod of the Department of Energy-Atmospheric Radiation Measurement Program s Mixed-Phase Arctic Cloud Experiment (M-PACE). The CRM is implemented with an advanced two-moment microphysics scheme, a state-of-the-art radiative transfer scheme, and a complicated third-order turbulence closure. Concurrent meteorological, aerosol, and ice nucleus measurements are used to initialize the CRM. The CRM is prescribed by time-varying large-scale advective tendencies of temperature and moisture and surface turbulent fluxes of sensible and latent heat. The CRM reproduces the occurrences of the single- and double-layer MPS clouds as revealed by the M-PACE observations. However, the simulated first cloud layer is lower and the second cloud layer thicker compared to observations. The magnitude of the simulated liquid water path agrees with that observed, but its temporal variation is more pronounced than that observed. As in an earlier study of single-layer cloud, the CRM also captures the major characteristics in the vertical distributions and temporal variations of liquid water content (LWC), total ice water content (IWC), droplet number concentration and ice crystal number concentration (nis) as suggested by the aircraft observations. However, the simulated mean values differ significantly from the observed. The magnitude of nis is especially underestimated by one order of magnitude. Sensitivity experiments suggest that the lower cloud layer is closely related to the surface fluxes of sensible and latent heat; the upper cloud layer is probably initialized by the large-scale advective cooling/moistening and maintained through the strong longwave (LW) radiative cooling near the cloud top which enhances the dynamical circulation; artificially turning off all ice-phase microphysical processes results in an increase in LWP by a factor of 3 due to interactions between the excessive LW radiative cooling and extra cloud water; heating caused by phase change of hydrometeors could affect the LWC and cloud top height by partially canceling out the LW radiative cooling. It is further shown that the resolved dynamical circulation appears to contribute more greatly to the evolution of the MPS cloud layers than the parameterized subgrid-scale circulation.

Luo, Yali↗

NASA NTRS2006

Thermal Dispersion Within a Porous Medium Near a Solid Wall

The regenerator is a key component to Stirling cycle machine efficiency. Typical regenerators are of sintered fine wires or layers of fine-wire screens. Such porous materials are contained within solid-waH casings. Thermal energy exchange between the regenerator and the casing is important to cycle performance for the matrix and casing would not have the same axial temperature profile in an actual machine. Exchange from one to the other may allow shunting of thermal energy, reducing cycle efficiency. In this paper, temperature profiles within the near-wall region of the matrix are measured and thermal energy transport, termed thermal dispersion, is inferred. The data show how the wall affects thermal transport. Transport normal to the mean flow direction is by conduction within the solid and fluid and by advective transport within the matrix. In the near-wall region, both may be interrupted from their normal in-core pattern. Solid conduction paths are broken and scales of advective transport are damped. An equation is presented which describes this change for a wire screen mesh. The near-wall layer typically acts as an insulating layer. This should be considered in design or analysis. Effective thermal conductivity within the core is uniform. In-core transverse thermal effective conductivity values are compared to direct and indirect measurements reported elsewhere and to 3D numerical simulation results, computed previously and reported elsewhere. The 3-D CFD model is composed of six cylinders in cross flow, staggered in arrangement to match the dimensions and porosity of the matrix used in the experiments. The commercial code FLUENT is used to obtain the flow and thermal fields. The thermal dispersion and effective thermal conductivities for the matrix are computed from the results.

Simon, T.↗

NASA NTRS2010

Modeling the Frozen-In Anticyclone in the 2005 Arctic Summer Stratosphere

Immediately following the breakup of the 2005 Arctic spring stratospheric vortex, a tropical air mass, characterized by low potential vorticity (PV) and high nitrous oxide (N2O), was advected poleward and became trapped in the easterly summer polar vortex. This feature, known as a "Frozen-In Anticyclone (FrIAC)", was observed in Earth Observing System (EOS) Aura Microwave Limb Sounder (MLS) data to span the potential temperature range from approximately 580 to 1100 K (approximately 25 to 40 km altitude) and to persist from late March to late August 2005. This study compares MLS N2O observations with simulations from the Global Modeling Initiative (GMI) chemistry and transport model, the GEOS-5/MERRA Replay model, and the VanLeer Icosahedral Triangular Advection isentropic transport model to elucidate the processes involved in the lifecycle of the FrIAC which is here divided into three distinct phases. During the "spin-up phase" (March to early April), strong poleward flow resulted in a tight isolated anticyclonic vortex at approximately 70-90 deg N, marked with elevated N2O. GMI, Replay, and VITA all reliably simulted the spin-up of the FrIAC, although the GMI and Replay peak N2O values were too low. The FrIAC became trapped in the developing summer easterly flow and circulated around the polar region during the "anticyclonic phase" (early April to the end of May). During this phase, the FrIAC crossed directly over the pole between the 7th and 14th of April. The VITA and Replay simulations transported the N2O anomaly intact during this crossing, in agreement with MLS, but unrealistic dispersion of the anomaly occurred in the GMI simulation due to excessive numerical mixing of the polar cap. The vortex associated with the FrIAC was apparently resistant to the weak vertical hear during the anticyclonic phase, and it thereby protected the embedded N20 anomaly from stretching. The vortex decayed in late May due to diabatic processes, leaving the N2O anomaly exposed to horizontal and vertical wind shears during the "shearing phase" (June to August). The observed lifetime of the FrIAC during this phase is consistent with time-scales calculated from the ambient horizontal and vertical wind shear. Replay maintained the horizontal structure of the N2O anomaly similar to NILS well into August. The VITA simulation also captured the horizontal structure of the FrIAC during this phase, but VITA eventually developed fine-scale N2O structure not observed in MLS data.

Allen, D. R.↗

NASA NTRS2011

Long-Term Evolution of the Aerosol Debris Cloud Produced by the 2009 Impact on Jupiter

We present a study of the long-term evolution of the cloud of aerosols produced in the atmosphere of Jupiter by the impact of an object on 19 July 2009. The work is based on images obtained during 5 months from the impact to 31 December 2009 taken in visible continuum wavelengths and from 20 July 2009 to 28 May 2010 taken in near-infrared deep hydrogen-methane absorption bands at 2.1-2.3 micron. The impact cloud expanded zonally from approximately 5000 km (July 19) to 225,000 km (29 October, about 180 deg in longitude), remaining meridionally localized within a latitude band from 53.5 deg S to 61.5 deg S planetographic latitude. During the first two months after its formation the site showed heterogeneous structure with 500-1000 km sized embedded spots. Later the reflectivity of the debris field became more homogeneous due to clump mergers. The cloud was mainly dispersed in longitude by the dominant zonal winds and their meridional shear, during the initial stages, localized motions may have been induced by thermal perturbation caused by the impact's energy deposition. The tracking of individual spots within the impact cloud shows that the westward jet at 56.5 deg S latitude increases its eastward velocity with altitude above the tropopause by 5- 10 m/s. The corresponding vertical wind shear is low, about 1 m/s per scale height in agreement with previous thermal wind estimations. We found evidence for discrete localized meridional motions with speeds of 1-2 m/s. Two numerical models are used to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the winds and their shears. The other uses the EPIC code, a nonlinear calculation of the evolution of the potential vorticity field generated by a heat pulse that simulates the impact. Both models reproduce the observed global structure of the cloud and the dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. The reflectivity of the impact cloud decreased exponentially with a characteristic timescale of 15 days; we can explain this behavior with a radiative transfer model of the cloud optical depth coupled to an advection model of the cloud dispersion by the wind shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) for the small aerosol particles forming the cloud is 45-200 days, thus aerosols were removed vertically over the long term following their zonal dispersion. No evidence of the cloud was detected 10 months after the impact.

Sanchez-Lavega, A.↗

NASA NTRS2012

Photospheric Magnetic Flux Transport - Supergranules Rule

Observations of the transport of magnetic flux in the Sun's photosphere show that active region magnetic flux is carried far from its origin by a combination of flows. These flows have previously been identified and modeled as separate axisymmetric processes: differential rotation, meridional flow, and supergranule diffusion. Experiments with a surface convective flow model reveal that the true nature of this transport is advection by the non-axisymmetric cellular flows themselves - supergranules. Magnetic elements are transported to the boundaries of the cells and then follow the evolving boundaries. The convective flows in supergranules have peak velocities near 500 m/s. These flows completely overpower the superimposed 20 m/s meridional flow and 100 m/s differential rotation. The magnetic elements remain pinned at the supergranule boundaries. Experiments with and without the superimposed axisymmetric photospheric flows show that the axisymmetric transport of magnetic flux is controlled by the advection of the cellular pattern by underlying flows representative of deeper layers. The magnetic elements follow the differential rotation and meridional flow associated with the convection cells themselves -- supergranules rule!

Hathaway, David H.↗

NASA NTRS2012

Closing the Seasonal Ocean Surface Temperature Balance in the Eastern Tropical Oceans from Remote Sensing and Model Reanalyses

The Eastern tropical ocean basins are regions of significant atmosphere-ocean interaction and are important to variability across subseasonal to decadal time scales. The numerous physical processes at play in these areas strain the abilities of coupled general circulation models to accurately reproduce observed upper ocean variability. Furthermore, limitations in the observing system of important terms in the surface temperature balance (e.g., turbulent and radiative heat fluxes, advection) introduce uncertainty into the analyses of processes controlling sea surface temperature variability. This study presents recent efforts to close the surface temperature balance through estimation of the terms in the mixed layer temperature budget using state-of-the-art remotely sensed and model-reanalysis derived products. A set of twelve net heat flux estimates constructed using combinations of radiative and turbulent heat flux products - including GEWEX-SRB, ISCCP-SRF, OAFlux, SeaFlux, among several others - are used with estimates of oceanic advection, entrainment, and mixed layer depth variability to investigate the seasonal variability of ocean surface temperatures. Particular emphasis is placed on how well the upper ocean temperature balance is, or is not, closed on these scales using the current generation of observational and model reanalysis products. That is, the magnitudes and spatial variability of residual imbalances are addressed. These residuals are placed into context within the current uncertainties of the surface net heat fluxes and the role of the mixed layer depth variability in scaling the impact of those uncertainties, particularly in the shallow mixed layers of the Eastern tropical ocean basins.

Roberts, J. Brent↗

NASA NTRS2014

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗

NASA NTRS2013

Natural Air-Sea Flux of CO2 in Simulations of the NASA-GISS Climate Model: Sensitivity to the Physical Ocean Model Formulation

Results from twin control simulations of the preindustrial CO2 gas exchange (natural flux of CO2) between the ocean and the atmosphere are presented here using the NASA-GISS climate model, in which the same atmospheric component (modelE2) is coupled to two different ocean models, the Russell ocean model and HYCOM. Both incarnations of the GISS climate model are also coupled to the same ocean biogeochemistry module (NOBM) which estimates prognostic distributions for biotic and abiotic fields that influence the air-sea flux of CO2. Model intercomparison is carried out at equilibrium conditions and model differences are contrasted with biases from present day climatologies. Although the models agree on the spatial patterns of the air-sea flux of CO2, they disagree on the strength of the North Atlantic and Southern Ocean sinks mainly because of kinematic (winds) and chemistry (pCO2) differences rather than thermodynamic (SST) ones. Biology/chemistry dissimilarities in the models stem from the different parameterizations of advective and diffusive processes, such as overturning, mixing and horizontal tracer advection and to a lesser degree from parameterizations of biogeochemical processes such as gravitational settling and sinking. The global meridional overturning circulation illustrates much of the different behavior of the biological pump in the two models, together with differences in mixed layer depth which are responsible for different SST, DIC and nutrient distributions in the two models and consequently different atmospheric feedbacks (in the wind, net heat and freshwater fluxes into the ocean).

Gas exchange↗