Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate boundary”

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

Algebraic grid generation using tensor product B-splines

Finite difference methods are more successful if the accompanying grid has lines which are smooth and nearly orthogonal. The development of an algorithm which produces such a grid when given the boundary description. Topological considerations in structuring the grid generation mapping are discussed. The concept of the degree of a mapping and how it can be used to determine what requirements are necessary if a mapping is to produce a suitable grid is examined. The grid generation algorithm uses a mapping composed of bicubic B-splines. Boundary coefficients are chosen so that the splines produce Schoenberg's variation diminishing spline approximation to the boundary. Interior coefficients are initially chosen to give a variation diminishing approximation to the transfinite bilinear interpolant of the function mapping the boundary of the unit square onto the boundary grid. The practicality of optimizing the grid by minimizing a functional involving the Jacobian of the grid generation mapping at each interior grid point and the dot product of vectors tangent to the grid lines is investigated. Grids generated by using the algorithm are presented.

Saunders, B. V.↗

Elastostatic stress analysis of orthotropic rectangular center-cracked plates

A mapping-collocation method was developed for the elastostatic stress analysis of finite, anisotropic plates with centrally located traction-free cracks. The method essentially consists of mapping the crack into the unit circle and satisfying the crack boundary conditions exactly with the help of Muskhelishvili's function extension concept. The conditions on the outer boundary are satisfied approximately by applying the method of least-squares boundary collocation. A parametric study of finite-plate stress intensity factors, employing this mapping-collocation method, is presented. It shows the effects of varying material properties, orientation angle, and crack-length-to-plate-width and plate-height-to-plate-width ratios for rectangular orthotropic plates under constant tensile and shear loads.

Gyekenyesi, G. S.↗

Some approximations in the linear dynamic equations of thin cylinders

Theoretical analysis is performed on the linear dynamic equations of thin cylindrical shells to find the error committed by making the Donnell assumption and the neglect of in-plane inertia. At first, the effect of these approximations is studied on a shell with classical simply supported boundary condition. The same approximations are then investigated for other boundary conditions from a consistent approximate solution of the eigenvalue problem. The Donnell assumption is valid at frequencies high compared with the ring frequencies, for finite length thin shells. The error in the eigenfrequencies from omitting tangential inertia is appreciable for modes with large circumferential and axial wavelengths, independent of shell thickness and boundary conditions.

El-Raheb, M.↗

Carbon Solubility in Silicon-Iron-Bearing Metals during Core Formation on Mercury

Recent results obtained from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft showed the surface of Mercury has high S abundances (approximately 4 wt%) and low Iron(II) Oxide abundances (less than 2 wt%). Based on these extreme values, the oxygen fugacity of Mercury's surface materials was estimated to be approximately 3 to 7 log(sub 10) units below the IW buffer (Delta IW-3 to Delta IW-7). This highly reducing nature of the planet has resulted in a large core and relatively thin mantle, extending to only approximately 420 km depth (corresponding to a core-mantle boundary pressure of approximately 4-7 GPa) within the planet. Furthermore, MESSENGER results have suggested the presence of carbon on the surface of the planet. Previous experimental results from have also suggested the possibility of a primary floatation crust on Mercury composed of graphite, produced after a global magma ocean event. With these exotic conditions of this compositional end-member planet, it begs the question, what is the core composition of Mercury? Although no definitive conclusion has been reached, previous studies have made advances towards answering this question. Riner et al. and Chen et al. looked at iron sulfide systems and implemented various crystallization and layered core scenarios to try and determine the composition and structure of Mercury's core. Malavergne et al. examined core crystallization scenarios in the presence of sulfur and silicon. Hauck et al. used the most recent geophysical constraints from the MESSENGER spacecraft to model the internal structure of Mercury, including the core, in a iron-sulfur-silicon system. More recently, Chabot et al. conducted a series of metal-silicate partitioning experiments in a iron-sulfur-silicon system. These results showed the core of Mercury has the potential to contain more than 15 wt% silicon. However, with the newest results from MESSENGER's low altitude campaign, carbon is another potential light element that could be incorporated into Mercury's core. The goal of this study is to determine the carbon concentration at graphite saturation in various silicon-iron bearing metals relevant to possible mercurian core compositions. Future experiments will include the addition of sulfur into these metals.

Vander Kaaden, Kathleen E.↗

Geologic Mapping of the Mawrth Vallis Region, Mars: MTM Quadrangles 25022, 25017, 25012, 20022, 20017, and 20012

Mawrth Vallis is a 15-25 km wide, 500 km long sinuous channel that winds through the highlands of Arabia Terra and debouches into the lowlands of Acidalia Planitia. The Mawrth Vallis region lies along the gradational zone between southern hemisphere thick crust and northern hemisphere thin crust, a topographically distinct portion of the Martian crustal dichotomy. The origin and age of the dichotomy boundary are controversial and are further complicated by the multi-stage and multi-process geologic history that has modified this approximately 6000 km section of the highland-lowland boundary (approximately 5 N, 330E to approximately 30 N, 80 E; herein referred as the Arabia Terra boundary). Furthermore, the Arabia Terra boundary has been subjected to many post-boundary processes such as outflow floods to the west, volcanism and tectonism to the east, and potential volatile deposition and glacial modification to the north. This study seeks to better understand the history of the Mawrth Vallis region by mapping six MTM quadrangles (17.5-27.5 N, 335-350 E) at 1:1M scale using traditional and modern digital geologic mapping techniques.

Chuang, F. C.↗

Hybrid finite element methods

