Finite element solution techniques for large-scale problems in computational fluid dynamics

Element-by-element approximate factorization, implicit-explicit and adaptive implicit-explicit approximation procedures are presented for the finite-element formulations of large-scale fluid dynamics problems. The element-by-element approximation scheme totally eliminates the need for formation, storage and inversion of large global matrices. Implicit-explicit schemes, which are approximations to implicit schemes, substantially reduce the computational burden associated with large global matrices. In the adaptive implicit-explicit scheme, the implicit elements are selected dynamically based on element level stability and accuracy considerations. This scheme provides implicit refinement where it is needed. The methods are applied to various problems governed by the convection-diffusion and incompressible Navier-Stokes equations. In all cases studied, the results obtained are indistinguishable from those obtained by the implicit formulations

A computational analysis of under-expanded jets in the hypersonic regime

Underexpanded axisymmetric jets are studied numerically using a full Navier-Stokes solver. Emphasis has been given to supersonic and hypersonic jets in supersonic and hypersonic ambient flows, a phenomenon previously overlooked. It is demonstrated that the shear layers and shock patterns in a jet plume can be captured without complicated viscous/inviscid and subsonic/supersonic coupling schemes. In addition, a supersonic pressure relief effect has been identified for underexpanded jets in supersonic ambient flows. While it is well known that an underexpanded jet in a quiescent ambience (or subsonic ambience) contains multiple shock cells, the present study shows that because of the supersonic pressure relief effect, an underexpanded jet in a supersonic or hypersonic ambience contains only one major shock cell

On the basic equations for the second-order modeling of compressible turbulence

Equations for the mean and turbulent quantities for compressible turbulent flows are derived. Both the conventional Reynolds average and the mass-weighted, Favre average were employed to decompose the flow variable into a mean and a turbulent quality. These equations are to be used later in developing second order Reynolds stress models for high speed compressible flows. A few recent advances in modeling some of the terms in the equations due to compressibility effects are also summarized

On the Basic Equations for the Second-order Modeling of Compressible Turbulence

Equations for the mean and the turbulence quantities of compressible turbulent flows are derived in this report. Both the conventional Reynolds average and the mass-weighted Favre average were employed to decompose the flow variable into mean and turbulent quantities. These equations are to be used later in developing second-order Reynolds stress models for high-speed compressible flows. A few recent advances in modeling some of the terms in the equation due to compressibility effects are also summarized

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. Part 2: Wall shear stress

An analysis is presented of the flow in the two inner layers, the Reynolds stress sublayer and the wall layer. Included is the calculation of the shear stress at the wall in the interaction region. The limit processes considered are those used for an inviscid flow

The Effect of Clouds on Air Showers Observation from Space

Issues relating to extensive air showers observation by a space-borne fluorescence detector and the effects of clouds on the observations are investigated using Monte Carlo simulation. The simulations assume the presence of clouds with varying altitudes and optical depths. Simulated events are reconstructed assuming a cloud-free atmosphere. While it is anticipated that auxiliary instruments, such as LIDAR (LIght Detection And Ranging), will be employed to measure the atmospheric conditions during actual observation, it is still possible that these instruments may fail to recognize the presence of a cloud in a particular shower observation. The purpose of this study is to investigate the effects on the reconstructed shower parameters in such cases. Reconstruction results are shown for both monocular and stereo detectors and for the two limiting cases of optically thin, and optically thick clouds.Comment: 31 pages, 13 figure

Interaction between a normal shock wave and a turbulent boundary layer at high transonic speeds. Part 1: Pressure distribution. Part 2: Wall shear stress. Part 3: Simplified formulas for the prediction of surface pressures and skin friction

An asymptotic description is derived for the interaction between a shock wave and a turbulent boundary layer in transonic flow, for a particular limiting case. The dimensionless difference between the external flow velocity and critical sound speed is taken to be much smaller than one, but large in comparison with the dimensionless friction velocity. The basic results are derived for a flat plate, and corrections for longitudinal wall curvature and for flow in a circular pipe are also shown. Solutions are given for the wall pressure distribution and the shape of the shock wave. Solutions for the wall shear stress are obtained, and a criterion for incipient separation is derived. Simplified solutions for both the wall pressure and skin friction distributions in the interaction region are given. These results are presented in a form suitable for use in computer programs

