Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boussinesq”

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 181 records · Page 10

Numerical Modeling of Crystal of ZnSe by Physical Vapor Transport - Towards a more Comprehensive Formulations

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of PVT growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided. The effect of gravity on the process from numerical computations shows subtle effects although experimental evidence from vertically and horizontally grown samples show dramatic evidence of gravitational effects. The shortcomings of the problem formulation will be discussed and a framework will be provided leading up towards a more comprehensive model of PVT systems.

Ramachandran, N.↗

Crystal Growth of ZnSe by Physical Vapor Transport: A Modeling Study

Crystal growth from the vapor phase has various advantages over melt growth. The main advantage is from a lower processing temperature which makes the process more amenable in instances where the melting temperature of the crystal is high. Other benefits stem from the inherent purification mechanism in the process due to differences in the vapor pressures of the native elements and impurities, and the enhanced interfacial morphological stability during the growth process. Further, the implementation of Physical Vapor Transport (PVT) growth in closed ampoules affords experimental simplicity with minimal needs for complex process control which makes it an ideal candidate for space investigations in systems where gravity tends to have undesirable effects on the growth process. Bulk growth of wide band gap II-VI semiconductors by physical vapor transport has been developed and refined over the past several years at NASA MSFC. Results from a modeling study of PVT crystal growth of ZnSe arc reported in this paper. The PVI process is numerically investigated using both two-dimensional and fully three-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and the compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The preliminary results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles will be provided in the final manuscript along with computed heat and mass transfer rates. Comparisons with the 1-D model will also be provided.

Ramachandran, Narayanan↗

The Topology of Three-Dimensional Symmetric Tensor Fields

We study the topology of 3-D symmetric tensor fields. The goal is to represent their complex structure by a simple set of carefully chosen points and lines analogous to vector field topology. The basic constituents of tensor topology are the degenerate points, or points where eigenvalues are equal to each other. First, we introduce a new method for locating 3-D degenerate points. We then extract the topological skeletons of the eigenvector fields and use them for a compact, comprehensive description of the tensor field. Finally, we demonstrate the use of tensor field topology for the interpretation of the two-force Boussinesq problem.

Lavin, Yingmei↗

Unique Aeolian Transport Mechanisms on Mars: Respective Roles of Percussive and Repercussive Grain Populations in the Sediment Load

Experiments show that when sand-size grains impact a sediment surface with energy levels commensurate for Mars, small craters are formed by the ejection of several hundred grains from the bed. The experiments were conducted with a modified crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism. Individual grains of sand could be fired at loose sand targets to observe ballistic effects unhindered by aerodynamic mobilization of the bed. Impact trajectories simulated the saltation process on dune surfaces. Impact craters were not elongated despite glancing (15 deg.) bed impact; the craters were very close to being circular. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 um-diameter grains into similar material. This behavior is explained by deposition of elastic energy in the bed by the "percussive" grain. Impact creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains to open packing and they consequently become forcefully ejected from the site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. A stress model based on repercussive bed dilatancy and interparticle adhesive forces (for smaller grains) predicts, to first order, the observed crater volumes for various impact conditions. On earth, only a few grains are mobilized by a percussive saltating grain; some grains are "knudged" along the ground, and some are partly expelled on short trajectories. These motions constitute reptation transport. On Mars, saltation and reptation become indistinct: secondary or "repercussive" trajectories have sufficient vertical impulse to create a dense saltation population of many tens or hundreds of grains for each single high-speed saltation percussion of the bed. Impact cascading will lead to near-surface distortion of the boundary layer, and choked flow formed by a dense "slurry" of sand, with the majority of grains mobilized by repercussive forces rather than by aerodynamic lift. This proceeds until a fully-matured transport layer imposes self-limitations as grain-population density constrains the free-path motion of individial grains.

Marshall, John R.↗

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗

Unique Aeolian Transport Mechanisms on Mars: Respective Roles of Percussive and Repercussive Grain Populations in the Sediment Load

