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 289 records · Page 16

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗

Zemax Modeling of a Self-aligned Focusing Schlieren System with Simulated Compressible Flow

A self-aligned focusing schlieren system is simulated in the ANSYS Zemax OpticStudio software environment. The model uses a custom scatter-defining dynamic link library on the retroreflective background in order to scatter rays in the direction of the incident ray rather than the specular ray. Simulated compressible flow refractive index fields are imported into Zemax by use of the grid gradient object and a workaround to the maximum point limit is introduced. The Ronchi ruling is built as a high resolution bitmap image in MATLAB and set in Zemax as a slide object allowing the bitmap image to either transmit or terminate rays at each line pair. The Zemax simulations are compared to computational numerical schlieren and experimental conventional, path-integrated schlieren results from a hypersonic stagnation point injection flow configuration and an experimental image of a shock wave from a laser-induced breakdown. In the stagnation point injection simulations, the Zemax schlieren images capture complex flow features such as the Mach disk, shear layers, contact surface, and bow shock well, but contains long, curved features (striations) between the model and the bow shock that are not present in the synthetic schlieren image from the CFD or the experimental image. The image from Zemax of the spark-induced shock wave captures the shape and curvature of the shock but appears blurry due to grid discretization in the computational data.

Nicholas A Mejia↗

Benchmarking a Soil Moisture Data Assimilation System for Agricultural Drought Monitoring

Despite considerable interest in the application of land surface data assimilation systems (LDAS) for agricultural drought applications, relatively little is known about the large-scale performance of such systems and, thus, the optimal methodological approach for implementing them. To address this need, this paper evaluates an LDAS for agricultural drought monitoring by benchmarking individual components of the system (i.e., a satellite soil moisture retrieval algorithm, a soil water balance model and a sequential data assimilation filter) against a series of linear models which perform the same function (i.e., have the same basic inputoutput structure) as the full system component. Benchmarking is based on the calculation of the lagged rank cross-correlation between the normalized difference vegetation index (NDVI) and soil moisture estimates acquired for various components of the system. Lagged soil moistureNDVI correlations obtained using individual LDAS components versus their linear analogs reveal the degree to which non-linearities andor complexities contained within each component actually contribute to the performance of the LDAS system as a whole. Here, a particular system based on surface soil moisture retrievals from the Land Parameter Retrieval Model (LPRM), a two-layer Palmer soil water balance model and an Ensemble Kalman filter (EnKF) is benchmarked. Results suggest significant room for improvement in each component of the system.

soil moisture↗

Lunar rock surfaces as detectors of solar processes

Lunar rock surfaces exposed at or just below the lunar surface are considered as detectors of the solar wind, solar flares and solar-derived magnetic fields through their interactions with galactic cosmic rays. The degradation of the solar detector capabilities of lunar surface rocks by meteoroid impact erosion, accreta deposition, loose dust, and sputtering, amorphous layer formation and accelerated diffusion due to solar particles and illumination is discussed, and it is noted that the complex interactions of factors affecting the outer micron of exposed surface material has so far prevented the development of a satisfactory model for a particle detector on the submicron scale. Methods for the determination of surface exposure ages based on the accumulation of light solar wind noble gases, Fe and Mg, impact craters, solar flare tracks, and cosmogenic Kr isotopes are examined, and the systematic variations in the ages determined by the various clocks are discussed. It is concluded that a means of obtaining satisfactory quantitative rate or flux data has not yet been established.

Hartung, J. B.↗

Effects of vegetation on passive microwave estimates of soil moisture

Data collected over plots using a microwave radiometer with a wavelength of 21 cm were used to evaluate a theoretical model of the effect of vegetative cover on the emission of soils. The theoretical model was used as a basis of a simplified model which predicts the emissivity of a soil-vegetation complex with a full canopy cover. This model, when transformed to predict soil moisture, requires measurements of the microwave brightness temperature, surface temperature and vegetation biomass. Since surface temperature can be determined using thermal infrared measurements and vegetation biomass could be estimated with visible and near infrared measurements, the simplified model could be applied using only remotely sensed data.

Jackson, T. J.↗

Rolling DICE to Advance Knowledge of Land-Atmosphere Interactions

The Diurnal Land-Atmosphere Coupling Experiment (DICE) aims to explore the complex interactions between the land surface and atmospheric boundary layer, which are generally not well understood and difficult to isolate in models. The project involves over 10 different models, combining expertise from both land surface and atmospheric boundary-layer modelling groups. A simple three-stage methodology is designed to assess land-atmosphere feedbacks. Stage 1: the individual components are assessed in isolation, driven and evaluated against observational data; stage 2: the impact of coupling is investigated; stage 3: the sensitivity of the stand-alone models to variations in driving data is explored. For this initial study, a three-day clear-sky period in the mid-west USA over, an assumed simple, predominantly grass surface was simulated using data from the CASES-99 field campaign. Key conclusions from the study include: (1) the memory of vegetation state within Land Surface Models (LSMs) needs attention; (2) the height of atmospheric forcing for LSMs is important, particularly for the nocturnal boundary layer, and this has implications for both observations and vertical resolution for atmospheric models; (3) land-atmosphere feedbacks reduce errors in simulated surface fluxes at the expense of the accuracy of the variables that the models are designed to simulate (e.g., temperature, humidity and wind speed); (4) problems remain in representing the stable boundary layer in atmospheric models; (5) the mixing of temperature and humidity within the boundary layer may need to be represented separately; (6) differences in daytime profiles of heat, moisture and momentum between models are mainly due to the way the models erode the inversion at the top of the boundary layer, rather than differences in the surface fluxes. Resultant variations in modelled boundary layer heights have a substantial impact on relative humidity and could partially explain variations in coupling strength between models in the GLACE experiment.