Remote sensing of tropical tropopause layer radiation balance using A-train measurements

Determining the level of zero net radiative heating (LZH) is critical to understanding parcel trajectory in the Tropical Tropopause Layer (TTL) and associated stratospheric hydration processes. Previous studies of the TTL radiative balance have focused on using radiosonde data, but remote sensing measurements from polar-orbiting satellites may provide the relevant horizontal and vertical information for assessing TTL solar heating and infrared cooling rates, especially across the Pacific Ocean. CloudSat provides a considerable amount of vertical information about the distribution of cloud properties relevant to heating rate analysis. The ability of CloudSat measurements and ancillary information to constrain LZH is explored. We employ formal error propagation analysis for derived heating rate uncertainty given the CloudSat cloud property retrieval algorithms. Estimation of the LZH to within approximately 0.5 to 1 km is achievable with CloudSat, but it has a low-altitude bias because the radar is unable to detect thin cirrus. This can be remedied with the proper utilization of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) lidar backscatter information. By utilizing an orbital simulation with the GISS data set, we explore the representativeness of non-cross-track scanning active sounders in terms of describing the LZH distribution. In order to supplement CloudSat, we explore the ability of Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Scanning Radiometer-EOS (AMSR-E) to constrain LZH and find that these passive sounders are useful where the cloud top height does not exceed 7 km. The spatiotemporal distributions of LZH derived from CloudSat and CALIPSO measurements are presented which suggest that thin cirrus have a limited effect on LZH mean values but affect LZH variability

Computer Simulation Of Three-dimensional Mechanical Assemblies: Part II - Computer Simulation

Computer simulation of the kinematic and dynamic behaviors of mechanical assemblies has become a very important tool in design and manufacturing, because the designer can foresee how a product is going to perform before the product is actually fabricated. However, up to now, the most current simulation modules are based on analysis from another kinematic or dynamic module by specifying the mating conditions between components, and then displaying the motion on the screen. This computer simulation actually performs similarly to a movie and can only provide visual checking. The drawback of this simulation approach is that designers are forced to use the available joint models and may lose their creativity. In part I of this paper, general mathematical modeling of the multi-body system is presented, while part n of this paper, a prototype convex-feature modeling system is presented with which a designer can interactively create an assembly of mechanical components ready for dynamic analysis. It can provide a state-of-the-art technology for real simulation of any mechanical systems, and act as a cost-effective test bed for concepts, final design, and control algorithms

Previously Unidentified Changes in Renal Cell Carcinoma Gene Expression Identified by Parametric Analysis of Microarray Data

BACKGROUND. Renal cell carcinoma is a common malignancy that often presents as a metastatic-disease for which there are no effective treatments. To gain insights into the mechanism of renal cell carcinogenesis, a number of genome-wide expression profiling studies have been performed. Surprisingly, there is very poor agreement among these studies as to which genes are differentially regulated. To better understand this lack of agreement we profiled renal cell tumor gene expression using genome-wide microarrays (45,000 probe sets) and compare our analysis to previous microarray studies. METHODS. We hybridized total RNA isolated from renal cell tumors and adjacent normal tissue to Affymetrix U133A and U133B arrays. We removed samples with technical defects and removed probesets that failed to exhibit sequence-specific hybridization in any of the samples. We detected differential gene expression in the resulting dataset with parametric methods and identified keywords that are overrepresented in the differentially expressed genes with the Fisher-exact test. RESULTS. We identify 1,234 genes that are more than three-fold changed in renal tumors by t-test, 800 of which have not been previously reported to be altered in renal cell tumors. Of the only 37 genes that have been identified as being differentially expressed in three or more of five previous microarray studies of renal tumor gene expression, our analysis finds 33 of these genes (89%). A key to the sensitivity and power of our analysis is filtering out defective samples and genes that are not reliably detected. CONCLUSIONS. The widespread use of sample-wise voting schemes for detecting differential expression that do not control for false positives likely account for the poor overlap among previous studies. Among the many genes we identified using parametric methods that were not previously reported as being differentially expressed in renal cell tumors are several oncogenes and tumor suppressor genes that likely play important roles in renal cell carcinogenesis. This highlights the need for rigorous statistical approaches in microarray studies.National Institutes of Healt