Experiments show that when sand-size grains impact a sediment surface with energy levels commensurate for Mars, small craters are formed by the ejection of several hundred grains from the bed. The experiments were conducted with a modified crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism. Individual grains of sand could be fired at loose sand targets to observe ballistic effects unhindered by aerodynamic mobilization of the bed. Impact trajectories simulated the saltation process on dune surfaces. Impact craters were not elongated despite glancing (15 deg.) bed impact; the craters were very close to being circular. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. This behavior is explained by deposition of elastic energy in the bed by the "percussive" grain. Impact creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains to open packing and they consequently become forcefully ejected from the site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. A stress model based on repercussive bed dilatancy and interparticle adhesive forces (for smaller grains) predicts, to first order, the observed crater volumes for various impact conditions. On earth, only a few grains are mobilized by a percussive saltating grain; some grains are "knudged" along the ground, and some are partly expelled on short trajectories. These motions constitute reptation transport. On Mars, saltation and reptation become indistinct: secondary or "repercussive" trajectories have sufficient vertical impulse to create a dense saltation population of many tens or hundreds of grains for each single high-speed saltation percussion of the bed. Impact cascading will lead to near-surface distortion of the boundary layer, and choked flow formed by a dense "slurry" of sand, with the majority of grains mobilized by repercussive forces rather than by aerodynamic lift. This proceeds until a fully-matured transport layer imposes self- limitations as grain-population density constrains the free-path motion of individual grains.

Marshall, John R.↗

Diagnostic Clast-Texture Criteria for Recognition of Impact Deposits

It is difficult to find definitive evidence for impact in the geological record because there are many endogenous geological processes that can produce diamictites similar to those generated by impact ejecta. The classic impact criteria of shock fabrics in certain minerals, and iridium layers, for example, may be either difficult to find, or long-since erased from the impact site (shock fabrics also anneal with time). It is important to be able to recognize impact-generated materials in order to understand earth's crustal development and biological evolution. In future exploration of Mars and other solar-system bodies, recognition of impact materials will be important for elucidating planetary evolution, planetary volatile inventories, and exobiological issues. The cobble depicted is typical of many that have been found in diamictite deposits in Belize generated by the Chicxulub K-T impact event. The pebbles are roughly-hewn in general shape with smoothed corners and edges. Surfaces are almost uniformly frosted (on both protuberances and hollows), but some asperities are glazed. Optical microscopy and thin-section petrographic microscopy reveal the frosting to be only a few microns thick, with a well-defined granular structure; grains are the same size as those composing the bulk of the limestone, but their clearer appearance may represent annealing. One or two adjacent pebble faces are often decorated with striated gouges and closely-spaced hemispherical depressions representing indentation hollows produced by well-rounded impacting clasts of up to 0.5 cm in diameter. Some of the impactors are still embedded in the cobble surface. Non-destructive x-ray diffraction techniques showed the impactors to be of the same mineralogy as the target cobble. We believe this unusual glazing and frosting to be related to the impact event, but this must be reconciled with its survival for over 60 my. since it is composed of one of the most alterable substances, CaCO3. We focus, however, on the non-fractured rounded pebbles that appear to have impacted the larger cobble as a swarm of clasts, probably encountered during ballistic flight. It cannot be defined whether the cobble swept through the pebble cloud, or the pebbles rained upon a slowermoving cobble. Two interesting questions arise: (1) Where did large numbers of such well rounded pebbles come from in the ejecta curtain?, and (2) How did they embed themselves in a nominally brittle rock without suffering damage? Are the well-rounded pebbles crystalline (devitrified) melt spherules? Further investigations are in progress. If the cobble was cold and brittle, impact of well-rounded pebbles would have produced Hertzian fracture patterns (a) in the virtually elastically isotropic cobble target. For penetration depths of about 0.5 of the impactor diameter, the pebbles would require a relative impact velocity sufficient to cause pebble fragmentation and crushing, and the development of a complex Hertzian-Boussinesq fracture field (b) involving deep fracturing and lateral surface spalling (impact velocities > 50 m/s). The existing relationship could only evolve by the impact of cold, hard pebbles into a soft, plastic cobble surface (c). For limestone to have been plastic, it must have been at elevated temperature, but to prevent calcification of the material, the ambient pressure must also have been elevated. This would be possible either in the impact's gas plume, or within the confines of a thin aerodynamically-produced shock bow generated by supersonic ballistic motion of the cobble. In the latter case, it is implied that the cobble swept through the pebbles rather than vice-versa. Although there are high T-P conditions associated with volcanism, such textures have not been reported on volcaniclastic materials, nor from other high T-P environments; e.g., metamorphic (as far as the authors are aware). We propose that this very easily recognizable embedding and indentation surface texture can be used as a diagnostic criterion for the recognition of impact ejecta. The term "peening texture" is suggested, because it is absolutely analogous to the plastic-deformation induced, metal-surface textures generated by ball-bearing bombardment used in engineering metallurgy to work-harden metal surfaces. Additional information is contained in the original (Figures).

