Engineering Papers⌕ Search

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.

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

At least 523 records · Page 29

Use of Maps, GTE and UARS-MLS Data in Understanding Tropospheric Processes Critical to Model Development

Two of the major problems in using meteorological models to explain observed tropospheric trace constituent distributions and thereby to understand the global budgets of the tracers are to properly define the vertical layered structure in the free atmosphere, and to understand the contribution of advection processes in generating horizontal inhomogeneities at all scales. We proposed to tackle these problems through the examination of an extensive collection of trace constituent data from research and commercial aircraft in conjunction with meteorological data from the European Center for Medium-Range Weather Forecasts. The physical mechanisms responsible for these advection and layering processes were explored and their implications for theories and models assessed. In addition, we calculated examples of how thin layers (not currently resolved by models) affect the radiative heating/cooling rates. We developed an improved algorithm for trace constituent layer detection, and used it to analyze over 100,000 km of ozone and humidity vertical profiles collected by instruments piggybacked on commercial aircraft. The same method was also used to examine ozone, humidity, carbon monoxide, and methane data from the NASA Pacific Exploratory Missions. The major conclusions from these studies were that tropospheric trace constituent layers are ubiquitous, and that their characteristics are remarkably universal. These results support the notion that ozone layers coming in from the stratosphere play a large role in creating tropospheric layers, whether by themselves or by capping buoyant pollution plumes. This, in turn, has important implications for the transport (especially in the vertical) of trace gases and how it is represented in models.

Newell, Reginald E.↗

A Study of Meteorological Conditions Associated With Noctilucent Clouds

Temperature measurements were obtained in the upper stratosphere and mesosphere between 50 and 95 km with passive inflatable falling spheres launched on small meteorological rockets as part of the DROPPS (Distribution and Role of Particles in the Polar Summer Mesosphere) program. Temperatures of the neutral atmosphere have been combined with similar measurements obtained during 1991 and 1993. Temperatures were found to change monatonically with altitude except during the Nocticulent Clouds (NLC) occurrences during DROPPS. The temperature lapse rate changed between 5 July 1999, 2313 UTC and 6 July 1999, 0209 UTC; this included a lowering of the altitude of minimum temperature by about 5 km. Furthermore, winds backed from a northeasterly direction to a northwesterly direction. Whether the change in temperature observed is a result of advection related to the changes of the wind field due to advection. Comparisons will also concentrate on the meteorological conditions during the NLC event during DROPPS and earlier 1991 and 1993 NLC'S.

Schmidlin, F. J.↗

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.↗

Thermocapillary Flows with Low Frequency g-Jitter

A thermocapillary parallel flow is established in a fluid filled slot with an applied temperature gradient. Low frequency jitter is imposed in arbitrary directions. Vertical jitter proves to be relatively uninteresting, merely augmenting or opposing the basic thermocapillary flow. Streamwise jitter still produces parallel flows, but these now exhibit boundary layers or layered cellular structures for large Rayleigh number as the applied stratification alternates between stable and unstable. Runaways are possible for unstable stratification and these correspond to resonant excitation of stationary long wave Rayleigh-Benard modes. Spanwise jitter produces fully three dimensional motion. A spanwise-streamwise circulation results for weak spanwise jitter, which advects the interfacial temperature establishing a subsidiary spanwise thermocapillary flow. This flow is strong at small Biot number when advected temperature is trapped in the slot, and has a counter-intuitive dependence on the spanwise-streamwise aspect ratio.

Grassia, P.↗

Satellite Retrieval of Atmospheric Water Budget over Gulf of Mexico- Caribbean Basin: Seasonal Variability

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 hourly-sampled monthly datasets derived from the GOES-8 5 Imager and the DMSP 7-channel passive microwave radiometer (SSM/I) have been acquired for the Gulf of Mexico-Caribbean Sea basin. Whereas the methodology is being 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 SSM/I passive microwave signals 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, we have sought 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 partly validated by first cross-checking all the algorithm components through multiple-algorithm retrieval intercomparisons. More fundamental validation is obtained by directly comparing water vapor transports into the targeted basin diagnosed from the satellite algorithm to those obtained observationally from a network of land-based upper air stations that nearly uniformly surround the basin. Total columnar atmospheric water budget results will be 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 will also be presented in the context of sensitivity testing to help understand the intrinsic uncertainties in the water budget terms.

Smith, Eric A.↗

A Study of the Response of Deep Tropical Clouds to Mesoscale Processes: Modeling Strategies and Simulations of TOGA-COARE Convective Systems - Part 1

