Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 343 records · Page 19

A Dynamic Nonlinear Subgrid-Scale Model for Large-Eddy Simulation of Complex Turbulent Flows

We present a new dynamic nonlinear subgrid-scale (SGS) model for large-eddy simulations (LES) and apply it to compute a flow involving pressure gradients, surface curvature and separation, for which data from a direct numerical simulation are available for comparison. The model, inspired by the triple model idea of Bardina et al. (“Improved Turbulence Models Based on Large Eddy Simulation of Homogeneous, Incompressible, Turbulent Flows,” Report No. TF-19, Thermosciences Division, Department of Mechanical Engineering, Stanford University, 1983), includes a Galilean-invariant term called the modified Leonard stress tensor, and two nonlinear terms comprised of the products of the strain-rate and rotation-rate tensors for an improved representation of the subgrid-scale dissipation, backscatter and anisotropy effects. The model does not employ any ad hoc averaging or clipping procedures, and does not require the specification of a characteristic length scale; hence, it naturally avoids the ambiguities associated with defining a proper length scale for anisotropic grids. Results from the wall-resolved LES of flow past a Gaussian bump using the new model demonstrate improved prediction of skin-friction, flow separation, mean flow profiles and turbulent quantities when compared to implicit LES as well as explicit LES using the Vreman SGS model on the same grid.

Turbulence↗

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Ground-Based Passive Microwave Remote Sensing Observations of Soil Moisture at S and L Band with Insight into Measurement Accuracy

Accurate estimates of spatially heterogeneous algorithm variables and parameters are required in determining the spatial distribution of soil moisture using radiometer data from aircraft and satellites. A ground-based experiment in passive microwave remote sensing of soil moisture was conducted in Huntsville, Alabama from July 1-14, 1996 to study retrieval algorithms and their sensitivity to variable and parameter specification. With high temporal frequency observations at S and L band, we were able to observe large scale moisture changes following irrigation and rainfall events, as well as diurnal behavior of surface moisture among three plots, one bare, one covered with short grass and another covered with alfalfa. The L band emitting depth was determined to be on the order of 0-3 or 0-5 cm below 0.30 cubic centimeter/cubic centimeter with an indication of a shallower emitting depth at higher moisture values. Surface moisture behavior was less apparent on the vegetated plots than it was on the bare plot because there was less moisture gradient and because of difficulty in determining vegetation water content and estimating the vegetation b parameter. Discrepancies between remotely sensed and gravimetric, soil moisture estimates on the vegetated plots point to an incomplete understanding of the requirements needed to correct for the effects of vegetation attenuation. Quantifying the uncertainty in moisture estimates is vital if applications are to utilize remotely-sensed soil moisture data. Computations based only on the real part of the complex dielectric constant and/or an alternative dielectric mixing model contribute a relatively insignificant amount of uncertainty to estimates of soil moisture. Rather, the retrieval algorithm is much more sensitive to soil properties, surface roughness and biomass.

Laymon, Charles A.↗

Measurement Via Optical Near-Nulling and Subaperture Stitching