Marshall, J. R.↗

Modeling of a Non-Dilute Alloy Solidification Under Terrestrial and Microgravity Conditions

A general finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe from the melt. 2-D and 3-D models with and without Oberbeck-Boussinesq (O-B) approximation were used. Verification of the 3-D model was based on experimental data. The location and the shape of the interface were determined using melting temperatures obtained from the HgCdTe liquidus curve. Double-diffusive convection in the melt is caused by radial temperature gradients and by rejection at the solid-melt interface of the heavier HgTe-rich solute. Because of these complicating factors, dimensional rather than non-dimensional modeling was performed. Models based on complete set of Navier-Stokes equations for compressible media are superior to the O-B approximation in microgravity, because they provide more accurate estimation of convection.

Bune, Andris V.↗

Measurement of the Surface Dilatational Viscosity of an Insoluble Surfactant Monolayer at the Air/Water Interface Using a Pendant Drop Apparatus

When a fluid interface with surfactants is at rest, the interfacial stress is isotropic (as given by the equilibrium interfacial tension), and is described by the equation of state which relates the surface tension to the surfactant surface concentration. When surfactants are subjected to shear and dilatational flows, flow induced interaction of the surfactants; can create interfacial stresses apart from the equilibrium surface tension. The simplest relationship between surface strain rate and surface stress is the Boussinesq-Scriven constitutive equation completely characterized by three coefficients: equilibrium interfacial tension, surface shear viscosity, and surface dilatational viscosity Equilibrium interfacial tension and surface shear viscosity measurements are very well established. On the other hand, surface dilatational viscosity measurements are difficult because a flow which change the surface area also changes the surfactant surface concentration creating changes in the equilibrium interfacial tension that must be also taken into account. Surface dilatational viscosity measurements of existing techniques differ by five orders of magnitude and use spatially damped surface waves and rapidly expanding bubbles. In this presentation we introduce a new technique for measuring the surface dilatational viscosity by contracting an aqueous pendant drop attached to a needle tip and having and insoluble surfactant monolayer at the air-water interface. The isotropic total tension on the surface consists of the equilibrium surface tension and the tension due to the dilation. Compression rates are undertaken slow enough so that bulk hydrodynamic stresses are small compared to the surface tension force. Under these conditions we show that the total tension is uniform along the surface and that the Young-Laplace equation governs the drop shape with the equilibrium surface tension replaced by the constant surface isotropic stress. We illustrate this technique using DPPC as the insoluble surfacant monolayer and measured for it a surface dilatational viscosity in the LE phase that is 20 surface poise.

Lorenzo, Jose↗

Incompressible Navier-Stokes Computations with Heat Transfer

The existing pseudocompressibility method for the system of incompressible Navier-Stokes equations is extended to heat transfer problems by including the energy equation. The solution method is based on the pseudo compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Current computations use one-equation Baldwin-Barth turbulence model which is derived from a simplified form of the standard k-epsilon model equations. Both forced and natural convection problems are examined. Numerical results from turbulent reattaching flow behind a backward-facing step will be compared against experimental measurements for the forced convection case. The validity of Boussinesq approximation to simplify the buoyancy force term will be investigated. The natural convective flow structure generated by heat transfer in a vertical rectangular cavity will be studied. The numerical results will be compared by experimental measurements by Morrison and Tran.

Kiris, Cetin↗

Experimental and Numerical Investigation of Buoyancy Driven Convection During PDAMNA Thin Film Growth

This paper presents results from numerical simulations as well as laboratory experiments of buoyancy driven convection in an ampoule under varying heating and gravitational acceleration loadings. The modeling effort in this work resolves the large scale natural convective motion that occurs in the fluid during photodeposition of polydiacetelene films which is due to energy absorbed by the growth solution from a UV source. Consequently, the growth kinetics of the film are ignored in the model discussed here, and also a much simplified ampoule geometry is considered. The objective of this work is to validate the numerical prediction on the strength and structure of buoyancy driven convection that could occur under terrestrial conditions during nonlinear optical film growth. The validation is used to enable a reliable predictive capability on the nature and strength of the convective motion under low gravity conditions. The ampoule geometry is in the form of a parallelepiped with rectangular faces. The numerical results obtained from the solution to the Boussinesq equations show that natural convection will occur regardless of the orientation of the UV source with respect to the gravity vector. The least strong convective motion occurred with the UV beam directed at the top face of the parallelepiped. The strength of the convective motion was found to be almost linearly proportional to the total power of the UV source. Also, it was found that the strength of the convective motion decreased linearly with the gravity due to acceleration. The pattern of the convective flow on the other hand, depended on the source location.