The purpose of this paper is to show how Lagrange multipliers can be used with finite elements to achieve a number of desirable properties in the underlying approximation. For elliptic boundary value problems, variational principles can be developed in which all boundary conditions are natural. In fluid flow problems, one can endow the approximations with physically essential conservation laws.

Fix, G. M.↗

AgRISTARS. Supporting research: Algorithms for scene modelling

The requirements for a comprehensive analysis of LANDSAT or other visual data scenes are defined. The development of a general model of a scene and a computer algorithm for finding the particular model for a given scene is discussed. The modelling system includes a boundary analysis subsystem, which detects all the boundaries and lines in the image and builds a boundary graph; a continuous variation analysis subsystem, which finds gradual variations not well approximated by a boundary structure; and a miscellaneous features analysis, which includes texture, line parallelism, etc. The noise reduction capabilities of this method and its use in image rectification and registration are discussed.

Rassbach, M. E.↗

Experimental Study of Bypass Transition in a Boundary Layer

A detailed investigation to compare the boundary layer transition process in a low intensity disturbance environment to that in an environment in which the disturbances are initially non-linear in amplitude was conducted using a flat plate model. The transition mechanism based on linear growth of Tollmien Schlichting (T-S) waves was associated with a freestream turbulence level of 0.3 percent; however, for a freestream turbulence intensity of 0.65 percent and higher, the bypass transition mechanism prevailed. The results of detailed measurements acquired to study and compare the two transition mechanisms indicate that there exists a critical value for the peak rms of the velocity fluctuations within the boundary layer of approximately 3 to 3.5 percent of the freestream velocity. Once the unsteadiness within the boundary layer reached this critical value, turbulent bursting initiated, regardless of the transition mechanism. The two point correlations and simultaneous time traces within the transition region illustrate the features of a turbulent burst and its effect on the surrounding flowfield.

Suder, Kenneth L.↗

A Von Karman integral approach to a two phase boundary layer problem

A Von Karman integral approximation of a two phase boundary layer is developed for bodies of arbitrary shape. The flow field considered is that of the injection of water through a porous airfoil. A solution for the special case of a flat plate is presented. The equations for the airfoil solution are developed and possible effects on airflow separation are discussed.

Henry, R.↗

Finite Element A Posteriori Error Estimation for Heat Conduction

This research investigates residual-based a posteriori error estimates for finite element approximations of heat conduction in single-layer and multi-layered materials. The finite element approximation, based upon hierarchical modelling combined with p-version finite elements, is described with specific application to a two-dimensional, steady state, heat-conduction problem. Element error indicators are determined by solving an element equation for the error with the element residual as a source, and a global error estimate in the energy norm is computed by collecting the element contributions. Numerical results of the performance of the error estimate are presented by comparisons to the actual error. Two methods are discussed and compared for approximating the element boundary flux. The equilibrated flux method provides more accurate results for estimating the error than the average flux method. The error estimation is applied to multi-layered materials with a modification to the equilibrated flux method to approximate the discontinuous flux along a boundary at the material interfaces. A directional error indicator is developed which distinguishes between the hierarchical modeling error and the finite element error. Numerical results are presented for single-layered materials which show that the directional indicators accurately determine which contribution to the total error dominates.

Lang, Christapher G.↗

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗

Analysis of separated boundary-layer flows

A method for predicting the strong interaction between the viscous and inviscid flows which occurs in flow separation is reviewed. An inverse boundary-layer procedure approximately accounts for normal pressure gradients that may be important in strongly interacting flows. Transformed boundary-layer equations are written in which the pressure gradient is set equal to the inviscid pressure gradient. As the boundary-layer edge is approached and the viscous shear and heat conduction terms vanish, the viscous flow solution is required to asymptotically approach the inviscid solution over the generalized displacement body. Attention is then focused on viscous-inviscid interacting flows with a first-order viscous formulation and constant pressure across the boundary layer. Results obtained with this procedure are presented for: (1) transitional separation bubbles near an airfoil leading edge, (2) subsonic boattail separated turbulent flow, and (3) transonic turbulent shock wave boundary-layer interaction on an axisymmetric bump configuration.

Carter, J. E.↗

Middle atmosphere project. A semi-spectral numerical model for the large-scale stratospheric circulation

The complete model is a semispectral model in which the longitudinal dependence is represented by expansion in zonal harmonics while the latitude and height dependencies are represented by a finite difference grid. The model is based on the primitive equations in the log pressure coordinate system. The lower boundary of the model domain is set at the 100 mb level (i.e., near the tropopause) and the effects of tropospheric forcing are included in the lower boundary condition. The upper boundary is at approximately 96 km, and the latitudinal extent is either global or hemispheric. The basic differential equations and boundary conditions are outlined. The finite difference equations are described. The initial conditions are discussed and a sample calculation is presented. The FORTRAN code is given in the appendix.

Holton, J. R.↗

Three-dimensional transonic nacelle/inlet flowfield computations using an efficient approximate factorization algorithm

A highly efficient computer analysis has been developed for predicting transonic nacelle/inlet flowfields. This algorithm can compute the three-dimensional transonic flowfield about axisymmetric or asymmetric nacelle/inlet configurations at zero or nonzero incidence. The flowfield is determined by solving the full-potential equation in conservative form on a body-fitted curvilinear computational mesh. The difference equations are solved using the AF2 approximate factorization scheme. The effects of boundary layer viscous entrainment are approximated in the inviscid algorithm by applying a surface transpiration velocity which is determined from the calculated boundary layer growth. Computed results and correlations with existing methods and experiment are presented to illustrate application of the analysis.

Vadyak, J.↗