A subaperture stitching interferometer system provides near-nulling of a subaperture wavefront reflected from an object of interest over a portion of a surface of the object. A variable optical element located in the radiation path adjustably provides near-nulling to facilitate stitching of subaperture interferograms, creating an interferogram representative of the entire surface of interest. This enables testing of aspheric surfaces without null optics customized for each surface prescription. The surface shapes of objects such as lenses and other precision components are often measured with interferometry. However, interferometers have a limited capture range, and thus the test wavefront cannot be too different from the reference or the interference cannot be analyzed. Furthermore, the performance of the interferometer is usually best when the test and reference wavefronts are nearly identical (referred to as a null condition). Thus, it is necessary when performing such measurements to correct for known variations in shape to ensure that unintended variations are within the capture range of the interferometer and accurately measured. This invention is a system for nearnulling within a subaperture stitching interferometer, although in principle, the concept can be employed by wavefront measuring gauges other than interferometers. The system employs a light source for providing coherent radiation of a subaperture extent. An object of interest is placed to modify the radiation (e.g., to reflect or pass the radiation), and a variable optical element is located to interact with, and nearly null, the affected radiation. A detector or imaging device is situated to obtain interference patterns in the modified radiation. Multiple subaperture interferograms are taken and are stitched, or joined, to provide an interferogram representative of the entire surface of the object of interest. The primary aspect of the invention is the use of adjustable corrective optics in the context of subaperture stitching near-nulling interferometry, wherein a complex surface is analyzed via multiple, separate, overlapping interferograms. For complex surfaces, the problem of managing the identification and placement of corrective optics becomes even more pronounced, to the extent that in most cases the null corrector optics are specific to the particular asphere prescription and no others (i.e. another asphere requires completely different null correction optics). In principle, the near-nulling technique does not require subaperture stitching at all. Building a near-null system that is practically useful relies on two key features: simplicity and universality. If the system is too complex, it will be difficult to calibrate and model its manufacturing errors, rendering it useless as a precision metrology tool and/or prohibitively expensive. If the system is not applicable to a wide range of test parts, then it does not provide significant value over conventional null-correction technology. Subaperture stitching enables simpler and more universal near-null systems to be effective, because a fraction of a surface is necessarily less complex than the whole surface (excepting the extreme case of a fractal surface description). The technique of near-nulling can significantly enhance aspheric subaperture stitching capability by allowing the interferometer to capture a wider range of aspheres. More over, subaperture stitching is essential to a truly effective near-nulling system, since looking at a fraction of the surface keeps the wavefront complexity within the capability of a relatively simple nearnull apparatus. Furthermore, by reducing the subaperture size, the complexity of the measured wavefront can be reduced until it is within the capability of the near-null design.

Forbes, Greg↗

The interaction of moderately strong shock waves with thick perforated walls of low porosity

A theoretical prediction is given of the flow through thick perforated walls of low porosity resulting from the impingement of a moderately strong traveling shock wave. The model was a flat plate positioned normal to the direction of the flow. Holes bored in the plate parallel to the direction of the flow provided nominal hole length-to-diameter ratios of 10:1 and an axial porosity of 25 percent of the flow channel cross section. The flow field behind the reflected shock wave was assumed to behave as a reservoir producing a quasi-steady duct flow through the model. Rayleigh and Fanno duct flow theoretical computations for each of three possible auxiliary wave patterns that can be associated with the transmitted shock (to satisfy contact surface compatibility) were used to provide bounding solutions as an alternative to the more complex influence coefficients method. Qualitative and quantitative behavior was verified in a 1.5- by 2.0-in. helium shock tube. High speed Schlieren photography, piezoelectric pressure-time histories, and electronic-counter wave speed measurements were used to assess the extent of correlation with the theoretical flow models. Reduced data indicated the adequacy of the bounding theory approach to predict wave phenomena and quantitative response.

Grant, D. J.↗

Meteorological and air pollution modeling for an urban airport

Results are presented of numerical experiments modeling meteorology, multiple pollutant sources, and nonlinear photochemical reactions for the case of an airport in a large urban area with complex terrain. A planetary boundary-layer model which predicts the mixing depth and generates wind, moisture, and temperature fields was used; it utilizes only surface and synoptic boundary conditions as input data. A version of the Hecht-Seinfeld-Dodge chemical kinetics model is integrated with a new, rapid numerical technique; both the San Francisco Bay Area Air Quality Management District source inventory and the San Jose Airport aircraft inventory are utilized. The air quality model results are presented in contour plots; the combined results illustrate that the highly nonlinear interactions which are present require that the chemistry and meteorology be considered simultaneously to make a valid assessment of the effects of individual sources on regional air quality.

Swan, P. R.↗

Thermoplastic Composites for Research-Model Components