Interactions between deep tropical clouds over the western Pacific warm pool and the larger-scale environment are key to understanding climate change. Cloud models are an extremely useful tool in simulating and providing statistical information on heat and moisture transfer processes between cloud systems and the environment, and can therefore be utilized to substantially improve cloud parameterizations in climate models. In this paper, the Goddard Cumulus Ensemble (GCE) cloud-resolving model is used in multi-day simulations of deep tropical convective activity over the Tropical Ocean-Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). Large-scale temperature and moisture advective tendencies, and horizontal momentum from the TOGA-COARE Intensive Flux Array (IFA) region, are applied to the GCE version which incorporates cyclical boundary conditions. Sensitivity experiments show that grid domain size produces the largest response to domain-mean temperature and moisture deviations, as well as cloudiness, when compared to grid horizontal or vertical resolution, and advection scheme. It is found that a minimum grid-domain size of 500 km is needed to adequately resolve the convective cloud features. The control experiment shows that the atmospheric heating and moistening is primarily a response to cloud latent processes of condensation/evaporation, and deposition/sublimation, and to a lesser extent, melting of ice particles. Air-sea exchange of heat and moisture is found to be significant, but of secondary importance, while the radiational response is small. The simulated rainfall and atmospheric heating and moistening, agrees well with observations, and performs favorably to other models simulating this case.

Johnson, Daniel E.↗

On Animating 2D Velocity Fields

A velocity field. even one that represents a steady state flow implies a dynamical system. Animated velocity fields is an important tool in understanding such complex phenomena. This paper looks at a number of techniques that animate velocity fields and propose two new alternatives, These are texture advection and streamline cycling. The common theme among these techniques is the use of advection on some texture to generate a realistic animation of the velocity field. Texture synthesis and selection for these methods are presented. Strengths and weaknesses of the techniques are also discussed in conjunction with several examples.

Kao, David↗

On Animating 2D Velocity Fields

A velocity field, even one that represents a steady state flow, implies a dynamical system. Animated velocity fields is an important tool in understanding such complex phenomena. This paper looks at a number of techniques that animate velocity fields and propose two new alternatives. These are texture advection and streamline cycling. The common theme among these techniques is the use of advection on some texture to generate a realistic animation of the velocity field. Texture synthesis and selection for these methods are presented. Strengths and weaknesses of the techniques are also discussed in conjunctions with several examples.

Kao, David↗

Algebraic Nonoverlapping Domain Decomposition Methods for Stabilized FEM and FV Discretizations

We consider preconditioning methods for convection dominated fluid flow problems based on a nonoverlapping Schur complement domain decomposition procedure for arbitrary triangulated domains. The triangulation is first partitioned into a number of subdomains and interfaces which induce a natural 2 x 2 partitioning of the p.d.e. discretization matrix. We view the Schur complement induced by this partitioning as an algebraically derived coarse space approximation. This avoids the known difficulties associated with the direct formation of an effective coarse discretization for advection dominated equations. By considering various approximations of the block factorization of the 2 x 2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 using MPI message passing protocol. A number of 2-D CFD calculations will be presented for both scalar advection-diffusion equations and the Euler equations discretized using stabilized finite element and finite volume methods. These results show very good scalability of the preconditioner for various discretizations as the number of processors is increased while the number of degrees of freedom per processor is fixed.

Barth, Timothy J.↗

Are Stronger North-Atlantic Southwesterlies the Forcing to the Late-Winter Warming in Europe?

We examine a possible mechanism leading to late-winter warming and thus to an early spring in Europe. From the NCEP Reanalysis, we extract for the years 1948-1999 ocean-surface winds over the eastern North Atlantic, and air temperatures at the surface, T(sub s), and at the 500 mb level, T(sub 500), in late-winter and spring. T(sub s) is extracted at six European locations, all at 50.5 N, ranging in longitude from 1.9 E (northeastern France) to 26.2 E (Ukraine). To quantify the advection of maritime air into Europe, we evaluate for 3-pentad groups the Index I(sub na) of the southwesterlies at 45 N; 20 W: I(sub na) is the average wind speed at this point if the direction is from the quadrant 180-270 deg (when the direction is different, the contribution counts as zero). In late winter correlations C(sub it) between the Index I(sub na) and the temperature T(sub s) are substantial, up to the 0.6 level, in western Europe (but weaker correlations for Poland and Ukraine). C(sub it) drops sharply by mid-March, taking occasionally negative values subsequently. This drop in C(sub it) indicates that maritime air advection is no longer associated closely with the surface-air warming, the role of immolation becomes important, and thus the drop in C(sub it) marks the arrival of spring. Correlations C(sub i delta) between I(sub na) and our lapse-rate parameter delta, the difference between T(sub s) and T(sub 500), indicate that the flow of warm maritime-air from the North Atlantic into this 'corridor' at 50.5 N is predominantly at low tropospheric level. By computing the best linear fit to I(sub na) and T(sub s), the trends for the period 1948-1999 are evaluated. The trends are appreciable in the second half of February and the first half of March. Our 3-pentad analysis points to the interval from mid-February to mid-March as the end-of-winter period in which the southwesterlies over the eastern North Atlantic become stronger and the surface-air temperatures in Europe rise markedly, the lapse rate becomes steeper, and concurrently the longitudinal temperature gradient between the Somme (France) and the Oder (Germany/Poland border) is reduced by 0.8 C, that is, by 20% of its 1948 value. Our thesis, that the observed late-winter warming and the corollary advancement of spring in Europe resulted at least in part from stronger southwesterlies over the North Atlantic, merits further investigations.

