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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: Methodology and application to high-order compact schemes

We present a systematic method for constructing boundary conditions (numerical and physical) of the required accuracy, for compact (Pade-like) high-order finite-difference schemes for hyperbolic systems. First, a roper summation-by-parts formula is found for the approximate derivative. A 'simultaneous approximation term' (SAT) is then introduced to treat the boundary conditions. This procedure leads to time-stable schemes even in the system case. An explicit construction of the fourth-order compact case is given. Numerical studies are presented to verify the efficacy of the approach.

Carpenter, Mark H.↗
DOE OSTI2021

Boundary Layer Climatology at ARM Southern Great Plains

Operational since 1992, the Atmospheric Radiation Measurement (ARM) southern great plains (SGP) site at Oklahoma, USA has become a reference research site for meteorological studies. Due to an open data policy the ARM data are used by researchers all over the world. In this report, we review the long-term climatology of the atmospheric boundary layer, SGP instrumentations, the site and some site-specific atmospheric conditions which potentially effect wind turbines within the region. As the atmospheric boundary layer is bounded and influenced by the surface, observations of surface radiation components and heat fluxes are crucial in understanding land-atmosphere interactions. The entrainment of air, updrafts, downdrafts and boundary layer height characteristics is needed for understanding the structure and growth of the atmospheric boundary layer. Therefore, measurements from both surface in-situ and remote sensing observations at SGP provide an overall climatology and their interactions from surface up to the boundary layer. Measurements from a 60 m meteorological tower, surface flux stations, disdrometers, soil temperature and moisture flux plates, coherent Doppler lidar, Raman lidar, radiosondes, and satellite data at SGP central facility were analyzed. All the measurements were generally within a few square kilometers of each other at the central facility. This report focuses on data from January 2010 to June 2020 at SGP central facility. The various sections describe the ARM SGP site and surrounding wind turbines; in-situ and remote sensing instrumentation used in the report; provides mathematical equations to analyze fluxes, turbulence and other boundary layer parameters; a climatological analysis of surface winds, fluxes and thermodynamic parameters for several years; an analysis of observed winds in the framework of Monin-Obukhov Similarity Theory; an analysis of the boundary layer winds and direction from a Doppler lidar; multi-year turbulence estimates through the boundary layer from a Doppler lidar; atmospheric boundary layer water vapor and relative humidity profiles from Raman lidar; cloud base height and boundary layer height from multiple sensors and satellite data; and finally site specific atmospheric conditions, such as nocturnal low-level jets. Diurnal, seasonal and yearly variations of surface, sub-surface and boundary layer quantities, such as wind speed, direction, temperature, atmospheric stability, soil temperature, and various atmospheric fluxes at SGP showed distinct trends useful for focused modeling studies. The applicability of surface similarity theory on ARM SGP data is also evaluated, which showed northerly flows are aligned with MO theory estimates compared to southerly flows. Boundary layer winds and direction profiles for several years from a Doppler lidar shows a consistent presence of a nocturnal low-level jet and predominant southerly wind directions through the boundary layer at SGP. The inter-annual variability at SGP is low (<3.5%), with a mean annual wind speed of approximately 7 m s -1 at 100 m above ground level. Boundary layer turbulence and moisture transport from Doppler and Raman lidars are evaluated, which provides evidence of increased water vapor mass flux into the great plains during nocturnal low-level jets. The moisture flux from nocturnal low-level jets is observed to be maximum during summer periods. A novel machine learning algorithm is implemented to accurately estimate the planetary boundary layer height, providing further insights into growth and destruction of the convective boundary layer height during various seasons and land-atmosphere conditions. The high frequency of low-level clouds during winter, spring and fall seasons is validated using the multi-sensor array and satellite estimates of cloud top height. Satellite vegetative fraction data provides insight into seasonal surface roughness and vegetation variability around SGP site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Application of Koopman operator for model-based control of fracture propagation and proppant transport in hydraulic fracturing operation

This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.

42 ENGINEERING↗
NASA NTRS1988

Cretaceous-Tertiary boundary in the Antarctic: Climatic cooling precedes biotic crisis