Oriented unidirectional prepreg tapes formed in ceramic molds. New technique developed at Langley Research Center, using ceramic mold for fabrication of models from graphite/thermoplastic materials. Upper surface of wings and fuselage of advanced airplane selected as shape to test fabrication technique. Technique well suited for other complex shapes as well. Thermoplastic composite easily workable with normal shop equipment and painted and repaired by standard methods.

Guenther, B. F.↗

Influence of airfoil geometry on delta wing leading-edge vortices and vortex-induced aerodynamics at supersonic speeds

An assessment of the influence of airfoil geometry on delta wing leading edge vortex flow and vortex induced aerodynamics at supersonic speeds is discussed. A series of delta wing wind tunnel models were tested over a Mach number range from 1.7 to 2.0. The model geometric variables included leading edge sweep and airfoil shape. Surface pressure data, vapor screen, and oil flow photograph data were taken to evaluate the complex structure of the vortices and shocks on the family of wings tested. The data show that airfoil shape has a significant impact on the wing upper surface flow structure and pressure distribution, but has a minimal impact on the integrated upper surface pressure increments.

Wood, Richard M.↗

Implementation and operation of three fractal measurement algorithms for analysis of remote-sensing data

Fractal geometry is increasingly becoming a useful tool for modeling natural phenomena. As an alternative to Euclidean concepts, fractals allow for a more accurate representation of the nature of complexity in natural boundaries and surfaces. The purpose of this paper is to introduce and implement three algorithms in C code for deriving fractal measurement from remotely sensed data. These three methods are: the line-divider method, the variogram method, and the triangular prism method. Remote-sensing data acquired by NASA's Calibrated Airborne Multispectral Scanner (CAMS) are used to compute the fractal dimension using each of the three methods. These data were obtained as a 30 m pixel spatial resolution over a portion of western Puerto Rico in January 1990. A description of the three methods, their implementation in PC-compatible environment, and some results of applying these algorithms to remotely sensed image data are presented.

Jaggi, S.↗

Multi-resolution processing for fractal analysis of airborne remotely sensed data

Fractal geometry is increasingly becoming a useful tool for modeling natural phenomenon. As an alternative to Euclidean concepts, fractals allow for a more accurate representation of the nature of complexity in natural boundaries and surfaces. Since they are characterized by self-similarity, an ideal fractal surface is scale-independent; i.e. at different scales a fractal surface looks the same. This is not exactly true for natural surfaces. When viewed at different spatial resolutions parts of natural surfaces look alike in a statistical manner and only for a limited range of scales. Images acquired by NASA's Thermal Infrared Multispectral Scanner are used to compute the fractal dimension as a function of spatial resolution. Three methods are used to determine the fractal dimension - Schelberg's line-divider method, the variogram method, and the triangular prism method. A description of these methods and the results of applying these methods to a remotely-sensed image is also presented. Five flights were flown in succession at altitudes of 2 km (low), 6 km (mid), 12 km (high), and then back again at 6 km and 2 km. The area selected was the Ross Barnett reservoir near Jackson, Mississippi. The mission was flown during the predawn hours of 1 Feb. 1992. Radiosonde data was collected for that duration to profile the characteristics of the atmosphere. This corresponds to 3 different pixel sizes - 5m, 15m, and 30m. After, simulating different spatial sampling intervals within the same image for each of the 3 image sets, the results are cross-correlated to compare the extent of detail and complexity that is obtained when data is taken at lower spatial intervals.

Jaggi, S.↗

Monthly atmospheric and oceanographic surface fields for the western North Atlantic: December 1986 - April 1989

Monthly atmospheric and oceanographic variables for the western North Atlantic Ocean from various sources are presented as contour or vector maps. These fields were assembled for a study of the upper ocean heat budget. Atmospheric fields include the net surface heat fluxes and wind stress derived from the 1000 mb winds from the European Center for Medium-range Weather Forecasting (ECMWF). Oceanographic fields include the sea surface height from the Geosat radar altimeter and sea surface temperature from the Advanced Very High Resolution Radiometer (AVHRR). An additional estimate of net surface heat flux is shown; this estimate was derived by assimilating winds, currents and ocean temperatures into a mixed layer model. The maps show a complex interplay of fluctuations in the winds and heat fluxes, and in the structure and temperature gradients of the Gulf Stream system. Some comments are offered on a comparison of the two heat flux estimates.