Otterman, J.↗

Lidar Observations of Tropospheric Aerosols Over Northeastern South Africa During the ARREX and SAFARI-2000 Dry Season Experiments

During the ARREX-1999 and SAFARI-2000 Dry Season experiments a micropulse lidar (523 nm) instrument was operated at the Skukuza Airport in northeastern South Africa. The Mar was collocated with a diverse array of passive radiometric equipment. For SAFARI-2000 the processed Mar data yields a daytime time-series of layer mean/derived aerosol optical properties, including extinction-to-backscatter ratios and vertical extinction cross-section profile. Combined with 523 run aerosol optical depth and spectral Angstrom exponent calculations from available CIMEL sun-photometer data and normalized broadband flux measurements the temporal evolution of the near surface aerosol layer optical properties is analyzed for climatological trends. For the densest smoke/haze events the extinction-to-backscatter ratio is found to be between 60-80/sr, and corresponding Angstrom exponent calculations near and above 1.75. The optical characteristics of an evolving smoke event from SAFARI-2000 are extensively detailed. The advecting smoke was embedded within two distinct stratified thermodynamic layers, causing the particulate mass to advect over the instrument array in an incoherent manner on the afternoon of its occurrence. Surface broadband flux forcing due to the smoke is calculated, as is the evolution in the vertical aerosol extinction profile as measured by the Han Finally, observations of persistent elevated aerosol during ARREX-1999 are presented and discussed. The lack of corroborating observations the following year makes these observation; both unique and noteworthy in the scope of regional aerosol transport over southern Africa.

Campbell, James R.↗

Anomalous Heat Budgets in the Interior Pacific Ocean on Seasonal- to -Timescales and Gyre Spacescales

This study quantifies uncertainties in closing the seasonal cycle of diabatic heat storage over the Pacific Ocean from 20 degrees S to 60 degrees N through the synthesis of World Ocean Circulation Experiment (WOCE) products over 7 years from 1993-1999. We utilize WOCE reanalysis products from the following sources: diabatic heat storage (DHS) from the Scripps Institution of Oceanography (SIO); near-surface geostrophic and Ekman currents from the Earth and Space Research (ESR); and air-sea heat fluxes from Comprehensive Ocean-Atmosphere Data Set (COADS), National Centers for Environmental Prediction (NCEP), and European Center for Mid-Range Weather Forecasts (ECMWF). We interpolate these products onto a common grid, allowing the seasonal cycle of DHS to be modeled for comparison with that observed. Everywhere latent heat flux residuals dominate sensible heat flux residuals and shortwave heat flux residuals dominate longwave heat flux residuals, both comparable in magnitude to the residual horizontal heat advection. We find the root-mean-square (RMS) of the differences between observed and model residual DHS tendencies to be less than 15 W per square meters everywhere except in the Kuroshio extension. Comparable COADS and NCEP products perform better than ECMWF products in the extra-tropics, while the NCEP product performs best in the tropics. Radiative and turbulent air-sea heat flux residuals computed from ship-born measurements perform better than those computed from satellite cloud and wind measurements. Since the RMS differences derive largely from biases in measured wind speed and cloud fraction, least-squares minimization is used to correct the residual Ekman heat advection and air-sea heat flux. Minimization reduces RMS differences less than 5 W per square meters except in the Kuroshio extension, suggesting how winds, clouds, and exchange coefficients in the NCEP, ECMWF, and ESR products can be improved.

White, Warren↗

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

CODA: Coordinated Observations of Dynamics in Aurora