Stable isotopic investigations were conducted on calcareous microfossils across two deep sea Cretaceous-Tertiary boundary sequences on Maud Rise, Weddell Sea, Antarctica. The boundary is taken at the level of massive extinctions in calcareous planktonic microfossils, and coincides with a sharp lithologic change from pure calcareous ooze to calcareous ooze with a large volcanic clay component. The uppermost Maestrichtian is marked by a long-term decrease in delta value of 0 to 18 which spans most of the lower and middle A. mayaroensis Zone and represents a warming trend which culminated in surface water temperatures of about 16 C. At approximately 3 meters below the K-T boundary this warming trend terminates abruptly and benthic and planktonic isotopic records exhibit a rapid increase in delta value of 0 to 18 that continues up to the K-T boundary. The trend towards cooler surface water temperatures stops abruptly at the K-T boundary and delta value of 0 to 18 values remain relatively stable through the Paleocene. Comparison of the Antarctic sequence with the previously documented deep sea records in the South Atlantic reveal shifts of similar magnitude in the latest Maestrichtian. It is indicated that the Southern Ocean underwent the most significant, and apparently permanent, climatic change. The latest Cretaceous oxygen isotopic shift recorded at Maud Rise and other deep sea sites is similar in magnitude to large positive delta valve of 0 to 18 shifts in the middle Eocene, at the Eocene/Oligocene boundary and in the middle Miocene that marked large scale climatic transitions which ultimately lead to cryospheric development of the Antarctic. The climatic shift at the end of the Cretaceous represents one of the most significant climatic transitions recorded in the latest Phanerozoic and had a profound effect on global climate as well as oceanic circulation.

Stott, Lowell D.↗
DOE OSTI2021

Natively Periodic Fast Multipole Method: Approximating the Optimal Green’s Function

The fast multipole method (FMM) obeys periodic boundary conditions "natively" if it uses a periodic Green's function for computing the multipole expansion in the interaction zone of each FMM oct-tree node. One can define the "optimal" Green's function for such a method that results in the numerical solution that converges to the equivalent particle-mesh (PM) solution in the limit of sufficiently high order of multipoles. A discrete functional equation for the optimal Green's function can be derived, but is not practically useful as methods for its solution are not known. Instead, this paper presents an approximation for the optimal Green's function that is accurate to better than 10 -3 in ${L}_{\mathrm{MAX}}$ norm and 10 -4 in L 2 norm for practically useful multipole counts. Such an approximately optimal Green's function offers a practical way for implementing the FMM with periodic boundary conditions natively, without the need to compute lattice sums or to rely on hybrid FMM-PM approaches.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗
NASA NTRS2007

Algorithm Plans Collision-Free Path for Robotic Manipulator

An algorithm has been developed to enable a computer aboard a robot to autonomously plan the path of the manipulator arm of the robot to avoid collisions between the arm and any obstacle, which could be another part of the robot or an external object in the vicinity of the robot. In simplified terms, the algorithm generates trial path segments and tests each segment for potential collisions in an iterative process that ends when a sequence of collision-free segments reaches from the starting point to the destination. The main advantage of this algorithm, relative to prior such algorithms, is computational efficiency: the algorithm is designed to make minimal demands upon the limited computational resources available aboard a robot. This path-planning algorithm utilizes a modified version of the collision-detection method described in "Improved Collision-Detection Method for Robotic Manipulator" (NPO-30356), NASA Tech Briefs, Vol. 27, No. 3 (June 2003), page 72. The method involves utilization of mathematical models of the robot constructed prior to operation and similar models of external objects constructed automatically from sensory data acquired during operation. This method incorporates a previously developed method, known in the art as the method of oriented bounding boxes (OBBs), in which an object is represented approximately, for computational purposes, by a box that encloses its outer boundary. Because many parts of a robotic manipulator are cylindrical, the OBB method has been extended in this method to enable the approximate representation of cylindrical parts by use of octagonal or other multiple-OBB assemblies denoted oriented bounding prisms (OBPs). A multiresolution OBB/OBP representation of the robot and its manipulator arm and a multiresolution OBB representation of external objects (including terrain) are constructed and used in a process in which collisions at successively finer resolutions are detected through computational detection of overlaps between the corresponding OBB and OBP models. For computational efficiency, the process is started at the coarsest resolution and stopped as soon as possible, preferably before reaching the finest resolution. At the coarsest resolution, there is a single OBB enclosing all relevant external objects and a single OBB enclosing the entire robot. At the next finer level of resolution, the coarsest-resolution OBB is divided into two OBBs, and so forth. If no collision is detected at the coarsest resolution, then there is no need for further computation to detect collisions. If a collision is detected at the coarsest resolution, then tests for collisions are performed at the next finer level of resolution. This process is continued to successively finer resolutions until either no more collisions are detected or the finest resolution is reached.