Caruso, Michael J.↗

Observations and Modeling of the Transient General Circulation of the North Pacific Basin

Because of recent progress in satellite altimetry and numerical modeling and the accumulation and archiving of long records of hydrographic and meteorological variables, it is becoming feasible to describe and understand the transient general circulation of the ocean (i.e., variations with spatial scales larger than a few hundred kilometers and time scales of seasonal and longer-beyond the mesoscale). We have carried out various studies in investigation of the transient general circulation of the Pacific Ocean from a coordinated analysis of satellite altimeter data, historical hydrographic gauge data, scatterometer wind observations, reanalyzed operational wind fields, and a variety of ocean circulation models. Broadly stated, our goal was to achieve a phenomenological catalogue of different possible types of large-scale, low-frequency variability, as a context for understanding the observational record. The approach is to identify the simplest possible model from which particular observed phenomena can be isolated and understood dynamically and then to determine how well these dynamical processes are represented in more complex Oceanic General Circulation Models (OGCMs). Research results have been obtained on Rossby wave propagation and transformation, oceanic intrinsic low-frequency variability, effects of surface gravity waves, pacific data analyses, OGCM formulation and developments, and OGCM simulations of forced variability.

McWilliams, James C.↗

Solid Lubrication Fundamentals and Applications: Abrasion: Plowing and Cutting - Chapter 5

Chapter 5 discusses abrasion, a common wear phenomenon of great economic importance. It has been estimated that 50% of the wear encountered in industry is due to abrasion. Also, it is the mechanism involved in the finishing of many surfaces. Experiments are described to help in understanding the complex abrasion process and in predicting friction and wear behavior in plowing and/or cutting. These experimental modelings and measurements used a single spherical pin (asperity) and a single wedge pin (asperity). Other two-body and three-body abrasion studies used hard abrasive particles.

Miyoshi, Kazuhisa↗

Simulation Tools Model Icing for Aircraft Design

Here s a simple science experiment to try: Place an unopened bottle of distilled water in your freezer. After 2-3 hours, if the water is pure enough, you will notice that it has not frozen. Carefully pour the water into a bowl with a piece of ice in it. When it strikes the ice, the water will instantly freeze. One of the most basic and commonly known scientific facts is that water freezes at around 32 F. But this is not always the case. Water lacking any impurities for ice crystals to form around can be supercooled to even lower temperatures without freezing. High in the atmosphere, water droplets can achieve this delicate, supercooled state. When a plane flies through clouds containing these droplets, the water can strike the airframe and, like the supercooled water hitting the ice in the experiment above, freeze instantly. The ice buildup alters the aerodynamics of the plane - reducing lift and increasing drag - affecting its performance and presenting a safety issue if the plane can no longer fly effectively. In certain circumstances, ice can form inside aircraft engines, another potential hazard. NASA has long studied ways of detecting and countering atmospheric icing conditions as part of the Agency s efforts to enhance aviation safety. To do this, the Icing Branch at Glenn Research Center utilizes a number of world-class tools, including the Center s Icing Research Tunnel and the NASA 607 icing research aircraft, a "flying laboratory" for studying icing conditions. The branch has also developed a suite of software programs to help aircraft and icing protection system designers understand the behavior of ice accumulation on various surfaces and in various conditions. One of these innovations is the LEWICE ice accretion simulation software. Initially developed in the 1980s (when Glenn was known as Lewis Research Center), LEWICE has become one of the most widely used tools in icing research and aircraft design and certification. LEWICE has been transformed over the years from strictly a research tool to one used routinely by industry and other government agencies. Glenn contractor William Wright has been the architect of this development, supported by a team of researchers investigating icing physics, creating validation data, and ensuring development according to standard software engineering practices. The program provides a virtual simulation environment for determining where water droplets strike an airfoil in flight, what kind of ice would result, and what shape that ice would take. Users can enter geometries for specific, two-dimensional cross sections of an airfoil or other airframe surface and then apply a range of inputs - different droplet sizes, temperatures, airspeeds, and more - to model how ice would build up on the surface in various conditions. The program s versatility, ease of use, and speed - LEWICE can run through complex icing simulations in only a few minutes - have contributed to it becoming a popular resource in the aviation industry.