Martin Best↗

Advances in Land Data Assimilation Systems

Assimilation of remotely sensed land surface observations into regional to global scale numerical models have the potential to significantly advance our ability, to assess, understand, and predict surface water, energy, and carbon cycles. This session seeks to assess the state-of-the-art in data assimilation methods for integrating land surface remote sensing and modeling, with a focus on practical applications and techniques. Assimilated land surface variables of interest include (but are not limited to, soil moisture, surface temperature, snowpack, streamflow, vegetation dynamics, and carbon storage. Contributions describing the development of practical land surface data assimilation methods, multivariate land surface data assimilation strategies, evaluation of the required accuracy and resolution of remote sensing observations, the effects of scale, process complexity, and uncertainty on data assimilation, and the optimal treatment of model and observation errors are encouraged.

Houser, Paul R.↗

A tentative unified sea model for scattering and emission

This work is an attempt to relate active and passive observations of the sea surface to each other, considering that more information can be obtained from a set of instruments when they are not just considered as a juxtaposition of separate and independent boxes, but as a single complex instrument. A careful review of several approaches to the solution of the scattering problem for a random rough surface has been performed. Inherent simplifications, and inconsistencies, are pointed out. To reduce the computation time to reasonable values, an approximate model is selected and some numerical results will be presented in order to illustrate the benefit that can be gained by a unified approach.

Guissard, A.↗

Theoretical evaluation of a V/STOL fighter model utilizing the PAN AIR code

The PAN AIR computer code was investigated as a tool for predicting closely coupled aerodynamic and propulsive flowfields of arbitrary configurations. The NASA/Ames V/STOL fighter model, a configuration of complex geometry, was analyzed with the PAN AIR code. A successful solution for this configuration was obtained when the nozzle exit was treated as an impermeable surface and no wakes were included around the nozzle exit. When separated flow was simulated from the end of the nacelle, requiring the use of wake networks emanating from the nozzle exit, a number of problems were encountered. A circular body nacelle model was used to investigate various techniques for simulating the exhaust plume in PAN AIR. Several approaches were tested and eliminated because they could not correctly simulate the interference effects. Only one plume modeling technique gave good results. A PAN AIR computation that used a plume shape and inflow velocities obtained from the Navier-Stokes solution for the plume produced results for the effects of power that compared well with experimental data.

Howell, G. A.↗

Radiation exchange in large space structures and frames

The analysis of radiation heat transfer is determined by the need to determine the geometrical quantity, F(ij), which is commonly called the view factor, and which represents the fraction of energy leaving surface i which strikes surface j. The present investigation shows that the use of a rotating plane to model a one-dimensional radiation element is both computationally efficient and accurate. Because of the reduced computational cost of this element, radiation view factors for complex structures can be calculated directly. Under some conditions, the characteristics of the one-dimensional element may lead to a loss in accuracy when it is used to model obstructing cylinders. In cases in which the element is too close to the radiating surfaces, curved surfaces must be represented by a collection of flat planes. However, it appears that as long as the one-dimensional element is more than 5 diameters from the other surfaces, its use is valid.

Emery, A. F.↗

Retrieval of Surface Snow Grain Size and Melt Water from AVIRIS Spectra

The Earth's energy balance and hydrology are affected by the distribution and characteristics of snow cover on the surface. Snow grain size and snow melt influence surface albedo and hydrology. A model of snow reflectance that depends on both grain size and surface melt water was developed to derive these parameters from remote spectral measurements. This reflectance model is based on a discrete ordinate radiative transfer approach that uses Mie calculations of snow optical properties, which are based on the complex refractive index of ice and water. This snow model was linked to an atmospheric radiative transfer code and a nonlinear least squares fitting algorithm. The resulting combined algorithm was applied to an AVIRIS snow data set acquired over Mammoth Mountain, California. Maps of grain size and surface snow melt were generated that are consistent with the expected ranges and distributions for conditions at the site.

Green Robert O.↗

Simulation studies for surfaces and materials strength

Computer simulation studies were conducted to investigate materials properties. All the calculations were carried out using atomic level simulation techniques which are based on semiempirical or model functions. During this project, functions with varying degrees of complexity were derived and employed in simulations. Investigations include covalently bonded materials as well as systems involving metal atoms. For small clusters calculated results provided information on various energy- and structure-related properties along with vibrational characteristics. Also, energy barriers for configurational transitions were calculated for selected cases. Simulation calculations for surfaces produced new results in areas related to surface energies, reconstructions and relaxations, surface defects, surface stresses as well as adsorption and nucleation processes. Simulations were also performed on bulk materials. Calculated results provided an atomic level understanding on energetics and structures of point defects, crystal stability, elastic properties, and materials strength for various systems. Calculations involving polymeric materials include studies of polar polymer melts and polymer/solid interfaces. Simulations employing specially developed codes provided significant information about energetics and conformational characteristics of different polymeric chains.

Halicioglu, Timur↗

The behavior of unsteady thermocapillary flows

The thermocapillary-driven flow of a liquid in a shallow slot provides a simple model to explore the effects of free-surface deflection and end walls on the stability of the flow. Such an investigation may be useful in understanding the complex flow seen in processes such as the float-zone refining of single crystals. A linear stability analysis of the viscously-dominated slot flow indicates stable basic states to both tow- and three-dimensional infinitesimal disturbances for capillary numbers below 43.3. Above this capillary number steady solutions do not exist. We also discuss current work involving the development of an inertially-dominated slot-flow model using singular perturbation methods. If these flows are unstable, inertial effects would be the only possible cause of the instability.

Vrane, David R.↗

A reattaching free shear layer in compressible turbulent flow - A comparison of numerical and experimental results

An investigation of a two-dimensional, free turbulent shear layer reattaching on an inclined surface at Mach 2.92 and at a high Reynolds number is described. The test geometry is specifically designed to isolate the reattachment process of a high-speed separated flow. A numerical solution of the time-dependent, Reynolds-averaged, Navier-Stokes equations for the entire flow field, employing a two-equation eddy viscosity turbulence model, is presented. Detailed comparisons of prediction and experiment are made in the free shear layer, at reattachment, and in the developing boundary layer downstream. These comparisons include mean surface quantities as well as mean and fluctuating flowfield quantities. Although the overall features of this complex flow field are predicted, there are several deficiencies in the numerical solution, particularly in the region downstream of reattachment. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗

Static Aeroelastic Analysis with an Inviscid Cartesian Method

An embedded-boundary, Cartesian-mesh flow solver is coupled with a three degree-of-freedom structural model to perform static, aeroelastic analysis of complex aircraft geometries. The approach solves a nonlinear, aerostructural system of equations using a loosely-coupled strategy. An open-source, 3-D discrete-geometry engine is utilized to deform a triangulated surface geometry according to the shape predicted by the structural model under the computed aerodynamic loads. The deformation scheme is capable of modeling large deflections and is applicable to the design of modern, very-flexible transport wings. The coupling interface is modular so that aerodynamic or structural analysis methods can be easily swapped or enhanced. After verifying the structural model with comparisons to Euler beam theory, two applications of the analysis method are presented as validation. The first is a relatively stiff, transport wing model which was a subject of a recent workshop on aeroelasticity. The second is a very flexible model recently tested in a low speed wind tunnel. Both cases show that the aeroelastic analysis method produces results in excellent agreement with experimental data.

Aeroelastic↗

Design, test, and evaluation of three active flutter suppression controllers

Three control law design techniques for flutter suppression are presented. Each technique uses multiple control surfaces and/or sensors. The first method uses traditional tools (such as pole/zero loci and Nyquist diagrams) for producing a controller that has minimal complexity and which is sufficiently robust to handle plant uncertainty. The second procedure uses linear combinations of several accelerometer signals and dynamic compensation to synthesize the model rate of the critical mode for feedback to the distributed control surfaces. The third technique starts with a minimum-energy linear quadratic Gaussian controller, iteratively modifies intensity matrices corresponding to input and output noise, and applies controller order reduction to achieve a low-order, robust controller. The resulting designs were implemented digitally and tested subsonically on the active flexible wing wind-tunnel model in the Langley Transonic Dynamics Tunnel. Only the traditional pole/zero loci design was sufficiently robust to errors in the nominal plant to successfully suppress flutter during the test. The traditional pole/zero loci design provided simultaneous suppression of symmetric and antisymmetric flutter with a 24-percent increase in attainable dynamic pressure. Posttest analyses are shown which illustrate the problems encountered with the other laws.

Adams, William M., Jr.↗

Analysis of ground reflection of jet noise obtained with various microphone arrays over an asphalt surface

Ground reflection effects on the propagation of jet noise over an asphalt surface are discussed for data obtained using a 33.02-cm diameter nozzle with microphones at several heights and distances from the nozzle axis. Ground reflection effects are analyzed using the concept of a reflected signal transfer function which represents the influence of both the reflecting surface and the atmosphere on the propagation of the reflected signal in a mathematical model. The mathematical model used as a basis for the computer program was successful in significantly reducing the ground reflection effects. The range of values of the single complex number used to define the reflected signal transfer function was larger than expected when determined only by the asphalt surface. This may indicate that the atmosphere is affecting the propagation of the reflected signal more than the asphalt surface. The selective placement of the reinforcements and cancellations in the design of an experiment to minimize ground reflection effects is also discussed.

Miles, J. H.↗