Antar, Basil N.↗

Suppression of Buoyancy in Gaseous Media at High Temperatures

Consider a rectangular box filled with a fluid having a heated bottom and a cold top surface, and insulated side-walls (Benard problem). As the temperature difference between the horizontal top and bottom surfaces increases, a critical condition, defined quantitatively by the Rayleigh number, is reached beyond which density stratification can no longer be sustained by conduction and the fluid disrupts from its stable, quiescent state into an unstable, convective mode in which lighter and heavier gas mix. This paper suggests that such a statement is not necessarily true for gaseous media under normalized temperature differences that are much larger than justifiable for the Boussinesq approximation! In fact, there may be situations where a system cannot ever be made unstable with respect to the onset on buoyant convection no matter how large the temperature (density) difference becomes at a given pressure even under normal gravity! This unexpected behavior is primarily attributed to highly temperature-sensitive kinematic viscosity which counteracts the tendency toward instability and dampens convection by making the gas more viscous at higher temperatures. This compensation of the buoyant force by the viscous force exhibits itself by the formation of a peak hot-surface temperature beyond which a system will tend to be more stable as the hot-surface temperature increases.

Gokoglu, Suleyman A.↗

Partially-Averaged Navier Stokes Model for Turbulence: Implementation and Validation

Partially-averaged Navier Stokes (PANS) is a suite of turbulence closure models of various modeled-to-resolved scale ratios ranging from Reynolds-averaged Navier Stokes (RANS) to Navier-Stokes (direct numerical simulations). The objective of PANS, like hybrid models, is to resolve large scale structures at reasonable computational expense. The modeled-to-resolved scale ratio or the level of physical resolution in PANS is quantified by two parameters: the unresolved-to-total ratios of kinetic energy (f(sub k)) and dissipation (f(sub epsilon)). The unresolved-scale stress is modeled with the Boussinesq approximation and modeled transport equations are solved for the unresolved kinetic energy and dissipation. In this paper, we first present a brief discussion of the PANS philosophy followed by a description of the implementation procedure and finally perform preliminary evaluation in benchmark problems.

Girimaji, Sharath S.↗

Antarctic Circumpolar Current Transport Variability during 2003-05 from GRACE

Gravity Recovery and Climate Experiment (GRACE) gravity data spanning January 2003-November 2005 are used as proxies for ocean bottom pressure (BP) averaged over 1 month, spherical Gaussian caps 500 km in radius, and along paths bracketing the Antarctic Circumpolar Current's various fronts. The GRACE BP signals are compared with those derived from the Estimating the Circulation and Climate of the Ocean (ECCO) ocean modeling-assimilation system, and to a non-Boussinesq version of the Regional Ocean Model System (ROMS). The discrepancy found between GRACE and the models is 1.7 cm H2O (1 cm H2O approx. 1 hPa), slightly lower than the 1.9 cmH2O estimated by the authors independently from propagation of GRACE errors. The northern signals are weak and uncorrelated among basins. The southern signals are strong, with a common seasonality. The seasonal cycle GRACE data observed in the Pacific and Indian Ocean sectors of the ACC are consistent, with annual and semiannual amplitudes of 3.6 and 0.6 cmH2O (1.1 and 0.6 cm H2O with ECCO), the average over the full southern path peaks (stronger ACC) in the southern winter, on days of year 197 and 97 for the annual and semiannual components, respectively; the Atlantic Ocean annual peak is 20 days earlier. An approximate conversion factor of 3.1 Sv (Sv equiv 10(exp 6)cu m/s) of barotropic transport variability per cm H2O of BP change is estimated. Wind stress data time series from the Quick Scatterometer (QuikSCAT), averaged monthly, zonally, and over the latitude band 40(deg)- 65(deg)S, are also constructed and subsampled at the same months as with the GRACE data. The annual and semiannual harmonics of the wind stress peak on days 198 and 82, respectively. A decreasing trend over the 3 yr is observed in the three data types.