Source record↗

Bias Correction of Hydrologic Projections Strongly Impacts Inferred Climate Vulnerabilities in Institutionally Complex Water Systems

Water-resources planners use regional water management models (WMMs) to identify vulnerabilities to climate change. Frequently, dynamically downscaled climate inputs are used in conjunction with land-surface models (LSMs) to provide hydrologic streamflow projections, which serve as critical inputs for WMMs. Here, we show how even modest projection errors can strongly affect assessments of water availability and financial stability for irrigation districts in California. Specifically, our results highlight that LSM errors in projections of flood and drought extremes are highly interactive across timescales, path-dependent, and can be amplified when modeling infrastructure systems (e.g., misrepresenting banked groundwater). Common strategies for reducing errors in deterministic LSM hydrologic projections (e.g., bias correction) can themselves strongly distort projected climate vulnerabilities and misrepresent their inferred financial consequences. Overall, our results indicate a need to move beyond standard deterministic climate projection and error management frameworks that are dependent on single simulated climate change scenario outcomes.

Keyvan Malek↗

Relationship between radiation-induced aberrations in individual chromosomes and their DNA content: effects of interaction distance

PURPOSE: To study the effect of the interaction distance on the frequency of inter- and intrachromosome exchanges in individual chromosomes with respect to their DNA content. Assumptions: Chromosome exchanges are formed by misrejoining of two DNA double-strand breaks (DSB) induced within an interaction distance, d. It is assumed that chromosomes in G(0)/G(1) phase of the cell cycle occupy a spherical domain in a cell nucleus, with no spatial overlap between individual chromosome domains. RESULTS: Formulae are derived for the probability of formation of inter-, as well as intra-, chromosome exchanges relating to the DNA content of the chromosome for a given interaction distance. For interaction distances <1 microm, the relative frequency of interchromosome exchanges predicted by the present model is similar to that by Cigarran et al. (1998) based on the assumption that the probability of interchromosome exchanges is proportional to the "surface area" of the chromosome territory. The "surface area" assumption is shown to be a limiting case of d-->0 in the present model. The present model also predicts that the probability of intrachromosome exchanges occurring in individual chromosomes is proportional to their DNA content with correction terms. CONCLUSION: When the interaction distance is small, the "surface area" distribution for chromosome participation in interchromosome exchanges has been expected. However, the present model shows that for the interaction distance as large as 1 microm, the predicted probability of interchromosome exchange formation is still close to the surface area distribution. Therefore, this distribution does not necessarily rule out the formation of complex chromosomal aberrations by long-range misrejoining of DSB.

Non-programmatic↗

Effects of roughness on heating and boundary-layer transition. Part 1: Effects of simulated panel joints on boundary-layer transition

The effects of distributed roughness, arising from space shuttle orbiter panel joints, on boundary layer transition are investigated. The North American Rockwell 134B delta wing shuttle configuration was used for the study. Results show: (1) Premature boundary layer transition occurred in models having simulated heat shield panels with rased joints. (2) Laminar flow was maintained with slot type panel joint models to a Reynolds number twice that at which transition occurred in raised joints. (3) Significant increases in peak surface temperature and the time during which turbulent flow occurs may result from distributed roughness of heat shield panel joints. (4) Laminar and turbulent heating levels were predicted within available theories. (5) A complex interference between the wing and fuselage flow was observed in the delta wing model.

Seegmiller, H. L.↗