Backes, Paul↗
NASA NTRS2003

Improved Collision-Detection Method for Robotic Manipulator

An improved method has been devised for the computational prediction of a collision between (1) a robotic manipulator and (2) another part of the robot or an external object in the vicinity of the robot. The method is intended to be used to test commanded manipulator trajectories in advance so that execution of the commands can be stopped before damage is done. The method involves utilization of both (1) mathematical models of the robot and its environment constructed manually prior to operation and (2) similar models constructed automatically from sensory data acquired during operation. The representation of objects in this method is simpler and more efficient (with respect to both computation time and computer memory), relative to the representations used in most prior methods. The present method was developed especially for use on a robotic land vehicle (rover) equipped with a manipulator arm and a vision system that includes stereoscopic electronic cameras. In this method, objects are represented and collisions detected by use of a previously developed technique known in the art as the method of oriented bounding boxes (OBBs). As the name of this technique indicates, an object is represented approximately, for computational purposes, by a box that encloses its outer boundary. Because many parts of a robotic manipulator are cylindrical, the OBB method has been extended in this method to enable the approximate representation of cylindrical parts by use of octagonal or other multiple-OBB assemblies denoted oriented bounding prisms (OBPs), as in the example of Figure 1. Unlike prior methods, the OBB/OBP method does not require any divisions or transcendental functions; this feature leads to greater robustness and numerical accuracy. The OBB/OBP method was selected for incorporation into the present method because it offers the best compromise between accuracy on the one hand and computational efficiency (and thus computational speed) on the other hand.

Leger, Chris↗
NASA NTRS1990

A preliminary analysis of LANDSAT and UK C130 data from Mission 197-A, on 16 July 1987

The co-ordinated LANDSAT-aircraft mission flown on 16 July to the west of San Nicolas Island shows a number of interesting features. There is considerable structure in the Band 4 (0.76 to 0.90 microns) reflectance on km to 10 km scales in otherwise complete cloud cover. In extremely light winds, the C130 performed a series of runs throughout the boundary layer and above for a period of approximately 90 mins on either side of the LANDSAT overpass time. Because of the light winds, the same features can be recognized on the various aircraft penetrations and composited. Furthermore, the accurate renavigation of the aircraft positions enables these to be unambiguously related to the features seen on the LANDSAT image. Early indications suggest that the reflectance features are related to local thickenings on the cloud layer, associated with slightly deeper convection rising into the overlying stratiform deck. If this is true, then LANDSAT images could possibly be used to remotely infer aspects of the convective regime within the boundary layer. In addition, the LANDSAT images also display a quasi-linear feature extending for some 50 km, although only a km or so wide. This feature was, in fact, photographed from the C130 flying just above cloud top and termed a 'cloud-cliff'. It is clearly a rapid and sizeable local change in cloud top height and is therefore presumably a propagating feature. An explanation is still outstanding, but it appears to resemble an undular bore. Further descriptions of this unexpected feature, and its possible origin will be discussed.

Nicholls, S.↗
NASA NTRS2001

Thermocouple Rakes for Measuring Boundary Layer Flows Extremely Close to Surface