Gravity Recovery and Climate Experiment (GRACE)↗

Antarctic Circumpolar Current Transport Variability during 2003-05 from GRACE

Gravity Recovery and Climate Experiment (GRACE) gravity data spanning January 2003 - November 2005 are used as proxies for ocean bottom pressure (BP) averaged over 1 month, spherical Gaussian caps 500 km in radius, and along paths bracketing the Antarctic Circumpolar Current's various fronts. The GRACE BP signals are compared with those derived from the Estimating the Circulation and Climate of the Ocean (ECCO) ocean modeling-assimilation system, and to a non-Boussinesq version of the Regional Ocean Model System (ROMS). The discrepancy found between GRACE and the models is 1.7 cm(sub H2O) (1 cm(sub H2O) similar to 1 hPa), slightly lower than the 1.9 cm(sub H2O) estimated by the authors independently from propagation of GRACE errors. The northern signals are weak and uncorrelated among basins. The southern signals are strong, with a common seasonality. The seasonal cycle GRACE data observed in the Pacific and Indian Ocean sectors of the ACC are consistent, with annual and semiannual amplitudes of 3.6 and 0.6 cm(sub H2O) (1.1 and 0.6 cm(sub H2O) with ECCO), the average over the full southern path peaks (stronger ACC) in the southern winter, on days of year 197 and 97 for the annual and semiannual components, respectively; the Atlantic Ocean annual peak is 20 days earlier. An approximate conversion factor of 3.1 Sv ( Sv equivalent to 10(exp 6) m(exp 3) s(exp -1)) of barotropic transport variability per cm(sub H2O) of BP change is estimated. Wind stress data time series from the Quick Scatterometer (QuikSCAT), averaged monthly, zonally, and over the latitude band 40 de - 65 deg S, are also constructed and subsampled at the same months as with the GRACE data. The annual and semiannual harmonics of the wind stress peak on days 198 and 82, respectively. A decreasing trend over the 3 yr is observed in the three data types.

oceanographic regions↗

Enhancement of USM3D Unstructured Flow Solver for High-Speed High-Temperature Shear Flows

Large temperature and pressure fluctuations have a profound effect on turbulence development in transonic and supersonic jets. For high-speed, high-temperature jet flows, standard turbulence models lack the ability to predict the observed mixing rate of a shear layer. Several proposals to address this deficiency have been advanced in the literature to modify the turbulence transport equations in a variety of ways. In the present study, some of the most proven and simple modifications to two-equation turbulence models have been selected and implemented in NASA's USM3D tetrahedral Navier-Stokes flow solver. The modifications include the addition of compressibility correction and pressure dilatation terms in the turbulence transport equations for high-speed flows, and the addition of a simple modification to the Boussinesq's closure model coefficient for high-temperature jets. The efficacy of the extended models is demonstrated by comparison with experimental data for two supersonic axisymmetric jet test cases at design pressure ratio.

Pandya, Mohagna J.↗

Lessons Learned from Numerical Simulations of the F-16XL Aircraft at Flight Conditions

Nine groups participating in the Cranked Arrow Wing Aerodynamics Project International (CAWAPI) project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, steady or unsteady, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗

What was Learned from Numerical Simulations of F-16XL (CAWAPI) at Flight Conditions

Nine groups participating in the CAWAPI project have contributed steady and unsteady viscous simulations of a full-scale, semi-span model of the F-16XL aircraft. Three different categories of flight Reynolds/Mach number combinations were computed and compared with flight-test measurements for the purpose of code validation and improved understanding of the flight physics. Steady-state simulations are done with several turbulence models of different complexity with no topology information required and which overcome Boussinesq-assumption problems in vortical flows. Detached-eddy simulation (DES) and its successor delayed detached-eddy simulation (DDES) have been used to compute the time accurate flow development. Common structured and unstructured grids as well as individually-adapted unstructured grids were used. Although discrepancies are observed in the comparisons, overall reasonable agreement is demonstrated for surface pressure distribution, local skin friction and boundary velocity profiles at subsonic speeds. The physical modeling, be it steady or unsteady flow, and the grid resolution both contribute to the discrepancies observed in the comparisons with flight data, but at this time it cannot be determined how much each part contributes to the whole. Overall it can be said that the technology readiness of CFD-simulation technology for the study of vehicle performance has matured since 2001 such that it can be used today with a reasonable level of confidence for complex configurations.

Rizzi, Arthur↗