The grant funds were provided to carry out chemical tracer wind and turbulence measurements in the mesosphere and lower thermosphere as part of the Coordinated Observations of Dynamics in Aurora (CODA) experiment for which Dr. Andy Christensen of the Aerospace Corp. was the Principal Investigator. Clemson University was a co-investigator institution on the project and was tasked with designing, building, and testing a total of four chemical tracer release payloads which were used to release trimethyl aluminum (TMA) to measure the winds and turbulent diffusion. All four payloads were flown from the rocket range at Poker Flat, Alaska, in support of instrumented payloads flown from the same location. The CODA experiment was designed to investigate the changes in the chemical composition in the lower thermosphere during an auroral substorm event. In particular, the objective was to determine the relative importance of two competing mechanisms. The first mechanism is associated with the advection of air from the polar cap into the auroral oval. A flow in that direction is typical in the postmidnight sector of the oval. As the air moves from the region of negligible forcing in the polar cap to the region of strong forcing and heating in the oval, changes in the composition are expected to occur. Alternatively, enhanced local mixing will also cause such changes. Since the auroral forcing accelerates the neutral winds in the lower E region, but not at lower altitudes, the forcing results in large shears that can be large enough to be unstable in the Richardson number sense. The unstable shears are expected to generate turbulence and therefore enhanced mixing. Such mixing will also create composition changes similar to what is observed. In the CODA experiment an instrumented payload designed and built by the Aerospace Corporation was used to measure the composition and plasma parameters needed to define the auroral forcing. In addition, three chemical tracer payloads were to be flown along three trajectories widely separated in azimuth to map out the winds and wind gradients. The measurements would thus provide information about the local composition changes above the launch site as well as detailed information about the advection and location of stable or unstable shears both locally and over a broader area in the vicinity of the launch site.

Larsen, M. F.↗

Various Numerical Applications on Tropical Convective Systems Using a Cloud Resolving Model

In recent years, increasing attention has been given to cloud resolving models (CRMs or cloud ensemble models-CEMs) for their ability to simulate the radiative-convective system, which plays a significant role in determining the regional heat and moisture budgets in the Tropics. The growing popularity of CRM usage can be credited to its inclusion of crucial and physically relatively realistic features such as explicit cloud-scale dynamics, sophisticated microphysical processes, and explicit cloud-radiation interaction. On the other hand, impacts of the environmental conditions (for example, the large-scale wind fields, heat and moisture advections as well as sea surface temperature) on the convective system can also be plausibly investigated using the CRMs with imposed explicit forcing. In this paper, by basically using a Goddard Cumulus Ensemble (GCE) model, three different studies on tropical convective systems are briefly presented. Each of these studies serves a different goal as well as uses a different approach. In the first study, which uses more of an idealized approach, the respective impacts of the large-scale horizontal wind shear and surface fluxes on the modeled tropical quasi-equilibrium states of temperature and water vapor are examined. In this 2-D study, the imposed large-scale horizontal wind shear is ideally either nudged (wind shear maintained strong) or mixed (wind shear weakened), while the minimum surface wind speed used for computing surface fluxes varies among various numerical experiments. For the second study, a handful of real tropical episodes (TRMM Kwajalein Experiment - KWAJEX, 1999; TRMM South China Sea Monsoon Experiment - SCSMEX, 1998) have been simulated such that several major atmospheric characteristics such as the rainfall amount and its associated stratiform contribution, the Qlheat and Q2/moisture budgets are investigated. In this study, the observed large-scale heat and moisture advections are continuously applied to the 2-D model. The modeled cloud generated from such an approach is termed continuously forced convection or continuous large-scale forced convection. A third study, which focuses on the respective impact of atmospheric components on upper Ocean heat and salt budgets, will be presented in the end. Unlike the two previous 2-D studies, this study employs the 3-D GCE-simulated diabatic source terms (using TOGA COARE observations) - radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the Ocean mixed-layer (OML) model.

Shie, C.-L.↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗

Nonlinear Gulf Stream Interaction with the Deep Western Boundary Current System: Observations and a Numerical Simulation

Gulf Stream (GS) separation near its observed Cape Hatteras (CH) separation location, and its ensuing path and dynamics, is a challenging ocean modeling problem. If a model GS separates much farther north than CH, then northward GS meanders, which pinch off warm core eddies (rings), are not possible or are strongly constrained by the Grand Banks shelfbreak. Cold core rings pinch off the southward GS meanders. The rings are often re-absorbed by the GS. The important warm core rings enhance heat exchange and, especially, affect the northern GS branch after GS bifurcation near the New England Seamount Chain. This northern branch gains heat by contact with the southern branch water upstream of bifurcation, and warms the Arctic Ocean and northern seas, thus playing a major role in ice dynamics, thermohaline circulation and possible global climate warming. These rings transport heat northward between the separated GS and shelf slope/Deep Western Boundary Current system (DWBC). This region has nearly level time mean isopycnals. The eddy heat transport convergence/divergence enhances the shelfbreak and GS front intensities and thus also increases watermass transformation. The fronts are maintained by warm advection by the Florida Current and cool advection by the DWBC. Thus, the GS interaction with the DWBC through the intermediate eddy field is climatologically important.

Dietrich, David E.↗

Computations of Wall Distances Based on Differential Equations

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Tucker, Paul G.↗