Of vital interest to aerodynamic researchers is precise knowledge of the flow velocity profile next to the surface. This information is needed for turbulence model development and the calculation of viscous shear force. Though many instruments can determine the flow velocity profile near the surface, none of them can make measurements closer than approximately 0.01 in. from the surface. The thermocouple boundary-layer rake can measure much closer to the surface than conventional instruments can, such as a total pressure boundary layer rake, hot wire, or hot film. By embedding the sensors (thermocouples) in the region where the velocity is equivalent to the velocity ahead of a constant thickness strut, the boundary-layer flow profile can be obtained. The present device fabricated at the NASA Glenn Research Center microsystem clean room has a heater made of platinum and thermocouples made of platinum and gold. Equal numbers of thermocouples are placed both upstream and downstream of the heater, so that the voltage generated by each pair at the same distance from the surface is indicative of the difference in temperature between the upstream and downstream thermocouple locations. This voltage differential is a function of the flow velocity, and like the conventional total pressure rake, it can provide the velocity profile. In order to measure flow extremely close to the surface, the strut is made of fused quartz with extremely low heat conductivity. A large size thermocouple boundary layer rake is shown in the following photo. The latest medium size sensors already provide smooth velocity profiles well into the boundary layer, as close as 0.0025 in. from the surface. This is about 4 times closer to the surface than the previously used total pressure rakes. This device also has the advantage of providing the flow profile of separated flow and also it is possible to measure simultaneous turbulence levels within the boundary layer.

Hwang, Danny P.↗
NASA NTRS2013

Trip-Induced Transition Measurements in a Hypersonic Boundary Layer Using Molecular Tagging Velocimetry

Measurements of mean streamwise velocity, fluctuating streamwise velocity, and instantaneous streamwise velocity profiles in a hypersonic boundary layer were obtained over a 10-degree half-angle wedge model. A laser-induced fluorescence-based molecular tagging velocimetry technique was used to make the measurements. The nominal edge Mach number was 4.2. Velocity profiles were measured both in an untripped boundary layer and in the wake of a 4-mm diameter cylindrical tripping element centered 75.4 mm downstream of the sharp leading edge. Three different trip heights were investigated: k = 0.53 mm, k = 1.0 mm and k = 2.0 mm. The laminar boundary layer thickness at the position of the measurements was approximately 1 mm, though the exact thickness was dependent on Reynolds number and wall temperature. All of the measurements were made starting from a streamwise location approximately 18 mm downstream of the tripping element. This measurement region continued approximately 30 mm in the streamwise direction. Additionally, measurements were made at several spanwise locations. An analysis of flow features show how the magnitude, spatial location, and spatial growth of streamwise velocity instabilities are affected by parameters such as the ratio of trip height to boundary layer thickness and roughness Reynolds number. The fluctuating component of streamwise velocity measured along the centerline of the model increased from approximately 75 m/s with no trip to +/-225 m/s with a 0.53-mm trip, and to +/-240 m/s with a 1-mm trip, while holding the freestream Reynolds number constant. These measurements were performed in the 31-inch Mach 10 Air Tunnel at the NASA Langley Research Center.

Bathel, Brett F.↗
NASA NTRS1951

Method of successive approximations for the solution of certain problems in aerodynamics

A method of successive approximations for the solution of problems in the fields of diffusion, boundary-layer flow, and heat-transfer is illustrated by solving problems in each of these fields. In most of the examples, the approximate solutions are compared with known accurate solutions and the agreement is shown to be good.

FLOW, TURBULENT - DIFFUSION↗
DOE OSTI2024

A streamline coordinate analysis of a turbulent boundary layer subject to pressure gradients and curvature on the windward side of a bump

Direct numerical simulation (DNS) of a turbulent boundary layer over the Gaussian (Boeing) bump is performed. This boundary layer exhibits a series of adverse and favourable pressure gradients and convex and concave curvature effects before separating. These effects on turbulent boundary layers are characterised and compared with a lower-Reynolds-number flow over the same geometry. The momentum budgets are analysed in the streamline-aligned coordinate system upstream of the separation region. These momentum budgets allow the simplification of equations to facilitate an integral analysis. Integral-analysis-based approximations for Reynolds stresses in the inner and outer regions of the boundary layer are also formulated. The shear and wall-normal Reynolds stress profiles normalised by these approximations exhibit a better collapse compared with friction velocity and Zagarola–Smits normalisations in the strong favourable pressure gradient region and in the mild adverse pressure region that precedes it in this flow. Simplification of these Reynolds stress approximations along with results from the DNS are used to obtain semi-empirical approximations that are able to provide stress closure in terms of wall solution fields for the turbulent boundary layer under consideration.

42 ENGINEERING↗
NASA NTRS1972

Compressible boundary layer with normal pressure gradients: Quasi-similarity equations - Their properties at the wall and at sharp and blunt leading edges.

The failure of most viscous-inviscid interaction methods at strong interactions is attributed to the presence of a normal pressure gradient. A new theory is proposed for supersonic laminar boundary layers that can generate normal pressure gradients. The Navier-Stokes equations are reexamined by an order of magnitude analysis and all first and second order terms are retained. The approximation is found to be dependent not only on the boundary layer thickness but also on the ratio of the dimensionless viscosity and density. The equations are transformed into two quasi-similar, nonlinear, third order, ordinary integro-differential equations for the velocity and pressure as functions of a single transverse variable. The properties of the equations at the boundaries are discussed.

Seginer, A.↗
NASA NTRS1995

Drag measurements of an axisymmetric nacelle mounted on a flat plate at supersonic speeds

An experimental investigation was conducted to determine the effect of diverter wedge half-angle and nacelle lip height on the drag characteristics of an assembly consisting of a nacelle fore cowl from a typical high-speed civil transport (HSCT) and a diverter mounted on a flat plate. Data were obtained for diverter wedge half-angles of 4.0 deg, 6.0 deg, and 8.0 deg and ratios of the nacelle lip height above a flat plate to the boundary-layer thickness (h(sub n)/delta) of approximately 0.87 to 2.45. Limited drag data were also obtained on a complete nacelle/diverter configuration that included fore and aft cowls. Although the nacelle/diverter drag data were not corrected for base pressures or internal flow drag, the data are useful for comparing the relative drag of the configuration tested. The tests were conducted in the Langley Unitary Plan Wind Tunnel at Mach numbers of 1.50, 1.80, 2.10, and 2.40 and Reynolds numbers ranging from 2.00 x 10(exp 6) to 5.00 x 10(exp 6) per foot. The results of this investigation showed that the nacelle/diverter drag essentially increased linearly with increasing h(sub n)/delta except near 1.0 where the data showed a nonlinear behavior. This nonlinear behavior was probably caused by the interaction of the shock waves from the nacelle/diverter configuration with the flat-plate boundary layer. At the lowest h(sub n)/delta tested, the diverter wedge half-angle had virtually no effect on the nacelle/diverter drag. However, as h(sub n)/delta increased, the nacelle/diverter drag increased as diverter wedge half-angle increased.

Flamm, Jeffrey D.↗
DOE OSTI2025

Effects of Strain and Strain Rate on Dynamic Grain Growth and Subgrain Evolution During Plastic Deformation of an Interstitial-Free Steel at 850 ° C

Here, the effects of strain and strain rate on dynamic grain growth (DGG) and subgrain evolution are reported for an interstitial-free steel deformed at 850 ° C. Microstructures produced during tension tests at true-strain rates of 10 -4 and to 10 -3 s -1 true strains ranging from 0.02 to 0.2 were preserved following deformation. These were characterized using electron backscatter diffraction (EBSD), including the application of spherical harmonic transform indexing to produce high-angular-resolution EBSD (HR-EBSD) data. HR-EBSD data resolved the small misorientation angles of subgrain boundaries while imaging much larger data fields than possible with previously available techniques. The resulting data confirmed that steady-state flow stress is inversely proportional to the average subgrain size and that subgrain boundary misorientation angle increases with strain. The following new observations are reported. The rate of DGG increased with respect to time but decreased with respect to strain as strain rate increased. This behavior is rationalized through a simple model using separate rate parameters for the effects of time and strain. Subgrain size was not constant during steady-state deformation, but decreased slowly with increasing strain. Subgrain size distributions and subgrain boundary misorientation angle distributions were measured, and both remained approximately log-normal during steady-state deformation. Subgrain evolution demonstrated no dependence on parent grain size, crystallographic orientation, or Taylor factor. These new data suggest that steady-state flow stress is more likely controlled by the dislocation density internal to subgrains than by the spacing between subgrain boundaries.

dynamic grain growth↗